Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
New PacBio system: Sequel SNPsaurus Pacific Biosciences 62 09-28-2017 04:33 AM
PacBio Sequel demultiplexing Ali May Pacific Biosciences 4 12-06-2016 09:43 PM
Sequel System NA12878 Structural Variant Call Set pacbio Pacific Biosciences 0 12-06-2016 03:38 PM
PacBio Sequel NA12878 Structural Variant Call Set amw Bioinformatics 0 10-26-2016 10:55 AM
Variant calling from Trios Data d_t_nguyen Bioinformatics 1 10-25-2013 01:33 PM

Thread Tools
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

pacbio, reference assembly, variant calling

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 12:59 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO