annotate package/delly.xml @ 1:5f241f2398c3 draft

Uploaded
author malbuquerque
date Tue, 20 Jan 2015 14:21:53 -0500
parents ace1798624b4
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
1 <tool id="delly" name="DELLY" version="0.6.1">
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
2
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
3 <description>structural variant discovery by integrated paired-end and split-read analysis</description>
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
4
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
5 <requirements>
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
6 <requirement type="package" version="0.6.1">delly</requirement>
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
7 <requirement type="set_environment">DELLY_DIR</requirement>
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
8 </requirements>
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
9
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
10
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
11 <command>
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
12
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
13 <!-- BAM and BAI linking, (1) link BAM to new BAM file & (2) link BAM metadata to new BAI file -->
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
14 #for $i, $s in enumerate( $repeatBam )
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
15 ln -s $s.sortedBam ./input$(i).bam;
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
16 ln -s $s.sortedBam.metadata.bam_index ./input$(i).bam.bai;
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
17 #end for
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
18
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
19 <!-- Sets args to a list of types selected -->
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
20 #if not isinstance( $type.value, list ):
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
21 #set $args = [ $type.value ]
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
22 #else:
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
23 #set $args = $type.value
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
24 #end if
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
25
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
26 <!-- Run Delly Jobs for each type selected -->
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
27 #for $option in $args
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
28 <!-- NEED TO FIX -->
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
29 <!-- Delly should be automatically installed into the galaxy instance and should not be running off
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
30 the computers specific tool set -->
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
31 \$DELLY_DIR/src/delly -t $option -o ./output.$(option).vcf -q $mapQual -s $madCutoff
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
32 #if $option == "DEL"
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
33 -m $minFlank
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
34 #end if
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
35 -u $genoQual -v $vcfgeno -g $genome
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
36
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
37 <!-- add each input bam to command -->
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
38 #for $i, $s in enumerate( $repeatBam )
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
39 ./input$(i).bam
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
40 #end for
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
41 ;
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
42 #end for
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
43
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
44 <!-- Combine VCF Files and Sort Lexographically -->
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
45 #set $option = $args[0]
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
46 grep ^\# output.$(option).vcf > $outfile;
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
47 grep ^\# -v output.$(option).vcf > variants.txt;
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
48
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
49 <!-- If we called more than a single variant type, concatenate all the other types variant output -->
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
50 #if isinstance( $type.value, list ):
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
51 #for $option in $args[1:]
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
52 grep ^\# -v output.$(option).vcf >> variants.txt;
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
53 #end for
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
54 #end if
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
55
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
56 <!-- Sort all variant output, assuming that it will sort lexographically by chromosome, then position, ID -->
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
57 <!-- In future, maybe develop a script to sort by bam header -->
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
58 sort -k1,1d -k2,2n -k3,3d variants.txt > sortedVariants.txt;
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
59
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
60 <!-- Filter Variants that have Passed Quality Checks -->
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
61 #if $filterCalls
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
62 awk '{if ($7 == "PASS") print $0;}' sortedVariants.txt >> $outfile;
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
63 #else
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
64 cat sortedVariants.txt >> $outfile;
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
65 #end if
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
66
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
67 </command>
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
68
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
69 <inputs>
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
70
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
71 <!-- General Options -->
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
72 <param name="type" type="select" multiple="true" display="checkboxes" label="Variant Types">
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
73 <option value="DEL" selected="true">Deletions</option>
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
74 <option value="DUP" selected="true">Tandem Duplications</option>
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
75 <option value="INV" selected="true">Inversions</option>
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
76 <option value="TRA" selected="true">Translocations</option>
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
77 </param>
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
78 <repeat name="repeatBam" title="Bam Alignment" min="1" default="1" >
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
79 <param format="bam" name="sortedBam" type="data" label="File" />
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
80 </repeat>
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
81 <param name="excludeFile" type="data" optional="true" label="Chromosomes to Exclude"/>
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
82 <param name="filterCalls" type="boolean" value="false" label="Filter Poor Variant Calls"/>
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
83
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
84 <!-- Paired End Options -->
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
85 <param name="mapQual" type="integer" value="0" min="0" max="255" label="PE - Minimum Mapping Quality" />
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
86 <param name="madCutoff" type="integer" value="9" min="0" max="255" label="PE - Insert Size Cutoff" />
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
87
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
88 <!-- SR Options -->
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
89 <param format="fasta" name="genome" type="data" optional="true" label="SR - Genome Fasta File" />
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
90 <param name="minFlank" type="integer" value="13" label="SR - Minimum Flanking Sequence" />
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
91
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
92 <!-- Genotyping Options -->
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
93 <param format="vcf" name="vcfgeno" type="data" optional="true" label="GT - Input VCF" />
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
94 <param name="genoQual" type="integer" value="20" min="0" max="255" label="GT - Minimum Mapping Quality" />
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
95
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
96 </inputs>
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
97
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
98 <outputs>
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
99 <data format="vcf" name="outfile" />
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
100 </outputs>
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
101
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
102 <help>
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
103
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
104 </help>
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
105
ace1798624b4 Uploaded
malbuquerque
parents:
diff changeset
106 </tool>