This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Blastp filter by Identity

I want to filter blast results by identity. I can find an option/parameter in blastn which is -perc_identity. but it is not available in blastp. I can get the blast result and the filter out results as a post process. but can't I do it by specifying a parameter like -perc_identity in blastp? why that option is missing?

BLAST 2.7.1+

sequence
• 4,560 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by sureshhewabi • sureshhewabi

2 answers

I think it's mainly because identity does not mean much on protein level (at least not as it does for blastn, nucleotide level), for protein comparison % positives is more informative. Nonetheless that's also not an option for blastP though.

In general it's not the best approach to let blast "filter" your results (look around on biostar and/or internet for reasons why) , so doing it in post processing is the preferred way , and easily done with a bit of awk for instance

• 1 views • link
written by lieven.sterck lieven.sterck

if you're using a xml output, using a xslt stylesheet, one can filter the Hsp on the identity.

xsltproc --novalid  biostar362495.xsl input.blastp.xml

• 1 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Thanks for the reply, BTW I am using TAB format

• 1 views • link
written by sureshhewabi • sureshhewabi

Log in to answer this question.

More posts like this