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

--no-trimming leads to significant drop in BUSCO scores #77

Open
JWDebler opened this issue Jun 5, 2023 · 0 comments
Open

--no-trimming leads to significant drop in BUSCO scores #77

JWDebler opened this issue Jun 5, 2023 · 0 comments

Comments

@JWDebler
Copy link

JWDebler commented Jun 5, 2023

Hi,
I just came accross something strange when comparing my genome after different rounds and settings of polishing.

I assembled Q20 Nanopore reads with Canu and then polished once with racon with or without the --no-trimming option, because in the past I have noticed that telomeres get removed without the --no-trimming option.
But this time the problem seems different. When the Canu assembly gets polished with the --no-trimming option, the assembly actually seems to be getting worse.

image

Going back over some 'old' assemblies from ONT R9 series flowcell and chemistry I see the same thing. Only happens when --no-trimming is used though. I always run Medaka after racon which seems to 'rescue' whatever racon was doing, but this makes me wonder if I should skip racon completely in the future.

My full command is:
minimap2 assembly.fasta reads.fastq.gz > minimap.racon.paf

racon -m 8 -x -6 -g -8 -w 500 -t 14 --no-trimming reads.fastq.gz minimap.racon.paf assembly.fasta > assembly_racon.fasta

I have tried this with the latest version 1.5.0 as well as the previous 1.4.21.

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