|
2
|
1 <?xml version="1.0" ?>
|
|
|
2 <tool id="qiime_feature-classifier_extract-reads" name="qiime feature-classifier extract-reads" version="2019.4">
|
|
|
3 <description> - Extract reads from reference</description>
|
|
|
4 <requirements>
|
|
|
5 <requirement type="package" version="2019.4">qiime2</requirement>
|
|
|
6 </requirements>
|
|
|
7 <command><![CDATA[
|
|
4
|
8
|
|
|
9
|
|
|
10 #if str( $input_sequences.selector ) == 'history'
|
|
|
11 #set $seq = $input_sequences.i_sequences
|
|
|
12
|
|
|
13 #else:
|
|
|
14 qiime tools import --type 'FeatureData[Sequence]' --input-path '$input_sequences.i_sequences.fields.path' --output-path database.qza;
|
|
|
15 #set $seq = 'database.qza'
|
|
|
16 #end if
|
|
|
17
|
|
|
18
|
|
2
|
19 qiime feature-classifier extract-reads
|
|
|
20
|
|
4
|
21 --i-sequences=$seq
|
|
|
22
|
|
2
|
23 --p-f-primer="$pfprimer"
|
|
|
24 --p-r-primer="$prprimer"
|
|
|
25
|
|
4
|
26 #if str($ptrunclen):
|
|
2
|
27 --p-trunc-len=$ptrunclen
|
|
|
28 #end if
|
|
|
29
|
|
4
|
30 #if str($ptrimleft):
|
|
2
|
31 --p-trim-left=$ptrimleft
|
|
|
32 #end if
|
|
|
33
|
|
4
|
34 #if str($pidentity):
|
|
2
|
35 --p-identity=$pidentity
|
|
|
36 #end if
|
|
|
37
|
|
4
|
38 #if str($pminlength):
|
|
2
|
39 --p-min-length=$pminlength
|
|
|
40 #end if
|
|
|
41
|
|
4
|
42 #if str($pmaxlength):
|
|
2
|
43 --p-max-length=$pmaxlength
|
|
|
44 #end if
|
|
|
45
|
|
|
46 --o-reads=oreads
|
|
|
47 ;
|
|
|
48 cp oreads.qza $oreads
|
|
|
49 ]]></command>
|
|
|
50 <inputs>
|
|
4
|
51 <conditional name="input_sequences" optional="False">
|
|
|
52 <param name="selector" type="select" label="--i-sequences: FeatureData[Sequence] [required] - Reference sequences to query">
|
|
|
53 <option value="cached">Public databases</option>
|
|
|
54 <option value="history">Databases from your history</option>
|
|
|
55 </param>
|
|
|
56 <when value="cached">
|
|
|
57 <param name="i_sequences" label="--i-sequences: FeatureData[Sequence] [required] - Reference sequences" type="select" optional="True">
|
|
|
58 <options from_data_table="qiime_rep_set" />
|
|
|
59 </param>
|
|
|
60 </when>
|
|
|
61 <when value="history">
|
|
|
62 <param name="i_sequences" type="data" format="qza,no_unzip.zip" label="Reference databases" optional="True" />
|
|
|
63 </when>
|
|
|
64 </conditional>
|
|
|
65
|
|
2
|
66 <param label="--p-f-primer: TEXT forward primer sequence [required]" name="pfprimer" optional="False" type="text"/>
|
|
|
67 <param label="--p-r-primer: TEXT reverse primer sequence [required]" name="prprimer" optional="False" type="text"/>
|
|
|
68 <param label="--p-trunc-len: INTEGER read is cut to trunc-len if trunc-len is positive. Applied before trim-left. [default: 0]" name="ptrunclen" optional="True" type="integer" value="0"/>
|
|
|
69 <param label="--p-trim-left: INTEGER trim-left nucleotides are removed from the 5' end if trim-left is positive. Applied after trunc-len. [default: 0]" name="ptrimleft" optional="True" type="integer" value="0"/>
|
|
|
70 <param label="--p-identity: NUMBER minimum combined primer match identity threshold. [default: 0.8]" name="pidentity" optional="True" type="float" value="0.8"/>
|
|
|
71 <param label="--p-min-length: INTEGER Minimum amplicon length. Shorter amplicons are Range(0, None) discarded. Applied after trimming and truncation, so be aware that trimming may impact sequence retention. Set to zero to disable min length filtering. [default: 50]" name="pminlength" optional="True" type="integer" min="0" value="50"/>
|
|
|
72 <param label="--p-max-length: INTEGER Maximum amplicon length. Longer amplicons are Range(0, None) discarded. Applied before trimming and truncation, so plan accordingly. Set to zero (default) to disable max length filtering. [default: 0]" name="pmaxlength" optional="True" type="integer" min="0" value="0"/>
|
|
|
73 </inputs>
|
|
|
74 <outputs>
|
|
|
75 <data format="qza" label="${tool.name} on ${on_string}: reads.qza" name="oreads"/>
|
|
|
76 </outputs>
|
|
|
77 <help><![CDATA[
|
|
|
78 Extract reads from reference
|
|
|
79 ############################
|
|
|
80
|
|
|
81 Extract sequencing-like reads from a reference database.
|
|
|
82
|
|
|
83 Parameters
|
|
|
84 ----------
|
|
|
85 sequences : FeatureData[Sequence]
|
|
|
86 \
|
|
|
87
|
|
|
88 f_primer : Str
|
|
|
89 forward primer sequence
|
|
|
90 r_primer : Str
|
|
|
91 reverse primer sequence
|
|
|
92 trunc_len : Int, optional
|
|
|
93 read is cut to trunc_len if trunc_len is positive. Applied before
|
|
|
94 trim_left.
|
|
|
95 trim_left : Int, optional
|
|
|
96 trim_left nucleotides are removed from the 5' end if trim_left is
|
|
|
97 positive. Applied after trunc_len.
|
|
|
98 identity : Float, optional
|
|
|
99 minimum combined primer match identity threshold.
|
|
|
100 min_length : Int % Range(0, None), optional
|
|
|
101 Minimum amplicon length. Shorter amplicons are discarded. Applied after
|
|
|
102 trimming and truncation, so be aware that trimming may impact sequence
|
|
|
103 retention. Set to zero to disable min length filtering.
|
|
|
104 max_length : Int % Range(0, None), optional
|
|
|
105 Maximum amplicon length. Longer amplicons are discarded. Applied before
|
|
|
106 trimming and truncation, so plan accordingly. Set to zero (default) to
|
|
|
107 disable max length filtering.
|
|
|
108
|
|
|
109 Returns
|
|
|
110 -------
|
|
|
111 reads : FeatureData[Sequence]
|
|
|
112 ]]></help>
|
|
|
113 <macros>
|
|
|
114 <import>qiime_citation.xml</import>
|
|
|
115 </macros>
|
|
|
116 <expand macro="qiime_citation"/>
|
|
|
117 </tool>
|
