Mercurial > repos > jjohnson > snpeff
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 |
