Question: Problems with breakdancer (sv caller) output.
gravatar for unique379
18 months ago by
unique37990 wrote:

Hi, I am trying to use results come from breakdancer and found multiple problems associated with the output vcf (converted using a python script by

  1. Start and End Coordinates are not always correct (some END has smaller than start)
  2. SVLEN is not correct as other tools does, meaning "END-STAR=SVTLEN"
  3. What is definition of SVTYPE=INS by breakdancer? I found start and end has big difference? Which should not be as per my understanding from different tool representation such as Delly. "START and END should have only "1" nucleotide difference. OR some other tool such as WHAMG represents "Insertions start and end at the same position.".

WHAMG : Insertions start and end at the same position

chr1 1442939 . T <ins> . . SVTYPE=INS;END=1442939;

Delly represents start and start+1=END


but BREAKDancer output has some different way of representation which I really don't understand.

chr1 24662379 . . . . PASS PROGRAM=breakdancer;SVTYPE=INS;SVLEN=-3686;SVEND=24721723

Original output from breakdancer

chr 1 24662379 59+59- chr 1 24721723 59+59- INS -3686 99 37 /projects/grs-lab/RupeshWork/ResultsMP/180920-110831_JGBF_LeJ_GT17-02246-MP_AGTCAA_S51_L008/mapped_bowtie2/JGBF_LeJ_GT17-02246-MP_AGTCAA_S51_L008_realigned.bam|37 -nan NA

Can anyone have some advice to overcome these problems? As while using SURVIVOR merge its not working due to conflict of END within the breakdancer vcf or should say native output format.

ADD COMMENTlink written 18 months ago by unique37990
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: 1520 users visited in the last hour