Question: Subset BAM file based on selected coverage, PacBio reads
gravatar for MyrNat
2.4 years ago by
MyrNat0 wrote:

Hi, I have a Canu assembly (Pacbio reads). I want to use arrow to polish my assembly but beforehand I want to subset my raw reads to have a coverage of around 80X. The output from pacbio was a single subreads.bam file with coverage of 150X. I read that polishing coverage above 100X does not really improve the quality, so I want to test different coverage values. My question is how to create a new BAM file from the original that has e.g 80X coverage and is compatible with pbalign/blasr. I found several posts saying that using samtools is not going to work as the BAM generated file does not contain all the info needed by blasr. Thanks, Myrsini

pacbio bam coverage pbalign arrow • 718 views
ADD COMMENTlink written 2.4 years ago by MyrNat0
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: 1789 users visited in the last hour