Does DADA2 can give me ASV results from pair end reads after denoising step even if they necessarily don't get merged? And if I'm getting the ASVs by this way(pair-ends not overlapping and hence not merging) are those reliable?
I didn't het this question anywhere before hence posting.
I'm running QIIME 2 2024.5 Amplicon Distribution using conda.
All reads that don't merge are removed! They don't appear as ASVs or as counts in the table.
pair-ends not overlapping and hence not merging) are those reliable?
Uh... yes, they are reliably absent, because DADA2 throws them away.
This is why it's so important to pick DADA2 trimming settings that maximize the number of reads that can be merged!
This process is covered in this part of the Atacama Soils tutorial.
Do your reads overlap enough for them to merge?
(I've edited your post so it has a short title and the full question is in the body of the post.)
Oh. Got you. And I will have to shift to other tools for the downstream analyses, am I right? So, the take home is I can't proceed with pair reads which are just "concatenated" in Qiime2?
Discontiguous amplicons break the assumptions made by many tools, including Qiiem2.
There is interest in supporting them, but new tools and methods need to be made.
It sounds like you are looking to use tools, not build new ones, which is fair!
So, one more question, is there is separate/other dedicated forum for query regarding, analysis using DADA2 in R and not in Qiime2? Or can I continue here?
For DADA2, lots of questions have been asked and answered on the GitHub:
I want to remind you that many of the questions you have while learning about a new method have already been asked by other folks who have learned before you. Checking to see if the question has already been asked and finding the following discussion is often faster then opening a new thread.
Basically, remember to do a lit review! Especially on the internet!