comparison snpEff.xml @ 23:924af057bbca draft default tip

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tool_collections/snpeff commit 93465571cf180554c3548338a68fa0f1604985dc-dirty
author jjohnson
date Sat, 09 Apr 2016 11:51:40 -0400
parents bb0797deab78
children
comparison
equal deleted inserted replaced
22:bb0797deab78 23:924af057bbca
6 <expand macro="requirements" /> 6 <expand macro="requirements" />
7 <expand macro="stdio" /> 7 <expand macro="stdio" />
8 <expand macro="version_command" /> 8 <expand macro="version_command" />
9 <command> 9 <command>
10 <![CDATA[ 10 <![CDATA[
11 java -Xmx6G -jar \$SNPEFF_JAR_PATH/snpEff.jar eff 11 java -Xmx6G -jar "\$SNPEFF_JAR_PATH/snpEff.jar" eff
12 -c \$SNPEFF_JAR_PATH/snpEff.config 12 -c "\$SNPEFF_JAR_PATH/snpEff.config"
13 -i $inputFormat -o ${outputConditional.outputFormat} -upDownStreamLen $udLength 13 -i $inputFormat -o ${outputConditional.outputFormat} -upDownStreamLen $udLength
14 #if $spliceSiteSize and $spliceSiteSize.__str__ != '': 14 #if $spliceSiteSize and str($spliceSiteSize) != '':
15 -spliceSiteSize $spliceSiteSize 15 -spliceSiteSize "$spliceSiteSize"
16 #end if 16 #end if
17 #if $spliceRegion.setSpliceRegions == 'yes': 17 #if $spliceRegion.setSpliceRegions == 'yes':
18 #if $spliceRegion.spliceRegionExonSize and $spliceRegion.spliceRegionExonSize.__str__ != '': 18 #if $spliceRegion.spliceRegionExonSize and str($spliceRegion.spliceRegionExonSize) != '':
19 -spliceRegionExonSize $spliceRegion.spliceRegionExonSize 19 -spliceRegionExonSize $spliceRegion.spliceRegionExonSize
20 #end if 20 #end if
21 #if $spliceRegion.spliceRegionIntronMin and $spliceRegion.spliceRegionIntronMin.__str__ != '': 21 #if $spliceRegion.spliceRegionIntronMin and str($spliceRegion.spliceRegionIntronMin) != '':
22 -spliceRegionIntronMin $spliceRegion.spliceRegionIntronMin 22 -spliceRegionIntronMin $spliceRegion.spliceRegionIntronMin
23 #end if 23 #end if
24 #if $spliceRegion.spliceRegionIntronMax and $spliceRegion.spliceRegionIntronMax.__str__ != '': 24 #if $spliceRegion.spliceRegionIntronMax and str($spliceRegion.spliceRegionIntronMax) != '':
25 -spliceRegionIntronMax $spliceRegion.spliceRegionIntronMax 25 -spliceRegionIntronMax $spliceRegion.spliceRegionIntronMax
26 #end if 26 #end if
27 #end if 27 #end if
28 #if $filterHomHet and $filterHomHet.__str__ != 'no_filter': 28 #if $annotations and str($annotations) != '':
29 $filterHomHet
30 #end if
31 #if $annotations and $annotations.__str__ != '':
32 #echo " " 29 #echo " "
33 #echo ' '.join($annotations.__str__.split(',')) 30 #echo ' '.join(str($annotations).split(','))
34 #end if 31 #end if
35 #if $filterOut and $filterOut.__str__ != '': 32 #if $filterOut and str($filterOut) != '':
36 #echo " " 33 #echo " "
37 #echo ' '.join($filterOut.__str__.split(',')) 34 #echo ' '.join(str($filterOut).split(','))
38 #end if 35 #end if
39 #if $filter.specificEffects == 'yes' and $filter.effects: 36 #if $filter.specificEffects == 'yes' and $filter.effects:
40 #for $eff in str($filter.effects).split(','): 37 #for $eff in str($filter.effects).split(','):
41 -no $eff 38 -no $eff
42 #end for 39 #end for
48 -interval $intervals 45 -interval $intervals
49 #end if 46 #end if
50 #if $statsFile: 47 #if $statsFile:
51 -stats $statsFile 48 -stats $statsFile
52 #end if 49 #end if
53 #if $offset.__str__ != 'default': 50 #if str($offset) != 'default':
54 ${offset} 51 ${offset}
55 #end if 52 #end if
56 #if $chr.__str__.strip() != '': 53 #if str($chr).strip() != '':
57 -chr "$chr" 54 -chr "$chr"
58 #end if 55 #end if
59 $noLog 56 $noLog
60 #if $snpDb.genomeSrc == 'cached': 57 #if $snpDb.genomeSrc == 'cached':
61 -dataDir ${snpDb.genomeVersion.fields.path} 58 -dataDir ${snpDb.genomeVersion.fields.path}
62 #if $snpDb.extra_annotations and $snpDb.extra_annotations.__str__ != '': 59 #if $snpDb.extra_annotations and str($snpDb.extra_annotations) != '':
63 #echo " " 60 #echo " "
64 #echo ' '.join($snpDb.extra_annotations.__str__.split(',')) 61 #echo ' '.join(str($snpDb.extra_annotations).split(','))
65 #end if 62 #end if
66 #if $snpDb.regulation and $snpDb.regulation.__str__ != '': 63 #if $snpDb.regulation and str($snpDb.regulation) != '':
67 -reg #echo ' -reg '.join($snpDb.regulation.__str__.split(','))# 64 -reg #echo ' -reg '.join(str($snpDb.regulation).split(','))#
68 #end if 65 #end if
69 $snpDb.genomeVersion 66 $snpDb.genomeVersion
70 #elif $snpDb.genomeSrc == 'history': 67 #elif $snpDb.genomeSrc == 'history':
71 -dataDir ${snpDb.snpeff_db.extra_files_path} 68 -dataDir ${snpDb.snpeff_db.extra_files_path}
72 #if $snpDb.extra_annotations and $snpDb.extra_annotations.__str__ != '': 69 #if $snpDb.extra_annotations and str($snpDb.extra_annotations) != '':
73 #set xannotations = [' '] + $snpDb.extra_annotations.__str__.split(',') 70 #set xannotations = [' '] + str($snpDb.extra_annotations).split(',')
74 #echo " " 71 #echo " "
75 #echo ' -'.join($xannotations) 72 #echo ' -'.join($xannotations)
76 #end if 73 #end if
77 #if $snpDb.regulation and $snpDb.regulation.__str__ != '': 74 #if $snpDb.regulation and str($snpDb.regulation) != '':
78 -reg #echo ' -reg '.join($snpDb.regulation.__str__.split(','))# 75 -reg #echo ' -reg '.join(str($snpDb.regulation).split(','))#
79 #end if 76 #end if
80 ${snpDb.snpeff_db.metadata.genome_version} 77 ${snpDb.snpeff_db.metadata.genome_version}
81 #else 78 #else
82 -download 79 -download
83 $snpDb.genome_version 80 $snpDb.genome_version
84 #end if 81 #end if
85 $input > $snpeff_output ; 82 "$input" > "$snpeff_output";
86 #if $statsFile: 83 #if $statsFile:
87 #import os 84 #import os
88 #set $genes_file = str($statsFile) + '.genes.txt' 85 #set $genes_file = str($statsFile) + '.genes.txt'
89 #set $genes_file_name = os.path.split($genes_file)[-1] 86 #set $genes_file_name = os.path.split($genes_file)[-1]
90 mkdir $statsFile.files_path; 87 mkdir $statsFile.files_path;
91 mv $genes_file #echo os.path.join($statsFile.files_path, $genes_file_name)#; 88 mv "$genes_file" #echo os.path.join($statsFile.files_path, $genes_file_name)#;
92 #end if 89 #end if
93 #if $outputConditional.outputFormat == 'gatk' and $outputConditional.gatk_v1 90 #if $outputConditional.outputFormat == 'gatk' and $outputConditional.gatk_v1
94 ## Replace real SnpEff version with 2.0.5 to prevent this GATK 1.x error: "The version of SnpEff used to generate the SnpEff input file (x.x) is not currently supported by the GATK. Supported versions are: [2.0.5]" 91 ## Replace real SnpEff version with 2.0.5 to prevent this GATK 1.x error: "The version of SnpEff used to generate the SnpEff input file (x.x) is not currently supported by the GATK. Supported versions are: [2.0.5]"
95 sed -i -e 's/^\#\#SnpEffVersion="\(\S*\s\)/\#\#SnpEffVersion="2.0.5 - real is \1/' $snpeff_output 92 sed -i -e 's/^\#\#SnpEffVersion="\(\S*\s\)/\#\#SnpEffVersion="2.0.5 - real is \1/' "$snpeff_output"
96 #end if 93 #end if
97 ]]> 94 ]]>
98 </command> 95 </command>
99 <inputs> 96 <inputs>
100 <param format="vcf,tabular,pileup,bed" name="input" type="data" label="Sequence changes (SNPs, MNPs, InDels)"/> 97 <param format="vcf,tabular,pileup,bed" name="input" type="data" label="Sequence changes (SNPs, MNPs, InDels)"/>
210 <param name="spliceRegionExonSize" type="integer" value="" min="1" max="10" optional="true" label="Set size for splice site region within exons. Default: 3 bases"/> 207 <param name="spliceRegionExonSize" type="integer" value="" min="1" max="10" optional="true" label="Set size for splice site region within exons. Default: 3 bases"/>
211 <param name="spliceRegionIntronMin" type="integer" value="" min="1" max="10" optional="true" label="Set minimum number of bases for splice site region within intron. Default: 3 bases"/> 208 <param name="spliceRegionIntronMin" type="integer" value="" min="1" max="10" optional="true" label="Set minimum number of bases for splice site region within intron. Default: 3 bases"/>
212 <param name="spliceRegionIntronMax" type="integer" value="" min="1" max="10" optional="true" label="Set maximum number of bases for splice site region within intron. Default: 8 bases"/> 209 <param name="spliceRegionIntronMax" type="integer" value="" min="1" max="10" optional="true" label="Set maximum number of bases for splice site region within intron. Default: 8 bases"/>
213 </when> 210 </when>
214 </conditional> 211 </conditional>
215
216 <param name="filterHomHet" type="select" display="radio" label="Filter homozygous / heterozygous changes">
217 <option value="no_filter" selected="true">No filter (analyze everything)</option>
218 <option value="-hom">Analyze homozygous sequence changes only</option>
219 <option value="-het">Analyze heterozygous sequence changes only</option>
220 </param>
221 212
222 <param name="annotations" type="select" display="checkboxes" multiple="true" label="Annotation options"> 213 <param name="annotations" type="select" display="checkboxes" multiple="true" label="Annotation options">
223 <option value="-cancer">Perform 'cancer' comparisons (somatic vs. germline)</option> 214 <option value="-cancer">Perform 'cancer' comparisons (somatic vs. germline)</option>
224 <option value="-canon">Only use canonical transcripts</option> 215 <option value="-canon">Only use canonical transcripts</option>
225 <option value="-geneId">Use gene ID instead of gene name (VCF output)</option> 216 <option value="-geneId">Use gene ID instead of gene name (VCF output)</option>
255 <param name="effects" type="select" display="checkboxes" multiple="true" label="Filter output: do not report these Effects"> 246 <param name="effects" type="select" display="checkboxes" multiple="true" label="Filter output: do not report these Effects">
256 <option value="CDS">CDS (coding_sequence_variant) The variant hits a CDS. MODIFIER</option> 247 <option value="CDS">CDS (coding_sequence_variant) The variant hits a CDS. MODIFIER</option>
257 <option value="CHROMOSOME_LARGE_DELETION">CHROMOSOME_LARGE_DELETION (chromosome) A large parte (over 1%) of the chromosome was deleted. HIGH</option> 248 <option value="CHROMOSOME_LARGE_DELETION">CHROMOSOME_LARGE_DELETION (chromosome) A large parte (over 1%) of the chromosome was deleted. HIGH</option>
258 <option value="CODON_CHANGE">CODON_CHANGE (coding_sequence_variant) One or many codons are changed e.g.: An MNP of size multiple of 3 MODERATE</option> 249 <option value="CODON_CHANGE">CODON_CHANGE (coding_sequence_variant) One or many codons are changed e.g.: An MNP of size multiple of 3 MODERATE</option>
259 <option value="CODON_INSERTION">CODON_INSERTION (inframe_insertion) One or many codons are inserted e.g.: An insert multiple of three in a codon boundary MODERATE</option> 250 <option value="CODON_INSERTION">CODON_INSERTION (inframe_insertion) One or many codons are inserted e.g.: An insert multiple of three in a codon boundary MODERATE</option>
260 <option value="CODON_CHANGE_PLUS CODON_INSERTION">CODON_CHANGE_PLUS CODON_INSERTION (disruptive_inframe_insertion) One codon is changed and one or many codons are inserted e.g.: An insert of size multiple of three, not at codon boundary MODERATE</option> 251 <option value="CODON_CHANGE_PLUS_CODON_INSERTION">CODON_CHANGE_PLUS_CODON_INSERTION (disruptive_inframe_insertion) One codon is changed and one or many codons are inserted e.g.: An insert of size multiple of three, not at codon boundary MODERATE</option>
261 <option value="CODON_DELETION">CODON_DELETION (inframe_deletion) One or many codons are deleted e.g.: A deletion multiple of three at codon boundary MODERATE</option> 252 <option value="CODON_DELETION">CODON_DELETION (inframe_deletion) One or many codons are deleted e.g.: A deletion multiple of three at codon boundary MODERATE</option>
262 <option value="CODON_CHANGE_PLUS CODON_DELETION">CODON_CHANGE_PLUS CODON_DELETION (disruptive_inframe_deletion) One codon is changed and one or more codons are deleted e.g.: A deletion of size multiple of three, not at codon boundary MODERATE</option> 253 <option value="CODON_CHANGE_PLUS_CODON_DELETION">CODON_CHANGE_PLUS_CODON_DELETION (disruptive_inframe_deletion) One codon is changed and one or more codons are deleted e.g.: A deletion of size multiple of three, not at codon boundary MODERATE</option>
263 <option value="DOWNSTREAM">DOWNSTREAM (downstream_gene_variant) Downstream of a gene (default length: 5K bases) MODIFIER</option> 254 <option value="DOWNSTREAM">DOWNSTREAM (downstream_gene_variant) Downstream of a gene (default length: 5K bases) MODIFIER</option>
264 <option value="EXON">EXON (exon_variant) The variant hits an exon (from a non-coding transcript) or a retained intron. MODIFIER</option> 255 <option value="EXON">EXON (exon_variant) The variant hits an exon (from a non-coding transcript) or a retained intron. MODIFIER</option>
265 <option value="EXON_DELETED">EXON_DELETED (exon_loss_variant) A deletion removes the whole exon. HIGH</option> 256 <option value="EXON_DELETED">EXON_DELETED (exon_loss_variant) A deletion removes the whole exon. HIGH</option>
266 <option value="FRAME_SHIFT">FRAME_SHIFT (frameshift_variant) Insertion or deletion causes a frame shift e.g.: An indel size is not multple of 3 HIGH</option> 257 <option value="FRAME_SHIFT">FRAME_SHIFT (frameshift_variant) Insertion or deletion causes a frame shift e.g.: An indel size is not multple of 3 HIGH</option>
267 <option value="GENE">GENE (gene_variant) The variant hits a gene. MODIFIER</option> 258 <option value="GENE">GENE (gene_variant) The variant hits a gene. MODIFIER</option>
311 You can prepend any string you want to the chromosome name. 302 You can prepend any string you want to the chromosome name.
312 </help> 303 </help>
313 <validator type="regex" message="No whitespace allowed">^\S*$</validator> 304 <validator type="regex" message="No whitespace allowed">^\S*$</validator>
314 </param> 305 </param>
315 <param name="generate_stats" type="boolean" truevalue="" falsevalue="-noStats" checked="true" label="Produce Summary Stats"/> 306 <param name="generate_stats" type="boolean" truevalue="" falsevalue="-noStats" checked="true" label="Produce Summary Stats"/>
316 <param name="noLog" type="boolean" truevalue="-noLog" falsevalue="" checked="true" label="Do not report usage statistics to server"/> 307 <param name="noLog" type="boolean" truevalue="-noLog" falsevalue="" checked="true" label="Suppress reporting usage statistics to server"/>
317 </inputs> 308 </inputs>
318 <outputs> 309 <outputs>
319 <data format="vcf" name="snpeff_output" > 310 <data format="vcf" name="snpeff_output" >
320 <change_format> 311 <change_format>
321 <when input="outputConditional.outputFormat" value="bed" format="bed" /> 312 <when input="outputConditional.outputFormat" value="bed" format="bed" />
342 <param name="inputFormat" value="vcf"/> 333 <param name="inputFormat" value="vcf"/>
343 <param name="outputFormat" value="vcf"/> 334 <param name="outputFormat" value="vcf"/>
344 <param name="genomeSrc" value="named"/> 335 <param name="genomeSrc" value="named"/>
345 <param name="genome_version" value="testCase"/> 336 <param name="genome_version" value="testCase"/>
346 <param name="udLength" value="0"/> 337 <param name="udLength" value="0"/>
347 <param name="filterHomHet" value="no_filter"/>
348 <param name="generate_stats" value="False"/> 338 <param name="generate_stats" value="False"/>
349 <param name="filterOut" value="+-no-upstream"/> 339 <param name="filterOut" value="+-no-upstream"/>
350 <output name="snpeff_output"> 340 <output name="snpeff_output">
351 <assert_contents> 341 <assert_contents>
352 <has_text text="EFF=" /> 342 <has_text text="EFF=" />
360 <param name="inputFormat" value="vcf"/> 350 <param name="inputFormat" value="vcf"/>
361 <param name="outputFormat" value="vcf"/> 351 <param name="outputFormat" value="vcf"/>
362 <param name="genomeSrc" value="named"/> 352 <param name="genomeSrc" value="named"/>
363 <param name="genome_version" value="testCase"/> 353 <param name="genome_version" value="testCase"/>
364 <param name="udLength" value="0"/> 354 <param name="udLength" value="0"/>
365 <param name="filterHomHet" value="+-het"/>
366 <!--
367 <param name="filterOut" value=""/>
368 -->
369 <param name="generate_stats" value="False"/>
370 <output name="snpeff_output">
371 <assert_contents>
372 <!-- Check that NO effects were added since -het is set -->
373 <not_has_text text="EFF=NON_SYNONYMOUS_CODING" />
374 </assert_contents>
375 </output>
376 </test>
377
378 <test>
379 <param name="input" ftype="vcf" value="vcf_homhet.vcf"/>
380 <param name="inputFormat" value="vcf"/>
381 <param name="outputFormat" value="vcf"/>
382 <param name="genomeSrc" value="named"/>
383 <param name="genome_version" value="testCase"/>
384 <param name="udLength" value="0"/>
385 <param name="filterHomHet" value="no_filter"/>
386 <!-- 355 <!--
387 <param name="filterOut" value=""/> 356 <param name="filterOut" value=""/>
388 --> 357 -->
389 <param name="generate_stats" value="False"/> 358 <param name="generate_stats" value="False"/>
390 <output name="snpeff_output"> 359 <output name="snpeff_output">
404 <param name="inputFormat" value="vcf"/> 373 <param name="inputFormat" value="vcf"/>
405 <param name="outputFormat" value="vcf"/> 374 <param name="outputFormat" value="vcf"/>
406 <param name="genomeSrc" value="named"/> 375 <param name="genomeSrc" value="named"/>
407 <param name="genome_version" value="testCase"/> 376 <param name="genome_version" value="testCase"/>
408 <param name="udLength" value="0"/> 377 <param name="udLength" value="0"/>
409 <param name="filterHomHet" value="no_filter"/>
410 <param name="filterOut" value="+-no-upstream"/> 378 <param name="filterOut" value="+-no-upstream"/>
411 <param name="generate_stats" value="False"/> 379 <param name="generate_stats" value="False"/>
412 <output name="snpeff_output"> 380 <output name="snpeff_output">
413 <assert_contents> 381 <assert_contents>
414 <not_has_text text="UPSTREAM" /> 382 <not_has_text text="UPSTREAM" />
422 390
423 This tool calculate the effect of variants (SNPs/MNPs/Insertions) and deletions. 391 This tool calculate the effect of variants (SNPs/MNPs/Insertions) and deletions.
424 392
425 @EXTERNAL_DOCUMENTATION@ 393 @EXTERNAL_DOCUMENTATION@
426 394
427 @CITATION_SECTION@
428
429 ]]> 395 ]]>
430 </help> 396 </help>
431 <expand macro="citations" /> 397 <expand macro="citations" />
432 </tool> 398 </tool>
433 399