|
0
|
1 <?xml version="1.0" ?>
|
|
|
2 <tool id="qiime_feature-classifier_classify-consensus-vsearch" name="qiime feature-classifier classify-consensus-vsearch" version="2019.4">
|
|
|
3 <description> - VSEARCH consensus taxonomy classifier</description>
|
|
|
4 <requirements>
|
|
|
5 <requirement type="package" version="2019.4">qiime2</requirement>
|
|
|
6 </requirements>
|
|
|
7 <command><![CDATA[
|
|
|
8 qiime feature-classifier classify-consensus-vsearch
|
|
|
9
|
|
|
10 --i-query=$iquery
|
|
|
11 --i-reference-reads=$ireferencereads
|
|
|
12
|
|
|
13
|
|
|
14 #if str( $id_to_taxonomy_fp.selector ) == 'history'
|
|
|
15 #set $tax = $id_to_taxonomy_fp.taxonomy_fp
|
|
|
16 --i-reference-taxonomy '$tax'
|
|
|
17 #else:
|
|
|
18 #set $tax = $id_to_taxonomy_fp.taxonomy_fp.fields.path
|
|
|
19 --i-reference-taxonomy '$tax'
|
|
|
20 #end if
|
|
|
21
|
|
|
22
|
|
4
|
23
|
|
6
|
24 #if str($pmaxaccepts):
|
|
0
|
25 --p-maxaccepts=$pmaxaccepts
|
|
|
26 #end if
|
|
|
27
|
|
6
|
28 #if str($ppercidentity):
|
|
0
|
29 --p-perc-identity=$ppercidentity
|
|
|
30 #end if
|
|
|
31
|
|
6
|
32 #if str($pquerycov):
|
|
0
|
33 --p-query-cov=$pquerycov
|
|
|
34 #end if
|
|
|
35
|
|
|
36 #if str($pstrand) != 'None':
|
|
|
37 --p-strand=$pstrand
|
|
|
38 #end if
|
|
|
39
|
|
6
|
40 #if str($pminconsensus):
|
|
0
|
41 --p-min-consensus=$pminconsensus
|
|
|
42 #end if
|
|
|
43
|
|
4
|
44
|
|
|
45
|
|
|
46
|
|
|
47
|
|
|
48 #if '__sq__' in str($punassignablelabel):
|
|
|
49 #set $punassignablelabel_temp = $punassignablelabel.replace('__sq__', "'")
|
|
|
50 #set $punassignablelabel = $punassignablelabel_temp
|
|
|
51 #end if
|
|
|
52
|
|
0
|
53 #if str($punassignablelabel):
|
|
|
54 --p-unassignable-label="$punassignablelabel"
|
|
|
55 #end if
|
|
|
56
|
|
4
|
57
|
|
|
58
|
|
|
59
|
|
0
|
60 #set $pthreads = '${GALAXY_SLOTS:-4}'
|
|
|
61 #if str($pthreads):
|
|
|
62 --p-threads="$pthreads"
|
|
|
63 #end if
|
|
|
64
|
|
|
65 --o-classification=oclassification
|
|
|
66 ;
|
|
|
67 cp oclassification.qza $oclassification
|
|
|
68 ]]></command>
|
|
|
69 <inputs>
|
|
|
70 <param format="qza,no_unzip.zip" label="--i-query: ARTIFACT FeatureData[Sequence] Sequences to classify taxonomically. [required]" name="iquery" optional="False" type="data"/>
|
|
|
71 <param format="qza,no_unzip.zip" label="--i-reference-reads: ARTIFACT FeatureData[Sequence] reference sequences. [required]" name="ireferencereads" optional="False" type="data"/>
|
|
|
72
|
|
|
73 <conditional name="id_to_taxonomy_fp" optional="True">
|
|
|
74 <param name="selector" type="select" label="Reference taxonomy to query">
|
|
|
75 <option value="cached">Public databases</option>
|
|
|
76 <option value="history">Databases from your history</option>
|
|
|
77 </param>
|
|
|
78 <when value="cached">
|
|
|
79 <param argument="--taxonomy_fp" label="Reference taxonomy" type="select" optional="True">
|
|
|
80 <options from_data_table="qiime_taxonomy" />
|
|
|
81 </param>
|
|
|
82 </when>
|
|
|
83 <when value="history">
|
|
|
84 <param argument="--taxonomy_fp" type="data" format="qza,no_unzip.zip" label="Reference databases" optional="True" />
|
|
|
85 </when>
|
|
|
86 </conditional>
|
|
|
87
|
|
|
88 <param label="--p-maxaccepts: INTEGER Maximum number of hits to keep for each query. Set Range(0, None) to 0 to keep all hits > perc-identity similarity. Must be in range [0, infinity]. [default: 10]" name="pmaxaccepts" optional="True" type="integer" min="0" value="10"/>
|
|
|
89 <param label="--p-perc-identity: PROPORTION Range(0.0, 1.0, inclusive_end=True) Reject match if percent identity to query is lower. Must be in range [0.0, 1.0]. [default: 0.8]" name="ppercidentity" optional="True" type="float" min="0" max="1" exclude_max="False" value="0.8"/>
|
|
|
90 <param label="--p-query-cov: PROPORTION Range(0.0, 1.0, inclusive_end=True) Reject match if query alignment coverage per high-scoring pair is lower. Must be in range [0.0, 1.0]. [default: 0.8]" name="pquerycov" optional="True" type="float" min="0" max="1" exclude_max="False" value="0.8"/>
|
|
|
91 <param label="--p-strand: Align against reference sequences in forward ('plus') or both directions ('both'). [default: 'both']" name="pstrand" optional="True" type="select">
|
|
|
92 <option value="None">Selection is Optional</option>
|
|
|
93 <option selected="True" value="both">both</option>
|
|
|
94 <option value="plus">plus</option>
|
|
|
95 </param>
|
|
|
96 <param label="--p-min-consensus: NUMBER Range(0.5, 1.0, inclusive_start=False, inclusive_end=True) Minimum fraction of assignments must match top hit to be accepted as consensus assignment. Must be in range (0.5, 1.0]. [default: 0.51]" name="pminconsensus" optional="True" type="float" min="0.5" max="1" exclude_min="True" exclude_max="False" value="0.51"/>
|
|
4
|
97 <param label="--p-unassignable-label: TEXT [default: 'Unassigned']" name="punassignablelabel" optional="True" type="text" value="Unassigned"/>
|
|
0
|
98 </inputs>
|
|
|
99 <outputs>
|
|
|
100 <data format="qza" label="${tool.name} on ${on_string}: classification.qza" name="oclassification"/>
|
|
|
101 </outputs>
|
|
|
102 <help><![CDATA[
|
|
|
103 VSEARCH consensus taxonomy classifier
|
|
|
104 #####################################
|
|
|
105
|
|
|
106 Assign taxonomy to query sequences using VSEARCH. Performs VSEARCH global
|
|
|
107 alignment between query and reference_reads, then assigns consensus
|
|
|
108 taxonomy to each query sequence from among maxaccepts top hits,
|
|
|
109 min_consensus of which share that taxonomic assignment. Unlike classify-
|
|
|
110 consensus-blast, this method searches the entire reference database before
|
|
|
111 choosing the top N hits, not the first N hits.
|
|
|
112
|
|
|
113 Parameters
|
|
|
114 ----------
|
|
|
115 query : FeatureData[Sequence]
|
|
|
116 Sequences to classify taxonomically.
|
|
|
117 reference_reads : FeatureData[Sequence]
|
|
|
118 reference sequences.
|
|
|
119 reference_taxonomy : FeatureData[Taxonomy]
|
|
|
120 reference taxonomy labels.
|
|
|
121 maxaccepts : Int % Range(0, None), optional
|
|
|
122 Maximum number of hits to keep for each query. Set to 0 to keep all
|
|
|
123 hits > perc_identity similarity. Must be in range [0, infinity].
|
|
|
124 perc_identity : Float % Range(0.0, 1.0, inclusive_end=True), optional
|
|
|
125 Reject match if percent identity to query is lower. Must be in range
|
|
|
126 [0.0, 1.0].
|
|
|
127 query_cov : Float % Range(0.0, 1.0, inclusive_end=True), optional
|
|
|
128 Reject match if query alignment coverage per high-scoring pair is
|
|
|
129 lower. Must be in range [0.0, 1.0].
|
|
|
130 strand : Str % Choices('both', 'plus'), optional
|
|
|
131 Align against reference sequences in forward ("plus") or both
|
|
|
132 directions ("both").
|
|
|
133 min_consensus : Float % Range(0.5, 1.0, inclusive_start=False, inclusive_end=True), optional
|
|
|
134 Minimum fraction of assignments must match top hit to be accepted as
|
|
|
135 consensus assignment. Must be in range (0.5, 1.0].
|
|
|
136 unassignable_label : Str, optional
|
|
|
137 \
|
|
|
138
|
|
|
139 Returns
|
|
|
140 -------
|
|
|
141 classification : FeatureData[Taxonomy]
|
|
|
142 The resulting taxonomy classifications.
|
|
|
143 ]]></help>
|
|
|
144 <macros>
|
|
|
145 <import>qiime_citation.xml</import>
|
|
|
146 </macros>
|
|
|
147 <expand macro="qiime_citation"/>
|
|
|
148 </tool>
|
