|
9
|
1 <tool id="breakdancer" name="breakdancer" version="1.4.4">
|
|
|
2 <requirements>
|
|
|
3 <requirement type="package" version="1.4.4">breakdancer</requirement>
|
|
|
4 </requirements>
|
|
|
5
|
|
|
6 <description>
|
|
|
7 BreakDancer, released under GPLv3, is a Cpp package that provides genome-wide detection of structural variants from next generation paired-end sequencing reads. It includes two complementary programs. BreakDancerMax predicts five types of structural variants: insertions, deletions, inversions, inter- and intra-chromosomal translocations from next-generation short paired-end sequencing reads using read pairs that are mapped with unexpected separation distances or orientation. BreakDancerMini focuses on detecting small indels (usually between 10bp and 100bp) using normally mapped read pairs. Please read our paper for detailed algorithmic description. http://www.nature.com/nmeth/journal/v6/n9/abs/nmeth.1363.html
|
|
|
8
|
|
|
9 The input is a set of map files produced by a front-end aligner such as MAQ, BWA, NovoAlign and Bfast
|
|
|
10
|
|
|
11 The output format
|
|
|
12 ----------------------
|
|
|
13 BreakDancer's output file consists of the following columns:
|
|
|
14
|
|
|
15 1. Chromosome 1
|
|
|
16 2. Position 1
|
|
|
17 3. Orientation 1
|
|
|
18 4. Chromosome 2
|
|
|
19 5. Position 2
|
|
|
20 6. Orientation 2
|
|
|
21 7. Type of a SV
|
|
|
22 8. Size of a SV
|
|
|
23 9. Confidence Score
|
|
|
24 10. Total number of supporting read pairs
|
|
|
25 11. Total number of supporting read pairs from each map file
|
|
|
26 12. Estimated allele frequency
|
|
|
27 13. Software version
|
|
|
28 14. The run parameters
|
|
|
29 </description>
|
|
|
30
|
|
|
31 <parallelism method="basic"></parallelism>
|
|
|
32
|
|
|
33 <command interpreter="python">
|
|
|
34 breakdancer.py
|
|
|
35 <!-- required -->
|
|
|
36 -i1 $inputBamFile
|
|
|
37 -o1 $outputRawFile
|
|
|
38 <!-- optional -->
|
|
|
39 -o2 $outputVcfFile
|
|
|
40 -o $chromosome
|
|
|
41 -s $minLength
|
|
|
42 -c $cutoff
|
|
|
43 -m $maxSvSize
|
|
|
44 -q $minMapQuality
|
|
|
45 -r $minReadDepth
|
|
|
46 -x $maxHaploidCov
|
|
|
47 -b $bufferSize
|
|
|
48 #if $onlyTrans:
|
|
|
49 -t
|
|
|
50 #end if
|
|
|
51 -d $prefix
|
|
|
52 #if $bedFormat:
|
|
|
53 -g
|
|
|
54 #end if
|
|
|
55 <!-- -l
|
|
|
56 -a
|
|
|
57 -h
|
|
|
58 -y -->
|
|
|
59 <!-- parser.add_argument('-l', dest='matePair', required=False, help='analyze Illumina long insert (mate-pair) library')
|
|
|
60 # parser.add_argument('-a', dest='sortByLibrary', action='store_true', help='print out copy number and support reads per library rather than per bam', default=False)
|
|
|
61 # parser.add_argument('-h', dest='', action='store_true', help='print out Allele Frequency column', default=False)
|
|
|
62 parser.add_argument('-y', dest='scoreFilter', type=int, required=False, help='output score filter', default=30) -->
|
|
|
63
|
|
|
64 </command>
|
|
|
65
|
|
|
66 <inputs>
|
|
|
67 <param name="inputBamFile" format="bam" type="data" label="input bam file" />
|
|
|
68 <!-- optional -->
|
|
|
69 <param name="chromosome" type="text" label="operate on a single chromosome" value="" />
|
|
|
70 <param name="chromosome" type="integer" label="minimum length of a region" value="7" />
|
|
|
71 <param name="cutoff" type="integer" label="cutoff in unit of standard deviation" value="3" />
|
|
|
72 <param name="maxSvSize" type="integer" label="maximum SV size" value="1000000000" />
|
|
|
73 <param name="minMapQuality" type="integer" label="minimum alternative mapping quality" value="35" />
|
|
|
74 <param name="minReadDepth" type="integer" label="minimum number of read pairs required to establish a connection" value="2" />
|
|
|
75 <param name="maxHaploidCov" type="integer" label="maximum threshold of haploid sequence coverage for regions to be ignored" value="1000" />
|
|
|
76 <param name="bufferSize" type="integer" label="buffer size for building connection" value="100" />
|
|
|
77 <param name="onlyTrans" type="boolean" label="only detect transchromosomal rearrangement" checked="False" />
|
|
|
78 <param name="bufferSize" type="text" label="prefix of fastq files that SV supporting reads will be saved by library" value="" />
|
|
|
79 <param name="bedFormat" type="boolean" label="dump SVs and supporting reads in BED format for GBrowse" checked="False" />
|
|
|
80 </inputs>
|
|
|
81
|
|
|
82 <outputs>
|
|
|
83 <data name="outputRawFile" format="tabular" />
|
|
|
84 <data name="outputVcfFile" format="vcf" />
|
|
|
85 </outputs>
|
|
|
86
|
|
|
87 <help>
|
|
|
88 help
|
|
|
89 </help>
|
|
|
90 </tool> |
