Mercurial > repos > artbio > gatk4

comparison macros_filtermutectcalls.xml @ 0:c51c08cc9fcc draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
"planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gatk4 commit 408454e8d10befcc76f38ab446091778537d4f31"
author artbio
date Wed, 29 Dec 2021 01:36:41 +0000
parents
children 646e6943bcd2
comparison
equal deleted inserted replaced
-1:000000000000 0:c51c08cc9fcc
1 <?xml version="1.0"?>
2 <macros>
3 <token name="@VERSION@">4.1.7.0</token>
4 <token name="@WRAPPER_VERSION@">@VERSION@+galaxy1</token>
5
6 <xml name="requirements">
7 <requirements>
8 <requirement type="package" version="@VERSION@">gatk4</requirement>
9 <requirement type="package" version="0.2.5">tabix</requirement>
10 <requirement type="package" version="1.9">samtools</requirement>
11 <yield />
12 </requirements>
13 </xml>
14
15 <!--Hacky way to determine GATK version, for display in tool info-->
16 <xml name="version_cmd">
17 <version_command>gatk SortSam --version 2>&amp;1 | grep Version | cut -d ':' -f 2</version_command>
18 </xml>
19
20 <!--Define sections that parameters could exist within.-->
21 <template name="set_sections">
22 #set global $sections = ['', 'optional.', 'advanced.', 'deprecated.']
23 </template>
24
25 <!--Reference genome handling-->
26 <!--One template each for the different reference genome parameter names.-->
27 <!--TODO: Can the reference parameters all be the same?-->
28 <xml name="ref_sel">
29 <conditional name="reference_source">
30 <param name="reference_source_selector" type="select" label="Choose the source for the reference list">
31 <option value="cached">Locally cached</option>
32 <option value="history">History</option>
33 <option value="no_ref" selected="true">Do not pass</option>
34 </param>
35 <when value="cached">
36 <param name="reference_sequence" type="select" label="Reference" help="Reference sequence file." >
37 <options from_data_table="all_fasta" >
38 <validator type="no_options" message="A built-in reference genome is not available for the build associated with the selected input file" />
39 </options>
40 </param>
41 </when>
42 <when value="history">
43 <param name="reference_sequence" type="data" format="fasta" label="Reference" help="Reference sequence file." />
44 </when>
45 <when value="no_ref" />
46 </conditional>
47 </xml>
48
49 <template name="ref_opts">
50 #for $sect in $sections
51 #if $varExists($sect + "reference_source.reference_source_selector")
52 #if $getVar($sect + "reference_source.reference_source_selector") != "no_ref"
53 #if $getVar($sect + "reference_source.reference_source_selector") != "history"
54 --reference $getVar($sect + "reference_source.reference_sequence.fields.path")
55 #else
56 --reference reference.fa
57 #end if
58 #end if
59 #end if
60 #end for
61 </template>
62
63 <template name="picard_ref_opts">
64 #for $sect in $sections
65 #if $varExists($sect + "reference_source.reference_source_selector")
66 #if $getVar($sect + "reference_source.reference_source_selector") != "no_ref"
67 #if $getVar($sect + "reference_source.reference_source_selector") != "history"
68 --REFERENCE_SEQUENCE $getVar($sect + "reference_source.reference_sequence.fields.path")
69 #else
70 --REFERENCE_SEQUENCE $getVar($sect + "reference_source.reference_sequence")
71 #end if
72 #end if
73 #end if
74 #end for
75 </template>
76
77 <template name="picard_ref_opts_plain">
78 #for $sect in $sections
79 #if $varExists($sect + "reference_source.reference_source_selector")
80 #if $getVar($sect + "reference_source.reference_source_selector") != "no_ref"
81 #if $getVar($sect + "reference_source.reference_source_selector") != "history"
82 --REFERENCE $getVar($sect + "reference_source.reference_sequence.fields.path")
83 #else
84 --REFERENCE $getVar($sect + "reference_source.reference_sequence")
85 #end if
86 #end if
87 #end if
88 #end for
89 </template>
90
91 <template name="ref_opts_input">
92 #for $sect in $sections
93 #if $varExists($sect + "reference_source.reference_source_selector")
94 #if $getVar($sect + "reference_source.reference_source_selector") != "no_ref"
95 #if $getVar($sect + "reference_source.reference_source_selector") != "history"
96 --input $getVar($sect + "reference_source.reference_sequence.fields.path")
97 #else
98 --input $getVar($sect + "reference_source.reference_sequence")
99 #end if
100 #end if
101 #end if
102 #end for
103 </template>
104
105
106 <!--Interval Macros-->
107 <template name="gatk_ints_chth">
108 #for $sect in $sections
109 #if $varExists($sect + "ival_type.ival_type_sel")
110 #if $getVar($sect + "ival_type.ival_type_sel") == "ival_file"
111 #if $varExists($sect + "ival_type.intervals")
112 #if $getVar($sect + "ival_type.intervals").is_of_type("gatk_interval")
113 --intervals intervals.interval_list
114 #end if
115 #if $getVar($sect + "ival_type.intervals").is_of_type("bed")
116 --intervals intervals.bed
117 #end if
118 #if $getVar($sect + "ival_type.intervals").is_of_type("vcf")
119 --intervals intervals.vcf
120 #end if
121 #end if
122 #else
123 #if $varExists($sect + "ival_type.intervals")
124 --intervals $getVar($sect + "ival_type.intervals")
125 #end if
126 #end if
127 #if $varExists($sect + "ival_type.interval_padding")
128 --interval-padding $getVar($sect + "ival_type.interval_padding")
129 #end if
130 #end if
131 #end for
132 </template>
133
134
135 <template name="pre_gatk_ints_chth"><![CDATA[
136 #for $sect in $sections
137 #if $varExists($sect + "ival_type.ival_type_sel")
138 #if $getVar($sect + "ival_type.ival_type_sel") == "ival_file"
139 #if $varExists($sect + "ival_type.intervals")
140 #if $getVar($sect + "ival_type.intervals").is_of_type("gatk_interval")
141 ln -s $getVar($sect + "ival_type.intervals") intervals.interval_list &&
142 #end if
143 #if $getVar($sect + "ival_type.intervals").is_of_type("bed")
144 ln -s $getVar($sect + "ival_type.intervals") intervals.bed &&
145 #end if
146 #if $getVar($sect + "ival_type.intervals").is_of_type("vcf")
147 ln -s $getVar($sect + "ival_type.intervals") intervals.vcf &&
148 #end if
149 #end if
150 #end if
151 #end if
152 #end for
153 ]]></template>
154
155 <xml name="gatk_ints">
156 <conditional name="ival_type">
157 <param name="ival_type_sel" type="select" label="Choose Genomic Interval Source">
158 <option value="ival_file" selected="true">Interval File</option>
159 <option value="ival_text" selected="false">Interval Text Input</option>
160 </param>
161 <when value="ival_file">
162 <param name="intervals" argument="--intervals" type="data" optional="true" format="bed,vcf,gatk_interval" label="Intervals File" help="One or more genomic intervals over which to operate"/>
163 <param name="interval_padding" argument="--interval-padding" type="integer" optional="true" value="0" label="Interval Padding" help="Amount of padding (in bp) to add to each interval you are including."/>
164 </when>
165 <when value="ival_text">
166 <param name="intervals" argument="--intervals" type="text" optional="true" label="Intervals Text Input" help="One or more genomic intervals over which to operate. Enter in chrom:start-stop format."/>
167 <param name="interval_padding" argument="--interval-padding" type="integer" optional="true" value="0" label="Interval Padding" help="Amount of padding (in bp) to add to each interval you are including."/>
168 </when>
169 </conditional>
170 </xml>
171
172
173 <!--Exclude Intervals-->
174 <xml name="gatk_excl_ints">
175 <conditional name="excl_ival_type">
176 <param name="excl_ival_type_sel" type="select" label="Choose Genomic Interval Exclusion Source">
177 <option value="excl_ival_file" selected="true">Exclude Interval File</option>
178 <option value="excl_ival_text" selected="false">Exclude Interval Text Input</option>
179 </param>
180 <when value="excl_ival_file">
181 <param name="exclude_intervals" argument="--exclude-intervals" type="data" optional="true" format="bed,vcf,gatk_interval" label="Exclude Intervals File" help="One or more genomic intervals to exclude from processing"/>
182 <param name="interval_exclusion_padding" argument="--interval-exclusion-padding" type="integer" optional="true" value="0" label="Interval Exclusion Padding" help="Amount of padding (in bp) to add to each interval you are excluding."/>
183 </when>
184 <when value="excl_ival_text">
185 <param name="exclude_intervals" argument="--exclude-intervals" type="text" optional="true" label="Exclude Intervals Text Input" help="One or more genomic intervals to exclude from processing. Enter in chrom:start-stop format."/>
186 <param name="interval_exclusion_padding" argument="--interval-exclusion-padding" type="integer" optional="true" value="0" label="Interval Exclusion Padding" help="Amount of padding (in bp) to add to each interval you are excluding."/>
187 </when>
188 </conditional>
189 </xml>
190
191 <template name="gatk_excl_ints_chth">
192 #for $sect in $sections
193 #if $varExists($sect + "excl_ival_type.excl_ival_type_sel")
194 #if $getVar($sect + "excl_ival_type.excl_ival_type_sel") == "ival_file"
195 #if $varExists($sect + "excl_ival_type.exclude_intervals")
196 #if $getVar($sect + "excl_ival_type.exclude_intervals").is_of_type("gatk_interval")
197 --exclude-intervals excl_intervals.interval_list
198 #end if
199 #if $getVar($sect + "excl_ival_type.exclude_intervals").is_of_type("bed")
200 --exclude-intervals excl_intervals.bed
201 #end if
202 #if $getVar($sect + "excl_ival_type.exclude_intervals").is_of_type("vcf")
203 --exclude-intervals excl_intervals.vcf
204 #end if
205 #end if
206 #elif $getVar($sect + "excl_ival_type.excl_ival_type_sel") == "excl_ival_text"
207 #if $varExists($sect + "excl_ival_type.exclude_intervals")
208 --exclude-intervals $getVar($sect + "excl_ival_type.exclude_intervals")
209 #end if
210 #else
211 #pass
212 #end if
213 #if $varExists($sect + "excl_ival_type.interval_exclusion_padding")
214 --interval-exclusion-padding $getVar($sect + "excl_ival_type.interval_exclusion_padding")
215 #end if
216 #end if
217 #end for
218 </template>
219
220 <template name="pre_gatk_excl_ints_chth"><![CDATA[
221 #for $sect in $sections
222 #if $varExists($sect + "excl_ival_type.excl_ival_type_sel")
223 #if $getVar($sect + "excl_ival_type.excl_ival_type_sel") == "excl_ival_file"
224 #if $varExists($sect + "excl_ival_type.exclude_intervals")
225 #if $getVar($sect + "excl_ival_type.exclude_intervals").is_of_type("gatk_interval")
226 ln -s $getVar($sect + "excl_ival_type.exclude_intervals") excl_intervals.interval_list &&
227 #end if
228 #if $getVar($sect + "excl_ival_type.exclude_intervals").is_of_type("bed")
229 ln -s $getVar($sect + "excl_ival_type.exclude_intervals") excl_intervals.bed &&
230 #end if
231 #if $getVar($sect + "excl_ival_type.exclude_intervals").is_of_type("vcf")
232 ln -s $getVar($sect + "excl_ival_type.exclude_intervals") excl_intervals.vcf &&
233 #end if
234 #end if
235 #end if
236 #end if
237 #end for
238 ]]></template>
239
240 <!--Sequence dictionary selection options for Picard type tools-->
241 <template name="picard_seqdict_opts">
242 #for $sect in $sections
243 #if $varExists($sect + "seqdict_source.seqdict_source_selector")
244 #if $getVar($sect + "seqdict_source.seqdict_source_selector") != "no_seq_dict"
245 #if $getVar($sect + "seqdict_source.seqdict_source_selector") != "history"
246 #set seq_dict_loc = ''.join($getVar($sect + seqdict_source.seqdict_sequence).fields.path.split('.')[:-1]) + '.dict'
247 --SEQUENCE_DICTIONARY $seq_dict_loc
248 #else
249 --SEQUENCE_DICTIONARY $getVar($sect + "seqdict_source.seqdict_sequence")
250 #end if
251 #end if
252 #end if
253 #end for
254 </template>
255
256 <template name="gatk_seqdict">
257 #for $sect in $sections
258 #if $varExists($sect + "seqdict_source.seqdict_source_selector")
259 #if $getVar($sect + "seqdict_source.seqdict_source_selector") != "no_seq_dict"
260 #if $getVar($sect + "seqdict_source.seqdict_source_selector") != "history"
261 #set $seq_dict_loc = ''.join($getVar($sect + "seqdict_source.seqdict_sequence").fields.path.split('.')[:-1]) + '.dict'
262 --sequence-dictionary $seq_dict_loc
263 #else
264 --sequence-dictionary $getVar($sect + "seqdict_source.seqdict_sequence")
265 #end if
266 #end if
267 #end if
268 #end for
269 </template>
270
271 <xml name="seq_dict_sel">
272 <conditional name="seqdict_source">
273 <param name="seqdict_source_selector" type="select" label="Choose the source for the sequence dictionary">
274 <option value="cached">Locally cached</option>
275 <option value="history">History</option>
276 <option value="no_seq_dict" selected="true">Do not pass</option>
277 </param>
278 <when value="cached">
279 <param name="seqdict_sequence" type="select" label="Sequence Dictionary" help="Sequence dictionary file." >
280 <options from_data_table="all_fasta" >
281 <validator type="no_options" message="A built-in reference genome is not available for the build associated with the selected input file" />
282 </options>
283 </param>
284 </when>
285 <when value="history">
286 <param name="seqdict_sequence" type="data" format="txt" label="Sequence Dictionary" help="Sequence dictionary file. Must be in dict format." />
287 </when>
288 <when value="no_seq_dict" />
289 </conditional>
290 </xml>
291
292 <!--Output goes to stdout, no output parameter exists.-->
293 <template name="stdout_to_output">
294 > output.txt
295 </template>
296
297 <xml name="stdout_to_output_params">
298 <data format="txt" name="output" label="${tool.name} on ${on_string}: txt" from_work_dir="output.txt" />
299 </xml>
300
301 <!--Multiple input files, true for List[File] types, and sometimes List[String] types -->
302 <template name="gatk_input_multi"><![CDATA[
303 #for $num, $file in enumerate($variant)
304 #if $file.is_of_type("vcf_bgzip")
305 --variant input${num}.vcf.gz
306 #elif $file.is_of_type("txt")
307 --variant input${num}.list
308 #else
309 --variant input${num}.vcf
310 #end if
311 #end for
312 ]]></template>
313
314 <!--Multiple input files, true for List[File] types, and sometimes List[String] types -->
315 <template name="gatk_input_single"><![CDATA[
316 #if $variant.is_of_type("vcf_bgzip")
317 --variant input${num}.vcf.gz
318 #elif $variant.is_of_type("txt")
319 --variant input${num}.list
320 #else
321 --variant input${num}.vcf
322 #end if
323 ]]></template>
324
325 <template name="gatk_tabix_multi"><![CDATA[
326 #for $num, $file in enumerate($variant)
327 #set datatype = $file.datatype
328 #if $file.is_of_type("vcf_bgzip")
329 ln -s $file input${num}.vcf.gz &&
330 tabix input${num}.vcf.gz &&
331 #elif $file.is_of_type("txt")
332 ln -s $file input${num}.list &&
333 #else
334 ln -s $file input${num}.vcf &&
335 #end if
336 #end for
337 ]]></template>
338
339 <xml name="vcf_input_params_multi">
340 <param name="input" type="data" multiple="true" format="vcf,vcf_bgzip" label="Input VCF file" help="Input VCF(s) to be sorted. Multiple inputs must have the same sample names (in order)"/>
341 </xml>
342
343
344 <!--Macros for multiple input tools. List[File] in GATK json.-->
345 <template name="vcf_tabix_multi"><![CDATA[
346 #for $num, $file in enumerate($input)
347 #set datatype = $file.datatype
348 #if $file.is_of_type("vcf_bgzip")
349 ln -s $file input${num}.vcf.gz &&
350 tabix input${num}.vcf.gz &&
351 #else
352 ln -s $file input${num}.vcf &&
353 #end if
354 #end for
355 ]]></template>
356
357 <template name="vcf_input_multi_picard"><![CDATA[
358 #for $num, $file in enumerate($input)
359 #if $file.is_of_type("vcf_bgzip")
360 --INPUT input${num}.vcf.gz
361 #else
362 --INPUT input${num}.vcf
363 #end if
364 #end for
365 ]]></template>
366
367 <template name="vcf_input_multi"><![CDATA[
368 #for $num, $file in enumerate($input)
369 #if $file.is_of_type("vcf_bgzip")
370 --input input${num}.vcf.gz
371 #else
372 --input input${num}.vcf
373 #end if
374 #end for
375 ]]></template>
376
377
378 <!--Picard single input tools-->
379 <template name="vcf_tabix"><![CDATA[
380 #set datatype = $input.datatype
381 #if $input.is_of_type("vcf_bgzip")
382 ln -s $input input.vcf.gz &&
383 tabix input.vcf.gz &&
384 #else
385 ln -s $input input.vcf &&
386 #end if
387 ]]></template>
388
389 <template name="gatk_tabix"><![CDATA[
390 #set datatype = $variant.datatype
391 #if $variant.is_of_type("vcf_bgzip")
392 ln -s $variant input.vcf.gz &&
393 tabix input.vcf.gz &&
394 #else
395 ln -s $variant input.vcf &&
396 gatk IndexFeatureFile -F input.vcf &&
397 #end if
398 ]]></template>
399
400 <template name="vcf_input_picard"><![CDATA[
401 #if $input.is_of_type("vcf_bgzip")
402 --INPUT input.vcf.gz
403 #else
404 --INPUT input.vcf
405 #end if
406 ]]></template>
407
408 <template name="vcf_input"><![CDATA[
409 #if $input.is_of_type("vcf_bgzip")
410 --input input.vcf.gz
411 #else
412 --input input.vcf
413 #end if
414 ]]></template>
415
416 <template name="gatk_input"><![CDATA[
417 #if $variant.is_of_type("vcf_bgzip")
418 --variant input.vcf.gz
419 #else
420 --variant input.vcf
421 #end if
422 ]]></template>
423
424 <template name="gatk_gvcf_tabix"><![CDATA[
425 #if $variant
426 ln -s $variant input.g.vcf &&
427 #end if
428 ]]></template>
429
430 <template name="gatk_gvcf_input"><![CDATA[
431 --variant input.g.vcf
432 ]]></template>
433
434 <xml name="gatk_gvcf_input_params">
435 <param name="variant" type="data" multiple="false" format="vcf" label="Input gVCF file" help=""/>
436 </xml>
437
438 <xml name="vcf_input_params">
439 <param name="input" type="data" multiple="false" format="vcf,vcf_bgzip" label="Input VCF file" help="Input VCF(s) to be sorted. Multiple inputs must have the same sample names (in order)"/>
440 </xml>
441
442 <xml name="gatk_vcf_input_params">
443 <param name="variant" type="data" multiple="false" format="vcf,vcf_bgzip" label="Input VCF file" help="A VCF file containing variants."/>
444 </xml>
445
446 <xml name="gatk_vcf_input_params_multi">
447 <param name="variant" type="data" multiple="true" format="vcf,vcf_bgzip,txt" label="Input VCF file(s)" help="A VCF file containing variants or a list of VCFs. Can be specified multiple times."/>
448 </xml>
449
450 <xml name="gatk_req_params">
451 <param name="input" type="data" format="bam,sam,cram" label="Input BAM/SAM/CRAM file" />
452 </xml>
453
454 <!--HDF5 Inputs-->
455
456 <xml name="hdf5_input">
457 <param name="input" type="data" format="h5,tabular" label="Input TSV or HDF5" help="Input TSV or HDF5 file containing integer read counts in genomic intervals for a single case sample (output of CollectReadCounts)." />
458 </xml>
459
460 <template name="hdf5_input_chth"><![CDATA[
461 --input "${input}"
462 ]]></template>
463
464 <template name="hdf5_output_chth">
465 --output "${output}"
466 </template>
467
468 <xml name="hdf5_output">
469 <data format="h5" name="output" label="${tool.name} on ${on_string}: HDF5" help="Output file for read counts." />
470 </xml>
471
472 <!--Output specific to ModelSegments. Files created based on prefix, so force that to be what we want, then pull important files with from_work_dir.-->
473 <!--${SAMPLE}.cr.seg-->
474 <!--${SAMPLE}.modelFinal.seg-->
475 <template name="modelsegments_chth"><![CDATA[
476 --output "."
477 --output-prefix "modelsegments"
478 ]]></template>
479
480 <xml name="modelsegments_output">
481 <data format="tabular" name="cr_seg" label="${tool.name} on ${on_string}: cr.seg" help="Copy-ratio segments." from_work_dir="modelsegments.cr.seg"/>
482 <data format="tabular" name="modelfinal_seg" label="${tool.name} on ${on_string}: modelFinal.seg" help="Modeled Segments" from_work_dir="modelsegments.modelFinal.seg"/>
483 </xml>
484
485 <!--deltaMAD.txt-->
486 <!--denoisedLimit4.png-->
487 <!--denoisedMAD.txt-->
488 <!--denoised.png-->
489 <!--scaledDeltaMAD.txt-->
490 <!--modeled.png-->
491 <!--standardizedMAD.txt-->
492
493 <template name="plotmodeledsegments_chth"><![CDATA[
494 --output "."
495 --output-prefix "plotmodeledsegments"
496 ]]></template>
497
498 <xml name="plotmodeledsegments_output">
499 <data format="png" name="modeled_png" label="${tool.name} on ${on_string}: modeled.png" help="Copy-Ratio Plot" from_work_dir="plotmodeledsegments.modeled.png"/>
500 </xml>
501
502 <!--Common Picard options-->
503 <template name="picard_opts">
504 #if $picard_adv.arguments_file
505 --arguments_file ${picard_adv.arguments_file}
506 #end if
507 --COMPRESSION_LEVEL ${picard_adv.COMPRESSION_LEVEL}
508 #if $picard_adv.GA4GH_CLIENT_SECRETS
509 --GA4GH_CLIENT_SECRETS ${picard_adv.GA4GH_CLIENT_SECRETS}
510 #end if
511 --MAX_RECORDS_IN_RAM ${picard_adv.MAX_RECORDS_IN_RAM}
512 --VALIDATION_STRINGENCY ${picard_adv.VALIDATION_STRINGENCY}
513 --VERBOSITY ${picard_adv.VERBOSITY}
514 ${picard_adv.CREATE_MD5_FILE}
515 ${picard_adv.USE_JDK_DEFLATER}
516 ${picard_adv.USE_JDK_INFLATER}
517 </template>
518
519 <xml name="picard_params">
520 <section name="picard_adv" title="Advanced Picard Options (Only change these if you know what you're doing.)" expanded="False">
521 <param argument="--arguments_file" type="data" optional="true" format="txt" label="Arguments File" help="read one or more arguments files and add them to the command line" />
522 <param argument="--COMPRESSION_LEVEL" type="integer" optional="true" value="5" min="1" max="9" label="Compression Level" help="Compression level for all compressed files created (e.g. BAM and VCF)." />
523 <param argument="--CREATE_MD5_FILE" truevalue="--CREATE_MD5_FILE" falsevalue="" type="boolean" optional="true" checked="false" label="Create MD5 File" help="Whether to create an MD5 digest for any BAM or FASTQ files created." />
524 <param argument="--GA4GH_CLIENT_SECRETS" type="data" format="json" optional="true" label="Ga4Gh Client Secrets" help="Google Genomics API client_secrets.json file path." />
525 <param argument="--MAX_RECORDS_IN_RAM" type="integer" optional="true" value="500000" label="Max Records In Ram" help="When writing files that need to be sorted, this will specify the number of records stored in RAM before spilling to disk. Increasing this number reduces the number of file handles needed to sort the file, and increases the amount of RAM needed." />
526 <param argument="--TMP_DIR" type="text" optional="true" label="Tmp Dir" help="One or more directories with space available to be used by this program for temporary storage of working files. Keep in mind, you must be able to access this directory from either your user, or from the Galaxy user, depending on your configuration." />
527 <param argument="--USE_JDK_DEFLATER" truevalue="--USE_JDK_DEFLATER" falsevalue="" type="boolean" optional="true" checked="false" label="Use Jdk Deflater" help="Use the JDK Deflater instead of the Intel Deflater for writing compressed output" />
528 <param argument="--USE_JDK_INFLATER" truevalue="--USE_JDK_INFLATER" falsevalue="" type="boolean" optional="true" checked="false" label="Use Jdk Inflater" help="Use the JDK Inflater instead of the Intel Inflater for reading compressed input" />
529 <param argument="--VALIDATION_STRINGENCY" type="select" optional="true" label="Validation Stringency" help="Validation stringency for all SAM files read by this program. Setting stringency to SILENT can improve performance when processing a BAM file in which variable-length data (read, qualities, tags) do not otherwise need to be decoded." >
530 <option value="STRICT" selected="true">STRICT</option>
531 <option value="LENIENT" selected="false">LENIENT</option>
532 <option value="SILENT" selected="false">SILENT</option>
533 </param>
534 <param argument="--VERBOSITY" type="select" optional="true" label="Verbosity" help="Control verbosity of logging." >
535 <option value="ERROR" selected="false">ERROR</option>
536 <option value="WARNING" selected="false">WARNING</option>
537 <option value="INFO" selected="true">INFO</option>
538 <option value="DEBUG" selected="false">DEBUG</option>
539 </param>
540 </section>
541 </xml>
542
543 <!--Provides option to create gzipped output for VCF files-->
544 <xml name="gzip_vcf_params">
545 <param name="gzipped_output" type="boolean" checked="true" label="GZIP Output?" help="If you would like gzipped output, check this box. In general, it would be preferable to do this, unless your downstream tool does not support handling of gzipped files." />
546 </xml>
547
548 <!--Output related Picard options-->
549 <xml name="gzip_vcf_output_params">
550 <data format="vcf" name="output_vcf" label="${tool.name} on ${on_string}: vcf" from_work_dir="output.vcf" >
551 <filter>not gzipped_output</filter>
552 </data>
553 <data format="vcf_bgzip" name="output_vcf_bgzip" label="${tool.name} on ${on_string}: vcf_bgzip" from_work_dir="output.vcf.gz" >
554 <filter>gzipped_output</filter>
555 </data>
556 </xml>
557
558
559 <!--These are the same, other than the capitalization of output, so maybe a better way to do this.-->
560 <template name="picard_vcf_output_opts">
561 #if $gzipped_output
562 --OUTPUT output.vcf.gz
563 #else
564 --OUTPUT output.vcf
565 #end if
566 </template>
567
568 <template name="vcf_output_opts">
569 #if $gzipped_output
570 --output output.vcf.gz
571 #else
572 --output output.vcf
573 #end if
574 </template>
575
576 <xml name="picard_output_params">
577 <data format="txt" name="output_md5" label="${tool.name} on ${on_string}: md5sum(txt)" from_work_dir="output.bam.md5" >
578 <filter>picard_adv['CREATE_MD5_FILE']</filter>
579 </data>
580 </xml>
581
582
583 <!--Citations-->
584 <xml name="citations">
585 <citation type="doi">10.1101/gr.107524.110</citation>
586 <citation type="doi">10.1038/ng.806</citation>
587 <citation type="doi">10.1002/0471250953.bi1110s43</citation>
588 <citation type="doi">10.1101/201178</citation>
589 <yield />
590 </xml>
591
592 </macros>