Entering edit mode
4.4 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
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