|
6
|
1 <tool id="gatk2_variant_eval" name="Eval Variants" version="@VERSION@.0">
|
|
0
|
2 <description></description>
|
|
|
3 <expand macro="requirements" />
|
|
6
|
4 <expand macro="version_command" />
|
|
0
|
5 <macros>
|
|
|
6 <import>gatk2_macros.xml</import>
|
|
|
7 </macros>
|
|
1
|
8 <command interpreter="python">
|
|
0
|
9 #from binascii import hexlify
|
|
1
|
10
|
|
|
11 gatk2_wrapper.py
|
|
0
|
12 --stdout "${output_log}"
|
|
1
|
13 #for $var_count, $variant in enumerate( $reference_source.variants ):
|
|
|
14 -d "--eval:input_${var_count},%(file_type)s" "${variant.input_variant}" "${variant.input_variant.ext}" "input_variants_${var_count}"
|
|
|
15 #end for
|
|
|
16 -p '
|
|
|
17 @JAR_PATH@
|
|
0
|
18 -T "VariantEval"
|
|
|
19 --out "${output_report}"
|
|
|
20 \$GATK2_SITE_OPTIONS
|
|
|
21
|
|
|
22 @THREADS@
|
|
|
23
|
|
|
24 #if $reference_source.reference_source_selector != "history":
|
|
|
25 -R "${reference_source.ref_file.fields.path}"
|
|
|
26 #end if
|
|
|
27 '
|
|
|
28
|
|
|
29 #for $rod_binding in $comp_rod_bind:
|
|
|
30 -d "--comp:${rod_binding.comp_rod_name},%(file_type)s" "${rod_binding.comp_input_rod}" "${rod_binding.comp_input_rod.ext}" "input_comp_${rod_binding.comp_rod_name}"
|
|
|
31 #if str( $rod_binding.comp_known_names ):
|
|
|
32 -p '--known_names "${rod_binding.comp_rod_name}"'
|
|
|
33 #end if
|
|
|
34 #end for
|
|
|
35
|
|
4
|
36 #if $dbsnp_rod_bind_type.dbsnp_rod_bind_type_selector == 'set_dbsnp'
|
|
0
|
37 -d "--dbsnp:${dbsnp_rod_bind_type.dbsnp_rod_name},%(file_type)s" "${dbsnp_rod_bind_type.dbsnp_input_rod}" "${dbsnp_rod_bind_type.dbsnp_input_rod.ext}" "input_dbsnp_${dbsnp_rod_bind_type.dbsnp_rod_name}"
|
|
4
|
38 #if $dbsnp_rod_bind_type.dbsnp_known_names
|
|
0
|
39 -p '--known_names "${dbsnp_rod_bind_type.dbsnp_rod_name}"'
|
|
|
40 #end if
|
|
|
41 #end if
|
|
|
42
|
|
|
43 #include source=$standard_gatk_options#
|
|
|
44
|
|
|
45
|
|
|
46 ##start analysis specific options
|
|
|
47 #if $analysis_param_type.analysis_param_type_selector == "advanced":
|
|
|
48 #for $stratification in $analysis_param_type.stratifications:
|
|
|
49 #set $select_string = "--select_exps '%s' --select_names '%s'" % ( str( $stratification.select_exps ), str( $stratification.select_name ) )
|
|
|
50 -o '${ hexlify( $select_string ) }'
|
|
|
51 #end for
|
|
|
52 -p '
|
|
|
53
|
|
|
54 #for $sample in $analysis_param_type.samples:
|
|
|
55 --sample "${sample.sample}"
|
|
|
56 #end for
|
|
|
57
|
|
|
58 #if str( $analysis_param_type.stratification_modules ) != "None":
|
|
|
59 #for $stratification_module in str( $analysis_param_type.stratification_modules).split( ',' ):
|
|
|
60 --stratificationModule "${stratification_module}"
|
|
|
61 #end for
|
|
|
62 #end if
|
|
|
63
|
|
|
64 ${analysis_param_type.do_not_use_all_standard_stratifications}
|
|
|
65
|
|
|
66 #for $variant_type in $analysis_param_type.only_variants_of_type:
|
|
|
67 --onlyVariantsOfType "${variant_type.variant_type}"
|
|
|
68 #end for
|
|
|
69
|
|
|
70 #if str( $analysis_param_type.eval_modules ) != "None":
|
|
|
71 #for $eval_module in str( $analysis_param_type.eval_modules).split( ',' ):
|
|
|
72 --evalModule "${eval_module}"
|
|
|
73 #end for
|
|
|
74 #end if
|
|
|
75
|
|
|
76 ${analysis_param_type.do_not_use_all_standard_modules}
|
|
|
77
|
|
|
78 #if str( $analysis_param_type.num_samples ) != "0":
|
|
|
79 --numSamples "${analysis_param_type.num_samples}"
|
|
|
80 #end if
|
|
|
81
|
|
|
82 --minPhaseQuality "${analysis_param_type.min_phase_quality}"
|
|
|
83
|
|
|
84 --mendelianViolationQualThreshold "${analysis_param_type.mendelian_violation_qual_threshold}"
|
|
|
85
|
|
|
86 #if str( $analysis_param_type.ancestral_alignments ) != "None":
|
|
|
87 --ancestralAlignments "${analysis_param_type.ancestral_alignments}"
|
|
|
88 #end if
|
|
|
89 '
|
|
|
90 #if str( $analysis_param_type.known_cnvs ) != "None":
|
|
|
91 -d "--knownCNVs" "${analysis_param_type.known_cnvs}" "${analysis_param_type.known_cnvs.ext}" "input_known_cnvs"
|
|
|
92 #end if
|
|
|
93
|
|
|
94 #if str( $analysis_param_type.strat_intervals ) != "None":
|
|
|
95 -d "--stratIntervals" "${analysis_param_type.strat_intervals}" "${analysis_param_type.strat_intervals.ext}" "input_strat_intervals"
|
|
|
96 #end if
|
|
|
97 #end if
|
|
|
98 </command>
|
|
|
99 <inputs>
|
|
|
100
|
|
|
101 <conditional name="reference_source">
|
|
|
102 <expand macro="reference_source_selector_param" />
|
|
|
103 <when value="cached">
|
|
|
104 <repeat name="variants" title="Variant" min="1" help="-eval,--eval &lt;eval&gt;">
|
|
|
105 <param name="input_variant" type="data" format="vcf" label="Input variant file" />
|
|
|
106 </repeat>
|
|
|
107 <param name="ref_file" type="select" label="Using reference genome" help="-R,--reference_sequence &lt;reference_sequence&gt;">
|
|
|
108 <options from_data_table="gatk2_picard_indexes">
|
|
|
109 <!-- <filter type="data_meta" key="dbkey" ref="input_variant" column="dbkey"/> -->
|
|
|
110 </options>
|
|
|
111 <validator type="no_options" message="A built-in reference genome is not available for the build associated with the selected input file"/>
|
|
|
112 </param>
|
|
|
113 </when>
|
|
|
114 <when value="history"> <!-- FIX ME!!!! -->
|
|
|
115 <repeat name="variants" title="Variant" min="1" help="-eval,--eval &lt;eval&gt;">
|
|
|
116 <param name="input_variant" type="data" format="vcf" label="Input variant file" />
|
|
|
117 </repeat>
|
|
|
118 <param name="ref_file" type="data" format="fasta" label="Using reference file" help="-R,--reference_sequence &lt;reference_sequence&gt;" />
|
|
|
119 </when>
|
|
|
120 </conditional>
|
|
|
121
|
|
4
|
122 <repeat name="comp_rod_bind" title="Comparison Reference-Ordered Data (ROD) file" help="-comp,--comp &lt;comp&gt;">
|
|
0
|
123 <param name="comp_input_rod" type="data" format="vcf" label="Comparison ROD file" />
|
|
4
|
124 <param name="comp_rod_name" type="text" value="" label="Comparison ROD name">
|
|
|
125 <validator type="regex" message="Value must be a not empty string composed by alphanumeric characters and underscores">^\w+$</validator>
|
|
|
126 </param>
|
|
|
127 <param name="comp_known_names" type="boolean" label="Use comparison ROD file as known_names" help="-knownName,--known_names &lt;known_names&gt;"/>
|
|
0
|
128 </repeat>
|
|
|
129 <conditional name="dbsnp_rod_bind_type">
|
|
4
|
130 <param name="dbsnp_rod_bind_type_selector" type="select" label="Provide a dbSNP Reference-Ordered Data (ROD) file" help="-D,--dbsnp &lt;dbsnp&gt;">
|
|
0
|
131 <option value="set_dbsnp" selected="True">Set dbSNP</option>
|
|
|
132 <option value="exclude_dbsnp">Don't set dbSNP</option>
|
|
|
133 </param>
|
|
4
|
134 <when value="exclude_dbsnp" />
|
|
0
|
135 <when value="set_dbsnp">
|
|
|
136 <param name="dbsnp_input_rod" type="data" format="vcf" label="dbSNP ROD file" />
|
|
4
|
137 <param name="dbsnp_rod_name" type="text" value="dbsnp" label="dbsnp ROD name">
|
|
|
138 <validator type="regex" message="Value must be a not empty string composed by alphanumeric characters and underscores">^\w+$</validator>
|
|
|
139 </param>
|
|
|
140 <param name="dbsnp_known_names" type="boolean" label="Use dbSNP ROD file as known_names" help="-knownName,--known_names &lt;known_names&gt;" />
|
|
0
|
141 </when>
|
|
|
142 </conditional>
|
|
|
143
|
|
|
144 <expand macro="gatk_param_type_conditional" />
|
|
|
145
|
|
|
146
|
|
|
147 <expand macro="analysis_type_conditional">
|
|
|
148 <repeat name="stratifications" title="Stratification">
|
|
|
149 <param name="select_exps" value="" type="text" label="Stratification Expression" help="-select,--select_exps &lt;select_exps&gt;">
|
|
|
150 <sanitizer>
|
|
|
151 <valid initial="string.printable">
|
|
|
152 <remove value="'"/>
|
|
|
153 </valid>
|
|
|
154 <mapping initial="none"/>
|
|
|
155 </sanitizer>
|
|
|
156 </param>
|
|
|
157 <param name="select_name" value="" type="text" label="Name" help="-selectName,--select_names &lt;select_names&gt;"/>
|
|
|
158 </repeat>
|
|
|
159
|
|
|
160 <repeat name="samples" title="Sample" help="-sn,--sample &lt;sample&gt;">
|
|
|
161 <param name="sample" value="" type="text" label="Derive eval and comp contexts using only these sample genotypes, when genotypes are available in the original context"/>
|
|
|
162 </repeat>
|
|
|
163
|
|
|
164 <param name="stratification_modules" type="select" multiple="True" display="checkboxes" label="Stratification modules to apply to the eval track(s)" help="-ST,--stratificationModule &lt;stratificationModule&gt;" >
|
|
6
|
165 <option value="AlleleCount" />
|
|
0
|
166 <option value="AlleleFrequency" />
|
|
|
167 <option value="CompRod" />
|
|
|
168 <option value="Contig" />
|
|
|
169 <option value="CpG" />
|
|
|
170 <option value="Degeneracy" />
|
|
|
171 <option value="EvalRod" />
|
|
|
172 <option value="Filter" />
|
|
|
173 <option value="FunctionalClass" />
|
|
6
|
174 <option value="IndelSize" />
|
|
|
175 <option value="IntervalStratification" />
|
|
0
|
176 <option value="JexlExpression" />
|
|
6
|
177 <option value="Novelty" />
|
|
|
178 <option value="OneBPIndel" />
|
|
0
|
179 <option value="Sample" />
|
|
6
|
180 <option value="SnpEffPositionModifier" />
|
|
|
181 <option value="TandemRepeat" />
|
|
|
182 <option value="VariantType" />
|
|
0
|
183 </param>
|
|
|
184 <param name="do_not_use_all_standard_stratifications" checked="false" type="boolean" truevalue="--doNotUseAllStandardStratifications" falsevalue="" label="Do not use the standard stratification modules by default" help="-noST,--doNotUseAllStandardStratifications" />
|
|
|
185
|
|
|
186 <repeat name="only_variants_of_type" title="only Variants Of Type" help="--onlyVariantsOfType">
|
|
|
187 <param name="variant_type" type="text" value="" label="only variants of these types will be considered during the evaluation"/>
|
|
|
188 </repeat>
|
|
|
189
|
|
|
190 <param name="eval_modules" type="select" multiple="True" display="checkboxes" label="Eval modules to apply to the eval track(s)" help="-EV,--evalModule &lt;evalModule&gt;" >
|
|
|
191 <option value="CompOverlap" />
|
|
|
192 <option value="CountVariants" />
|
|
6
|
193 <option value="IndelLengthHistogram" />
|
|
|
194 <option value="IndelSummary" />
|
|
0
|
195 <option value="MendelianViolationEvaluator" />
|
|
6
|
196 <option value="MultiallelicSummary" />
|
|
0
|
197 <option value="PrintMissingComp" />
|
|
|
198 <option value="ThetaVariantEvaluator" />
|
|
|
199 <option value="TiTvVariantEvaluator" />
|
|
6
|
200 <option value="ValidationReport" />
|
|
|
201 <option value="VariantSummary" />
|
|
0
|
202 </param>
|
|
|
203 <param name="do_not_use_all_standard_modules" checked="false" type="boolean" truevalue="--doNotUseAllStandardModules" falsevalue="" label="Do not use the standard eval modules by default" help="-noEV,--doNotUseAllStandardModules" />
|
|
|
204
|
|
|
205 <param name="num_samples" type="integer" label="Number of samples (used if no samples are available in the VCF file" value="0" help="-ns,--numSamples &lt;numSamples&gt;"/>
|
|
|
206 <param name="min_phase_quality" type="float" label="Minimum phasing quality " value="10.0" help="-mpq,--minPhaseQuality &lt;minPhaseQuality&gt;"/>
|
|
|
207 <param name="mendelian_violation_qual_threshold" type="integer" label="Minimum genotype QUAL score for each trio member required to accept a site as a violation" value="50" help="-mvq,--mendelianViolationQualThreshold &lt;mendelianViolationQualThreshold&gt;"/>
|
|
|
208 <param name="ancestral_alignments" type="data" format="fasta" optional="True" label="Fasta file with ancestral alleles" help="-aa,--ancestralAlignments &lt;ancestralAlignments&gt;" />
|
|
|
209 <param name="known_cnvs" type="data" format="bed,gatk_interval,picard_interval_list" optional="True" label="File containing tribble-readable features describing a known list of copy number variants" help="-knownCNVs,--knownCNVs &lt;knownCNVs&gt;" />
|
|
|
210 <param name="strat_intervals" type="data" format="bed,gatk_interval,picard_interval_list" optional="True" label="File containing tribble-readable features for the IntervalStratificiation" help="-stratIntervals,--stratIntervals &lt;stratIntervals&gt;" />
|
|
|
211
|
|
|
212 </expand>
|
|
|
213
|
|
|
214
|
|
|
215 </inputs>
|
|
|
216 <outputs>
|
|
|
217 <data format="gatk_report" name="output_report" label="${tool.name} on ${on_string} (report)" />
|
|
|
218 <data format="txt" name="output_log" label="${tool.name} on ${on_string} (log)" />
|
|
|
219 </outputs>
|
|
|
220 <tests>
|
|
|
221 <test>
|
|
|
222 <param name="reference_source_selector" value="history" />
|
|
|
223 <param name="ref_file" value="phiX.fasta" ftype="fasta" />
|
|
|
224 <param name="input_variant" value="gatk/gatk_variant_annotator/gatk_variant_annotator_out_1.vcf" ftype="vcf" />
|
|
|
225 <param name="dbsnp_rod_bind_type_selector" value="set_dbsnp" />
|
|
|
226 <param name="dbsnp_input_rod" value="gatk/fake_phiX_variant_locations.vcf" ftype="vcf" />
|
|
4
|
227 <param name="dbsnp_rod_name" value="dbsnp" />
|
|
0
|
228 <param name="dbsnp_known_names" value="True"/>
|
|
|
229 <param name="comp_rod_bind" value="0" />
|
|
|
230 <param name="gatk_param_type_selector" value="basic" />
|
|
|
231 <param name="analysis_param_type_selector" value="basic" />
|
|
|
232 <output name="output_report" file="gatk/gatk_variant_eval/gatk_variant_eval_out_1.gatk_report" />
|
|
|
233 <output name="output_log" file="gatk/gatk_variant_eval/gatk_variant_eval_out_1.log.contains" compare="contains" />
|
|
|
234 </test>
|
|
|
235 </tests>
|
|
|
236 <help>
|
|
|
237 **What it does**
|
|
|
238
|
|
|
239 General-purpose tool for variant evaluation (% in dbSNP, genotype concordance, Ti/Tv ratios, and a lot more)
|
|
|
240
|
|
|
241 For more information on using the VariantEval module, see this `tool specific page <http://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_sting_gatk_walkers_varianteval_VariantEval.html>`_.
|
|
|
242
|
|
|
243 To learn about best practices for variant detection using GATK, see this `overview <http://www.broadinstitute.org/gatk/guide/topic?name=best-practices>`_.
|
|
|
244
|
|
|
245 If you encounter errors, please view the `GATK FAQ <http://www.broadinstitute.org/gatk/guide/topic?name=faqs>`_.
|
|
|
246
|
|
|
247 ------
|
|
|
248
|
|
|
249 **Inputs**
|
|
|
250
|
|
|
251 GenomeAnalysisTK: VariantEval accepts variant files as input.
|
|
|
252
|
|
|
253
|
|
|
254 **Outputs**
|
|
|
255
|
|
|
256 The output is a table of variant evaluation.
|
|
|
257
|
|
|
258
|
|
|
259 Go `here <http://www.broadinstitute.org/gatk/guide/topic?name=intro>`_ for details on GATK file formats.
|
|
|
260
|
|
|
261 -------
|
|
|
262
|
|
|
263 **Settings**::
|
|
|
264
|
|
|
265 out An output file presented to the walker. Will overwrite contents if file exists.
|
|
|
266 list List the available eval modules and exit
|
|
|
267 select_exps One or more stratifications to use when evaluating the data
|
|
|
268 select_names Names to use for the list of stratifications (must be a 1-to-1 mapping)
|
|
|
269 sample Derive eval and comp contexts using only these sample genotypes, when genotypes are available in the original context
|
|
|
270 known_names Name of ROD bindings containing variant sites that should be treated as known when splitting eval rods into known and novel subsets
|
|
|
271 stratificationModule One or more specific stratification modules to apply to the eval track(s) (in addition to the standard stratifications, unless -noS is specified)
|
|
|
272 doNotUseAllStandardStratifications Do not use the standard stratification modules by default (instead, only those that are specified with the -S option)
|
|
|
273 onlyVariantsOfType If provided, only variants of these types will be considered during the evaluation, in
|
|
|
274 evalModule One or more specific eval modules to apply to the eval track(s) (in addition to the standard modules, unless -noE is specified)
|
|
|
275 doNotUseAllStandardModules Do not use the standard modules by default (instead, only those that are specified with the -E option)
|
|
|
276 numSamples Number of samples (used if no samples are available in the VCF file
|
|
|
277 minPhaseQuality Minimum phasing quality
|
|
|
278 mendelianViolationQualThreshold Minimum genotype QUAL score for each trio member required to accept a site as a violation
|
|
|
279 ancestralAlignments Fasta file with ancestral alleles
|
|
|
280
|
|
|
281 @CITATION_SECTION@
|
|
|
282 </help>
|
|
6
|
283 <expand macro="citations" />
|
|
0
|
284 </tool>
|
