|
1
|
1 <tool id="classify_blast" name="Classify BLAST" version="1.0.0">
|
|
|
2 <description>results using KronaTools</description>
|
|
|
3
|
|
|
4 <requirements>
|
|
|
5 <requirement name="package" version="2.8.1">krona</requirement>
|
|
|
6 </requirements>
|
|
|
7
|
|
|
8 <command>
|
|
|
9 ktClassifyBLAST
|
|
|
10 -o $output
|
|
|
11 -e $factor
|
|
|
12 $random
|
|
|
13 #if $score == 'percentIdentity':
|
|
|
14 -p
|
|
|
15 #elif $score == 'bitScore':
|
|
|
16 -b
|
|
|
17 #end if
|
|
|
18 $summarize
|
|
|
19 #for $input in $inputs
|
|
|
20 ${input.input}
|
|
|
21 #end for
|
|
|
22 </command>
|
|
|
23
|
|
|
24 <inputs>
|
|
|
25 <repeat name="inputs" title="BLAST result">
|
|
|
26 <param name="input" format="tabular" type="data" label="BLAST result"/>
|
|
|
27 </repeat>
|
|
|
28 <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 (or picking randomly)."/>
|
|
|
29 <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"/>
|
|
|
30 <param name="score" type="select" label="Score by">
|
|
|
31 <option value="evalue">E-value</option>
|
|
|
32 <option value="percentIdentity">Percent identity</option>
|
|
|
33 <option value="bitScore">Bit score</option>
|
|
|
34 </param>
|
|
|
35 <param name="summarize" type="boolean" truevalue="-s" falsevalue="" label="Summarize" help="Sum results by taxonomy ID and list totals in the first column instead of query IDs"/>
|
|
|
36 </inputs>
|
|
|
37
|
|
|
38 <outputs>
|
|
|
39 <data format="tabular" name="output"/>
|
|
|
40 </outputs>
|
|
|
41
|
|
|
42 <help>
|
|
|
43
|
|
|
44 **What it does**
|
|
|
45
|
|
|
46 Creates taxonomic classifications for each query in BLAST results by finding the lowest common ancestor (or by picking randomly, if specified) of "best" hits (as determined by e-value factor above). Requires a KronaTools_ installation (2.0 or higher) with taxonomy downloaded.
|
|
|
47
|
|
|
48 .. _KronaTools: http://krona.sourceforge.net
|
|
|
49
|
|
|
50 **Input**
|
|
|
51
|
|
|
52 Tabular BLAST results containing gi numbers. If comment lines are present, queries with no hits will be assigned to taxon -1 (consistent with MEGAN).
|
|
|
53
|
|
|
54 Example input (BLAST)::
|
|
|
55
|
|
|
56 1 2 3 4 5 6 7 8 9 10 11 12
|
|
|
57 read1 gi|211853080|gb|EQ846228.1| 97.33 150 1 3 6 154 9518999 9519146 5e-64 252
|
|
|
58 read2 gi|211853080|gb|EQ846228.1| 96.69 151 0 5 6 154 9496147 9496294 2e-62 246
|
|
|
59
|
|
|
60 **Output**
|
|
|
61
|
|
|
62 The output is a list of query IDs (or counts, if **summarize** is specified) followed by NCBI taxonomy IDs and scores. This can be used to create a Krona chart by specifying "Taxonomy ID list" as the **input type**. If the output is summarized, the **summarized** option must also be specified when creating the chart.
|
|
|
63
|
|
|
64 Example output (default)::
|
|
|
65
|
|
|
66 1 (query ID) 2 (tax ID) 3 (score)
|
|
|
67 read1 9606 94.34
|
|
|
68 read2 9616 32.17
|
|
|
69
|
|
|
70 Example output (summarized)::
|
|
|
71
|
|
|
72 1 (count) 2 (tax ID) 3 (score)
|
|
|
73 243 9606 94.34
|
|
|
74 17 9616 32.17
|
|
|
75
|
|
|
76 </help>
|
|
|
77 <requirements>
|
|
|
78 <requirement type="binary">ktClassifyBLAST</requirement>
|
|
|
79 </requirements>
|
|
|
80 </tool>
|
