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

about plotQualityProfile and derepFastq function #2003

Open
dstrgcv opened this issue Aug 22, 2024 · 1 comment
Open

about plotQualityProfile and derepFastq function #2003

dstrgcv opened this issue Aug 22, 2024 · 1 comment

Comments

@dstrgcv
Copy link

dstrgcv commented Aug 22, 2024

Q1: truncLen parameter in filterAndTrim function is determined by plotQualityProfile. however, I think its objective to see the green line to select the values. is there any standard such as the threshold value of quality score to select the best truncLen parameter?

Q2: I see the current tutorial omitted the derepFastq procedure, but the previous one have derepFastq ofter filterAndTrim

@benjjneb
Copy link
Owner

Q2: I see the current tutorial omitted the derepFastq procedure, but the previous one have derepFastq ofter filterAndTrim

derepFastq now happens "on the fly" within the learnErrors and dada functions. This reduces memory requirements dramatically on large datasets, as only one sample needs to be loaded into memory at a time instead of all at once as in the previous tutorial version.

Q1: truncLen parameter in filterAndTrim function is determined by plotQualityProfile. however, I think its objective to see the green line to select the values. is there any standard such as the threshold value of quality score to select the best truncLen parameter?

I wouldn't worry too much about trying to find the optimum truncation length. The goal is to cut off the tails of sequences after a "quality crash", as is often seen towards the end of Illumina reverse reads, while maintaining enough overlap between paired reads to allow them to merge after truncation.

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

2 participants