0
|
1 <tool id="velvet" name="velvet" version="0.1.0">
|
|
2 <requirements>
|
|
3 <requirement type="package" version="1.2.10">velvet</requirement>
|
|
4 </requirements>
|
|
5 <stdio>
|
|
6 <exit_code range="1:" />
|
|
7 </stdio>
|
|
8
|
|
9 <command><![CDATA[
|
|
10 velveth
|
|
11 out
|
|
12 $kmer
|
|
13 #for $i in $files:
|
|
14 -${i.filetype}
|
|
15
|
|
16 #if str($i.paired_type.paired_type_selector) == "paired":
|
|
17
|
|
18 #if $i.__index__ > 1:
|
|
19 -shortPaired$i.__index__
|
|
20 #else
|
|
21 -shortPaired
|
|
22 #end if
|
|
23
|
|
24 -separate $i.paired_type.input1 $i.paired_type.input2
|
|
25 #end if
|
|
26
|
|
27 #if str($i.paired_type.paired_type_selector) == "single"
|
|
28
|
|
29 #if $i.__index__ > 1
|
|
30 -short$i.__index__
|
|
31 #else
|
|
32 -short
|
|
33 #end if
|
|
34
|
|
35 $i.paired_type.input1
|
|
36 #end if
|
|
37
|
|
38 #if str($i.paired_type.paired_type_selector) == "paired_collection"
|
|
39 #if $i.__index__ > 1
|
|
40 -shortPaired$i.__index__
|
|
41 #else
|
|
42 -shortPaired
|
|
43 #end if
|
|
44 -separate $i.paired_type.input1.forward $i.paired_type.input1.reverse
|
|
45 #end if
|
|
46 #if str($i.paired_type.paired_type_selector) == "paired_iv"
|
|
47 #if $i.__index__ > 1
|
|
48 -shortPaired$i.__index__
|
|
49 #else
|
|
50 -shortPaired
|
|
51 #end if
|
|
52 $i.paired_type.input1
|
|
53 #end if
|
|
54 #end for
|
|
55
|
|
56 &&
|
|
57
|
|
58 velvetg out > $log
|
|
59 ]]></command>
|
|
60 <inputs>
|
|
61 <param name="kmer" type="integer" value="29" label="K-mer" help="The kmer value to use for velvet - must be an odd integer between 19 and 191" />
|
|
62 <repeat name="files" title="Input files" min="1">
|
|
63 <param name="filetype" label="Input file type" type="select" help="Input file type">
|
|
64 <option value="fastq">Fastq</option>
|
|
65 <option value="fasta">Fasta</option>
|
|
66 <option value="bam">Bam</option>
|
|
67 </param>
|
|
68 <conditional name="paired_type">
|
|
69 <param name="paired_type_selector" type="select" label="Single or paired end reads" help="Select between paired and single end data">
|
|
70 <option value="paired">Paired</option>
|
|
71 <option value="single">Single</option>
|
|
72 <option value="paired_collection">Paired Collection</option>
|
|
73 <option value="paired_iv">Paired Interleaved</option>
|
|
74 </param>
|
|
75 <when value="paired">
|
|
76 <param name="input1" type="data" format="fastqsanger,fasta,bam" label="Select first set of reads" help="Specify dataset with forward reads"/>
|
|
77 <param name="input2" type="data" format="fastqsanger,fasta,bam" label="Select second set of reads" help="Specify dataset with reverse reads"/>
|
|
78 </when>
|
|
79 <when value="single">
|
|
80 <param name="input1" type="data" format="fastqsanger,fasta,bam" label="Select the reads" help="Specify dataset with the reads"/>
|
|
81 </when>
|
|
82 <when value="paired_collection">
|
|
83 <param name="input1" format="fastqsanger,fasta" type="data_collection" collection_type="paired" label="Select a paired collection" help="See help section for an explanation of dataset collections"/>
|
|
84 </when>
|
|
85 <when value="paired_iv">
|
|
86 <param name="input1" type="data" format="fastqsanger,fasta,bam" label="Select set of interleaved reads" help="Specify dataset with interleaved paired reads"/>
|
|
87 </when>
|
|
88 </conditional>
|
|
89 </repeat>
|
|
90
|
|
91 </inputs>
|
|
92 <outputs>
|
|
93 <data format="fasta" name="contigs" label="${tool.name} on ${on_string}: Contigs" from_work_dir="out/contigs.fa"/>
|
|
94 <data format="tabular" name="stats" label="${tool.name} on ${on_string}: Contig Stats" from_work_dir="out/stats.txt"/>
|
|
95 <data format="tabular" name="graph" label="${tool.name} on ${on_string}: Last Graph" from_work_dir="out/LastGraph"/>
|
|
96 <data format="tabular" name="log" label="${tool.name} on ${on_string}: log" />
|
|
97 </outputs>
|
|
98 <tests>
|
|
99 <test>
|
|
100 <param name="kmer" value="29" />
|
|
101 <param name="files_0|filetype" value="fasta" />
|
|
102 <param name="files_0|paired_type|paired_type_selector" value="single" />
|
|
103 <param name="files_0|paired_type|input1" value="test_reads.fa" ftype="fasta" />
|
|
104 <output name="contigs" ftype="fasta" file="contigs.fa" />
|
|
105 </test>
|
|
106 </tests>
|
|
107 <help><![CDATA[
|
|
108 The Velvet wrapper, wraps both velveth and velvetg for single kmer value with no parameters supplied to velvetg. It is mainly for teaching how kmer size affects the assembly.
|
|
109
|
|
110 It returns the contigs.fa file, the contig_stats.txt file, the LastGraph file and captures the screen output to the log.
|
|
111
|
|
112 This tool is mainly wrapped for assembly teaching purposes. It should not be used for actual assemblies. Please use the Velvet Optimiser for this purpose.
|
|
113
|
|
114 The kmer value supplied needs to be an odd integer.
|
|
115 ]]></help>
|
|
116 <citations>
|
|
117 <citation type="doi">10.1101/gr.074492.107</citation>
|
|
118 <citation type="doi">10.1371/journal.pone.0008407</citation>
|
|
119 </citations>
|
|
120 </tool> |