BreakDancer Pipeline


  • Step 1. Create a configuration file using bam2cfg.pl

    e.g.,
    myDirectory/bam2cfg.pl -g -h tumor.bam normal.bam > BRC6.cfg

    bam2cfg now only has the perl version.

    Manually view the insert size and flag distribution results in BRC6.cfg to see if there are any data quality issue. Usually std/mean should be < 0.2 or 0.3 at most. The flag 32(x%), represents percent of chimeric insert, this number (x%) should usually be smaller than 3%.

    View png files for the insert size distribution. You should usually see a normal distribution, a bimodal distribution is undesirable and it is not recommended to continue BreakDancerMax step with this situation existing.

  • Step 2. Detect inter-chromosomal translocations

    e.g.,
    myDirectory/BreakDancerMax.pl -t -q 10 -f -d BRC6.ctx BRC6.cfg > BRC6.ctx
          for perl version;
    myDirectory/breakdancer_max -t -q 10 -d BRC6.ctx BRC6.cfg > BRC6.ctx
          for cpp version.

    The -d option dumps CTX supporting read pairs into fastq files (in this case BRC6.ctx) by library.

    This step normally takes 12 hours or so for three bam files, 8 hours or so for two bam files for cpp version, around three days for perl version.