This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to compute the number of unique proteins hit found in the blastx search

Dear all,

Could anybody please let me know how to calculate the number of unique proteins hit found in the blastx search or the number of unique contigs hit by proteins in the tblastn search and create a unique sequence list? Sorry if you find it's so basic question, but I have to ask for making sure. Thanks

blast rna-seq alignment
• 2,405 views • link
updated by Ram Ram • written by seta seta

How to parse the blast results will depend on the output format.

• 0 views • link
written by h.mon h.mon

The output format is tabular (-outfmt 6)

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

1 answer

In case you used tabular output this is as simple as applying sort..

E.g., how many of my query seqs got a hit?

cut -f1 outputFile | sort -u | wc -l

How many subjects did my query seqs hit?

cut -f2 outputFile | sort -u | wc -l

The first commands works when each query ID is unique, which is not the case in at least blastx. For that, you could also cut the fields that specify where the alignment starts and ends, but it's complicated because you would still likely have multiple hits that differ slightly in their start and end coordinates. I circumvent such problems by using actual protein prediction algorithms (like Prodigal) and applying blastp. It's way faster this way too..

• 0 views • link
updated by Ram Ram • written by 5heikki 5heikki

Thanks friend, however I would like to have the number and list of sequences with unique hit for further analysis. Could you please help me out to access them?

• 0 views • link
updated by Ram Ram • written by seta seta
If you mean ids of contigs that your queries hit, just drop the last pipe.
• 0 views • link
written by 5heikki 5heikki

Log in to answer this question.

More posts like this