This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Taxidlist output question (BLAST+)

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

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.

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

0 answers

No answers yet.

Log in to answer this question.