Mercurial > repos > iuc > bcftools_plugin_split_vep
diff bcftools_plugin_split_vep.xml @ 0:3b6cd8086498 draft
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bcftools commit c45135e52ae5039e09272ac6f504d0ceb574aa70
author | iuc |
---|---|
date | Sat, 23 Jul 2022 13:49:21 +0000 |
parents | |
children | 70276425d001 |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/bcftools_plugin_split_vep.xml Sat Jul 23 13:49:21 2022 +0000 @@ -0,0 +1,112 @@ +<?xml version='1.0' encoding='utf-8'?> +<tool name="bcftools @EXECUTABLE@" id="bcftools_plugin_@PLUGIN_ID@" version="@TOOL_VERSION@"> + <description>plugin Extracts fields from structured annotations such as INFO/CSQ</description> + <macros> + <token name="@EXECUTABLE@">split-vep</token> + <token name="@PLUGIN_ID@">split_vep</token> + <import>macros.xml</import> + </macros> + <expand macro="bio_tools" /> + <expand macro="requirements" /> + <expand macro="version_command" /> + <command detect_errors="aggressive"><![CDATA[ +@PREPARE_ENV@ +@PREPARE_INPUT_FILE@ +#set $section = $sec_restrict +@PREPARE_TARGETS_FILE@ +@PREPARE_REGIONS_FILE@ + +bcftools plugin @EXECUTABLE@ +## VCF input section +#set $section = $sec_restrict +@INCLUDE@ +@EXCLUDE@ +@REGIONS@ +@TARGETS@ + +@OUTPUT_TYPE@ + +## Primary Input/Outputs +@INPUT_FILE@ +## Plugin options +## ToDo: Parameters not wrapped yet: -A, -f (not for BCF/VCF output!), -l, -S, -x +-a '$a' +-c '$c' +$d +#if $p: + -p '$p' +#end if +#if $s: + -s '$s' +#end if +> '$output_file' +]]> + </command> + <inputs> + <expand macro="macro_input" /> + <param argument="-a" type="text" value="CSQ" label="Select INFO annotation to parse" help="Usually, annotations are located in the CSQ, ANN or EFF INFO field, though any arbitrary tag may be entered here." /> + <param argument="-c" type="text" label="Enter fields to be extracted, listed either as indexes or names" help="The default type of the new annotation is String but can be also set to Integer/Int or Float/Real by adding ':type' to the field name or index. Examples can be found in the tool help." /> + <param argument="-d" type="boolean" label="Enable Output per transcript/allele consequences on a new line rather rather than as comma-separated fields on a single line" truevalue="-d" falsevalue="" /> + <param argument="-p" type="text" optional="true" label="Enter prefix of INFO annotations to be created after splitting the CSQ string (optional)" /> + <param argument="-s" type="text" optional="true" label="Select transcripts to extract by type and/or consequence (optional)" help="Argument has the form TR:CSQ. TR = transcript: worst,primary,all (default: all). CSQ = consequence: any,missense,missense+,etc (default: any). Examples can be found in the tool help." /> + + <section name="sec_restrict" expanded="false" title="Restrict to"> + <expand macro="macro_restrict" /> + <expand macro="macro_restrict" type="target" label_type="Target" /> + <expand macro="macro_include" /> + <expand macro="macro_exclude" /> + </section> + <expand macro="macro_select_output_type" /> + </inputs> + <outputs> + <expand macro="macro_vcf_output" /> + </outputs> + <tests> + <test> + <param name="input_file" ftype="vcf" value="plugin2.vcf" /> + <param name="a" value="ANN" /> + <param name="c" value="IMPACT,gnomAD_AF:Float" /> + <param name="p" value="NW_" /> + <param name="s" value="worst" /> + <param name="output_type" value="v" /> + <output name="output_file" ftype="vcf"> + <assert_contents> + <has_line_matching expression="##INFO=<ID=NW_IMPACT,.+" /> + <has_line_matching expression="##INFO=<ID=NW_gnomAD_AF,.+" /> + <has_line_matching expression=".+NW_IMPACT=MODERATE;NW_gnomAD_AF=0\.08.+" /> + <has_line_matching expression=".+NW_IMPACT=LOW;NW_gnomAD_AF=0\.9443.+" /> + <has_line_matching expression=".+NW_IMPACT=MODIFIER;NW_gnomAD_AF=\..+" /> + </assert_contents> + </output> + </test> + </tests> + <help><![CDATA[ +===================================== + bcftools @EXECUTABLE@ plugin +===================================== + +This plugin allows to extract fields from structured annotations such as ``INFO/CSQ`` created by bcftools/csq or VEP. These fields are then added to the VCF as new ``INFO`` fields. + +Examples:: + + # Extract Consequence, IMPACT and gene SYMBOL of the most severe consequence into + # INFO annotations starting with the prefix "vep". For brevity, the columns can + # be given also as 0-based indexes + bcftools +split-vep -c Consequence,IMPACT,SYMBOL -s worst -p vep file.vcf.gz + bcftools +split-vep -c 1-3 -s worst -p vep file.vcf.gz + + # Extract gnomAD_AF subfield into a new INFO/gnomAD_AF annotation of Type=Float so that + # numeric filtering can be used. + bcftools +split-vep -c gnomAD_AF:Float file.vcf.gz -i 'gnomAD_AF<0.001' + + # Similar to above, but add the annotation only if the consequence severity is missense + # or equivalent. + bcftools +split-vep -c gnomAD_AF:Float -s :missense file.vcf.gz + +@REGIONS_HELP@ +@TARGETS_HELP@ +@EXPRESSIONS_HELP@ + + ]]></help> + <expand macro="citations" /> +</tool>