My sequence count reduced immensely in all sample, For instance, Sample 12 lacks nearly 129,806 reads after dada analysis.
My forwards reads quality was good rather than reverse reads, So I truncated low quality reads from reverse reads. I kept phred score value as 20 for truncation. After truncation, My sequence count was reduced even more worse than previous dada analysis result ( without truncation ).
Hey there,
ITS processing before dereplication is a little bit different than 16S, and so a “regular pipeline”. Dada2 provides an ITS tutorial where you can have a look. https://benjjneb.github.io/dada2/ITS_workflow.html
Yes — we have a basic ITS tutorial that you can follow here (though note that as with 16S analysis there are many different possible workflows with QIIME 2 for analyzing ITS and you may want to check out plugins like q2-itsxpress and others to help with your ITS analysis):
Well maybe yes and maybe no. You need to check out the dada2 stats file to assess.
If you are losing these sequences at the filtering stage or any stage prior to read merging (joining) in the dada2 pipeline, then the loss should be random and will not impact taxonomic composition profiles or diversity estimates unless if loss is extreme.
If you are losing too many sequences at the merging step, then this loss will be non-random and will severely skew taxonomic profiles. ITS is hypervariable in length and sequences lost at the merging stage will bias for taxa with longer ITS that cannot merge. If this is the case, you either need to adjust your trimming settings or quite possibly only process your ITS data as single-end reads (e.g., only use the forward reads in your analysis).
