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

Major drop quality score with --trim adapters #53

Open
felipebatalini opened this issue Nov 13, 2024 · 0 comments
Open

Major drop quality score with --trim adapters #53

felipebatalini opened this issue Nov 13, 2024 · 0 comments
Labels
question Further information is requested

Comments

@felipebatalini
Copy link

felipebatalini commented Nov 13, 2024

Why are my q scores dropping so much with --trim adapters?

We are using FLO-MIN114 and R10 chemistry for a cDNA library derived from human RNA.
We noticed a high percentage (>50%) of unusable reads detected by pychopper with the wf-transcriptome workflow, and then identified that it could be improved to <10% if we trimmed the adapters (therefore keeping primers).

However, I was surprised to see a significant drop in the quality scores when I turn on --trim adapters:
nextflow run epi2me-labs/wf-basecalling \ -profile singularity \ --sample_name $sample_name \ --input $pod5_dir \ --dorado_ext pod5 \ --basecaller_cfg [email protected] \ --qscore_filter 10 \ --basecaller_args "--trim adapters" \ --output_fmt fastq \ --out_dir $results_folder
While it makes sense for pychopper to work better with the primers present, I can't understand while the basecalling quality drops do much. In the example below, I demonstrate the different q scores from the same sample.
BC10_without_trimming_parameters
BC10_with---trim adapters on

I appreciate any help to understand this!
Felipe

@felipebatalini felipebatalini added the question Further information is requested label Nov 13, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
question Further information is requested
Development

No branches or pull requests

1 participant