Thanks to visit codestin.com
Credit goes to github.com

Skip to content

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

@JWDebler

Description

@JWDebler

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.

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions