|
0
|
1 <tool id="krona" name="Krona chart" version="1.1.0">
|
|
|
2 <description>of taxonomic abundance</description>
|
|
|
3 <command>
|
|
|
4 #if $type.program == 'galaxy':
|
|
|
5 ktImportGalaxy
|
|
|
6 #else if $type.program == 'blast':
|
|
|
7 ktImportBLAST
|
|
|
8 -e ${type.factor}
|
|
|
9 ${type.random}
|
|
|
10 #else
|
|
|
11 ktImportTaxonomy
|
|
|
12 ${type.summary}
|
|
|
13 #end if
|
|
|
14 -o $output
|
|
|
15 ${type.include}
|
|
|
16
|
|
|
17 -d $depth
|
|
|
18
|
|
|
19 ## uncomment for isolated intranets (see README)
|
|
|
20 ##
|
|
|
21 ##-u /static/krona
|
|
|
22
|
|
|
23 #if $datasets.multiple == 'single':
|
|
|
24 $input
|
|
|
25 #else
|
|
|
26 $datasets.combine
|
|
|
27 #for $input in $datasets.inputs
|
|
|
28 ${input.file},'${input.name}'
|
|
|
29 #end for
|
|
|
30 #end if
|
|
|
31 </command>
|
|
|
32
|
|
|
33 <inputs>
|
|
|
34 <conditional name="type">
|
|
|
35 <param name="program" label="Input type" type="select">
|
|
|
36 <option value="galaxy">Galaxy taxonomic representation</option>
|
|
|
37 <option value="blast">Tabular BLAST results</option>
|
|
|
38 <option value="taxonomy">Taxonomy ID list</option>
|
|
|
39 </param>
|
|
|
40 <when value="galaxy">
|
|
|
41 <param name="factor" type="hidden" value=""/>
|
|
|
42 <param name="include" type="hidden" value=""/>
|
|
|
43 </when>
|
|
|
44 <when value="blast">
|
|
|
45 <param name="factor" type="float" value="10" label="E-value factor" help="For each query, hits with e-values within this factor of the best hit's e-value will be included when computing the lowest common ancestor."/>
|
|
|
46 <param name="include" type="boolean" truevalue="-i" falsevalue="" label="Include reads with no hits (comment lines must be present in BLAST results)"/>
|
|
|
47 <param name="random" type="boolean" truevalue="-r" falsevalue="" label="Random" help="Pick from hits within E-value factor randomly instead of finding the lowest common ancestor"/>
|
|
|
48 </when>
|
|
|
49 <when value="taxonomy">
|
|
|
50 <param name="factor" type="hidden" value=""/>
|
|
|
51 <param name="summary" type="boolean" truevalue="-m 1" falsevalue="" label="Summarized" help="The first column is counts rather than query IDs."/>
|
|
|
52 <param name="include" type="boolean" truevalue="-i" falsevalue="" label="Include reads with no hits"/>
|
|
|
53 </when>
|
|
|
54 </conditional>
|
|
|
55 <param name="depth" label="Depth" type="integer" help="The maximum depth to show in the chart (0 for unlimited)." value="0"/>
|
|
|
56 <conditional name="datasets">
|
|
|
57 <param name="multiple" type="select" label="Number of datasets">
|
|
|
58 <option value="single">Single</option>
|
|
|
59 <option value="multiple">Multiple</option>
|
|
|
60 </param>
|
|
|
61 <when value="single">
|
|
|
62 <param name="input" type="data" format="tabular" label="Dataset"/>
|
|
|
63 </when>
|
|
|
64 <when value="multiple">
|
|
|
65 <param
|
|
|
66 name="combine"
|
|
|
67 type="boolean"
|
|
|
68 truevalue="-c"
|
|
|
69 falsevalue=""
|
|
|
70 label="Combine data sets"
|
|
|
71 help="By default, data sets can be stepped through to see varying abundances. Combining them will sum the data sets at each node to represent total abundances."
|
|
|
72 />
|
|
|
73 <repeat name="inputs" title="Data set">
|
|
|
74 <param
|
|
|
75 name="file"
|
|
|
76 format="tabular"
|
|
|
77 type="data"
|
|
|
78 label="Data"
|
|
|
79 />
|
|
|
80 <param name="name" label="Name" type="text">
|
|
|
81 <sanitizer>
|
|
|
82 <valid initial="string.printable">
|
|
|
83 <remove value="""/>
|
|
|
84 </valid>
|
|
|
85 <mapping initial="none">
|
|
|
86 <add source=""" target="&quot;"/>
|
|
|
87 </mapping>
|
|
|
88 </sanitizer>
|
|
|
89 </param>
|
|
|
90 </repeat>
|
|
|
91 </when>
|
|
|
92 </conditional>
|
|
|
93 </inputs>
|
|
|
94
|
|
|
95 <outputs>
|
|
|
96 <data format="html" name="output"/>
|
|
|
97 </outputs>
|
|
|
98
|
|
|
99 <help>
|
|
|
100
|
|
|
101 **What it does**
|
|
|
102
|
|
|
103 Creates an interactive Krona_ chart of taxonomic abundance. Requires a KronaTools installation (2.0 or higher) with taxonomy downloaded.
|
|
|
104
|
|
|
105 .. _Krona: http://krona.sourceforge.net
|
|
|
106
|
|
|
107 **Inputs**
|
|
|
108
|
|
|
109 The inputs can be Galaxy taxonomic representations, for example::
|
|
|
110
|
|
|
111 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
|
|
|
112 1L_EYKX4VC01BXWX1_265 9606 root Eukaryota Metazoa n n Chordata Craniata Gnathostomata Mammalia n Euarchontoglires Primates Haplorrhini Hominoidea Hominidae n n n Homo n Homo sapiens n 1430919
|
|
|
113
|
|
|
114 ...or, they can be tabular BLAST results containing gi numbers, which will be used to look up taxonomy IDs and find lowest common ancestors for best hits (comment lines are optional but required for showing queries with no hits in the chart)::
|
|
|
115
|
|
|
116 1 2 3 4 5 6 7 8 9 10 11 12
|
|
|
117 read1 gi|211853080|gb|EQ846228.1| 97.33 150 1 3 6 154 9518999 9519146 5e-64 252
|
|
|
118 read2 gi|211853080|gb|EQ846228.1| 96.69 151 0 5 6 154 9496147 9496294 2e-62 246
|
|
|
119
|
|
|
120 ...or, they can be lists of taxonomy IDs preceded by query IDs or counts (these formats can be exported by ktClassifyBLAST or MEGAN), but **Summarized** must be checked if they are counts::
|
|
|
121
|
|
|
122 1 2
|
|
|
123 read1 9606
|
|
|
124 read2 9616
|
|
|
125
|
|
|
126 Summarized::
|
|
|
127
|
|
|
128 1 2
|
|
|
129 243 9606
|
|
|
130 17 9616
|
|
|
131
|
|
|
132 Both forms can optionally have scores::
|
|
|
133
|
|
|
134 1 2 3
|
|
|
135 read1 9606 94.34
|
|
|
136 read2 9616 32.17
|
|
|
137
|
|
|
138 Summarized::
|
|
|
139
|
|
|
140 1 2 3
|
|
|
141 243 9606 94.34
|
|
|
142 17 9616 32.17
|
|
|
143
|
|
|
144 </help>
|
|
|
145 <requirements>
|
|
|
146 <requirement type="binary">ktImportGalaxy</requirement>
|
|
|
147 <requirement type="binary">ktImportBLAST</requirement>
|
|
|
148 <requirement type="binary">ktImportTaxonomy</requirement>
|
|
|
149 </requirements>
|
|
|
150 </tool>
|
