Question: mpileup for PacBio reads
gravatar for Biomonika (Noolean)
3.3 years ago by
State College, PA, USA
Biomonika (Noolean)3.1k wrote:

I am interested in calling errors (especially indels) in my PacBio alignments. I want ALL errors to be reported, not only high frequency ones. Thus, the best for me would an alternative to mpileup for Illumina reads. I am considering two options here:

1) using PacBio-specific software, something along the lines of: --algorithm=quiver chr22_P6.cmp.h5 -r hg19.fa -o variants.gff 
-o consensus.fasta -o consensus.fastq

2) Tweaking mpileup by turning of all the probabilistic realignments, something along the lines of:

samtools mpileup sam/forward/chr21_P6.cmp.h5.bam.forward.bam -f reference.fa 
-l feature.bed -uv --no-BAQ --open-prob 15 -Q 0 -D -t INFO/DPR

I have both .cmp.h5 and bam files with alignments.

ADD COMMENTlink modified 2.7 years ago by tjduncan260 • written 3.3 years ago by Biomonika (Noolean)3.1k
gravatar for tjduncan
2.7 years ago by
Indianapolis, IN
tjduncan260 wrote:

I think that taking a look at the recent biorxiv paper may answer some questions... but significantly complicate things...

Unfortunately I don't know of a great method reliably call "true" errors in your PacBio alignments that are mapped back to the currently used reference genomes (NA12878, PlatGen, GIAB, HG19, HG38 ect). This is because the variants that were annotated as "true" in these reference genomes were generated using primarily short read assemblers and the consensus of many different short-read variant callers. Thus these reference genomes are biased towards the error profile of what Illumina short reads can characterize and will not play well with novel variant calling methods and sequencing technologies (PacBio and ONT) that have different error profiles than short read sequencing.

ADD COMMENTlink modified 2.7 years ago • written 2.7 years ago by tjduncan260
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: 1491 users visited in the last hour