Error correcting PacBio reads prior to alignment with bwa mem
I have PacBio reads of my organism of interest that I would like to align to a reference sequence pulled from GenBank. Does it make sense to error correct and trim the PacBio reads, with a tool like canu and then align with bwa mem?
• 2,073 views
•
link
updated
by
Charles Warden
•
written
by
aindap
•
written
by
aindap
1 answer
I would say it is OK to use CCS reads, where it might be a good idea to impose 5x or 10x cutoffs (the default is 3x).
I'm not as certain about error-corrected subreads if the library size is too large to define CCS reads. I guess it should be OK, but I can't really provide you with evidence from 1st hand experience either way. I think of that as mostly an internal step for de novo assembly.
• 0 views
•
link
written
by
Charles Warden
Log in to answer this question.
More posts like this
-
Attempting to Use circlator running into issues with CANU
written by dlera.lozano •
-
Mapping a small number of contigs to a reference subsequence
written by alexandru.bologa.marian •
-
mapping long-reads to a reference library
written by mthm •
-
Align pacbio subreads against assembled pacbio reads
written by wes •
-
Canu diploid assembly with nanopore reads
written by Irsan
-
Canu assembling with multiple pacbio reads files at a single run
written by ranjeet_maurya •
-
slurm script for canu assembler
written by bio_d •
-
denovo pacbio assembly
written by 1234anjalianjali1234 •
-
Canu: failed with exit code 134
written by biomagician
-
finding repeats in de-novo assembled contig from PacBio reads
written by aindap
