|
2
|
1 <tool id="t-vs-vnormal_junctions" name="Virtual Normal Correction SVs" version="1.8">
|
|
0
|
2 <description> Filter SVs based on presence in VN set </description>
|
|
|
3
|
|
|
4 <requirements>
|
|
|
5 <requirement type="package" version="1.7">cgatools</requirement>
|
|
|
6 </requirements>
|
|
|
7
|
|
|
8 <command interpreter="bash">
|
|
|
9 JunctionDiff-vs-background.sh
|
|
|
10 --variants $variants
|
|
|
11 --reference ${reference.fields.reference_crr_cgatools}
|
|
2
|
12 --VN_junctions ${reference.fields.VN_genomes_junctionfile_list}${VNset}
|
|
0
|
13 --cgatools_binary cgatools
|
|
|
14 --outputfile_filtered $output_filtered
|
|
|
15 --scoreThresholdA $scoreThresholdA
|
|
|
16 --scoreThresholdB $scoreThresholdB
|
|
|
17 --distance $distance
|
|
2
|
18 --minlength $minlength
|
|
0
|
19 </command>
|
|
|
20
|
|
|
21 <inputs>
|
|
2
|
22 <param name="variants" type="data" format="tabular" label="CG Junctions file"/>
|
|
|
23
|
|
0
|
24 <!--select build-->
|
|
|
25 <param name="reference" type="select" label="Select Build">
|
|
|
26 <options from_data_table="virtual_normal_correction" />
|
|
|
27 </param>
|
|
2
|
28
|
|
|
29 <param name="VNset" type="select" label="Select Virtual Normal set to use" help="1000Genomes set can only be used for hg19 samples, for hg18 54 genomes will be used.">
|
|
|
30 <option value="46_diversity.txt" > CG Diversity Panel (46 Genomes) </option>
|
|
|
31 <option value="433_1000g.txt" > CG 1000G project genomes (433 Genomes) (hg19 only) </option>
|
|
|
32 <option value="479_diversity_1000g.txt" > Diversity and 1000G (479 genomes) (hg19 only) </option>
|
|
|
33 <option value="10_tutorial.txt" > Small VN for tutorial (10 Genomes) </option>
|
|
0
|
34 </param>
|
|
2
|
35
|
|
|
36
|
|
0
|
37 <param name="scoreThresholdA" type="text" value="10" label="scoreThreshold" help="The minimum number of discordant mate pair alignments supporting the junction from input genome"/>
|
|
|
38 <param name="scoreThresholdB" type="text" value="10" label="scoreThreshold" help="The minimum number of discordant mate pair alignments supporting the junction from background genomes"/>
|
|
|
39 <param name="distance" type="text" value="200" label="Maximum distance between coordinates of potentially compatible junctions."/>
|
|
|
40 <param name="minlength" type="text" value="500" label="Minimum deletion junctions length to be included into the difference file."/>
|
|
|
41 <param name="report" type="select" label="Generate report file?">
|
|
|
42 <option value="N" selected="true"> No </option>
|
|
|
43 <option value="Y"> Yes </option>
|
|
|
44 </param>
|
|
|
45 <param name="fname" type="text" value="" label="Prefix for your output file" help="Optional"/>
|
|
|
46 </inputs>
|
|
|
47
|
|
|
48 <outputs>
|
|
|
49 <data format="tabular" name="output_filtered" label="${fname} Filtered junctions for ${tool.name} on ${on_string}"/>
|
|
2
|
50 <data format="tabular" name="output_report" from_work_dir= "output_reports.tsv" label="${fname} report for ${tool.name} on ${on_string}">
|
|
|
51 <filter> report == "Y" </filter>
|
|
|
52 </data>
|
|
0
|
53 </outputs>
|
|
|
54
|
|
|
55 <help>
|
|
|
56 **What it does**
|
|
|
57
|
|
|
58
|
|
|
59
|
|
|
60 **Input Files**
|
|
|
61 Complete Genomics Junctions file
|
|
|
62
|
|
|
63 **Output Files**
|
|
|
64 Junctions remaining after filtering
|
|
|
65
|
|
|
66
|
|
|
67 </help>
|
|
|
68
|
|
|
69 </tool>
|
|
|
70
|
|
|
71
|
