0
|
1 <tool id="vcfs2fasta" name="VCFs to fasta" version="1.0">
|
|
2 <description>Takes a set of VCF files and outputs a multi fasta file with only the variant positions.</description>
|
|
3 <requirements>
|
|
4 <requirement type="package" version="2.7.10">python</requirement>
|
|
5 <requirement type="package" version="0.6.8dev">pyvcf</requirement>
|
|
6 <requirement type="package" version="3.11">pyyaml</requirement>
|
|
7 <requirement type="package" version="2.0.2">bintrees</requirement>
|
|
8 <requirement type="package" version="1.66">biopython</requirement>
|
|
9 <requirement type="package" version="1.4">matplotlib</requirement>
|
|
10 <requirement type="package" version="0.23.1">cython</requirement>
|
|
11 <requirement type="package" version="1.0">phephenix</requirement>
|
|
12 </requirements>
|
|
13 <stdio>
|
|
14 <!-- Assume anything other than zero is an error -->
|
|
15 <exit_code range="1:" />
|
|
16 <exit_code range=":-1" />
|
|
17 </stdio>
|
|
18 <command interpreter="bash">
|
|
19 vcfs2fasta.sh
|
8
|
20 $output
|
|
21 #if str($mix_cond.mix) == "yes":
|
0
|
22 $mix_cond.mix_value
|
|
23 #else
|
|
24 NOTTHERE
|
|
25 #end if
|
8
|
26 #if str($coln_cond.coln) == "yes":
|
|
27 $coln_cond.column_ns
|
0
|
28 #else
|
|
29 NOTTHERE
|
|
30 #end if
|
8
|
31 #if str($colgap_cond.colgap) == "yes":
|
9
|
32 $colgap_cond.column_gaps
|
8
|
33 #else
|
|
34 NOTTHERE
|
|
35 #end if
|
|
36 #if str($samplen_cond.samplen) == "yes":
|
|
37 $samplen_cond.sample_ns
|
|
38 #else
|
|
39 NOTTHERE
|
|
40 #end if
|
|
41 #if str($samplegap_cond.samplegap) == "yes":
|
|
42 $samplegap_cond.sample_gaps
|
0
|
43 #else
|
|
44 NOTTHERE
|
|
45 #end if
|
8
|
46 #if str($reference_cond.reference) == "yes":
|
|
47 $reference_cond.ref_fa
|
|
48 #else
|
|
49 NOTTHERE
|
|
50 #end if
|
|
51 #if str($include_cond.include) == "yes":
|
|
52 $include_cond.in_bed
|
|
53 #else
|
0
|
54 NOTTHERE
|
|
55 #end if
|
8
|
56 #if str($exclude_cond.exclude) == "yes":
|
|
57 $exclude_cond.ex_bed
|
|
58 #else
|
0
|
59 NOTTHERE
|
|
60 #end if
|
18
|
61 #if $withstats:
|
|
62 $statsfile
|
|
63 #else
|
|
64 NOTTHERE
|
|
65 #end if
|
8
|
66 #for $i, $input_vcf in enumerate( $input_vcfs ):
|
0
|
67 "${input_vcf}"
|
|
68 #end for
|
|
69 </command>
|
|
70
|
|
71 <inputs>
|
8
|
72 <param name="input_vcfs" type="data" multiple="true" format="vcf" label="Input VCF file(s)" />
|
|
73
|
|
74 <conditional name="mix_cond">
|
0
|
75 <param name="mix" type="select" label="With Mixtures" help="Specify this option with a threshold to output mixtures above this threshold.">
|
|
76 <option value="yes">Specify</option>
|
|
77 <option value="no" selected="true">Do not specify</option>
|
|
78 </param>
|
|
79 <when value="yes">
|
8
|
80 <param name="mix_value" type="float" value="0.5" label="Mixture value" />
|
0
|
81 </when>
|
|
82 </conditional>
|
8
|
83
|
|
84 <conditional name="coln_cond">
|
|
85 <param name="coln" type="select" label="Column Ns" help="Keeps columns with fraction of Ns below specified threshold.">
|
|
86 <option value="yes">Specify</option>
|
|
87 <option value="no" selected="true">Do not specify</option>
|
|
88 </param>
|
|
89 <when value="yes">
|
|
90 <param name="column_ns" type="float" value="0.5" label="Column Ns value" />
|
|
91 </when>
|
|
92 </conditional>
|
|
93
|
|
94 <conditional name="colgap_cond">
|
|
95 <param name="colgap" type="select" label="Column Gaps" help="Keeps columns with fraction of gaps below specified threshold.">
|
0
|
96 <option value="yes">Specify</option>
|
|
97 <option value="no" selected="true">Do not specify</option>
|
|
98 </param>
|
|
99 <when value="yes">
|
8
|
100 <param name="column_gaps" type="float" value="0.5" label="Column gaps value" />
|
0
|
101 </when>
|
|
102 </conditional>
|
8
|
103
|
|
104 <conditional name="samplen_cond">
|
|
105 <param name="samplen" type="select" label="Sample Ns" help="Keeps samples with fraction of Ns below specified threshold.">
|
0
|
106 <option value="yes">Specify</option>
|
|
107 <option value="no" selected="true">Do not specify</option>
|
|
108 </param>
|
|
109 <when value="yes">
|
8
|
110 <param name="sample_ns" type="float" value="0.5" label="Sample Ns value" />
|
0
|
111 </when>
|
|
112 </conditional>
|
8
|
113
|
|
114 <conditional name="samplegap_cond">
|
|
115 <param name="samplegap" type="select" label="Sample gaps" help="Keeps samples with fraction of gaps below specified threshold.">
|
|
116 <option value="yes">Specify</option>
|
|
117 <option value="no" selected="true">Do not specify</option>
|
|
118 </param>
|
|
119 <when value="yes">
|
|
120 <param name="sample_gaps" type="float" value="0.5" label="Sample gaps value" />
|
|
121 </when>
|
|
122 </conditional>
|
|
123
|
|
124 <conditional name="reference_cond">
|
0
|
125 <param name="reference" type="select" label="Reference genome file" help="If path to reference specified, then whole genome will be outputted.">
|
|
126 <option value="yes">Specify</option>
|
|
127 <option value="no" selected="true">Do not specify</option>
|
|
128 </param>
|
|
129 <when value="yes">
|
8
|
130 <param name="ref_fa" type="data" format="fasta" label="Reference fasta file" help="Fasta format"/>
|
0
|
131 </when>
|
|
132 </conditional>
|
8
|
133
|
|
134 <conditional name="include_cond">
|
0
|
135 <param name="include" type="select" label="Include region" help="Specify regions to include in a bed file.">
|
|
136 <option value="yes">Specify</option>
|
|
137 <option value="no" selected="true">Do not specify</option>
|
|
138 </param>
|
|
139 <when value="yes">
|
8
|
140 <param name="in_bed" type="data" format="bed" label="Include regions bed file" help="bed format"/>
|
0
|
141 </when>
|
|
142 </conditional>
|
8
|
143
|
|
144 <conditional name="exclude_cond">
|
0
|
145 <param name="exclude" type="select" label="Exclude region" help="Specify regions to exclude in a bed file.">
|
|
146 <option value="yes">Specify</option>
|
|
147 <option value="no" selected="true">Do not specify</option>
|
|
148 </param>
|
|
149 <when value="yes">
|
8
|
150 <param name="ex_bed" type="data" format="bed" label="Exclude regions bed file" help="bed format"/>
|
0
|
151 </when>
|
|
152 </conditional>
|
8
|
153
|
18
|
154 <param name="withstats" type="boolean" truevalue="true" falsevalue="false" checked="False" label="Write stats to separate file?" />
|
|
155
|
0
|
156 </inputs>
|
|
157
|
|
158 <outputs>
|
|
159 <data format="fasta" name="output" label="${tool.name} on ${on_string}: FASTA file" />
|
|
160 </outputs>
|
|
161 <test>
|
|
162 <param name="input_vcfs" value="1_short.vcf" ftype="vcf" />
|
|
163 <param name="input_vcfs" value="2_short.vcf" ftype="vcf" />
|
|
164 <output name="output" file="testresult.fa" ftype="fasta" />
|
|
165 </test>
|
|
166 <help>
|
|
167
|
|
168 </help>
|
|
169 </tool>
|