|
0
|
1 <?xml version="1.0"?>
|
|
|
2 <tool id="mitoflow" name="MITOflow" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="18.01">
|
|
|
3 <description>de novo assembler and annotation for fungi mitogenomes</description>
|
|
|
4 <macros>
|
|
|
5 <xml name="content" token_min="" token_max="">
|
|
|
6 <param name="genome_range_min" type="text" value="@MIN@" label="Set lower limit for the expected genome size range" help="(Genome Range)"/>
|
|
|
7 <param name="genome_range_max" type="text" value="@MAX@" label="Set upper limit for the expected genome size range" help="(Genome Range)"/>
|
|
|
8 </xml>
|
|
|
9 </macros>
|
|
|
10 <requirements>
|
|
|
11 <requirement type="package" version="4.2.0">novoplasty</requirement>
|
|
|
12 <requirement type="package" version="2.0.6">mitos</requirement>
|
|
|
13 <requirement type="package">zip</requirement>
|
|
|
14 </requirements>
|
|
|
15 <version_command>1.0</version_command>
|
|
|
16 <command detect_errors="exit_code"><![CDATA[
|
|
|
17 ## initialize
|
|
|
18 ## samples require extension .fasta or .fasta.gz
|
|
|
19 #if $reads_cond.reads_sel == 'combined'
|
|
|
20 ln -s '$reads_cond.combined_reads' 'combined.${reads_cond.combined_reads.datatype.file_ext}' &&
|
|
|
21 #else if $reads_cond.reads_sel == 'separate'
|
|
|
22 ln -s '$reads_cond.forward_reads' 'forward.${reads_cond.forward_reads.datatype.file_ext}' &&
|
|
|
23 ln -s '$reads_cond.reverse_reads' 'reverse.${reads_cond.reverse_reads.datatype.file_ext}' &&
|
|
|
24 #end if
|
|
|
25
|
|
|
26 ## run
|
|
|
27 NOVOPlasty.pl -c '$config'
|
|
|
28
|
|
|
29 ## postprocessing
|
|
|
30 && cp '$config' 'config.txt'
|
|
|
31 ]]></command>
|
|
|
32 <configfiles>
|
|
|
33 <configfile name="config"><![CDATA[
|
|
|
34 Project:
|
|
|
35 -----------------------
|
|
|
36 Project name = result
|
|
|
37 Type = ${type_cond.type_sel}
|
|
|
38 Genome Range = ${type_cond.genome_range_min}-${type_cond.genome_range_max}
|
|
|
39 K-mer = ${kmer}
|
|
|
40 Max memory =
|
|
|
41 Extended log = 1
|
|
|
42 Save assembled reads = ${save_assembled_reads}
|
|
|
43 Seed Input = ${seed_input}
|
|
|
44 Extend seed directly = ${extend_seed_directly}
|
|
|
45 Reference sequence = #if $reference then $reference else ''#
|
|
|
46 Variance detection = #if $platform_cond.platform_sel == 'illumina' and $platform_cond.mode_cond.mode_sel == 'variance' then 'yes' else ''#
|
|
|
47 Heteroplasmy =
|
|
|
48 MAF = #if $platform_cond.platform_sel == 'illumina' and $platform_cond.mode_cond.mode_sel == 'heteroplasmy' then $platform_cond.mode_cond.maf else ''#
|
|
|
49 HP exclude list =
|
|
|
50 Chloroplast sequence = #if $type_cond.type_sel == 'mito_plant' and $type_cond.chloroplast_sequence then $type_cond.chloroplast_sequence else ''#
|
|
|
51
|
|
|
52 Dataset 1:
|
|
|
53 -----------------------
|
|
|
54 Read Length = ${read_length}
|
|
|
55 Insert size = ${insert_size}
|
|
|
56 Platform = ${$platform_cond.platform_sel}
|
|
|
57 Single/Paired = ${single_paired}
|
|
|
58 Combined reads = #if $reads_cond.reads_sel == 'combined' then 'combined.'+$reads_cond.combined_reads.datatype.file_ext else ''#
|
|
|
59 Forward reads = #if $reads_cond.reads_sel == 'separate' then 'forward.'+$reads_cond.forward_reads.datatype.file_ext else ''#
|
|
|
60 Reverse reads = #if $reads_cond.reads_sel == 'separate' then 'reverse.'+$reads_cond.reverse_reads.datatype.file_ext else ''#
|
|
|
61
|
|
|
62 Optional:
|
|
|
63 -----------------------
|
|
|
64 Insert size auto = ${insert_size_auto}
|
|
|
65 Use Quality Scores = ${use_quality_scores}
|
|
|
66 ]]></configfile>
|
|
|
67 </configfiles>
|
|
|
68 <inputs>
|
|
|
69 <conditional name="reads_cond">
|
|
|
70 <param name="reads_sel" type="select" label="Select read file type" help="(Combined reads/Forward reads/Reverse reads)">
|
|
|
71 <option value="combined">Combined reads</option>
|
|
|
72 <option value="separate">Separate reads</option>
|
|
|
73 </param>
|
|
|
74 <when value="combined">
|
|
|
75 <param name="combined_reads" type="data" format="fasta,fasta.gz,fastq,fastq.gz" label="Select file with combined reads" help="(Combined reads)"/>
|
|
|
76 </when>
|
|
|
77 <when value="separate">
|
|
|
78 <param name="forward_reads" type="data" format="fasta,fasta.gz,fastq,fastq.gz" label="Select file with forward reads" help="(Forward reads)"/>
|
|
|
79 <param name="reverse_reads" type="data" format="fasta,fasta.gz,fastq,fastq.gz" label="Select file with reverse reads" help="(Reverse reads)"/>
|
|
|
80 </when>
|
|
|
81 </conditional>t
|
|
|
82 <!-- seed_input doesn't support fasta.gz, throws error INVALID SEED, PLEASE TRY AGAIN WITH A NEW ONE -->
|
|
|
83 <param name="seed_input" type="data" format="fasta" label="Select file with seed sequence" help="(Seed Input)"/>
|
|
|
84 <param name="reference" type="data" format="fasta,fasta.gz" optional="true" label="Select reference file" help="(Reference sequence)"/>
|
|
|
85 <conditional name="platform_cond">
|
|
|
86 <param name="platform_sel" type="select" label="Select platform" help="(Platform)">
|
|
|
87 <option value="illumina" selected="true">Illumina</option>
|
|
|
88 <option value="ion">ION</option>
|
|
|
89 </param>
|
|
|
90 <when value="illumina">
|
|
|
91 <conditional name="mode_cond">
|
|
|
92 <param name="mode_sel" type="select" label="Select mode" help="(Variance detection, Heteroplasmy)">
|
|
|
93 <option value="none" selected="true">None</option>
|
|
|
94 <option value="variance">Variance detection</option>
|
|
|
95 <option value="heteroplasmy">Heteroplasmy</option>
|
|
|
96 </param>
|
|
|
97 <when value="none"/>
|
|
|
98 <when value="variance"/>
|
|
|
99 <when value="heteroplasmy">
|
|
|
100 <param name="maf" type="float" value="0.007" min="0.007" max="0.49" label="Set minimum minor allele frequency" help="If you want to detect heteroplasmy, first assemble the genome without this option. Then give the resulting sequence as a reference and as a seed input. (Heteroplasmy)"/>
|
|
|
101 </when>
|
|
|
102 </conditional>
|
|
|
103 </when>
|
|
|
104 <when value="ion"/>
|
|
|
105 </conditional>
|
|
|
106 <conditional name="type_cond">
|
|
|
107 <param name="type_sel" type="select" label="Select assembly type" help="(Type)">
|
|
|
108 <option value="mito">Mitochondrial assembly (mito)</option>
|
|
|
109 </param>
|
|
|
110 <when value="mito">
|
|
|
111 <expand macro="content" min="12000" max="20000"/>
|
|
|
112 </when>
|
|
|
113 </conditional>
|
|
|
114 <param name="read_length" type="integer" value="151" label="Set read length" help="(Read Length)"/>
|
|
|
115 <param name="insert_size" type="integer" value="300" label="Set total insert size of your paired end reads" help="It doesn't have to be accurate but should be close enough. (Insert size)"/>
|
|
|
116 <param name="single_paired" type="select" label="Select read type" help="(Single/Paired)">
|
|
|
117 <option value="PE" selected="true">Paired</option>
|
|
|
118 <option value="SE">Single</option>
|
|
|
119 </param>
|
|
|
120 <param name="extend_seed_directly" type="boolean" truevalue="yes" falsevalue="no" checked="true" label="Extend seeds instead of finding matching reads" help="Only use this when your seed originates from the same sample and there are no possible mismatches. (Extend seed directly)"/>
|
|
|
121 <param name="kmer" type="integer" value="39" min="1" label="Set length of overlap between matching reads" help="If reads are shorter then 90 bp or you have low coverage data, this value should be decreased down to 23. For reads longer then 101 bp, this value can be increased, but this is not necessary."/>
|
|
|
122 <!-- Assembled reads are saved anyway. If this option is selected, the parameter is set to "2", otherwise "yes" -->
|
|
|
123 <param name="save_assembled_reads" type="boolean" truevalue="2" falsevalue="yes" label="Retain original IDs when saving assemled reads?"/>
|
|
|
124 <param name="insert_size_auto" type="boolean" truevalue="yes" falsevalue="no" checked="true" label="Finetune your insert size automatically?" help="(Insert size auto)"/>
|
|
|
125 <param name="use_quality_scores" type="boolean" truevalue="yes" falsevalue="no" label="Take quality scores into account?" help="Only use this when reads have low quality, like with the 300 bp reads of Illumina. (Use Quality Scores)"/>
|
|
|
126 <param name="out" type="select" multiple="true" optional="false" label="Select output file(s)" help="Only shown in history if selected here and generated by the specific NOVOplasty run.">
|
|
|
127 <!-- Assembly -->
|
|
|
128 <option value="ar" selected="true">Assembled reads</option>
|
|
|
129 <option value="c" selected="true">Contigs</option>
|
|
|
130 <option value="ca" selected="true">Circularized assembly</option>
|
|
|
131 <option value="cc" selected="true">Contig combinations</option>
|
|
|
132 <option value="ct" selected="true">Contigs Tmp</option>
|
|
|
133 <option value="mc" selected="true">Merged contigs</option>
|
|
|
134 <option value="ua" selected="true">Uncircularized assemblies</option>
|
|
|
135 <!-- Variance -->
|
|
|
136 <option value="v" selected="true">Variance</option>
|
|
|
137 <!-- Heteroplasmy -->
|
|
|
138 <option value="h" selected="true">Heteroplasmy</option>
|
|
|
139 <option value="ha" selected="true">Heteroplasmy assemblies</option>
|
|
|
140 <option value="lth" selected="true">Linkage table heteroplasmy</option>
|
|
|
141 <option value="pn" selected="true">Possible NUMTs</option>
|
|
|
142 <option value="pna" selected="true">Possible NUMTs assemblies</option>
|
|
|
143 <option value="ltn" selected="true">Linkage table NUMTs</option>
|
|
|
144 <option value="cm" selected="true">Circos mutations</option>
|
|
|
145 <option value="cl" selected="true">Circos links</option>
|
|
|
146 <!-- Others -->
|
|
|
147 <option value="config">Config file</option>
|
|
|
148 <option value="el">Extended log</option>
|
|
|
149 <option value="l">Log</option>
|
|
|
150 </param>
|
|
|
151 </inputs>
|
|
|
152 <outputs>
|
|
|
153 <!-- Assembly -->
|
|
|
154 <data name="out_c" format="fasta" from_work_dir="Contigs_1_result.fasta" label="${tool.name} on ${on_string}: Contigs">
|
|
|
155 <filter>'c' in out</filter>
|
|
|
156 </data>
|
|
|
157 <data name="out_ct" format="txt" from_work_dir="contigs_tmp_result.txt" label="${tool.name} on ${on_string}: Contigs Tmp">
|
|
|
158 <filter>'ct' in out</filter>
|
|
|
159 </data>
|
|
|
160 <data name="out_mc" format="txt" from_work_dir="Merged_contigs_result.txt" label="${tool.name} on ${on_string}: Merged Contigs">
|
|
|
161 <filter>'mc' in out</filter>
|
|
|
162 </data>
|
|
|
163 <data name="out_cc" format="txt">
|
|
|
164 <discover_datasets pattern="Option\_(?P<designation>.+)\_result\.txt" format="tabular" directory="" visible="true"/>
|
|
|
165 <filter>'cc' in out</filter>
|
|
|
166 </data>
|
|
|
167 <data name="out_ar1" format="fasta" from_work_dir="Assembled_reads_result_R1.fasta" label="${tool.name} on ${on_string}: Assembled reads R1">
|
|
|
168 <filter>'ar' in out</filter>
|
|
|
169 </data>
|
|
|
170 <data name="out_ar2" format="fasta" from_work_dir="Assembled_reads_result_R2.fasta" label="${tool.name} on ${on_string}: Assembled reads R2">
|
|
|
171 <filter>'ar' in out</filter>
|
|
|
172 </data>
|
|
|
173 <data name="out_ua" format="fasta" from_work_dir="Uncircularized_assemblies_1_result.fasta" label="${tool.name} on ${on_string}: Uncircularized assemblies">
|
|
|
174 <filter>'ua' in out</filter>
|
|
|
175 </data>
|
|
|
176 <data name="out_ca" format="fasta" from_work_dir="Circularized_assembly_1_result.fasta" label="${tool.name} on ${on_string}: Circularized assembly">
|
|
|
177 <filter>'ca' in out</filter>
|
|
|
178 </data>
|
|
|
179 <!-- Variance -->
|
|
|
180 <data name="out_v" format="vcf" from_work_dir="Variance_result.vcf" label="${tool.name} on ${on_string}: Variance">
|
|
|
181 <filter>'illumina' == platform_cond['platform_sel'] and 'variance' == platform_cond['mode_cond']['mode_sel'] and 'v' in out</filter>
|
|
|
182 </data>
|
|
|
183 <!-- Heteroplasmy -->
|
|
|
184 <data name="out_h" format="vcf" from_work_dir="Heteroplasmy_result.vcf" label="${tool.name} on ${on_string}: Heteroplasmy">
|
|
|
185 <filter>'heteroplasmy' == platform_cond['mode_cond']['mode_sel'] and 'h' in out</filter>
|
|
|
186 </data>
|
|
|
187 <data name="out_ha" format="fasta" from_work_dir="Heteroplasmy_assemblies_result.fasta" label="${tool.name} on ${on_string}: Heteroplasmy assemblies">
|
|
|
188 <filter>'heteroplasmy' == platform_cond['mode_cond']['mode_sel'] and 'ha' in out</filter>
|
|
|
189 </data>
|
|
|
190 <data name="out_lth" format="tabular" from_work_dir="Linkage_table_heteroplasmy_result.txt" label="${tool.name} on ${on_string}: Linkage table heteroplasmy">
|
|
|
191 <filter>'heteroplasmy' == platform_cond['mode_cond']['mode_sel'] and 'lth' in out</filter>
|
|
|
192 </data>
|
|
|
193 <data name="out_pn" format="vcf" from_work_dir="Possible_NUMTs_result.vcf" label="${tool.name} on ${on_string}: Possible NUMTs">
|
|
|
194 <filter>'heteroplasmy' == platform_cond['mode_cond']['mode_sel'] and 'pn' in out</filter>
|
|
|
195 </data>
|
|
|
196 <data name="out_pna" format="fasta" from_work_dir="Possible_NUMTs_assemblies_result.fasta" label="${tool.name} on ${on_string}: Possible_NUMTs_assemblies_project.fasta">
|
|
|
197 <filter>'heteroplasmy' == platform_cond['mode_cond']['mode_sel'] and 'pna' in out</filter>
|
|
|
198 </data>
|
|
|
199 <data name="out_ltn" format="tabular" from_work_dir="Linkage_table_NUMTs_result.txt" label="${tool.name} on ${on_string}: Linkage table NUMTs">
|
|
|
200 <filter>'heteroplasmy' == platform_cond['mode_cond']['mode_sel'] and 'ltn' in out</filter>
|
|
|
201 </data>
|
|
|
202 <data name="out_cm" format="txt" from_work_dir="Circos_mutations_result.txt" label="${tool.name} on ${on_string}: Circos mutations">
|
|
|
203 <filter>'heteroplasmy' == platform_cond['mode_cond']['mode_sel'] and 'cm' in out</filter>
|
|
|
204 </data>
|
|
|
205 <data name="out_cl" format="txt" from_work_dir="Circos_links_result.txt" label="${tool.name} on ${on_string}: Circos links">
|
|
|
206 <filter>'heteroplasmy' == platform_cond['mode_cond']['mode_sel'] and 'cl' in out</filter>
|
|
|
207 </data>
|
|
|
208 <!-- Others -->
|
|
|
209 <data name="out_l" format="txt" from_work_dir="log_result.txt" label="${tool.name} on ${on_string}: Log">
|
|
|
210 <filter>'l' in out</filter>
|
|
|
211 </data>
|
|
|
212 <data name="out_el" format="txt" from_work_dir="log_extended_result.txt" label="${tool.name} on ${on_string}: Extended log">
|
|
|
213 <filter>'el' in out</filter>
|
|
|
214 </data>
|
|
|
215 <data name="out_config" format="txt" from_work_dir="config.txt" label="${tool.name} on ${on_string}: Config file">
|
|
|
216 <filter>'config' in out</filter>
|
|
|
217 </data>
|
|
|
218 </outputs>
|
|
|
219 <tests>
|
|
|
220 <!-- #1 mito, default, separate, compressed inputs -->
|
|
|
221 <test expect_num_outputs="8">
|
|
|
222 <conditional name="reads_cond">
|
|
|
223 <param name="reads_sel" value="separate"/>
|
|
|
224 <param name="forward_reads" value="forward.fasta.gz"/>
|
|
|
225 <param name="reverse_reads" value="reverse.fasta.gz"/>
|
|
|
226 </conditional>
|
|
|
227 <conditional name="type_cond">
|
|
|
228 <param name="type_sel" value="mito"/>
|
|
|
229 </conditional>
|
|
|
230 <!-- fasta.gz not supported -->
|
|
|
231 <param name="seed_input" value="seed.fasta"/>
|
|
|
232 <!-- Assembly -->
|
|
|
233 <output name="out_c">
|
|
|
234 <assert_contents>
|
|
|
235 <has_n_lines n="110"/>
|
|
|
236 <has_line line=">Contig01+seed"/>
|
|
|
237 </assert_contents>
|
|
|
238 </output>
|
|
|
239 <output name="out_ct">
|
|
|
240 <assert_contents>
|
|
|
241 <has_n_lines n="2"/>
|
|
|
242 <has_line line=">seed"/>
|
|
|
243 </assert_contents>
|
|
|
244 </output>
|
|
|
245 <output name="out_ar1">
|
|
|
246 <assert_contents>
|
|
|
247 <has_n_lines n="10"/>
|
|
|
248 <has_line line=">676/1"/>
|
|
|
249 </assert_contents>
|
|
|
250 </output>
|
|
|
251 <output name="out_ar2">
|
|
|
252 <assert_contents>
|
|
|
253 <has_n_lines n="10"/>
|
|
|
254 <has_line line=">676/2"/>
|
|
|
255 </assert_contents>
|
|
|
256 </output>
|
|
|
257 <output name="out_ua">
|
|
|
258 <assert_contents>
|
|
|
259 <has_n_lines n="0"/>
|
|
|
260 </assert_contents>
|
|
|
261 </output>
|
|
|
262 </test>
|
|
|
263 <!-- #2 mito, combined, compressed inputs -->
|
|
|
264 <test expect_num_outputs="11">
|
|
|
265 <conditional name="reads_cond">
|
|
|
266 <param name="reads_sel" value="combined"/>
|
|
|
267 <param name="combined_reads" value="combined.fasta.gz"/>
|
|
|
268 </conditional>
|
|
|
269 <param name="seed_input" value="seed.fasta"/>
|
|
|
270 <param name="reference" value="reference.fasta.gz"/>
|
|
|
271 <conditional name="type_cond">
|
|
|
272 <param name="type_sel" value="mito"/>
|
|
|
273 <param name="genome_range_min" value="12001"/>
|
|
|
274 <param name="genome_range_max" value="20001"/>
|
|
|
275 </conditional>
|
|
|
276 <param name="read_length" value="150"/>
|
|
|
277 <param name="insert_size" value="301"/>
|
|
|
278 <param name="extend_seed_directly" value="no"/>
|
|
|
279 <param name="kmer" value="40"/>
|
|
|
280 <param name="save_assembled_reads" value="2"/>
|
|
|
281 <param name="insert_size_auto" value="yes"/>
|
|
|
282 <param name="use_quality_scores" value="yes"/>
|
|
|
283 <param name="out" value="c,ca,ua,mc,cc,ct,ar,l,el,config"/>
|
|
|
284 <!-- Assembly -->
|
|
|
285 <output name="out_c">
|
|
|
286 <assert_contents>
|
|
|
287 <has_n_lines n="0"/>
|
|
|
288 </assert_contents>
|
|
|
289 </output>
|
|
|
290 <output name="out_ct">
|
|
|
291 <assert_contents>
|
|
|
292 <has_n_lines n="2"/>
|
|
|
293 <has_line line=">4191"/>
|
|
|
294 </assert_contents>
|
|
|
295 </output>
|
|
|
296 <output name="out_ar1">
|
|
|
297 <assert_contents>
|
|
|
298 <has_n_lines n="142"/>
|
|
|
299 <has_line line=">38326//1"/>
|
|
|
300 </assert_contents>
|
|
|
301 </output>
|
|
|
302 <output name="out_ar2">
|
|
|
303 <assert_contents>
|
|
|
304 <has_n_lines n="142"/>
|
|
|
305 <has_line line=">38326//2"/>
|
|
|
306 </assert_contents>
|
|
|
307 </output>
|
|
|
308 <output name="out_ua">
|
|
|
309 <assert_contents>
|
|
|
310 <has_n_lines n="13"/>
|
|
|
311 <has_text_matching expression=">Contig01.+"/>
|
|
|
312 </assert_contents>
|
|
|
313 </output>
|
|
|
314 <!-- Others -->
|
|
|
315 <output name="out_l">
|
|
|
316 <assert_contents>
|
|
|
317 <has_text_matching expression=".+Assembly 1 finished.+"/>
|
|
|
318 </assert_contents>
|
|
|
319 </output>
|
|
|
320 <output name="out_el">
|
|
|
321 <assert_contents>
|
|
|
322 <has_text_matching expression=".+Assembly 1 finished.+"/>
|
|
|
323 </assert_contents>
|
|
|
324 </output>
|
|
|
325 <output name="out_config">
|
|
|
326 <assert_contents>
|
|
|
327 <has_n_lines n="34"/>
|
|
|
328 <has_line line="Project:"/>
|
|
|
329 <has_text_matching expression=".+"/>
|
|
|
330 </assert_contents>
|
|
|
331 </output>
|
|
|
332 </test>
|
|
|
333 <!-- #3 mito, variance, decompressed inputs -->
|
|
|
334 <test expect_num_outputs="12">
|
|
|
335 <conditional name="reads_cond">
|
|
|
336 <param name="reads_sel" value="combined"/>
|
|
|
337 <param name="combined_reads" value="combined.fasta"/>
|
|
|
338 </conditional>
|
|
|
339 <param name="seed_input" value="seed.fasta"/>
|
|
|
340 <param name="reference" value="reference.fasta"/>
|
|
|
341 <conditional name="platform_cond">
|
|
|
342 <param name="platform_sel" value="illumina"/>
|
|
|
343 <conditional name="mode_cond">
|
|
|
344 <param name="mode_sel" value="variance"/>
|
|
|
345 </conditional>
|
|
|
346 </conditional>
|
|
|
347 <conditional name="type_cond">
|
|
|
348 <param name="type_sel" value="mito"/>
|
|
|
349 </conditional>
|
|
|
350 <param name="out" value="c,ca,ua,mc,cc,ct,ar,v,l,el,config"/>
|
|
|
351 <!-- Assembly -->
|
|
|
352 <output name="out_c">
|
|
|
353 <assert_contents>
|
|
|
354 <has_n_lines n="0"/>
|
|
|
355 </assert_contents>
|
|
|
356 </output>
|
|
|
357 <output name="out_ct">
|
|
|
358 <assert_contents>
|
|
|
359 <has_n_lines n="2"/>
|
|
|
360 <has_line line=">seed"/>
|
|
|
361 </assert_contents>
|
|
|
362 </output>
|
|
|
363 <output name="out_ar1">
|
|
|
364 <assert_contents>
|
|
|
365 <has_n_lines n="10"/>
|
|
|
366 <has_line line=">20/1"/>
|
|
|
367 </assert_contents>
|
|
|
368 </output>
|
|
|
369 <output name="out_ar2">
|
|
|
370 <assert_contents>
|
|
|
371 <has_n_lines n="10"/>
|
|
|
372 <has_line line=">20/2"/>
|
|
|
373 </assert_contents>
|
|
|
374 </output>
|
|
|
375 <output name="out_ua">
|
|
|
376 <assert_contents>
|
|
|
377 <has_n_lines n="110"/>
|
|
|
378 <has_text_matching expression=">Contig01.+"/>
|
|
|
379 </assert_contents>
|
|
|
380 </output>
|
|
|
381 <!-- Variance -->
|
|
|
382 <output name="out_v">
|
|
|
383 <assert_contents>
|
|
|
384 <has_n_lines n="6"/>
|
|
|
385 <has_line line="##fileformat=VCFv4.0"/>
|
|
|
386 <has_text_matching expression="#CHROM	POS	ID	REF	ALT	QUAL	FILTER	INFO"/>
|
|
|
387 </assert_contents>
|
|
|
388 </output>
|
|
|
389 <!-- Others -->
|
|
|
390 <output name="out_l">
|
|
|
391 <assert_contents>
|
|
|
392 <has_text_matching expression=".+Assembly 1 finished.+"/>
|
|
|
393 </assert_contents>
|
|
|
394 </output>
|
|
|
395 <output name="out_el">
|
|
|
396 <assert_contents>
|
|
|
397 <has_text_matching expression=".+Assembly 1 finished.+"/>
|
|
|
398 </assert_contents>
|
|
|
399 </output>
|
|
|
400 <output name="out_config">
|
|
|
401 <assert_contents>
|
|
|
402 <has_n_lines n="34"/>
|
|
|
403 <has_line line="Project:"/>
|
|
|
404 <has_text_matching expression=".+"/>
|
|
|
405 </assert_contents>
|
|
|
406 </output>
|
|
|
407 </test>
|
|
|
408 <!-- #4 mito, untested parameters, compressed inputs -->
|
|
|
409 <test expect_num_outputs="16">
|
|
|
410 <conditional name="reads_cond">
|
|
|
411 <param name="reads_sel" value="separate"/>
|
|
|
412 <param name="forward_reads" value="forward.fasta.gz"/>
|
|
|
413 <param name="reverse_reads" value="reverse.fasta.gz"/>
|
|
|
414 </conditional>
|
|
|
415 <param name="seed_input" value="seed.fasta"/>
|
|
|
416 <conditional name="platform_cond">
|
|
|
417 <param name="platform_sel" value="ion"/>
|
|
|
418 </conditional>
|
|
|
419 <conditional name="type_cond">
|
|
|
420 <param name="type_sel" value="mito"/>
|
|
|
421 </conditional>
|
|
|
422 <param name="single_paired" value="SE"/>
|
|
|
423 <!-- Assembly -->
|
|
|
424 <output name="out_c">
|
|
|
425 <assert_contents>
|
|
|
426 <has_n_lines n="110"/>
|
|
|
427 <has_line line=">Contig01+seed"/>
|
|
|
428 </assert_contents>
|
|
|
429 </output>
|
|
|
430 <output name="out_ct">
|
|
|
431 <assert_contents>
|
|
|
432 <has_n_lines n="2"/>
|
|
|
433 <has_line line=">seed"/>
|
|
|
434 </assert_contents>
|
|
|
435 </output>
|
|
|
436 <output name="out_ar1">
|
|
|
437 <assert_contents>
|
|
|
438 <has_n_lines n="10"/>
|
|
|
439 <has_line line=">676/1"/>
|
|
|
440 </assert_contents>
|
|
|
441 </output>
|
|
|
442 <output name="out_ar2">
|
|
|
443 <assert_contents>
|
|
|
444 <has_n_lines n="10"/>
|
|
|
445 <has_line line=">676/2"/>
|
|
|
446 </assert_contents>
|
|
|
447 </output>
|
|
|
448 <output name="out_ua">
|
|
|
449 <assert_contents>
|
|
|
450 <has_n_lines n="0"/>
|
|
|
451 </assert_contents>
|
|
|
452 </output>
|
|
|
453 </test>
|
|
|
454 <!-- #5 mito, heteroplasmy -->
|
|
|
455 <test expect_num_outputs="16">
|
|
|
456 <conditional name="reads_cond">
|
|
|
457 <param name="reads_sel" value="combined"/>
|
|
|
458 <param name="combined_reads" value="combined.fasta"/>
|
|
|
459 </conditional>
|
|
|
460 <param name="seed_input" value="seed.fasta"/>
|
|
|
461 <param name="reference" value="reference.fasta"/>
|
|
|
462 <conditional name="platform_cond">
|
|
|
463 <param name="platform_sel" value="illumina"/>
|
|
|
464 <conditional name="mode_cond">
|
|
|
465 <param name="mode_sel" value="heteroplasmy"/>
|
|
|
466 <param name="maf" value="0.008"/>
|
|
|
467 </conditional>
|
|
|
468 </conditional>
|
|
|
469 <conditional name="type_cond">
|
|
|
470 <param name="type_sel" value="mito"/>
|
|
|
471 </conditional>
|
|
|
472 <!-- Assembly -->
|
|
|
473 <output name="out_c">
|
|
|
474 <assert_contents>
|
|
|
475 <has_n_lines n="0"/>
|
|
|
476 </assert_contents>
|
|
|
477 </output>
|
|
|
478 <output name="out_ct">
|
|
|
479 <assert_contents>
|
|
|
480 <has_n_lines n="0"/>
|
|
|
481 </assert_contents>
|
|
|
482 </output>
|
|
|
483 <output name="out_ar1">
|
|
|
484 <assert_contents>
|
|
|
485 <has_n_lines n="2"/>
|
|
|
486 <has_line line=">332/1"/>
|
|
|
487 </assert_contents>
|
|
|
488 </output>
|
|
|
489 <output name="out_ar2">
|
|
|
490 <assert_contents>
|
|
|
491 <has_n_lines n="2"/>
|
|
|
492 <has_line line=">332/2"/>
|
|
|
493 </assert_contents>
|
|
|
494 </output>
|
|
|
495 <output name="out_ua">
|
|
|
496 <assert_contents>
|
|
|
497 <has_n_lines n="0"/>
|
|
|
498 </assert_contents>
|
|
|
499 </output>
|
|
|
500 <!-- Heteroplasmy -->
|
|
|
501 <output name="out_h">
|
|
|
502 <assert_contents>
|
|
|
503 <has_n_lines n="8"/>
|
|
|
504 <has_text_matching expression="#CHROM	POS	ID	REF	ALT	QUAL	FILTER	INFO"/>
|
|
|
505 </assert_contents>
|
|
|
506 </output>
|
|
|
507 <output name="out_ha">
|
|
|
508 <assert_contents>
|
|
|
509 <has_n_lines n="0"/>
|
|
|
510 </assert_contents>
|
|
|
511 </output>
|
|
|
512 <output name="out_lth">
|
|
|
513 <assert_contents>
|
|
|
514 <has_n_lines n="0"/>
|
|
|
515 </assert_contents>
|
|
|
516 </output>
|
|
|
517 <output name="out_pn">
|
|
|
518 <assert_contents>
|
|
|
519 <has_n_lines n="0"/>
|
|
|
520 </assert_contents>
|
|
|
521 </output>
|
|
|
522 <output name="out_pna">
|
|
|
523 <assert_contents>
|
|
|
524 <has_n_lines n="0"/>
|
|
|
525 </assert_contents>
|
|
|
526 </output>
|
|
|
527 <output name="out_ltn">
|
|
|
528 <assert_contents>
|
|
|
529 <has_n_lines n="0"/>
|
|
|
530 </assert_contents>
|
|
|
531 </output>
|
|
|
532 <output name="out_cm">
|
|
|
533 <assert_contents>
|
|
|
534 <has_n_lines n="0"/>
|
|
|
535 </assert_contents>
|
|
|
536 </output>
|
|
|
537 <output name="out_cl">
|
|
|
538 <assert_contents>
|
|
|
539 <has_n_lines n="0"/>
|
|
|
540 </assert_contents>
|
|
|
541 </output>
|
|
|
542 </test>
|
|
|
543 </tests>
|
|
|
544 <help><![CDATA[
|
|
|
545 MitoFlow
|
|
|
546
|
|
|
547 ]]></help>
|
|
|
548 <citations>
|
|
|
549 <citation type="doi">10.1093/nar/gkw955</citation>
|
|
|
550 <citation type="doi">10.1093/nargab/lqz011</citation>
|
|
|
551 </citations>
|
|
|
552 </tool>
|
