2
|
1 <tool id="kraken_combine" name="kraken_combine">
|
|
2 <description>extract unique and common sequences from two kraken runs</description>
|
|
3 <requirements>
|
|
4 <requirement type="package" version="2.7">python</requirement>
|
|
5 </requirements>
|
|
6 <command interpreter="python">kraken_combine.py
|
|
7 --in1="${dataFormat.input1}"
|
|
8 --set1-output-left=${set1OutputLeft}
|
|
9 --set2-output-left=${set2OutputLeft}
|
|
10 --unclassified-output-left=${unclassifiedOutputLeft}
|
|
11 --intersection-output-left=${intersectionOutputLeft}
|
|
12 --kraken-results1="${classificationSet1}"
|
|
13 --kraken-results2="${classificationSet2}"
|
|
14
|
|
15 #if $dataFormat.inputFormat == "pairedFASTQ" or $dataFormat.inputFormat == "pairedFASTA":
|
|
16 --in2="${dataFormat.input2}"
|
|
17 --set1-output-right=${set1OutputRight}
|
|
18 --set2-output-right=${set2OutputRight}
|
|
19 --unclassified-output-right=${unclassifiedOutputRight}
|
|
20 --intersection-output-right=${intersectionOutputRight}
|
|
21 #end if
|
|
22
|
|
23 #if $dataFormat.inputFormat == "singleFASTQ" or $dataFormat.inputFormat == "pairedFASTQ":
|
|
24 --input-format="fq"
|
|
25 #else:
|
|
26 --input-format="fa"
|
|
27 #end if
|
|
28
|
|
29
|
|
30
|
|
31 </command>
|
|
32 <inputs>
|
|
33 <param name="classificationSet1" type="data" format="tabular" label="Kraken classification output set1" />
|
|
34 <param name="classificationSet2" type="data" format="tabular" label="Kraken classification output set2" />
|
|
35
|
|
36 <conditional name="dataFormat">
|
|
37 <param name="inputFormat" type="select" label="Please select input file type and library type.">
|
|
38 <option value="singleFASTQ">Single-end FASTQ</option>
|
|
39 <option value="pairedFASTQ">Paired-end FASTQ</option>
|
|
40 <option value="singleFASTA">Single-end FASTA</option>
|
|
41 <option value="pairedFASTA">Paired-end FASTA</option>
|
|
42 </param>
|
|
43 <when value="singleFASTQ">
|
|
44 <param name="input1" type="data" format="fastq,fq.gz,fq" label="FASTQ file" />
|
|
45 </when>
|
|
46 <when value="pairedFASTQ">
|
|
47 <param name="input1" type="data" format="fastq,fq.gz,fq" label="Forward/Left FASTQ file" />
|
|
48 <param name="input2" type="data" format="fastq,fq.gz,fq" label="Reverse/Right FASTQ file" />
|
|
49 </when>
|
|
50 <when value="singleFASTA">
|
|
51 <param name="input1" type="data" format="fasta,fasta.gz,fasta.bz2,fa.gz" label="FASTA file" />
|
|
52 </when>
|
|
53 <when value="pairedFASTA">
|
|
54 <param name="input1" type="data" format="fasta,fasta.gz,fasta.bz2,fa.gz" label="Forward/Left FASTA file" />
|
|
55 <param name="input2" type="data" format="fasta,fasta.gz,fasta.bz2,fa.gz" label="Reverse/Right FASTA file" />
|
|
56 </when>
|
|
57 </conditional>
|
|
58 </inputs>
|
|
59
|
|
60 <outputs>
|
|
61 <data format="input1" name="set1OutputLeft" label="Taxonomy1-unique sequences (R1) of ${on_string}" />
|
|
62 <data format="input1" name="set2OutputLeft" label="Taxonomy2-unique sequences (R1) of ${on_string}" />
|
|
63 <data format="input1" name="intersectionOutputLeft" label="Common sequences (R1) of ${on_string}" />
|
|
64 <data format="input1" name="unclassifiedOutputLeft" label="Unclassified sequences (R1) of ${on_string}" />
|
|
65
|
|
66 <data format="input2" name="set1OutputRight" label="Taxonomy1-unique sequences (R2) of ${on_string}">
|
|
67 <filter>inputFormat == "pairedFASTQ" or inputFormat == "pairedFASTA"</filter>
|
|
68 </data>
|
|
69 <data format="input2" name="set2OutputRight" label="Taxonomy2-unique sequences (R2) of ${on_string}">
|
|
70 <filter>inputFormat == "pairedFASTQ" or inputFormat == "pairedFASTA"</filter>
|
|
71 </data>
|
|
72 <data format="input2" name="intersectionOutputRight" label="Common sequences (R2) of ${on_string}">
|
|
73 <filter>inputFormat == "pairedFASTQ" or inputFormat == "pairedFASTA"</filter>
|
|
74 </data>
|
|
75 <data format="input2" name="unclassifiedOutputRight" label="Unclassified sequences (R2) of ${on_string}">
|
|
76 <filter>inputFormat == "pairedFASTQ" or inputFormat == "pairedFASTA"</filter>
|
|
77 </data>
|
|
78 </outputs>
|
|
79
|
|
80 <help>
|
|
81 This tool compares the classification results on a dataset using two different kraken-databases.
|
|
82 </help>
|
|
83 </tool>
|