# Taxonomy results of mothur do not match with NCBI blast results **URL:** \ **Category:** Uncategorized\ **Created:** [September 23, 2022, 10:45am UTC](https://forum.mothur.org/t/taxonomy-results-of-mothur-do-not-match-with-ncbi-blast-results/21562 "2022-09-23T10:45:29Z")\ **Posts on this page:** 6\ **Page:** 1
**Author:** ![Hira](https://avatars.discourse-cdn.com/v4/letter/h/7993a0/32.png) [@Hira](https://forum.mothur.org/u/Hira)\ **Post date:** [September 23, 2022, 10:45am UTC](https://forum.mothur.org/t/taxonomy-results-of-mothur-do-not-match-with-ncbi-blast-results/21562/1 "2022-09-23T10:45:29Z")
I have a question related the cons.taxonomy table generated via mothur. I followed the MiSeq SOP ([MiSeq SOP](https://mothur.org/wiki/miseq_sop/)). After removing the mock groups from the data, the final: fasta file, count table and taxonomy file remained for further analysis. I clustered the sequences into otus using agc method, then created the otu table (shared file) and finally classified the otus. mothur \> cluster(fasta=final.ng.fasta, count=final.count\_table, vsearch=/users/desktop/mothur/vsearch, method=agc) mothur \> make.shared(list=final.ng.agc.list, count=final.count\_table, label=0.03) mothur \> classify.otu(list=final.ng.agc.list, count=final.count\_table, taxonomy=final.taxonomy, label=0.03) I selected OTU 2 from list file: Otu00002 M06339\_192\_000000000-KBLJK\_1\_2104\_20594\_11589 I searched for this tag “ **M06339\_192\_000000000-KBLJK\_1\_2104\_20594\_11589** ” in my final.fasta file. It was as follows: > M06339\_192\_000000000-KBLJK\_1\_2104\_20594\_11589 > TAC–GT-AG-GGT----GCG-A-G–C–G–T-T–AA-T-CGG-AA–TT-A-C-T–GG-GC–GT-A–AA-GC-GC-GC----G-TA-G-G-T-G-------G–T-TT-G-G-T-------AA----G-A-T-G--------G-A-T–G–TG–A-AA-TC–C-C-CG-G-G--------------CT-C-AA------------C-C-T-G-G-G-A–A-C----T-G–C-A–T–C—C–AT-A-A—C------T–G-C–CT–G-A-C------------------------------------------------------------------------------T-A-G-A-G-T–A----C-GG----TA-G-A----G-G-G-T—GG-T---------GG–A–ATT----T-C-C-T-GT–GT-A-G-CG-GT–G-A-A-A----TG-C-GT-AG–AT-A-TA-----------G-G-----A-A–G-G-A-AC-A-CC-------------AG–T–G–GC-GAA-G–G-C—G—A–C-C-A-C—CTG–G–GC-T-C------------------------A-T------A-C-T–GA–CA–C–T-G–A-GG–T-G-CG-A–AA-G-C-----G-TG–GG-G–AG-C-A-AA–CAGG I copied the sequence and did standard nucleotide blast as follows: ![Screenshot 2022-09-23 at 3.39.44 PM](https://canada1.discourse-cdn.com/flex036/uploads/mothur/original/1X/7d651fc928f50bddd03dfaffcf714db0b1e953ea.png) The taxonomy of OTU 2 is not the same as in the taxonomy file generated using classify.otu command: ![Screenshot 2022-09-23 at 3.42.33 PM](https://canada1.discourse-cdn.com/flex036/uploads/mothur/original/1X/b92ba85a1b761aae9b3f937cbee1767e7b15046a.png) Can you please explain me that why are there such mismatches in data? And am i checking this correctly? If not then what can i do to check weather the sequences in my fasta file with the taxonomy generated by mothur match with blast results as well? Regards Hira ---
**Author:** ![Alexandre\_Thibodeau](https://avatars.discourse-cdn.com/v4/letter/a/e480ec/32.png) [@Alexandre\_Thibodeau](https://forum.mothur.org/u/Alexandre_Thibodeau)\ **Post date:** [September 26, 2022, 1:24pm UTC](https://forum.mothur.org/t/taxonomy-results-of-mothur-do-not-match-with-ncbi-blast-results/21562/2 "2022-09-26T13:24:41Z")
Hello! What database did you use at what cutoff? Latest version of Silva using cutoff of 80 is usually making it right for me. Ran into that kind of problems with older versions of taxonomy references database where, using my mock community control. For example, Pseudomonas was ok with RDp but not with Silva or vice versa I do not remember correctly. Had to lower the cutoff to 75 to get this peculiar genus right. Anyway, how is your positive control looking? ---
**Author:** ![pschloss](https://yyz2.discourse-cdn.com/flex036/user_avatar/forum.mothur.org/pschloss/32/4_2.png) [@pschloss](https://forum.mothur.org/u/pschloss)\ **Post date:** [September 30, 2022, 12:52pm UTC](https://forum.mothur.org/t/taxonomy-results-of-mothur-do-not-match-with-ncbi-blast-results/21562/3 "2022-09-30T12:52:29Z")
BLAST and the RDP call the same sequence different things. Also, you are look at the names that the scientist who deposited the sequence gave the sequence, which is not necessarily correct or as robust as a true taxonomy. Finally, the algorithms you are using are fundamentally different. The naive bayesian classifier used in classify.seqs is far more robust than blasting a sequence and taking the top match. Pat ---
**Author:** ![Hira](https://avatars.discourse-cdn.com/v4/letter/h/7993a0/32.png) [@Hira](https://forum.mothur.org/u/Hira)\ **Post date:** [October 6, 2022, 5:55pm UTC](https://forum.mothur.org/t/taxonomy-results-of-mothur-do-not-match-with-ncbi-blast-results/21562/4 "2022-10-06T17:55:30Z")
I used reference files from [Silva reference files](https://mothur.org/wiki/silva_reference_files/). Release 138.1. \* [Full-length sequences and taxonomy references](https://mothur.s3.us-east-2.amazonaws.com/wiki/silva.nr_v138_1.tgz)(128884 bacteria, 2846 archaea, and 14871 eukarya sequences). This reference could be customized for alignments, but could also be used for classification. The uncompressed version is ~6.8 GB and the compressed version is 241 MB. I did not use any cutoff parameter in classify.otu command. The default is 51. Moreover, I did another thing. I found out the representative sequence of each OTU using the following: mothur \> get.oturep(fasta=final.ng.fasta, count=final.count\_table, method=abundance, list=final.ng.agc.list) You did not provide a label, using 0.03. 0.03 17032 Output File Names: /users/desktop/mothur/untitled folder/final.ng.agc.0.03.rep.count\_table /users/desktop/mothur/untitled folder/final.ng.agc.0.03.rep.fasta I did ncbi blast of sequences of otus from final.ng.agc.0.03.rep.fasta file. The blast results were similar to the taxonomy given in my final.ng.agc.0.03.cons.taxonomy file. Is this method correct? ---
**Author:** ![Alexandre\_Thibodeau](https://avatars.discourse-cdn.com/v4/letter/a/e480ec/32.png) [@Alexandre\_Thibodeau](https://forum.mothur.org/u/Alexandre_Thibodeau)\ **Post date:** [October 12, 2022, 1:14pm UTC](https://forum.mothur.org/t/taxonomy-results-of-mothur-do-not-match-with-ncbi-blast-results/21562/5 "2022-10-12T13:14:04Z")
Almost, use cutoff of 0.8 for the classify.otu. Follow the online SOP and you should be fine until you get familiar with Mothur and start tailoring the pipeline to your liking. ---
**Author:** ![pschloss](https://yyz2.discourse-cdn.com/flex036/user_avatar/forum.mothur.org/pschloss/32/4_2.png) [@pschloss](https://forum.mothur.org/u/pschloss)\ **Post date:** [October 14, 2022, 3:17pm UTC](https://forum.mothur.org/t/taxonomy-results-of-mothur-do-not-match-with-ncbi-blast-results/21562/6 "2022-10-14T15:17:15Z")
Also, if you want a taxonomy for an OTU, I strongly encourage people to use `classify.otu` rather than `get.oturep` and classifying the output. You can see how we do this in the [SOP](https://mothur.org/wiki/miseq_sop/) Pat