Taxidlist output question (BLAST+)
0
0
Entering edit mode
3.8 years ago
mbrav005 • 0

hello all, I am wondering how when I use a Taxidlist of only 169 unique taxids in my blastp search, When I use this taxidlist in my next blastp search, i get 177 uniq taxids in the output blast search? How am I gaining 8 new unique taxids when I only started with 169? Shouldn't I only get less or atleast the same amount of Taxids if im restricting my blastp search to only 169?

example:

$ wc -l S5_rpoB.taxids
169 S5_rpoB.taxids

$ uniq S5_rpoB.taxids | wc -l
169

  $  blastp -db swissprot -query Ecoli_rpoB.fasta -taxidlist S5_rpoB.taxids -outfmt  '6 staxid' | uniq | wc -l
    177

Where are these 8 new taxids coming from and why?

Thank you

BLAST • 926 views
ADD COMMENT
0
Entering edit mode

TaxID's are hierarchical and it is possible that you are getting something a level higher or level lower (especially if there has been some reassignments of ID's). Something to check into.

ADD REPLY
0
Entering edit mode

Ok, thats also what I was thinking along the lines of, thank you

ADD REPLY

Login before adding your answer.

Traffic: 2417 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6