Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Failing to re-alligning sequences with same truncLen parameters #2089

Open
Debora6991 opened this issue Mar 4, 2025 · 0 comments
Open

Failing to re-alligning sequences with same truncLen parameters #2089

Debora6991 opened this issue Mar 4, 2025 · 0 comments

Comments

@Debora6991
Copy link

Hello,
I've been an intensive user of your dada2 pipeline until a few years ago. After 4 years, I need to re-run the same data in the pipeline to map the reads agains an updated database (2x300 paired end, V3-V4 region).
My samples present with these quality profiles after removing the primers with Cutadapt:

Image

I've tried to fine-tune the truncLen parameters from 200 to 250 for forward reads and 160 to 190 for reverse reads in all possible combinations, but I always get this output (or very similar) after the merging step

Image

The only combination that raises the merged reads a bit (15%) is 200 forward/170 reverse. Above and below these values, the merged reads drop again to 1.9%.

I am wondering what I am doing wrong, as a few years back I successfully merged reads with an average output abundance close to 50%. I compared the script used back then and now, and nothing significant has changed.

Could you please provide some guidance on how to resolve this issue? Thank you!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant