This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Filtering fasta file based on identifier

Hi, I have a fasta file with many segments and I want to filter out all the segments that have a "P" in the identifier of the segment. Is there a conventional way to do so? Thanks.

fasta filter
• 5,804 views • link
updated by Ram Ram • written by jyu429 jyu429

Thank you!

• 0 views • link
updated by Ram Ram • written by jyu429 jyu429

2 answers 

bioawk -c fastx '$name ~ /P/ { print ">"$name; print $seq }' <sequences.fa

If you wanna take all except those with a "P",

bioawk -c fastx '$name ! /P/ { print ">"$name; print $seq }' <sequences.fa

bioawk here

• 0 views • link
written by Ram Ram

Neat. I've not found a use for bioawk before but this seems perfect.

• 0 views • link
written by Matt Shirley Matt Shirley

It kinda clicked out of the blue for me yesterday. Now I'm gonna add this to my arsenal of regular-use tools :-)

• 0 views • link
written by Ram Ram

how can give the transcriptome.fasta and headerlist.txt in this command?

• 0 views • link
written by Shahzad • Shahzad

What are those two files?

• 0 views • link
written by Ram Ram

just awk

awk '/^>/{N=0} /^>P/{N=1} {if(N)print}' *.fa
• 0 views • link
updated by Ram Ram • written by Pierre Lindenbaum Pierre Lindenbaum

Maybe

/^>\S*P\S*/

To match identifiers (up to the first space) that contain P rather than just identifiers that start with P.

• 0 views • link
updated by Ram Ram • written by Rob Syme Rob Syme

Would this not print only headers, Pierre?

• 0 views • link
written by Ram Ram

no, if there is no 'next' statement, awk continues to scan all the patterns.

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Oops, I read it wrong. I read it as the if(N) print being in the same {} as the N=1. My bad!

• 0 views • link
written by Ram Ram

But where does the ouput go? Sorry for my ignorance.

• 0 views • link
written by jahn.davik • jahn.davik

"standard out" or "stdout". You can redirect this to a file like:

awk '/^>/{N=0} /^>P/{N=1} {if(N)print}' *.fa > out.fa
• 0 views • link
written by Matt Shirley Matt Shirley

Log in to answer this question.

More posts like this