Question: Error correcting PacBio reads prior to alignment with bwa mem
gravatar for aindap
2.0 years ago by
United States
aindap110 wrote:

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?

pacbio bwa • 841 views
ADD COMMENTlink modified 2.0 years ago by Charles Warden6.5k • written 2.0 years ago by aindap110
gravatar for Charles Warden
2.0 years ago by
Charles Warden6.5k
Duarte, CA
Charles Warden6.5k wrote:

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.

ADD COMMENTlink written 2.0 years ago by Charles Warden6.5k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1373 users visited in the last hour