Mercurial > repos > p.lucas > krona_with_conda
comparison classify_blast.xml @ 1:2d0679d0d8f7 draft
Uploaded classify blast
| author | p.lucas |
|---|---|
| date | Thu, 12 Jan 2023 07:59:14 +0000 |
| parents | |
| children | 0396eaee2bd2 |
comparison
equal
deleted
inserted
replaced
| 0:a0cb48398130 | 1:2d0679d0d8f7 |
|---|---|
| 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> |