This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Remove specific sequence group in fasta file

I have several files with the sequence of the organism or species and its reference sequence (CDS) and I would like to eliminate the reference sequences from them leaving only the sequence of the organism.

sequence genome r gene software error
• 1,305 views • link
written by USER • USER

1 answer

One way would be to linearize the fasta sequences (courtesy of @Pierre's gist which can be easily found by search for linearize fasta). Then grep "^gb" to keep the sequences you want and reformat back to fasta.

awk '/^>/ {printf("%s%s\t",(N>0?"\n":""),$0);N++;next;} {printf("%s",$0);} END {printf("\n");}' < input.fa  | grep "^>gb" | tr "\t" "\n" > final.fa
• 0 views • link
written by GenoMax GenoMax

I would like to remove all that start with lcl. The reference CDS. If so, does it work too?

• 0 views • link
written by USER • USER

Try the above. It will only keep sequences that start with >gb.

• 0 views • link
written by GenoMax GenoMax

Log in to answer this question.

More posts like this