0
|
1 <tool id="vt_@BINARY@" name="VT @BINARY@" version="@VERSION@.0">
|
|
2 <description>normalizes variants in a VCF file</description>
|
|
3 <expand macro="requirements" />
|
|
4 <expand macro="version_command" />
|
|
5 <macros>
|
|
6 <import>vt_macros.xml</import>
|
|
7 <token name="@BINARY@">normalize</token>
|
|
8 </macros>
|
|
9 <command>
|
|
10 <![CDATA[
|
|
11 vt @BINARY@
|
3
|
12 #if str($output_format) == 'bcf':
|
2
|
13 -o normalised.bcf
|
|
14 #else:
|
|
15 -o normalised.vcf
|
|
16 #end if
|
0
|
17 ##-q do not print options and summary []
|
1
|
18 -w $window
|
0
|
19
|
|
20 #if str( $intervals_file ) != 'None':
|
|
21 -I "${intervals_file}"
|
|
22 #end if
|
|
23
|
|
24 ##-i intervals []
|
|
25
|
|
26 #if $reference_source.reference_source_selector != "history":
|
1
|
27 -r "${reference_source.reference_genome.fields.path}"
|
0
|
28 #else:
|
|
29 -r "${reference_source.reference_genome}"
|
|
30 #end if
|
|
31
|
3
|
32 "${ infile }"
|
2
|
33
|
3
|
34 &&
|
|
35
|
|
36 #if str($output_format) == 'bcf':
|
2
|
37 mv normalised.bcf "${ outfile }";
|
|
38 #else:
|
|
39 mv normalised.vcf "${ outfile }";
|
|
40 #end if
|
|
41
|
0
|
42 ]]>
|
|
43 </command>
|
|
44 <expand macro="stdio" />
|
|
45 <inputs>
|
|
46 <param name="infile" type="data" format="vcf" label="VCF file to be normalised" />
|
|
47
|
|
48
|
|
49 <conditional name="reference_source">
|
|
50 <param name="reference_source_selector" type="select" label="Choose the source for the reference list">
|
|
51 <option value="cached" selected="True">Locally cached</option>
|
|
52 <option value="history">History</option>
|
|
53 </param>
|
|
54 <when value="cached">
|
|
55 <param name="reference_genome" type="select" label="Using reference genome">
|
|
56 <options from_data_table="fasta_indexes">
|
|
57 <filter type="data_meta" ref="infile" key="dbkey" column="1" />
|
|
58 <validator type="no_options" message="No reference genome is available for the build associated with the selected input dataset" />
|
|
59 </options>
|
|
60 </param>
|
|
61 </when>
|
|
62 <when value="history">
|
|
63 <param name="reference_genome" type="data" format="fasta" label="Using reference file" />
|
|
64 </when>
|
|
65 </conditional>
|
|
66
|
|
67 <param name="intervals_file" type="data" format="bed" optional="True" label="File containing list of intervals" />
|
|
68 <param name="window" type="integer" value="10000" size="10" label="Window size for local sorting of variants"
|
|
69 help="(-w)">
|
|
70 <validator type="in_range" min="0"/>
|
|
71 </param>
|
|
72
|
2
|
73 <param name="output_format" type="select" label="Choose the output format" help="">
|
|
74 <option value="bcf">BCF</option>
|
|
75 <option value="vcf" selected="true">VCF</option>
|
|
76 </param>
|
|
77
|
0
|
78 </inputs>
|
|
79 <outputs>
|
2
|
80 <data name="outfile" format="vcf" label="${tool.name} on ${on_string}">
|
|
81 <change_format>
|
|
82 <when input="output_format" value="bcf" format="bcf" />
|
|
83 </change_format>
|
|
84 </data>
|
0
|
85 </outputs>
|
|
86 <tests>
|
|
87 <test>
|
|
88 </test>
|
|
89 </tests>
|
|
90 <help>
|
|
91 **What it does**
|
|
92
|
|
93 Normalize variants in a VCF file. Normalized variants may have their positions changed; in such cases,
|
|
94 the normalized variants are reordered and output in an ordered fashion. The local reordering takes place over a window of 10000 base pairs.
|
|
95
|
|
96 @CITATION@
|
|
97 </help>
|
|
98 <!--expand macro="citations"/-->
|
|
99 </tool>
|