Patg13
(Patrick Gagné)
October 23, 2025, 7:02pm
1
Hello,
To give you more context, I do metabarcoding (amplicon 16S, ITS, COI) on environmental samples. Today,
I receive new dataset from Illumina new sequencing system (i100) which uses a new way to call quality per base (Reference: https://knowledge.illumina.com/instrumentation/miseq-i100-series/instrumentation-miseq-i100-series-reference_material-list/000009540 )
Now, only 4 quality are possible (2,12,24 and 38), so I got a high percentage of my reads that have 38 phred score (G) for every bases. What would be the best course of action to work with these reads. Normally, I use Dada2 to trim and reassemble my forward and reverse read to obtain my amplicon, but since Dada2 need to generate an error model, will it still work correctly when there is almost no quality variation ?
I don’t know if I’m clear enough, I can provide more info if needed.
Thank you in advance.
1 Like
Hello Pat,
Good news! It's a known issue and Ben is investigating!
opened 02:26PM - 09 Sep 25 UTC
I am working with a shorter dataset from the new MiSeq i100 and I am running int… o problems with the learnErrors function in dada2. As this is the first time that I have used this pipeline for the new MiSeq I chose a few samples (11 samples) to use as a test case. When I run learErrors I receive a message stating that the error rates can not be estimated. Looking back at the plots of the quality scores it seems like there is not much variation in the quality scores and I was wondering if this might be the problem? I have included pictures of the quality score plots and the error code that I received running default settings. I also increased the amount of bases used by learnErrors and received the same error.
I know that this step is essential in this pipeline as the error rate matrix is used in following steps. I am not sure if there is already a fix for this type of data/problem or if there is a work around that I can use to move forward with this data?
Thanks for all of your help!
Schyler
<img width="699" height="67" alt="Image" src="https://github.com/user-attachments/assets/b75ce52a-d4b5-48c3-8a46-de81726497cc" />
<img width="307" height="205" alt="Image" src="https://github.com/user-attachments/assets/2d5c6de5-21bd-413e-8dfd-0db285fea724" />
<img width="956" height="878" alt="Image" src="https://github.com/user-attachments/assets/936abbfe-c257-4f26-8867-245984aa5bdb" />
<img width="956" height="878" alt="Image" src="https://github.com/user-attachments/assets/75fdb29e-714b-4054-8cdc-568e2507c30b" />
The Illumina MiSeq 100 is another one of Illumina's "two-color" instruments: 2-Channel SBS Technology | Faster sequencing and data acquisition
Interesting!
Bad news, there's no fix for DADA2 yet, and once that's working we still need to get this merged into Qiime2.
Thank you for your patience.
1 Like