|
0
|
1 <tool id="vcftools_merge" name="Merge" version="0.1">
|
|
|
2 <description>multiple VCF datasets</description>
|
|
|
3
|
|
|
4 <requirements>
|
|
|
5 <requirement type="package">vcftools</requirement>
|
|
|
6 </requirements>
|
|
|
7
|
|
|
8 <command>
|
|
|
9 ## Preprocessing for each dataset.
|
|
|
10 #set dataset_names = []
|
|
|
11 #for i, $input in enumerate( $inputs ):
|
|
|
12 ## Sort file.
|
|
|
13 vcf-sort ${input.file} > ${i}.vcf.sorted ;
|
|
|
14
|
|
|
15 ## Compress.
|
|
|
16 bgzip ${i}.vcf.sorted ;
|
|
|
17
|
|
|
18 ## Index.
|
|
|
19 tabix -p vcf ${i}.vcf.sorted.gz ;
|
|
|
20
|
|
|
21 #silent dataset_names.append( str($i) + '.vcf.sorted.gz' )
|
|
|
22 #end for
|
|
|
23
|
|
|
24 ## Merge.
|
|
|
25 vcf-merge
|
|
|
26 #echo ' '.join( dataset_names ) # > ${output}
|
|
|
27
|
|
|
28 </command>
|
|
|
29 <inputs>
|
|
|
30 <repeat name="inputs" title="Datasets to Merge" min="2">
|
|
|
31 <param name="file" label="Dataset" type="data" format="vcf"/>
|
|
|
32 </repeat>
|
|
|
33 </inputs>
|
|
|
34
|
|
|
35 <outputs>
|
|
|
36 <data name="output" format="vcf"/>
|
|
|
37 </outputs>
|
|
|
38
|
|
|
39 <stdio>
|
|
|
40 <regex match=".*" source="both" level="log" description="tool progress"/>
|
|
|
41 </stdio>
|
|
|
42
|
|
|
43 <tests>
|
|
|
44 <!-- Cannot specify multiple repeats in test framework right now, so no tests possible. -->
|
|
|
45 </tests>
|
|
|
46
|
|
|
47 <help>
|
|
|
48 Please see the VCFtools `documentation`__ for help and further information.
|
|
|
49
|
|
|
50 .. __: http://vcftools.sourceforge.net/docs.html
|
|
|
51 </help>
|
|
|
52 </tool>
|
