View Single Post
Old 03-25-2017, 06:20 AM   #1
senior molecular biologist
Location: Belgium

Join Date: Jun 2009
Posts: 31
Default variant calling from pacBio Sequel data

I am struggling with variant calling in 'Re-sequencing' mode using SMRT v4 and so far found that the VCF format generate by SMRTv4 is not valid.
Their VCF v3.3 is a non-documented weird version which does not validate and includes strange annotations (I&D for INS and DEL and not reporting the last conserved base in INS calls - among probably other weirdnesses).

The PB GFF output is not very helpful either as I need to compare variants with other public and platform results (I need VCF4)

* Does anyone map Sequel reads to a reference assembly using external software and leading to VCF 4 results ?

* Does anyone have success calling variants from PacBio re-sequencing:BAM alignments and corresponding reference fasta assembly using other tools?

If yes please let us know which.
Thanks in advance.


Last edited by splaisan; 03-25-2017 at 06:29 AM.
splaisan is offline   Reply With Quote