Question: PacBio Variant Caller options
gravatar for Mariela Cortés-López
4 months ago by
Mariela Cortés-López10 wrote:


I'm trying to find the way to call variants for a PacBio library. All the reads come from a similar construct (~1300 bp) on which we introduced several random mutations, each mutation set is associated with a particular barcode, for this reason, I have demultiplexed my reads applying a custom script over the circular consensus sequence, so now I want to get the mutational set associated with the barcodes.

I have been searching for ways to call the variants for PacBio Reads, but I have run into problems with some of the tools. For example, I aligned my demultiplexed reads with blastr and then I saved the output in xml format, but when attempted to use quiver it returns "Parse error" that seems to be related with the xml. I tried also running quiver with a bam file but this also returned an invalid "PacBio BAM file" error, I think quiver needs the cmp.h5 file, but I don't have that file and I also don't know if that could be useful to the kind of information that I want to get.

I've tried MICADo but for that algorithm, I have a problem running the example data set that they provide. I also read that using GATK for long reads would not be ideal. Is there any other tools that work with xml/bam/sam/fastq files from PacBio which could be useful for my plasmid library?


variants pacbio blasr quiver • 194 views
ADD COMMENTlink written 4 months ago by Mariela Cortés-López10
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: 1886 users visited in the last hour