Blastp filter by Identity

Question

[sureshhewabi]

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+

Answer 1

[lieven.sterck]

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

[sureshhewabi]

Thanks. I will check what is % positives

Answer 2

[Pierre Lindenbaum]

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

[sureshhewabi]

Thanks for the reply, BTW I am using TAB format