Same set of analysis was done using mothur pipeline by company. When comparing with mothur results, I got more number of taxonomic classification with more abundance.
Note: We used control (buffer)rep-seqs.qza (10.2 KB) table.qza (10.8 KB) taxonomy_silva_200.qza (25.5 KB)
data for this analysis. We expect very less or no taxonomical classification.
I want to make sure I understand the problem. You ran the same pipeline in QIIME 2 and Mothur. So, you did your denoising via DADA2 in QIIME2. From my understanding, you also ran de novo OTU picking in Mothur? How did you perform taxonomic assignment in Mothur? What algorithm and what reference database?
Is your total number of features different between the Dada2 table and the Mothur table? Iād expect Dada2 to result in more features overall (since it handles ASVs and is able to give higher resolution than OTU-based clustering at 99%), but possibly fewer counts. (I am, to be honest, less familiar with the QC protocols under the hood of each algorithm.)
If you look at the over all taxonomy plots, how do they compare? If you use a metric like Bray-Curtis distance or Jaccard distance, and perform a procrustes, does the relationship between the two methods remain similar, or do you have large scale differences in the inter-sample community structure?