Mercurial > repos > devteam > vcftools_annotate
annotate vcftools_annotate.xml @ 6:a6e0edfcfdae draft default tip
planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
| author | devteam |
|---|---|
| date | Tue, 13 Oct 2015 12:59:18 -0400 |
| parents | 765f626ebd5a |
| children |
| rev | line source |
|---|---|
| 0 | 1 <tool id="vcftools_annotate" name="Annotate" version="0.1"> |
|
6
a6e0edfcfdae
planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents:
0
diff
changeset
|
2 <description>a VCF dataset with custom filters</description> |
| 0 | 3 |
| 4 <requirements> | |
| 5 <requirement type="binary">echo</requirement> | |
| 6 <requirement type="package" version="0.1.11">vcftools</requirement> | |
| 7 </requirements> | |
| 8 | |
| 9 <command> | |
| 10 ## Generate filter file. | |
| 11 echo "{ tag => '${tag}', name => '${filter}', desc => '${description}', test => sub { my @t = split('%', @\\$MATCH[0]); return @t[0] $condition ? \\$PASS : \\$FAIL }, }," > f.txt ; | |
| 12 | |
| 13 ## Annotate. | |
| 14 vcf-annotate -f f.txt ${input} > ${output} | |
| 15 </command> | |
| 16 <inputs> | |
| 17 <param name="input" label="Input" type="data" format="vcf"/> | |
| 18 <param name="tag" label="Tag name" type="text"/> | |
| 19 <param name="description" label="Tag description" type="text"> | |
| 20 <sanitizer sanitize="False"/> | |
| 21 </param> | |
| 22 <param name="filter" label="Filter name" type="text"/> | |
| 23 <param name="condition" label="Filter condition" type="text"> | |
| 24 <sanitizer sanitize="False"/> | |
| 25 </param> | |
| 26 </inputs> | |
| 27 | |
| 28 <outputs> | |
| 29 <data name="output" format="vcf"/> | |
| 30 </outputs> | |
| 31 | |
| 32 <stdio> | |
| 33 <regex match=".*" source="both" level="log" description="tool progress"/> | |
| 34 </stdio> | |
| 35 | |
| 36 <tests> | |
| 37 <test> | |
| 38 <param name="input" value="test_in1.vcf" /> | |
| 39 <param name="tag" value="FORMAT/FREQ" /> | |
| 40 <param name="description" value="MinAF [7]"/> | |
| 41 <param name="filter" value="MinAF" /> | |
| 42 <param name="condition" value=">= 7" /> | |
| 43 <!-- 2 lines diff because command line with full file path is included in output VCF, and | |
| 44 it not possible to match full file path. --> | |
| 45 <output name="output" file="test_out1.vcf" lines_diff="2" /> | |
| 46 </test> | |
| 47 </tests> | |
| 48 | |
| 49 <help> | |
|
6
a6e0edfcfdae
planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents:
0
diff
changeset
|
50 Annotates VCF dataset with custom annotations. For example, if this format tag is used for allele frequency: |
|
a6e0edfcfdae
planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents:
0
diff
changeset
|
51 |
|
a6e0edfcfdae
planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents:
0
diff
changeset
|
52 ##FORMAT=<ID=FREQ,Number=1,Type=String,Description="Variant allele frequency"> |
| 0 | 53 |
|
6
a6e0edfcfdae
planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents:
0
diff
changeset
|
54 you can add a filter for allele frequency using "FORMAT/FREQ" as the tag name and the condition ">= [desired allele freq]" |
|
a6e0edfcfdae
planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents:
0
diff
changeset
|
55 |
|
a6e0edfcfdae
planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents:
0
diff
changeset
|
56 Please see the VCFtools `documentation`__ for help and further information. |
|
a6e0edfcfdae
planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents:
0
diff
changeset
|
57 |
|
a6e0edfcfdae
planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents:
0
diff
changeset
|
58 .. __: http://vcftools.sourceforge.net/perl_module.html#vcf-annotate |
| 0 | 59 </help> |
| 60 </tool> |
