Mercurial > repos > pcingola > snpeff

view snpEff_2_1a/snpEff_2_1a/scripts/data/c_briggsae_WS230/create.sh @ 3:c052639fa666 default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author pcingola
date Fri, 20 Apr 2012 11:22:59 -0400
parents
children
line wrap: on
line source

#!/bin/sh -e

zcat c_briggsae.WS230.annotations.gff3.gz \
	| grep -v SNP \
	| grep -v translated_nucleotide_match \
	| grep -v repeat_region \
	| grep -v inverted_repeat \
	| grep -v tandem_repeat \
	| grep -v nucleotide_match \
	> genes.gff

# Add fasta sequence
echo "###" >> genes.gff
echo "##FASTA" >> genes.gff
zcat c_briggsae.WS230.genomic.fa.gz >> genes.gff