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> |