Mercurial > repos > rnateam > sortmerna

comparison sortmerna.xml @ 10:10b84b577117 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/sortmerna commit 7dca92baf6a91f10bced28b883f14d3f4d7630c5
author rnateam
date Sat, 21 Oct 2023 09:07:21 +0000
parents eb35257d2e29
children
comparison
equal deleted inserted replaced
9:eb35257d2e29 10:10b84b577117
1 <tool id="bg_sortmerna" name="Filter with SortMeRNA" version="@VERSION@.6"> 1 <tool id="bg_sortmerna" name="Filter with SortMeRNA" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="22.05">
2 <description>of ribosomal RNAs in metatranscriptomic data</description> 2 <description>of ribosomal RNAs in metatranscriptomic data</description>
3 <macros> 3 <macros>
4 <token name="@VERSION@">2.1b</token> 4 <import>macros.xml</import>
5 <xml name="db_prep">
6 <param name="seed_length" type="integer" min="0" max="100" value="18" label="Seed length for database indexing" help="(-L)"/>
7 <param name="max_pos" type="integer" min="0" max="100000" value="10000" label="Maximum number of positions to store for each k-mer for database indexing" help="With 0, all positions are stored (--max_pos)"/>
8 </xml>
9 <xml name="output_alignments">
10 <param name="print_all_reads" type="boolean" checked="false" truevalue="--print_all_reads" falsevalue="" label="Output null alignment strings for non-aligned reads"/>
11 <conditional name="blast">
12 <param name="blast_output" type="select" label="Output BLAST report?">
13 <option value="True">Yes</option>
14 <option value="False" selected="True">No</option>
15 </param>
16 <when value="True">
17 <param name="blast_format" type="select" label="BLAST-like format?">
18 <option value="0">pairwise (--blast '0')</option>
19 <option value="1">tabular BLAST -m 8 format (--blast '1')</option>
20 <option value="1 cigar">tabular + column for CIGAR (--blast '1 cigar')</option>
21 <option value="1 cigar qcov">tabular + columns for CIGAR and query coverage (--blast '1 cigar qcov')</option>
22 <option value="1 cigar qcov qstrand">tabular + columns for CIGAR, query coverage and strand (--blast '1 cigar qcov qstrand')</option>
23 </param>
24 </when>
25 <when value="False"/>
26 </conditional>
27 </xml>
28 <token name="@ALIGNMENTS@">
29 $report.print_all_reads
30 --sam
31 --SQ
32 #if $report.blast.blast_output == 'True'
33 --blast '$report.blast.blast_format'
34 #end if
35 </token>
36 </macros> 5 </macros>
37 <requirements> 6 <expand macro="requirements" />
38 <requirement type="package" version="@VERSION@">sortmerna</requirement>
39 <requirement type="package" version="1.5">samtools</requirement>
40 </requirements>
41 <stdio> 7 <stdio>
42 <regex match="This program builds a Burst trie on an input rRNA database" 8 <regex match="This program builds a Burst trie on an input rRNA database"
43 source="both" 9 source="both"
44 level="fatal" 10 level="fatal"
45 description="Buildtrie program failed to execute." /> 11 description="Buildtrie program failed to execute." />
57 sortmerna --version 2>&1|grep 'SortMeRNA version' 23 sortmerna --version 2>&1|grep 'SortMeRNA version'
58 ]]> 24 ]]>
59 </version_command> 25 </version_command>
60 <command> 26 <command>
61 <![CDATA[ 27 <![CDATA[
62 #set $ref = '' 28
63 #set $sep='' 29 @PREPROCESSING@
64 #if str( $databases_type.databases_selector ) == 'history' 30
65 #for $db in $databases_type.database_name 31 mkdir -p './aligned' &&
66 #set $ref += $sep + str($db) + ',' + $os.path.splitext($os.path.basename(str($db)))[0] 32 mkdir -p './kvdb_folder' &&
67 #set $sep = ':' 33 mkdir -p './readb_folder' &&
34 mkdir -p './idx_folder' &&
35
36 sortmerna
37 #for $reference in $ref
38 --ref '$reference'
68 #end for 39 #end for
69 #else if str( $databases_type.databases_selector ) == 'cached_to_index' 40 -L '$databases_type.seed_length'
70 #for $db in $databases_type.input_databases.fields.path.split(",") 41 --max_pos '$databases_type.max_pos'
71 #set $ref += $sep + $db + ',' + $os.path.splitext($db)[0] + '-reindexed'
72 #set $sep = ':'
73 #end for
74 #else
75 #for $db in $databases_type.input_databases.fields.path.split(",")
76 #set $ref += $sep + $db + ',' + $os.path.splitext($db)[0]
77 #set $sep = ':'
78 #end for
79 #end if
80
81 #if str( $databases_type.databases_selector ) != 'cached'
82 indexdb_rna
83 --ref '$ref'
84 -L '$databases_type.seed_length'
85 --max_pos '$databases_type.max_pos'
86 &&
87 #end if
88
89 #if str( $sequencing_type.sequencing_type_selector ) == 'paired'
90 merge-paired-reads.sh
91 '$sequencing_type.forward_reads'
92 '$sequencing_type.reverse_reads'
93 merged-reads
94 &&
95 #end if
96
97 sortmerna
98 --ref '$ref'
99 --aligned 'aligned' 42 --aligned 'aligned'
100 #if str( $sequencing_type.sequencing_type_selector ) == 'paired' 43 #if $sequencing_type.sequencing_type_selector == 'paired'
101 --reads 'merged-reads' 44 --reads '$sequencing_type.forward_reads'
45 --reads '$sequencing_type.reverse_reads'
46 $sequencing_type.paired_type
47 #elif $sequencing_type.sequencing_type_selector == 'interleaved'
48 --reads '$sequencing_type.reads'
49 --paired
50 $sequencing_type.paired_type
51 #else if $sequencing_type.sequencing_type_selector == 'paired_collection'
52 --reads '${$sequencing_type.reads.forward}'
53 --reads '${$sequencing_type.reads.reverse}'
102 $sequencing_type.paired_type 54 $sequencing_type.paired_type
103 #else 55 #else
104 --reads '$sequencing_type.reads' 56 --reads '$sequencing_type.reads'
105 #end if 57 #end if
106 $strand_search 58 $strand_search
107 $log
108 $aligned_fastx.aligned_fastx_selector 59 $aligned_fastx.aligned_fastx_selector
109 #if $aligned_fastx.aligned_fastx_selector == '--fastx' 60 #if $aligned_fastx.aligned_fastx_selector == '--fastx'
110 #if $aligned_fastx.other 61 #if $aligned_fastx.other
111 --other 'unaligned' 62 --other 'unaligned'
112 #end if 63 #end if
113 #end if 64 #end if
114 #if $report.report_type == 'best' 65 #if $report.report_type == 'number_alignments'
115 @ALIGNMENTS@ 66 @ALIGNMENTS@
116 #if $report.otu.otu_map == 'True' 67 @OTU_PICKING@
117 --otu_map 68
118 --id '$report.otu.id' 69 #if $report.report_num_alignments.output_alignments == 'all'
119 --coverage '$report.otu.coverage' 70 --num_alignments '0'
120 $report.otu.de_novo_otu 71 #else
72 --num_alignments $report.report_num_alignments.num_alignments
121 #end if 73 #end if
122 #if $report.report_best.report_best_type == '1' 74 $report.no_best
123 --best 1 75 #elif $report.report_type == 'min_lis'
124 --min_lis '$report.report_best.report_best_min_lis'
125 #else
126 --best '$report.report_best.report_best_value'
127 --min_lis '$report.report_best.report_best_min_lis'
128 #end if
129 #elif $report.report_type == 'num_alignments'
130 @ALIGNMENTS@ 76 @ALIGNMENTS@
131 #if $report.report_num_alignments.report_num_alignments_type == 'other_value' 77 @OTU_PICKING@
132 --num_alignments '$report.report_num_alignments.report_num_alignments_value' 78 --min_lis $report.min_lis
133 #else
134 --num_alignments '$report.report_num_alignments.report_num_alignments_type'
135 #end if
136 #end if 79 #end if
137 -e '$e_value' 80 -e '$e_value'
138 --match '$match' 81 --match '$match'
139 --mismatch '$mismatch' 82 --mismatch '$mismatch'
140 --gap_open '$gap_open' 83 --gap_open '$gap_open'
141 --gap_ext '$gap_ext' 84 --gap_ext '$gap_ext'
142 -N '$ambiguous_letter' 85 --threads \${GALAXY_SLOTS:-1}
143 -a \${GALAXY_SLOTS:-1} 86 -m \${GALAXY_MEMORY_MB:-8192}
87 --kvdb 'kvdb_folder'
88 --idx-dir 'idx_folder'
89 --readb 'readb_folder'
90 #if $aligned_fastx.aligned_fastx_selector == '--fastx' and str($sequencing_type.sequencing_type_selector) == 'paired'
91 --out2
92 #end if
144 #if $report.report_type != 'None' 93 #if $report.report_type != 'None'
145 && 94 && samtools sort aligned.sam -@ "\${GALAXY_SLOTS:-4}" -T tmp -O bam -o '$output_bam'
146 samtools view -@ "\${GALAXY_SLOTS:-4}" -u aligned.sam | samtools sort -@ "\${GALAXY_SLOTS:-4}" -T tmp -O bam -o '$output_bam'
147 #end if 95 #end if
148 96 && ls -lah
149 #if $aligned_fastx.aligned_fastx_selector == '--fastx' and str($sequencing_type.sequencing_type_selector) == 'paired'
150 #if str($sequencing_type.paired_type) != ''
151 &&
152 unmerge-paired-reads.sh
153 aligned.fast*
154 '$aligned_forward'
155 '$aligned_reverse'
156 #if $aligned_fastx.other
157 &&
158 unmerge-paired-reads.sh
159 unaligned.fast*
160 '$unaligned_forward'
161 '$unaligned_reverse'
162 #end if
163 #else
164 &&
165 mv aligned.fast* '$aligned_paired'
166 #if $aligned_fastx.other
167 &&
168 mv unaligned.fast* '$unaligned_paired'
169 #end if
170 #end if
171 #end if
172 ]]> 97 ]]>
173 </command> 98 </command>
174 <inputs> 99 <inputs>
175 <conditional name="sequencing_type"> 100 <conditional name="sequencing_type">
176 <param name="sequencing_type_selector" type="select" label="Sequencing type"> 101 <param name="sequencing_type_selector" type="select" label="Sequencing type">
177 <option value="not_paired">Reads are not paired</option> 102 <option value="not_paired">Single-end reads</option>
178 <option value="paired">Reads are paired</option> 103 <option value="paired">Paired-end reads</option>
104 <option value="interleaved">Interleaved paired-end reads</option>
105 <option value="paired_collection">Paired collection</option>
179 </param> 106 </param>
180 <when value="not_paired"> 107 <when value="not_paired">
181 <param argument="--reads" type="data" format="fasta,fastq" label="Querying sequences"/> 108 <expand macro="reads_macro"/>
182 </when> 109 </when>
183 <when value="paired"> 110 <when value="paired">
184 <param name="forward_reads" type="data" format="fasta,fastq" label="Forward reads"/> 111 <param name="forward_reads" type="data" format="fasta,fastq,fastq.gz,fasta.gz" label="Forward reads"/>
185 <param name="reverse_reads" type="data" format="fasta,fastq" label="Reverse reads"/> 112 <param name="reverse_reads" type="data" format="fasta,fastq,fastq.gz,fasta.gz" label="Reverse reads"/>
186 <param name="paired_type" type="select" display="radio" label="If one of the paired-end reads aligns and the other one does not"> 113 <expand macro="paired_type_macro"/>
187 <option value="">Leave the reads split between aligned and rejected files</option> 114 </when>
188 <option value="--paired_in">Output both reads to aligned file (--paired_in)</option> 115 <when value="interleaved">
189 <option value="--paired_out">Output both reads to rejected file (--paired_out)</option> 116 <expand macro="reads_macro"/>
190 </param> 117 <expand macro="paired_type_macro"/>
118 </when>
119 <when value="paired_collection">
120 <param name="reads" format="fasta,fastq,fasta.gz,fastq.gz" type="data_collection" collection_type="paired" label="Select a paired collection" help="See help section for an explanation of dataset collections"/>
121 <expand macro="paired_type_macro"/>
191 </when> 122 </when>
192 </conditional> 123 </conditional>
193 <param name="strand_search" type="select" label="Which strands to search"> 124 <param name="strand_search" type="select" label="Which strands to search">
194 <option value="">Search both strands</option> 125 <option value="">Search both strands</option>
195 <option value="-F">Search only the forward strand (-F)</option> 126 <option value="-F">Search only the forward strand (-F)</option>
196 <option value="-R">Search only the reverse-complementary strand (-R)</option> 127 <option value="-R">Search only the reverse-complementary strand (-R)</option>
197 </param> 128 </param>
198 <conditional name="databases_type"> 129 <conditional name="databases_type">
199 <param name="databases_selector" type="select" label="Databases to query" 130 <param name="databases_selector" type="select" label="Databases to query"
200 help="Public rRNA databases provided with SortMeRNA have been indexed. On the contrary, personal databases must be indexed each time SortMeRNA is launched. Please be patient, this may take some time depending on the size of the given database."> 131 help="Public rRNA databases provided with SortMeRNA have been indexed. On the contrary, personal databases must be indexed each time SortMeRNA is launched. Please be patient, this may take some time depending on the size of the given database.">
201 <option value="cached" selected="true">Public pre-indexed ribosomal databases</option> 132 <option value="cached" selected="true">Public ribosomal databases</option>
202 <option value="cached_to_index">Public ribosomal databases to index with non default parameters</option>
203 <option value="history">Databases from your history</option> 133 <option value="history">Databases from your history</option>
204 </param> 134 </param>
205 <when value="cached"> 135 <when value="cached">
206 <param name="input_databases" label="rRNA databases" type="select" display="checkboxes" multiple="true"> 136 <param name="input_databases" label="rRNA databases" type="select" optional="false" multiple="true">
207 <options from_data_table="rRNA_databases" /> 137 <options from_data_table="rRNA_databases">
208 <validator type="no_options" message="Select at least one database"/> 138 <column name="name" index="1"/>
209 </param> 139 <column name="value" index="2"/>
210 </when> 140 </options>
211 <when value="cached_to_index"> 141 <validator type="no_options" message="No options available. Contact your Galaxy administrator."/>
212 <param name="input_databases" label="rRNA databases" type="select" display="checkboxes" multiple="true">
213 <options from_data_table="rRNA_databases" />
214 <validator type="no_options" message="Select at least one database"/>
215 </param> 142 </param>
216 <expand macro="db_prep"/> 143 <expand macro="db_prep"/>
217 </when> 144 </when>
218 <when value="history"> 145 <when value="history">
219 <param name="database_name" type="data" format="fasta" multiple="true" label="rRNA databases" help="Your databases will be indexed first, which may take up to several minutes."/> 146 <param name="database_name" type="data" format="fasta,fasta.gz" multiple="true" label="rRNA databases" help="Your databases will be indexed first, which may take up to several minutes."/>
220 <expand macro="db_prep"/> 147 <expand macro="db_prep"/>
221 </when> 148 </when>
222 </conditional> 149 </conditional>
223 <!-- Outputs --> 150 <!-- Outputs -->
224 <conditional name="aligned_fastx"> 151 <conditional name="aligned_fastx">
229 <when value="--fastx"> 156 <when value="--fastx">
230 <param argument="--other" type="boolean" truevalue="True" falsevalue="False" label="Include rejected reads file?"/> 157 <param argument="--other" type="boolean" truevalue="True" falsevalue="False" label="Include rejected reads file?"/>
231 </when> 158 </when>
232 <when value=""/> 159 <when value=""/>
233 </conditional> 160 </conditional>
234 <param argument="--log" type="boolean" checked="false" truevalue="--log" falsevalue="" label="Generate statistics file" help="Generates statistics for the rRNA content of reads, as well as rRNA subunit distribution">
235 </param>
236 <conditional name="report"> 161 <conditional name="report">
237 <param name="report_type" type="select" label="Alignment report"> 162 <param name="report_type" type="select" label="Alignment report">
238 <option value="None">Do not report alignments</option> 163 <option value="None">Do not report alignments</option>
239 <option value="best">Report best alignments per read reaching E-value</option> 164 <option value="number_alignments">Report N number of alignments reaching E-value</option>
240 <option value="num_alignments">Report first alignements per read reaching E-value</option> 165 <option value="min_lis">Report alignments that have the LIS of at least N seeds long reaching E-value</option>
241 </param> 166 </param>
242 <when value="None"/> 167 <when value="None"/>
243 <when value="best"> 168 <when value="number_alignments">
169 <conditional name="report_num_alignments">
170 <param argument="output_alignments" type="select" label="Number of output alignments" help="Report all alignments can be time consuming; this option is not suggested for high similarity rRNA databases.">
171 <option value="all">All alignments reaching the E-value threshold are reported</option>
172 <option value="custom">Custom number of alignments</option>
173 </param>
174 <when value="all" />
175 <when value="custom">
176 <param argument="--num_alignments" type="integer" min="1" max="100" value="1" label="Number of alignments to be reported"/>
177 </when>
178 </conditional>
179 <param argument="--no-best" type="boolean" truevalue="--no-best" falsevalue="" checked="false" label="Disable best alignments search" help="The 'best' alignment is the highest
180 scoring alignment out of All alignments of a read, and the read can potentially be aligned (reaching E-value threshold) to multiple reference sequences. By default the
181 program searches for best alignments i.e. performs an exhaustive search over all references. Using '-no-best' will make the program to search just the first N alignments." />
244 <expand macro="output_alignments"/> 182 <expand macro="output_alignments"/>
245 <conditional name="otu"> 183 <expand macro="otu_picking"/>
246 <param name="otu_map" type="select" label="Pick OTUs?"> 184 </when>
247 <option value="True">Yes</option> 185 <when value="min_lis">
248 <option value="False" selected="true">No</option> 186 <param argument="--min_lis" type="integer" min="0" value="" label="Minimum Longest Increasing Subsequence (LIS)" help="It is computed using seeds, which are k-mers common to
249 </param> 187 the read and the reference sequence. Sorted sequences of such seeds are used to filter the candidate references prior performing the Smith-Waterman alignment." />
250 <when value="True">
251 <param argument="--id" type="float" value="0.97" label="Percentage id similarity threshold" help="The alignment must still pass the E-value threshold" />
252 <param argument="--coverage" type="float" value="0.97" label="Percentage query coverage threshold" help="The alignment must still pass the E-value threshold" />
253 <param name="de_novo_otu" type="boolean" truevalue="--de_novo_otu" falsevalue="" label="FASTA/FASTQ file for reads matching database below percentage id" help="--de_novo_otu" />
254 </when>
255 <when value="False"/>
256 </conditional>
257 <conditional name="report_best">
258 <param argument="report_best_type" type="select" label="Number of searched alignments" help="Only the best alignment is reported (--best)">
259 <option value="1" selected="true">Only one high-candidate reference sequence is searched for alignments (fast). The high-candidate sequences are determined heuristically using a LIS of seed matches)</option>
260 <option value="other_value">A custom number of reference sequences are searched for alignments (speed decrease for high value)</option>
261 </param>
262 <when value="1">
263 <param name="report_best_min_lis" type="integer" min="0" max="100" value="2" label="Number of longest LIS an alignement needs to be searched" help="The alignements having the first INT longest LIS. LIS stands for Longest Increasing Subsequence, it is computed using seeds' positions to expand hits into longer matches prior to Smith-Waterman alignment. (--min_lis)"/>
264 </when>
265 <when value="other_value">
266 <param name="report_best_value" type="integer" min="2" max="100" value="2" label="Number of alignments to be made" help="Only the best one is reported. The computation speed decrease with high value"/>
267 <param name="report_best_min_lis" type="integer" min="0" max="100" value="2" label="Number of longest LIS an alignement needs to be searched" help="The alignements having the first INT longest LIS. LIS stands for Longest Increasing Subsequence, it is computed using seeds' positions to expand hits into longer matches prior to Smith-Waterman alignment. (--min_lis)"/>
268 </when>
269 </conditional>
270 </when>
271 <when value="num_alignments">
272 <expand macro="output_alignments"/> 188 <expand macro="output_alignments"/>
273 <conditional name="report_num_alignments"> 189 <expand macro="otu_picking"/>
274 <param name="report_num_alignments_type" type="select" label="Number of output alignments" help="(--num_alignments)">
275 <option value="0">All alignments reaching the E-value threshold are reported (very slow, this option is not suggested for high similarity rRNA databases)</option>
276 <option value="1" selected="true">The first alignment passing E-value threshold are reported (very fast, best choice if only filtering is needed)</option>
277 <option value="other_value">A custom number of alignments are made and reported (speed decrease for high value)</option>
278 </param>
279 <when value="0" />
280 <when value="1" />
281 <when value="other_value">
282 <param name="report_num_alignments_value" type="integer" min="0" max="100" value="1" label="Number of alignments to be made and reported" help=""/>
283 </when>
284 </conditional>
285 </when> 190 </when>
286 </conditional> 191 </conditional>
287 <param name="e_value" type="float" min="0" max="10" value="1" label="E-value threshold" help="(-e)"/>
288 <param argument="--match" type="integer" min="0" max="10" value="2" label="SW score for a match"/> 192 <param argument="--match" type="integer" min="0" max="10" value="2" label="SW score for a match"/>
289 <param argument="--mismatch" type="integer" min="-10" max="0" value="-3" label="SW penalty for a mismatch"/> 193 <param argument="--mismatch" type="integer" min="-10" max="0" value="-3" label="SW penalty for a mismatch"/>
290 <param name="gap_open" type="integer" min="0" max="10" value="5" label="SW penalty for introducing a gap" help="(--gap_open)"/> 194 <param name="gap_open" type="integer" min="0" max="10" value="5" label="SW penalty for introducing a gap" help="(--gap_open)"/>
291 <param name="gap_ext" type="integer" min="0" max="10" value="2" label="SW penalty for extending a gap" help="(--gap_ext)"/> 195 <param name="gap_ext" type="integer" min="0" max="10" value="2" label="SW penalty for extending a gap" help="(--gap_ext)"/>
292 <param name="ambiguous_letter" type="integer" min="-10" max="0" value="-3" label="SW penalty for ambiguous letters (N's)" help="(-N)"/> 196 <param name="e_value" argument="-e" type="float" min="0" max="10" value="1" label="E-value threshold" help="Defines the 'statistical significance' of a local alignment. Exponentially
197 correllates with the Minimal Alignment score. Higher E-values (100, 1000, ...) cause more reads to pass the alignment threshold."/>
198 <param argument="-N" type="boolean" truevalue="True" falsevalue="False" checked="false" label="SW penalty for ambiguous letters (N's)" help="Scored as --mismatch" />
293 </inputs> 199 </inputs>
294 <outputs> 200 <outputs>
295 <data name="output_fastx" format_source="reads" from_work_dir="aligned.dat" label="${tool.name} on ${on_string}: Aligned reads"> 201 <data name="aligned" format_source="reads" from_work_dir="aligned.f*" label="${tool.name} on ${on_string}: Aligned reads">
296 <filter>aligned_fastx['aligned_fastx_selector'] != '' and sequencing_type['sequencing_type_selector'] != 'paired'</filter> 202 <filter>aligned_fastx['aligned_fastx_selector'] != ''</filter>
297 </data> 203 <filter>sequencing_type['sequencing_type_selector'] != 'paired'</filter>
298 <data name="aligned_paired" format_source="forward_reads" label="${tool.name} on ${on_string}: Aligned reads"> 204 </data>
299 <filter>aligned_fastx['aligned_fastx_selector'] != '' and sequencing_type['sequencing_type_selector'] == 'paired' and sequencing_type['paired_type'] == ''</filter> 205 <data name="aligned_forward" format_source="forward_reads" from_work_dir="aligned_fwd*" label="${tool.name} on ${on_string}: Aligned forward reads">
300 </data> 206 <filter>aligned_fastx['aligned_fastx_selector'] != '' and sequencing_type['sequencing_type_selector'] == 'paired'</filter>
301 <data name="aligned_forward" format_source="forward_reads" label="${tool.name} on ${on_string}: Aligned forward reads"> 207 </data>
302 <filter>aligned_fastx['aligned_fastx_selector'] != '' and sequencing_type['sequencing_type_selector'] == 'paired' and sequencing_type['paired_type'] != ''</filter> 208 <data name="aligned_reverse" format_source="reverse_reads" from_work_dir="aligned_rev*" label="${tool.name} on ${on_string}: Aligned reverse reads">
303 </data> 209 <filter>aligned_fastx['aligned_fastx_selector'] != '' and sequencing_type['sequencing_type_selector'] == 'paired'</filter>
304 <data name="aligned_reverse" format_source="reverse_reads" label="${tool.name} on ${on_string}: Aligned reverse reads"> 210 </data>
305 <filter>aligned_fastx['aligned_fastx_selector'] != '' and sequencing_type['sequencing_type_selector'] == 'paired' and sequencing_type['paired_type'] != ''</filter> 211 <data name="aligned_forward_singleton" format_source="forward_reads" from_work_dir="aligned_singleton_fwd*" label="${tool.name} on ${on_string}: Aligned forward singleton reads">
306 </data> 212 <filter>aligned_fastx['aligned_fastx_selector'] != ''</filter>
307 <data name="output_other" format_source="reads" from_work_dir="unaligned.dat" label="${tool.name} on ${on_string}: Unaligned reads"> 213 <filter>sequencing_type['sequencing_type_selector'] == 'paired'</filter>
308 <filter>aligned_fastx['aligned_fastx_selector'] != '' and aligned_fastx['other'] == True and sequencing_type['sequencing_type_selector'] != 'paired'</filter> 214 <filter>sequencing_type['paired_type'] == '--sout'</filter>
309 </data> 215 </data>
310 <data name="unaligned_paired" format_source="forward_reads" label="${tool.name} on ${on_string}: Unaligned reads"> 216 <data name="aligned_reverse_singleton" format_source="reverse_reads" from_work_dir="aligned_singleton_rev*" label="${tool.name} on ${on_string}: Aligned reverse singleton reads">
311 <filter>aligned_fastx['aligned_fastx_selector'] != '' and aligned_fastx['other'] == True and sequencing_type['sequencing_type_selector'] == 'paired' and sequencing_type['paired_type'] == ''</filter> 217 <filter>aligned_fastx['aligned_fastx_selector'] != ''</filter>
312 </data> 218 <filter>sequencing_type['sequencing_type_selector'] == 'paired'</filter>
313 <data name="unaligned_forward" format_source="forward_reads" label="${tool.name} on ${on_string}: Unaligned forward reads"> 219 <filter>sequencing_type['paired_type'] == '--sout'</filter>
314 <filter>aligned_fastx['aligned_fastx_selector'] != '' and aligned_fastx['other'] == True and sequencing_type['sequencing_type_selector'] == 'paired' and sequencing_type['paired_type'] != ''</filter> 220 </data>
315 </data> 221 <data name="unaligned" format_source="reads" from_work_dir="aligned.f*" label="${tool.name} on ${on_string}: Unaligned reads">
316 <data name="unaligned_reverse" format_source="reverse_reads" label="${tool.name} on ${on_string}: Unaligned reverse reads"> 222 <filter>aligned_fastx['aligned_fastx_selector'] != ''</filter>
317 <filter>aligned_fastx['aligned_fastx_selector'] != '' and aligned_fastx['other'] == True and sequencing_type['sequencing_type_selector'] == 'paired' and sequencing_type['paired_type'] != ''</filter> 223 <filter>aligned_fastx['other'] == True</filter>
318 </data> 224 <filter>sequencing_type['sequencing_type_selector'] != 'paired'</filter>
319 <data name="output_bam" format="bam" label="${tool.name} on ${on_string}: Alignments"> 225 </data>
226 <data name="unaligned_forward" format_source="forward_reads" from_work_dir="unaligned_fwd*" label="${tool.name} on ${on_string}: Unaligned forward reads">
227 <filter>aligned_fastx['aligned_fastx_selector'] != ''</filter>
228 <filter>aligned_fastx['other'] == True</filter>
229 <filter>sequencing_type['sequencing_type_selector'] == 'paired'</filter>
230 </data>
231 <data name="unaligned_reverse" format_source="reverse_reads" from_work_dir="unaligned_rev*" label="${tool.name} on ${on_string}: Unaligned reverse reads">
232 <filter>aligned_fastx['aligned_fastx_selector'] != ''</filter>
233 <filter>aligned_fastx['other'] == True</filter>
234 <filter>sequencing_type['sequencing_type_selector'] == 'paired'</filter>
235 </data>
236 <data name="unaligned_forward_singleton" format_source="forward_reads" from_work_dir="unaligned_fwd*" label="${tool.name} on ${on_string}: Unaligned forward singleton reads">
237 <filter>aligned_fastx['aligned_fastx_selector'] != ''</filter>
238 <filter>aligned_fastx['other'] == True</filter>
239 <filter>sequencing_type['sequencing_type_selector'] == 'paired'</filter>
240 <filter>sequencing_type['paired_type'] == '--sout'</filter>
241 </data>
242 <data name="unaligned_reverse_singleton" format_source="reverse_reads" from_work_dir="unaligned_rev*" label="${tool.name} on ${on_string}: Unaligned reverse singleton reads">
243 <filter>aligned_fastx['aligned_fastx_selector'] != ''</filter>
244 <filter>aligned_fastx['other'] == True</filter>
245 <filter>sequencing_type['sequencing_type_selector'] == 'paired'</filter>
246 <filter>sequencing_type['paired_type'] == '--sout'</filter>
247 </data>
248 <data name="output_bam" format="bam" label="${tool.name} on ${on_string}: Alignments (BAM)">
320 <filter>report['report_type'] != 'None'</filter> 249 <filter>report['report_type'] != 'None'</filter>
321 </data> 250 </data>
322 <data name="output_blast" format="tabular" from_work_dir="aligned.blast" label="${tool.name} on ${on_string}: BLAST report"> 251 <data name="output_blast" format="tabular" from_work_dir="aligned.blast" label="${tool.name} on ${on_string}: BLAST report">
323 <filter>report['report_type'] != 'None' and report['blast']['blast_output'] == 'True'</filter> 252 <filter>report['report_type'] != 'None' and report['blast']['blast_output'] == 'True'</filter>
324 <change_format> 253 <change_format>
325 <when input="aligned_blast.aligned_blast_format" value="0" format="txt" /> 254 <when input="aligned_blast.aligned_blast_format" value="0" format="txt" />
326 </change_format> 255 </change_format>
327 </data> 256 </data>
328 <data name="output_biom" format="txt" from_work_dir="aligned_otus.txt" label="${tool.name} on ${on_string}: OTU map"> 257 <data name="output_biom" format="txt" from_work_dir="aligned/otu_map.txt" label="${tool.name} on ${on_string}: OTU map">
329 <filter>report['report_type'] != 'None' and report['report_type'] == 'best' and report['otu']['otu_map'] == 'True'</filter> 258 <filter>report['report_type'] != 'None' and report['otu']['otu_map'] == 'True'</filter>
330 </data> 259 </data>
331 <data name="output_de_novo" format_source="reads" from_work_dir="aligned_denovo.dat" label="${tool.name} on ${on_string}: De novo reads matching database"> 260
332 <filter>report['report_type'] != 'None' and report['report_type'] == 'best' and report['otu']['otu_map'] == 'True' and report['otu']['de_novo_otu'] == True</filter> 261 <data name="output_de_novo" format_source="reads" from_work_dir="aligned_denovo*" label="${tool.name} on ${on_string}: De novo reads matching database">
333 </data> 262 <filter>report['report_type'] != 'None' and report['otu']['otu_map'] == 'True' and report['otu']['de_novo_otu'] == True</filter>
334 <data name="output_log" format="txt" label="${tool.name} on ${on_string}: Log" from_work_dir="aligned.log">
335 <filter>log == True</filter>
336 </data> 263 </data>
337 </outputs> 264 </outputs>
338 <tests> 265 <tests>
339 <test> 266 <test expect_num_outputs="4">
340 <conditional name="sequencing_type"> 267 <conditional name="sequencing_type">
341 <param name="sequencing_type_selector" value="not_paired" /> 268 <param name="sequencing_type_selector" value="not_paired" />
342 <param name="reads" value="read_small.fastq" /> 269 <param name="reads" value="read_small.fastq" />
343 </conditional> 270 </conditional>
344 <param name="strand_search" value="" /> 271 <param name="strand_search" value="" />
348 </conditional> 275 </conditional>
349 <conditional name="aligned_fastx"> 276 <conditional name="aligned_fastx">
350 <param name="aligned_fastx_selector" value="--fastx" /> 277 <param name="aligned_fastx_selector" value="--fastx" />
351 <param name="other" value="True" /> 278 <param name="other" value="True" />
352 </conditional> 279 </conditional>
353 <param name="log" value="False" /> 280 <conditional name="report">
354 <conditional name="report"> 281 <param name="report_type" value="number_alignments" />
355 <param name="report_type" value="num_alignments" /> 282 <conditional name="report_num_alignments">
356 <param name="report_best_type" value="1" /> 283 <param name="output_alignments" value="custom"/>
357 <param name="print_all_reads" value="False" /> 284 <param name="num_alignments" value="1"/>
285 </conditional>
358 <conditional name="blast"> 286 <conditional name="blast">
359 <param name="blast_output" value="True"/> 287 <param name="blast_output" value="True"/>
360 <param name="blast_format" value="1 cigar qcov"/> 288 <param name="blast_format" value="1 cigar qcov"/>
361 </conditional> 289 </conditional>
362 <conditional name="otu"> 290 </conditional>
363 <param name="otu_map" value="False"/> 291 <param name="e_value" value="1"/>
364 </conditional> 292 <param name="match" value="2"/>
365 </conditional> 293 <param name="mismatch" value="-3" />
366 <param name="e_value" value="1"/> 294 <param name="gap_open" value="5"/>
367 <param name="match" value="2"/> 295 <param name="gap_ext" value="2"/>
368 <param name="mismatch" value="-3" /> 296 <param name="N" value="True"/>
369 <param name="gap_open" value="5"/> 297 <output name="aligned" file="test1_aligned.fastq" />
370 <param name="gap_ext" value="2"/> 298 <output name="unaligned" file="test1_unaligned.fastq" />
371 <param name="ambiguous_letter" value="-3"/> 299 <output name="output_bam" file="test1_bam.bam" lines_diff="6" ftype="bam"/>
372 <output name="output_fastx" file="test1_aligned.fastq" />
373 <output name="output_other" file="test1_other.fastq" />
374 <output name="output_bam" file="test1_bam.bam" compare="sim_size" delta="200" />
375 <output name="output_blast" file="test1_blast.tabular"/> 300 <output name="output_blast" file="test1_blast.tabular"/>
376 </test> 301 </test>
377 <test> 302 <!-- test cached reference data -->
303 <test expect_num_outputs="4">
378 <conditional name="sequencing_type"> 304 <conditional name="sequencing_type">
379 <param name="sequencing_type_selector" value="not_paired" /> 305 <param name="sequencing_type_selector" value="not_paired" />
380 <param name="reads" value="read_small.fasta" /> 306 <param name="reads" value="read_small.fastq" />
307 </conditional>
308 <param name="strand_search" value="" />
309 <conditional name="databases_type">
310 <param name="databases_selector" value="cached" />
311 <param name="input_databases" value="ref_small" />
312 </conditional>
313 <conditional name="aligned_fastx">
314 <param name="aligned_fastx_selector" value="--fastx" />
315 <param name="other" value="True" />
316 </conditional>
317 <conditional name="report">
318 <param name="report_type" value="number_alignments" />
319 <conditional name="report_num_alignments">
320 <param name="output_alignments" value="custom"/>
321 <param name="num_alignments" value="1"/>
322 </conditional>
323 <conditional name="blast">
324 <param name="blast_output" value="True"/>
325 <param name="blast_format" value="1 cigar qcov"/>
326 </conditional>
327 </conditional>
328 <param name="e_value" value="1"/>
329 <param name="match" value="2"/>
330 <param name="mismatch" value="-3" />
331 <param name="gap_open" value="5"/>
332 <param name="gap_ext" value="2"/>
333 <param name="N" value="True"/>
334 <output name="aligned" file="test1_aligned.fastq" />
335 <output name="unaligned" file="test1_unaligned.fastq" />
336 <output name="output_bam" file="test1_bam.bam" lines_diff="6" ftype="bam"/>
337 <output name="output_blast" file="test1_blast.tabular"/>
338 </test>
339 <!-- testing 2 reference files from history .. apart from this same as previous test -->
340 <test expect_num_outputs="4">
341 <conditional name="sequencing_type">
342 <param name="sequencing_type_selector" value="not_paired" />
343 <param name="reads" value="read_small.fastq" />
344 </conditional>
345 <param name="strand_search" value="" />
346 <conditional name="databases_type">
347 <param name="databases_selector" value="history" />
348 <param name="database_name" value="ref_small.fasta,ref_small_copy.fasta" />
349 </conditional>
350 <conditional name="aligned_fastx">
351 <param name="aligned_fastx_selector" value="--fastx" />
352 <param name="other" value="True" />
353 </conditional>
354 <conditional name="report">
355 <param name="report_type" value="number_alignments" />
356 <conditional name="report_num_alignments">
357 <param name="output_alignments" value="custom"/>
358 <param name="num_alignments" value="1"/>
359 </conditional>
360 <conditional name="blast">
361 <param name="blast_output" value="True"/>
362 <param name="blast_format" value="1 cigar qcov"/>
363 </conditional>
364 </conditional>
365 <param name="e_value" value="1"/>
366 <param name="match" value="2"/>
367 <param name="mismatch" value="-3" />
368 <param name="gap_open" value="5"/>
369 <param name="gap_ext" value="2"/>
370 <param name="N" value="True"/>
371 <output name="aligned" file="test1_aligned.fastq" />
372 <output name="unaligned" file="test1_unaligned.fastq" />
373 <output name="output_bam" file="test1_bam.bam" lines_diff="6" ftype="bam"/>
374 <output name="output_blast" file="test1_blast.tabular"/>
375 <assert_command>
376 <has_text text="--ref" n="2"/>
377 </assert_command>
378 <assert_stdout>
379 <has_text text="Processing reference [2] out of total [2] references"/>
380 </assert_stdout>
381 </test>
382 <!-- testing 2 cached reference -->
383 <test expect_num_outputs="4">
384 <conditional name="sequencing_type">
385 <param name="sequencing_type_selector" value="not_paired" />
386 <param name="reads" value="read_small.fastq" />
387 </conditional>
388 <param name="strand_search" value="" />
389 <conditional name="databases_type">
390 <param name="databases_selector" value="cached" />
391 <param name="input_databases" value="ref_small,ref_small_copy" />
392 </conditional>
393 <conditional name="aligned_fastx">
394 <param name="aligned_fastx_selector" value="--fastx" />
395 <param name="other" value="True" />
396 </conditional>
397 <conditional name="report">
398 <param name="report_type" value="number_alignments" />
399 <conditional name="report_num_alignments">
400 <param name="output_alignments" value="custom"/>
401 <param name="num_alignments" value="1"/>
402 </conditional>
403 <conditional name="blast">
404 <param name="blast_output" value="True"/>
405 <param name="blast_format" value="1 cigar qcov"/>
406 </conditional>
407 </conditional>
408 <param name="e_value" value="1"/>
409 <param name="match" value="2"/>
410 <param name="mismatch" value="-3" />
411 <param name="gap_open" value="5"/>
412 <param name="gap_ext" value="2"/>
413 <param name="N" value="True"/>
414 <output name="aligned" file="test1_aligned.fastq" />
415 <output name="unaligned" file="test1_unaligned.fastq" />
416 <output name="output_bam" file="test1_bam.bam" lines_diff="6" ftype="bam"/>
417 <output name="output_blast" file="test1_blast.tabular"/>
418 <assert_command>
419 <has_text text="--ref" n="2"/>
420 </assert_command>
421 <assert_stdout>
422 <has_text text="Processing reference [2] out of total [2] references"/>
423 </assert_stdout>
424 </test>
425 <test expect_num_outputs="1">
426 <conditional name="sequencing_type">
427 <param name="sequencing_type_selector" value="not_paired" />
428 <param name="reads" value="read_small.fastq" />
381 </conditional> 429 </conditional>
382 <param name="strand_search" value="" /> 430 <param name="strand_search" value="" />
383 <conditional name="databases_type"> 431 <conditional name="databases_type">
384 <param name="databases_selector" value="history" /> 432 <param name="databases_selector" value="history" />
385 <param name="database_name" value="ref_small.fasta" /> 433 <param name="database_name" value="ref_small.fasta" />
386 </conditional> 434 </conditional>
387 <conditional name="aligned_fastx"> 435 <conditional name="aligned_fastx">
388 <param name="aligned_fastx_selector" value="--fastx" /> 436 <param name="aligned_fastx_selector" value="--fastx" />
389 <param name="other" value="False" /> 437 <param name="other" value="False" />
390 </conditional> 438 </conditional>
391 <param name="log" value="True" />
392 <conditional name="report"> 439 <conditional name="report">
393 <param name="report_type" value="None" /> 440 <param name="report_type" value="None" />
394 </conditional> 441 </conditional>
395 <param name="e_value" value="1"/> 442 <param name="e_value" value="1"/>
396 <param name="match" value="2"/> 443 <param name="match" value="2"/>
397 <param name="mismatch" value="-3" /> 444 <param name="mismatch" value="-3" />
398 <param name="gap_open" value="5"/> 445 <param name="gap_open" value="5"/>
399 <param name="gap_ext" value="2"/> 446 <param name="gap_ext" value="2"/>
400 <param name="ambiguous_letter" value="-3"/> 447 <param name="N" value="True"/>
401 <output name="output_fastx" file="test2_aligned.fasta" /> 448 <output name="aligned" file="test2_aligned.fasta" />
402 <output name="output_log" file="test2_log.txt" compare="sim_size" /> 449 </test>
403 </test> 450 <test expect_num_outputs="5">
404 <test>
405 <conditional name="sequencing_type"> 451 <conditional name="sequencing_type">
406 <param name="sequencing_type_selector" value="paired" /> 452 <param name="sequencing_type_selector" value="paired" />
407 <param name="forward_reads" value="forward_reads.fastq" /> 453 <param name="forward_reads" value="forward_reads.fastq" />
408 <param name="reverse_reads" value="reverse_reads.fastq" /> 454 <param name="reverse_reads" value="reverse_reads.fastq" />
409 <param name="paired_type" value="--paired_in"/> 455 <param name="paired_type" value="--paired_in"/>
417 </conditional> 463 </conditional>
418 <conditional name="aligned_fastx"> 464 <conditional name="aligned_fastx">
419 <param name="aligned_fastx_selector" value="--fastx" /> 465 <param name="aligned_fastx_selector" value="--fastx" />
420 <param name="other" value="True" /> 466 <param name="other" value="True" />
421 </conditional> 467 </conditional>
422 <param name="log" value="False" /> 468 <conditional name="report">
423 <conditional name="report"> 469 <param name="report_type" value="number_alignments" />
424 <param name="report_type" value="best" /> 470 <conditional name="report_num_alignments">
425 <param name="report_num_alignments_type" value="1"/> 471 <param name="output_alignments" value="all"/>
472 </conditional>
426 <param name="print_all_reads" value="False" /> 473 <param name="print_all_reads" value="False" />
427 <conditional name="blast"> 474 <conditional name="blast">
428 <param name="blast_output" value="False"/> 475 <param name="blast_output" value="False"/>
429 </conditional> 476 </conditional>
430 <conditional name="otu"> 477 <conditional name="otu">
434 <param name="e_value" value="1"/> 481 <param name="e_value" value="1"/>
435 <param name="match" value="2"/> 482 <param name="match" value="2"/>
436 <param name="mismatch" value="-3" /> 483 <param name="mismatch" value="-3" />
437 <param name="gap_open" value="5"/> 484 <param name="gap_open" value="5"/>
438 <param name="gap_ext" value="2"/> 485 <param name="gap_ext" value="2"/>
439 <param name="ambiguous_letter" value="-3"/> 486 <param name="N" value="True"/>
440 <output name="aligned_forward" file="test3_aligned_forward.fastq" /> 487 <output name="aligned_forward" file="test3_aligned_forward.fastq" />
441 <output name="aligned_reverse" file="test3_aligned_reverse.fastq" /> 488 <output name="aligned_reverse" file="test3_aligned_reverse.fastq" />
442 <output name="unaligned_forward" file="test3_unaligned_forward.fastq" /> 489 <output name="unaligned_forward" file="test3_unaligned_forward.fastq" />
443 <output name="unaligned_reverse" file="test3_unaligned_reverse.fastq" /> 490 <output name="unaligned_reverse" file="test3_unaligned_reverse.fastq" />
444 <output name="output_bam" file="test3_bam.bam" compare="sim_size" delta="200" /> 491 <output name="output_bam" file="test3_bam.bam" lines_diff="6" ftype="bam"/>
445 </test> 492 </test>
446 <test> 493 <test expect_num_outputs="3">
447 <conditional name="sequencing_type"> 494 <conditional name="sequencing_type">
448 <param name="sequencing_type_selector" value="not_paired" /> 495 <param name="sequencing_type_selector" value="not_paired" />
449 <param name="reads" value="test4_input.fasta" /> 496 <param name="reads" value="test4_input.fasta" />
450 </conditional> 497 </conditional>
451 <param name="strand_search" value="" /> 498 <param name="strand_search" value="" />
456 <param name="max_pos" value="100000"/> 503 <param name="max_pos" value="100000"/>
457 </conditional> 504 </conditional>
458 <conditional name="aligned_fastx"> 505 <conditional name="aligned_fastx">
459 <param name="aligned_fastx_selector" value="" /> 506 <param name="aligned_fastx_selector" value="" />
460 </conditional> 507 </conditional>
461 <param name="log" value="False" /> 508 <conditional name="report">
462 <conditional name="report"> 509 <param name="report_type" value="number_alignments" />
463 <param name="report_type" value="best" /> 510 <conditional name="report_num_alignments">
464 <param name="report_num_alignments_type" value="1"/> 511 <param name="output_alignments" value="all"/>
512 </conditional>
465 <param name="print_all_reads" value="False" /> 513 <param name="print_all_reads" value="False" />
466 <conditional name="blast"> 514 <conditional name="blast">
467 <param name="blast_output" value="False"/> 515 <param name="blast_output" value="False"/>
468 </conditional> 516 </conditional>
469 <conditional name="otu"> 517 <conditional name="otu">
476 <param name="e_value" value="1"/> 524 <param name="e_value" value="1"/>
477 <param name="match" value="2"/> 525 <param name="match" value="2"/>
478 <param name="mismatch" value="-3" /> 526 <param name="mismatch" value="-3" />
479 <param name="gap_open" value="5"/> 527 <param name="gap_open" value="5"/>
480 <param name="gap_ext" value="2"/> 528 <param name="gap_ext" value="2"/>
481 <param name="ambiguous_letter" value="-3"/> 529 <param name="N" value="True"/>
482 <output name="output_bam" file="test4_bam.bam" compare="sim_size" delta="200" /> 530 <output name="output_bam" file="test4_bam.bam" compare="sim_size" delta="200" />
483 <output name="output_biom" file="test4_biom.txt"/> 531 <output name="output_biom" file="test4_biom.txt"/>
484 <output name="output_de_novo" file="test4_de_novo.fasta"/> 532 <output name="output_de_novo" file="test4_de_novo.fasta"/>
485 </test> 533 </test>
486 <test> 534 <test expect_num_outputs="5">
487 <conditional name="sequencing_type"> 535 <conditional name="sequencing_type">
488 <param name="sequencing_type_selector" value="paired" /> 536 <param name="sequencing_type_selector" value="paired" />
489 <param name="forward_reads" value="forward_reads.fasta" /> 537 <param name="forward_reads" value="forward_reads.fasta" />
490 <param name="reverse_reads" value="reverse_reads.fasta" /> 538 <param name="reverse_reads" value="reverse_reads.fasta" />
491 <param name="paired_type" value=""/> 539 <param name="paired_type" value=""/>
499 </conditional> 547 </conditional>
500 <conditional name="aligned_fastx"> 548 <conditional name="aligned_fastx">
501 <param name="aligned_fastx_selector" value="--fastx" /> 549 <param name="aligned_fastx_selector" value="--fastx" />
502 <param name="other" value="True" /> 550 <param name="other" value="True" />
503 </conditional> 551 </conditional>
504 <param name="log" value="False" /> 552 <conditional name="report">
505 <conditional name="report"> 553 <param name="report_type" value="number_alignments" />
506 <param name="report_type" value="best" />
507 <param name="report_num_alignments_type" value="1"/>
508 <param name="print_all_reads" value="False" /> 554 <param name="print_all_reads" value="False" />
509 <conditional name="blast"> 555 <conditional name="blast">
510 <param name="blast_output" value="False"/> 556 <param name="blast_output" value="False"/>
511 </conditional> 557 </conditional>
512 <conditional name="otu"> 558 <conditional name="otu">
516 <param name="e_value" value="1"/> 562 <param name="e_value" value="1"/>
517 <param name="match" value="2"/> 563 <param name="match" value="2"/>
518 <param name="mismatch" value="-3" /> 564 <param name="mismatch" value="-3" />
519 <param name="gap_open" value="5"/> 565 <param name="gap_open" value="5"/>
520 <param name="gap_ext" value="2"/> 566 <param name="gap_ext" value="2"/>
521 <param name="ambiguous_letter" value="-3"/> 567 <param name="N" value="True"/>
522 <output name="aligned_paired" file="test5_aligned.fasta" /> 568 <output name="aligned_forward" file="test5_aligned_forward.fasta" />
523 <output name="unaligned_paired" file="test5_unaligned.fasta" /> 569 <output name="aligned_reverse" file="test5_aligned_reverse.fasta" />
524 </test> 570 <output name="output_bam" file="test5_bam.bam" lines_diff="6" ftype="bam"/>
525 <test> 571 </test>
572 <test expect_num_outputs="5">
526 <conditional name="sequencing_type"> 573 <conditional name="sequencing_type">
527 <param name="sequencing_type_selector" value="paired" /> 574 <param name="sequencing_type_selector" value="paired" />
528 <param name="forward_reads" value="forward_reads.fasta" /> 575 <param name="forward_reads" value="forward_reads.fasta" />
529 <param name="reverse_reads" value="reverse_reads.fasta" /> 576 <param name="reverse_reads" value="reverse_reads.fasta" />
530 <param name="paired_type" value="--paired_out"/> 577 <param name="paired_type" value="--paired_out"/>
538 </conditional> 585 </conditional>
539 <conditional name="aligned_fastx"> 586 <conditional name="aligned_fastx">
540 <param name="aligned_fastx_selector" value="--fastx" /> 587 <param name="aligned_fastx_selector" value="--fastx" />
541 <param name="other" value="True" /> 588 <param name="other" value="True" />
542 </conditional> 589 </conditional>
543 <param name="log" value="False" /> 590 <conditional name="report">
544 <conditional name="report"> 591 <param name="report_type" value="number_alignments" />
545 <param name="report_type" value="best" /> 592 <conditional name="report_num_alignments">
546 <param name="report_num_alignments_type" value="1"/> 593 <param name="output_alignments" value="all"/>
594 </conditional>
547 <param name="print_all_reads" value="False" /> 595 <param name="print_all_reads" value="False" />
548 <conditional name="blast"> 596 <conditional name="blast">
549 <param name="blast_output" value="False"/> 597 <param name="blast_output" value="False"/>
550 </conditional> 598 </conditional>
551 <conditional name="otu"> 599 <conditional name="otu">
555 <param name="e_value" value="1"/> 603 <param name="e_value" value="1"/>
556 <param name="match" value="2"/> 604 <param name="match" value="2"/>
557 <param name="mismatch" value="-3" /> 605 <param name="mismatch" value="-3" />
558 <param name="gap_open" value="5"/> 606 <param name="gap_open" value="5"/>
559 <param name="gap_ext" value="2"/> 607 <param name="gap_ext" value="2"/>
560 <param name="ambiguous_letter" value="-3"/> 608 <param name="N" value="True"/>
561 <output name="aligned_forward" file="test6_aligned_forward.fasta" /> 609 <output name="aligned_forward" file="test6_aligned_forward.fasta" />
562 <output name="aligned_reverse" file="test6_aligned_reverse.fasta" /> 610 <output name="aligned_reverse" file="test6_aligned_reverse.fasta" />
563 <output name="unaligned_forward" file="test6_unaligned_forward.fasta" /> 611 <output name="unaligned_forward" file="test6_unaligned_forward.fasta" />
564 <output name="unaligned_reverse" file="test6_unaligned_reverse.fasta" /> 612 <output name="unaligned_reverse" file="test6_unaligned_reverse.fasta" />
613 <output name="output_bam" file="test6_bam.bam" lines_diff="6" ftype="bam"/>
614 </test>
615 <test expect_num_outputs="1">
616 <conditional name="sequencing_type">
617 <param name="sequencing_type_selector" value="paired" />
618 <param name="forward_reads" value="forward_reads.fasta" />
619 <param name="reverse_reads" value="reverse_reads.fasta" />
620 </conditional>
621 <conditional name="databases_type">
622 <param name="databases_selector" value="history" />
623 <param name="database_name" value="ref_small.fasta" />
624 </conditional>
625 <conditional name="aligned_fastx">
626 <param name="aligned_fastx_selector" value="" />
627 </conditional>
628 <conditional name="report">
629 <param name="report_type" value="min_lis" />
630 <param name="min_lis" value="1"/>
631 </conditional>
632 <output name="output_bam" file="test7_bam.bam" lines_diff="6" ftype="bam"/>
633 </test>
634 <test expect_num_outputs="9">
635 <conditional name="sequencing_type">
636 <param name="sequencing_type_selector" value="paired" />
637 <param name="forward_reads" value="forward_reads.fastq" />
638 <param name="reverse_reads" value="reverse_reads.fastq" />
639 <param name="paired_type" value="--sout"/>
640 </conditional>
641 <param name="strand_search" value="" />
642 <conditional name="databases_type">
643 <param name="databases_selector" value="history" />
644 <param name="database_name" value="ref_small.fasta" />
645 <param name="seed_length" value="18" />
646 <param name="max_pos" value="100000"/>
647 </conditional>
648 <conditional name="aligned_fastx">
649 <param name="aligned_fastx_selector" value="--fastx" />
650 <param name="other" value="True" />
651 </conditional>
652 <conditional name="report">
653 <param name="report_type" value="number_alignments" />
654 <conditional name="report_num_alignments">
655 <param name="output_alignments" value="all"/>
656 </conditional>
657 <param name="print_all_reads" value="False" />
658 <conditional name="blast">
659 <param name="blast_output" value="False"/>
660 </conditional>
661 <conditional name="otu">
662 <param name="otu_map" value="False"/>
663 </conditional>
664 </conditional>
665 <param name="e_value" value="1"/>
666 <param name="match" value="2"/>
667 <param name="mismatch" value="-3" />
668 <param name="gap_open" value="5"/>
669 <param name="gap_ext" value="2"/>
670 <param name="N" value="True"/>
671 <output name="aligned_forward_singleton" file="test8_aligned_forward_singleton.fastq" />
672 <output name="aligned_reverse_singleton" file="test8_aligned_reverse_singleton.fastq" />
673 <output name="output_bam" file="test8_bam.bam" lines_diff="6" ftype="bam"/>
674 </test>
675 <test expect_num_outputs="1">
676 <conditional name="sequencing_type">
677 <param name="sequencing_type_selector" value="interleaved" />
678 <param name="reads" value="interlaced_reads.fastq.gz" />
679 </conditional>
680 <conditional name="databases_type">
681 <param name="databases_selector" value="history" />
682 <param name="database_name" value="ref_small.fasta" />
683 </conditional>
684 <output name="aligned" file="test9_aligned.fastq.gz" compare="sim_size"/>
685 </test>
686 <test expect_num_outputs="1">
687 <conditional name="sequencing_type">
688 <param name="sequencing_type_selector" value="paired_collection" />
689 <param name="reads">
690 <collection type="paired">
691 <element name="forward" value="forward_reads.fastq" />
692 <element name="reverse" value="reverse_reads.fastq" />
693 </collection>
694 </param>
695 </conditional>
696 <conditional name="databases_type">
697 <param name="databases_selector" value="history" />
698 <param name="database_name" value="ref_small.fasta" />
699 </conditional>
700 <conditional name="aligned_fastx">
701 <param name="aligned_fastx_selector" value="" />
702 </conditional>
703 <conditional name="report">
704 <param name="report_type" value="number_alignments" />
705 </conditional>
706 <output name="output_bam" file="test10_bam.bam" lines_diff="8" ftype="bam"/>
565 </test> 707 </test>
566 </tests> 708 </tests>
567 <help> 709 <help>
568 <![CDATA[ 710 <![CDATA[
569 **What it does** 711 **What it does**
632 databases are indexed when added, but they can be re-indexed with non-default indexing 774 databases are indexed when added, but they can be re-indexed with non-default indexing
633 parameters. The indexing may take some time depending on the size of the given database. 775 parameters. The indexing may take some time depending on the size of the given database.
634 776
635 ]]> 777 ]]>
636 </help> 778 </help>
637 779 <expand macro="citations" />
638 <citations>
639 <citation type="doi">10.1093/bioinformatics/bts611</citation>
640 <citation type="doi">10.1093/nar/gks1219</citation>
641 <citation type="doi">10.1093/nar/gks1005</citation>
642 <citation type="doi">10.1093/bioinformatics/btq461</citation>
643 <citation type="doi">10.1038/nbt.2198</citation>
644 </citations>
645 </tool> 780 </tool>