Dear all,
After processing the data with DADA2, I found that the number of sequences in the sample was 3W-40W
But the sequencing company said it was normal.
Then I filter the dataset by bacteria
qiime taxa filter-table \
--i-table ../table-dada2-3.qza \
--i-taxonomy ../taxonomy-3.qza \
--p-include d__Bacteria \
--o-filtered-table Bacteria-dada2-table.qza \
# filter sequence
qiime taxa filter-seqs \
--i-sequences ../rep-seq-dada2-3.qza \
--i-taxonomy ../taxonomy-3.qza \
--p-include d__Bacteria \
--o-filtered-sequences Bacteria-sequences.qza \
qiime feature-table summarize \
--i-table Bacteria-dada2-table.qza \
--o-visualization Bacteria-dada2-table.qzv \
--m-sample-metadata-file ../metadata_qiime2.txt
Bacteria-dada2-table.qzv
Results:
Align to tree
qiime phylogeny align-to-tree-mafft-fasttree \
--i-sequences Bacteria-sequences.qza \
--p-n-threads auto \
--o-alignment aligned-Bacteria-rep-seqs.qza \
--o-masked-alignment masked-aligned-Bacteria-rep-seqs.qza \
--o-tree Bacteria-unrooted-tree.qza \
--o-rooted-tree Bacteria-rooted-tree.qza
According to the search, this error may be due to too many ASVs.
So, what do I need to do with my data after DADA2?
Rarefaction or filter out the lower frequence before rarefaction to reduce the ASV number?





