|
6
|
1 <?xml version="1.0"?>
|
|
|
2 <tool id="novoplasty" name="NOVOplasty" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="18.01">
|
|
|
3 <description>de novo assembler for short circular genomes</description>
|
|
|
4 <macros>
|
|
|
5 <token name="@TOOL_VERSION@">4.2</token>
|
|
|
6 <token name="@VERSION_SUFFIX@">0</token>
|
|
|
7 <xml name="content" token_min="" token_max="">
|
|
|
8 <param name="genome_range_min" type="text" value="@MIN@" label="Set lower limit for the expected genome size range" help="(Genome Range)"/>
|
|
|
9 <param name="genome_range_max" type="text" value="@MAX@" label="Set upper limit for the expected genome size range" help="(Genome Range)"/>
|
|
|
10 </xml>
|
|
|
11 </macros>
|
|
|
12 <requirements>
|
|
|
13 <requirement type="package" version="@TOOL_VERSION@">novoplasty</requirement>
|
|
|
14 </requirements>
|
|
|
15 <version_command><![CDATA[NOVOPlasty.pl 2>&1 | grep 'Version' | awk -F' ' '{print $2}']]></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>
|
|
|
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 <section name="advanced" title="Advanced options">
|
|
|
86 <conditional name="platform_cond">
|
|
|
87 <param name="platform_sel" type="select" label="Select platform" help="(Platform)">
|
|
|
88 <option value="illumina" selected="true">Illumina</option>
|
|
|
89 <option value="ion">ION</option>
|
|
|
90 </param>
|
|
|
91 <when value="illumina">
|
|
|
92 <conditional name="mode_cond">
|
|
|
93 <param name="mode_sel" type="select" label="Select mode" help="(Variance detection, Heteroplasmy)">
|
|
|
94 <option value="none" selected="true">None</option>
|
|
|
95 <option value="variance">Variance detection</option>
|
|
|
96 <option value="heteroplasmy">Heteroplasmy</option>
|
|
|
97 </param>
|
|
|
98 <when value="none"/>
|
|
|
99 <when value="variance"/>
|
|
|
100 <when value="heteroplasmy">
|
|
|
101 <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)"/>
|
|
|
102 </when>
|
|
|
103 </conditional>
|
|
|
104 </when>
|
|
|
105 <when value="ion"/>
|
|
|
106 </conditional>
|
|
|
107 <conditional name="type_cond">
|
|
|
108 <param name="type_sel" type="select" label="Select assembly type" help="(Type)">
|
|
|
109 <option value="chloro" selected="true">Chloroplast assembly (chloro)</option>
|
|
|
110 <option value="mito">Mitochondrial assembly (mito)</option>
|
|
|
111 <option value="mito_plant">Mitochondrial assembly in plants (mito_plant)</option>
|
|
|
112 </param>
|
|
|
113 <when value="chloro">
|
|
|
114 <expand macro="content" min="120000" max="200000"/>
|
|
|
115 </when>
|
|
|
116 <when value="mito">
|
|
|
117 <expand macro="content" min="12000" max="20000"/>
|
|
|
118 </when>
|
|
|
119 <when value="mito_plant">
|
|
|
120 <expand macro="content" min="12000" max="20000"/>
|
|
|
121 <param name="chloroplast_sequence" type="data" format="fasta,fasta.gz" optional="true" label="Select chloroplast sequence" help="Chloroplast needs to be assembled before mitochondria. (Chloroplast sequence)"/>
|
|
|
122 </when>
|
|
|
123 </conditional>
|
|
|
124 <param name="read_length" type="integer" value="151" label="Set read length" help="(Read Length)"/>
|
|
|
125 <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)"/>
|
|
|
126 <param name="single_paired" type="select" label="Select read type" help="(Single/Paired)">
|
|
|
127 <option value="PE" selected="true">Paired</option>
|
|
|
128 <option value="SE">Single</option>
|
|
|
129 </param>
|
|
|
130 <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)"/>
|
|
|
131 <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."/>
|
|
|
132 <!-- Assembled reads are saved anyway. If this option is selected, the parameter is set to "2", otherwise "yes" -->
|
|
|
133 <param name="save_assembled_reads" type="boolean" truevalue="2" falsevalue="yes" label="Retain original IDs when saving assemled reads?"/>
|
|
|
134 <param name="insert_size_auto" type="boolean" truevalue="yes" falsevalue="no" checked="true" label="Finetune your insert size automatically?" help="(Insert size auto)"/>
|
|
|
135 <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)"/>
|
|
|
136 </section>
|
|
|
137 <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.">
|
|
|
138 <!-- Assembly -->
|
|
|
139 <option value="ar" selected="true">Assembled reads</option>
|
|
|
140 <option value="c" selected="true">Contigs</option>
|
|
|
141 <option value="ca" selected="true">Circularized assembly</option>
|
|
|
142 <option value="cc" selected="true">Contig combinations</option>
|
|
|
143 <option value="ct" selected="true">Contigs Tmp</option>
|
|
|
144 <option value="mc" selected="true">Merged contigs</option>
|
|
|
145 <option value="ua" selected="true">Uncircularized assemblies</option>
|
|
|
146 <!-- Variance -->
|
|
|
147 <option value="v" selected="true">Variance</option>
|
|
|
148 <!-- Heteroplasmy -->
|
|
|
149 <option value="h" selected="true">Heteroplasmy</option>
|
|
|
150 <option value="ha" selected="true">Heteroplasmy assemblies</option>
|
|
|
151 <option value="lth" selected="true">Linkage table heteroplasmy</option>
|
|
|
152 <option value="pn" selected="true">Possible NUMTs</option>
|
|
|
153 <option value="pna" selected="true">Possible NUMTs assemblies</option>
|
|
|
154 <option value="ltn" selected="true">Linkage table NUMTs</option>
|
|
|
155 <option value="cm" selected="true">Circos mutations</option>
|
|
|
156 <option value="cl" selected="true">Circos links</option>
|
|
|
157 <!-- Others -->
|
|
|
158 <option value="config">Config file</option>
|
|
|
159 <option value="el">Extended log</option>
|
|
|
160 <option value="l">Log</option>
|
|
|
161 </param>
|
|
|
162 </inputs>
|
|
|
163 <outputs>
|
|
|
164 <!-- Assembly -->
|
|
|
165 <data name="out_c" format="fasta" from_work_dir="Contigs_1_result.fasta" label="${tool.name} on ${on_string}: Contigs">
|
|
|
166 <filter>'c' in out</filter>
|
|
|
167 </data>
|
|
|
168 <data name="out_ct" format="txt" from_work_dir="contigs_tmp_result.txt" label="${tool.name} on ${on_string}: Contigs Tmp">
|
|
|
169 <filter>'ct' in out</filter>
|
|
|
170 </data>
|
|
|
171 <data name="out_mc" format="txt" from_work_dir="Merged_contigs_result.txt" label="${tool.name} on ${on_string}: Merged Contigs">
|
|
|
172 <filter>'mc' in out</filter>
|
|
|
173 </data>
|
|
|
174 <data name="out_cc" format="txt">
|
|
|
175 <discover_datasets pattern="Option\_(?P<designation>.+)\_result\.txt" format="tabular" directory="" visible="true"/>
|
|
|
176 <filter>'cc' in out</filter>
|
|
|
177 </data>
|
|
|
178 <data name="out_ar1" format="fasta" from_work_dir="Assembled_reads_result_R1.fasta" label="${tool.name} on ${on_string}: Assembled reads R1">
|
|
|
179 <filter>'ar' in out</filter>
|
|
|
180 </data>
|
|
|
181 <data name="out_ar2" format="fasta" from_work_dir="Assembled_reads_result_R2.fasta" label="${tool.name} on ${on_string}: Assembled reads R2">
|
|
|
182 <filter>'ar' in out</filter>
|
|
|
183 </data>
|
|
|
184 <data name="out_ua" format="fasta" from_work_dir="Uncircularized_assemblies_1_result.fasta" label="${tool.name} on ${on_string}: Uncircularized assemblies">
|
|
|
185 <filter>'ua' in out</filter>
|
|
|
186 </data>
|
|
|
187 <data name="out_ca" format="fasta" from_work_dir="Circularized_assembly_1_result.fasta" label="${tool.name} on ${on_string}: Circularized assembly">
|
|
|
188 <filter>'ca' in out</filter>
|
|
|
189 </data>
|
|
|
190 <!-- Variance -->
|
|
|
191 <data name="out_v" format="vcf" from_work_dir="Variance_result.vcf" label="${tool.name} on ${on_string}: Variance">
|
|
|
192 <filter>'illumina' == platform_cond['platform_sel'] and 'variance' == platform_cond['mode_cond']['mode_sel'] and 'v' in out</filter>
|
|
|
193 </data>
|
|
|
194 <!-- Heteroplasmy -->
|
|
|
195 <data name="out_h" format="vcf" from_work_dir="Heteroplasmy_result.vcf" label="${tool.name} on ${on_string}: Heteroplasmy">
|
|
|
196 <filter>'heteroplasmy' == platform_cond['mode_cond']['mode_sel'] and 'h' in out</filter>
|
|
|
197 </data>
|
|
|
198 <data name="out_ha" format="fasta" from_work_dir="Heteroplasmy_assemblies_result.fasta" label="${tool.name} on ${on_string}: Heteroplasmy assemblies">
|
|
|
199 <filter>'heteroplasmy' == platform_cond['mode_cond']['mode_sel'] and 'ha' in out</filter>
|
|
|
200 </data>
|
|
|
201 <data name="out_lth" format="tabular" from_work_dir="Linkage_table_heteroplasmy_result.txt" label="${tool.name} on ${on_string}: Linkage table heteroplasmy">
|
|
|
202 <filter>'heteroplasmy' == platform_cond['mode_cond']['mode_sel'] and 'lth' in out</filter>
|
|
|
203 </data>
|
|
|
204 <data name="out_pn" format="vcf" from_work_dir="Possible_NUMTs_result.vcf" label="${tool.name} on ${on_string}: Possible NUMTs">
|
|
|
205 <filter>'heteroplasmy' == platform_cond['mode_cond']['mode_sel'] and 'pn' in out</filter>
|
|
|
206 </data>
|
|
|
207 <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">
|
|
|
208 <filter>'heteroplasmy' == platform_cond['mode_cond']['mode_sel'] and 'pna' in out</filter>
|
|
|
209 </data>
|
|
|
210 <data name="out_ltn" format="tabular" from_work_dir="Linkage_table_NUMTs_result.txt" label="${tool.name} on ${on_string}: Linkage table NUMTs">
|
|
|
211 <filter>'heteroplasmy' == platform_cond['mode_cond']['mode_sel'] and 'ltn' in out</filter>
|
|
|
212 </data>
|
|
|
213 <data name="out_cm" format="txt" from_work_dir="Circos_mutations_result.txt" label="${tool.name} on ${on_string}: Circos mutations">
|
|
|
214 <filter>'heteroplasmy' == platform_cond['mode_cond']['mode_sel'] and 'cm' in out</filter>
|
|
|
215 </data>
|
|
|
216 <data name="out_cl" format="txt" from_work_dir="Circos_links_result.txt" label="${tool.name} on ${on_string}: Circos links">
|
|
|
217 <filter>'heteroplasmy' == platform_cond['mode_cond']['mode_sel'] and 'cl' in out</filter>
|
|
|
218 </data>
|
|
|
219 <!-- Others -->
|
|
|
220 <data name="out_l" format="txt" from_work_dir="log_result.txt" label="${tool.name} on ${on_string}: Log">
|
|
|
221 <filter>'l' in out</filter>
|
|
|
222 </data>
|
|
|
223 <data name="out_el" format="txt" from_work_dir="log_extended_result.txt" label="${tool.name} on ${on_string}: Extended log">
|
|
|
224 <filter>'el' in out</filter>
|
|
|
225 </data>
|
|
|
226 <data name="out_config" format="txt" from_work_dir="config.txt" label="${tool.name} on ${on_string}: Config file">
|
|
|
227 <filter>'config' in out</filter>
|
|
|
228 </data>
|
|
|
229 </outputs>
|
|
|
230 <tests>
|
|
|
231 <!-- #1 mito, default, separate, compressed inputs -->
|
|
|
232 <test expect_num_outputs="8">
|
|
|
233 <conditional name="reads_cond">
|
|
|
234 <param name="reads_sel" value="separate"/>
|
|
|
235 <param name="forward_reads" value="forward.fasta.gz"/>
|
|
|
236 <param name="reverse_reads" value="reverse.fasta.gz"/>
|
|
|
237 </conditional>
|
|
|
238 <conditional name="type_cond">
|
|
|
239 <param name="type_sel" value="mito"/>
|
|
|
240 </conditional>
|
|
|
241 <!-- fasta.gz not supported -->
|
|
|
242 <param name="seed_input" value="seed.fasta"/>
|
|
|
243 <!-- Assembly -->
|
|
|
244 <output name="out_c">
|
|
|
245 <assert_contents>
|
|
|
246 <has_n_lines n="110"/>
|
|
|
247 <has_line line=">Contig01+seed"/>
|
|
|
248 </assert_contents>
|
|
|
249 </output>
|
|
|
250 <output name="out_ct">
|
|
|
251 <assert_contents>
|
|
|
252 <has_n_lines n="2"/>
|
|
|
253 <has_line line=">seed"/>
|
|
|
254 </assert_contents>
|
|
|
255 </output>
|
|
|
256 <output name="out_ar1">
|
|
|
257 <assert_contents>
|
|
|
258 <has_n_lines n="10"/>
|
|
|
259 <has_line line=">676/1"/>
|
|
|
260 </assert_contents>
|
|
|
261 </output>
|
|
|
262 <output name="out_ar2">
|
|
|
263 <assert_contents>
|
|
|
264 <has_n_lines n="10"/>
|
|
|
265 <has_line line=">676/2"/>
|
|
|
266 </assert_contents>
|
|
|
267 </output>
|
|
|
268 <output name="out_ua">
|
|
|
269 <assert_contents>
|
|
|
270 <has_n_lines n="0"/>
|
|
|
271 </assert_contents>
|
|
|
272 </output>
|
|
|
273 </test>
|
|
|
274 <!-- #2 mito, combined, compressed inputs -->
|
|
|
275 <test expect_num_outputs="11">
|
|
|
276 <conditional name="reads_cond">
|
|
|
277 <param name="reads_sel" value="combined"/>
|
|
|
278 <param name="combined_reads" value="combined.fasta.gz"/>
|
|
|
279 </conditional>
|
|
|
280 <param name="seed_input" value="seed.fasta"/>
|
|
|
281 <param name="reference" value="reference.fasta.gz"/>
|
|
|
282 <conditional name="type_cond">
|
|
|
283 <param name="type_sel" value="mito"/>
|
|
|
284 <param name="genome_range_min" value="12001"/>
|
|
|
285 <param name="genome_range_max" value="20001"/>
|
|
|
286 </conditional>
|
|
|
287 <param name="read_length" value="150"/>
|
|
|
288 <param name="insert_size" value="301"/>
|
|
|
289 <param name="extend_seed_directly" value="no"/>
|
|
|
290 <param name="kmer" value="40"/>
|
|
|
291 <param name="save_assembled_reads" value="2"/>
|
|
|
292 <param name="insert_size_auto" value="yes"/>
|
|
|
293 <param name="use_quality_scores" value="yes"/>
|
|
|
294 <param name="out" value="c,ca,ua,mc,cc,ct,ar,l,el,config"/>
|
|
|
295 <!-- Assembly -->
|
|
|
296 <output name="out_c">
|
|
|
297 <assert_contents>
|
|
|
298 <has_n_lines n="0"/>
|
|
|
299 </assert_contents>
|
|
|
300 </output>
|
|
|
301 <output name="out_ct">
|
|
|
302 <assert_contents>
|
|
|
303 <has_n_lines n="2"/>
|
|
|
304 <has_line line=">4191"/>
|
|
|
305 </assert_contents>
|
|
|
306 </output>
|
|
|
307 <output name="out_ar1">
|
|
|
308 <assert_contents>
|
|
|
309 <has_n_lines n="142"/>
|
|
|
310 <has_line line=">38326//1"/>
|
|
|
311 </assert_contents>
|
|
|
312 </output>
|
|
|
313 <output name="out_ar2">
|
|
|
314 <assert_contents>
|
|
|
315 <has_n_lines n="142"/>
|
|
|
316 <has_line line=">38326//2"/>
|
|
|
317 </assert_contents>
|
|
|
318 </output>
|
|
|
319 <output name="out_ua">
|
|
|
320 <assert_contents>
|
|
|
321 <has_n_lines n="13"/>
|
|
|
322 <has_text_matching expression=">Contig01.+"/>
|
|
|
323 </assert_contents>
|
|
|
324 </output>
|
|
|
325 <!-- Others -->
|
|
|
326 <output name="out_l">
|
|
|
327 <assert_contents>
|
|
|
328 <has_text_matching expression=".+Assembly 1 finished.+"/>
|
|
|
329 </assert_contents>
|
|
|
330 </output>
|
|
|
331 <output name="out_el">
|
|
|
332 <assert_contents>
|
|
|
333 <has_text_matching expression=".+Assembly 1 finished.+"/>
|
|
|
334 </assert_contents>
|
|
|
335 </output>
|
|
|
336 <output name="out_config">
|
|
|
337 <assert_contents>
|
|
|
338 <has_n_lines n="34"/>
|
|
|
339 <has_line line="Project:"/>
|
|
|
340 <has_text_matching expression=".+"/>
|
|
|
341 </assert_contents>
|
|
|
342 </output>
|
|
|
343 </test>
|
|
|
344 <!-- #3 mito, variance, decompressed inputs -->
|
|
|
345 <test expect_num_outputs="12">
|
|
|
346 <conditional name="reads_cond">
|
|
|
347 <param name="reads_sel" value="combined"/>
|
|
|
348 <param name="combined_reads" value="combined.fasta"/>
|
|
|
349 </conditional>
|
|
|
350 <param name="seed_input" value="seed.fasta"/>
|
|
|
351 <param name="reference" value="reference.fasta"/>
|
|
|
352 <conditional name="platform_cond">
|
|
|
353 <param name="platform_sel" value="illumina"/>
|
|
|
354 <conditional name="mode_cond">
|
|
|
355 <param name="mode_sel" value="variance"/>
|
|
|
356 </conditional>
|
|
|
357 </conditional>
|
|
|
358 <conditional name="type_cond">
|
|
|
359 <param name="type_sel" value="mito"/>
|
|
|
360 </conditional>
|
|
|
361 <param name="out" value="c,ca,ua,mc,cc,ct,ar,v,l,el,config"/>
|
|
|
362 <!-- Assembly -->
|
|
|
363 <output name="out_c">
|
|
|
364 <assert_contents>
|
|
|
365 <has_n_lines n="0"/>
|
|
|
366 </assert_contents>
|
|
|
367 </output>
|
|
|
368 <output name="out_ct">
|
|
|
369 <assert_contents>
|
|
|
370 <has_n_lines n="2"/>
|
|
|
371 <has_line line=">seed"/>
|
|
|
372 </assert_contents>
|
|
|
373 </output>
|
|
|
374 <output name="out_ar1">
|
|
|
375 <assert_contents>
|
|
|
376 <has_n_lines n="10"/>
|
|
|
377 <has_line line=">20/1"/>
|
|
|
378 </assert_contents>
|
|
|
379 </output>
|
|
|
380 <output name="out_ar2">
|
|
|
381 <assert_contents>
|
|
|
382 <has_n_lines n="10"/>
|
|
|
383 <has_line line=">20/2"/>
|
|
|
384 </assert_contents>
|
|
|
385 </output>
|
|
|
386 <output name="out_ua">
|
|
|
387 <assert_contents>
|
|
|
388 <has_n_lines n="110"/>
|
|
|
389 <has_text_matching expression=">Contig01.+"/>
|
|
|
390 </assert_contents>
|
|
|
391 </output>
|
|
|
392 <!-- Variance -->
|
|
|
393 <output name="out_v">
|
|
|
394 <assert_contents>
|
|
|
395 <has_n_lines n="6"/>
|
|
|
396 <has_line line="##fileformat=VCFv4.0"/>
|
|
|
397 <has_text_matching expression="#CHROM	POS	ID	REF	ALT	QUAL	FILTER	INFO"/>
|
|
|
398 </assert_contents>
|
|
|
399 </output>
|
|
|
400 <!-- Others -->
|
|
|
401 <output name="out_l">
|
|
|
402 <assert_contents>
|
|
|
403 <has_text_matching expression=".+Assembly 1 finished.+"/>
|
|
|
404 </assert_contents>
|
|
|
405 </output>
|
|
|
406 <output name="out_el">
|
|
|
407 <assert_contents>
|
|
|
408 <has_text_matching expression=".+Assembly 1 finished.+"/>
|
|
|
409 </assert_contents>
|
|
|
410 </output>
|
|
|
411 <output name="out_config">
|
|
|
412 <assert_contents>
|
|
|
413 <has_n_lines n="34"/>
|
|
|
414 <has_line line="Project:"/>
|
|
|
415 <has_text_matching expression=".+"/>
|
|
|
416 </assert_contents>
|
|
|
417 </output>
|
|
|
418 </test>
|
|
|
419 <!-- #4 mito, untested parameters, compressed inputs -->
|
|
|
420 <test expect_num_outputs="16">
|
|
|
421 <conditional name="reads_cond">
|
|
|
422 <param name="reads_sel" value="separate"/>
|
|
|
423 <param name="forward_reads" value="forward.fasta.gz"/>
|
|
|
424 <param name="reverse_reads" value="reverse.fasta.gz"/>
|
|
|
425 </conditional>
|
|
|
426 <param name="seed_input" value="seed.fasta"/>
|
|
|
427 <conditional name="platform_cond">
|
|
|
428 <param name="platform_sel" value="ion"/>
|
|
|
429 </conditional>
|
|
|
430 <conditional name="type_cond">
|
|
|
431 <param name="type_sel" value="mito"/>
|
|
|
432 </conditional>
|
|
|
433 <param name="single_paired" value="SE"/>
|
|
|
434 <!-- Assembly -->
|
|
|
435 <output name="out_c">
|
|
|
436 <assert_contents>
|
|
|
437 <has_n_lines n="110"/>
|
|
|
438 <has_line line=">Contig01+seed"/>
|
|
|
439 </assert_contents>
|
|
|
440 </output>
|
|
|
441 <output name="out_ct">
|
|
|
442 <assert_contents>
|
|
|
443 <has_n_lines n="2"/>
|
|
|
444 <has_line line=">seed"/>
|
|
|
445 </assert_contents>
|
|
|
446 </output>
|
|
|
447 <output name="out_ar1">
|
|
|
448 <assert_contents>
|
|
|
449 <has_n_lines n="10"/>
|
|
|
450 <has_line line=">676/1"/>
|
|
|
451 </assert_contents>
|
|
|
452 </output>
|
|
|
453 <output name="out_ar2">
|
|
|
454 <assert_contents>
|
|
|
455 <has_n_lines n="10"/>
|
|
|
456 <has_line line=">676/2"/>
|
|
|
457 </assert_contents>
|
|
|
458 </output>
|
|
|
459 <output name="out_ua">
|
|
|
460 <assert_contents>
|
|
|
461 <has_n_lines n="0"/>
|
|
|
462 </assert_contents>
|
|
|
463 </output>
|
|
|
464 </test>
|
|
|
465 <!-- #5 mito, heteroplasmy -->
|
|
|
466 <test expect_num_outputs="16">
|
|
|
467 <conditional name="reads_cond">
|
|
|
468 <param name="reads_sel" value="combined"/>
|
|
|
469 <param name="combined_reads" value="combined.fasta"/>
|
|
|
470 </conditional>
|
|
|
471 <param name="seed_input" value="seed.fasta"/>
|
|
|
472 <param name="reference" value="reference.fasta"/>
|
|
|
473 <conditional name="platform_cond">
|
|
|
474 <param name="platform_sel" value="illumina"/>
|
|
|
475 <conditional name="mode_cond">
|
|
|
476 <param name="mode_sel" value="heteroplasmy"/>
|
|
|
477 <param name="maf" value="0.008"/>
|
|
|
478 </conditional>
|
|
|
479 </conditional>
|
|
|
480 <conditional name="type_cond">
|
|
|
481 <param name="type_sel" value="mito"/>
|
|
|
482 </conditional>
|
|
|
483 <!-- Assembly -->
|
|
|
484 <output name="out_c">
|
|
|
485 <assert_contents>
|
|
|
486 <has_n_lines n="0"/>
|
|
|
487 </assert_contents>
|
|
|
488 </output>
|
|
|
489 <output name="out_ct">
|
|
|
490 <assert_contents>
|
|
|
491 <has_n_lines n="0"/>
|
|
|
492 </assert_contents>
|
|
|
493 </output>
|
|
|
494 <output name="out_ar1">
|
|
|
495 <assert_contents>
|
|
|
496 <has_n_lines n="2"/>
|
|
|
497 <has_line line=">332/1"/>
|
|
|
498 </assert_contents>
|
|
|
499 </output>
|
|
|
500 <output name="out_ar2">
|
|
|
501 <assert_contents>
|
|
|
502 <has_n_lines n="2"/>
|
|
|
503 <has_line line=">332/2"/>
|
|
|
504 </assert_contents>
|
|
|
505 </output>
|
|
|
506 <output name="out_ua">
|
|
|
507 <assert_contents>
|
|
|
508 <has_n_lines n="0"/>
|
|
|
509 </assert_contents>
|
|
|
510 </output>
|
|
|
511 <!-- Heteroplasmy -->
|
|
|
512 <output name="out_h">
|
|
|
513 <assert_contents>
|
|
|
514 <has_n_lines n="8"/>
|
|
|
515 <has_text_matching expression="#CHROM	POS	ID	REF	ALT	QUAL	FILTER	INFO"/>
|
|
|
516 </assert_contents>
|
|
|
517 </output>
|
|
|
518 <output name="out_ha">
|
|
|
519 <assert_contents>
|
|
|
520 <has_n_lines n="0"/>
|
|
|
521 </assert_contents>
|
|
|
522 </output>
|
|
|
523 <output name="out_lth">
|
|
|
524 <assert_contents>
|
|
|
525 <has_n_lines n="0"/>
|
|
|
526 </assert_contents>
|
|
|
527 </output>
|
|
|
528 <output name="out_pn">
|
|
|
529 <assert_contents>
|
|
|
530 <has_n_lines n="0"/>
|
|
|
531 </assert_contents>
|
|
|
532 </output>
|
|
|
533 <output name="out_pna">
|
|
|
534 <assert_contents>
|
|
|
535 <has_n_lines n="0"/>
|
|
|
536 </assert_contents>
|
|
|
537 </output>
|
|
|
538 <output name="out_ltn">
|
|
|
539 <assert_contents>
|
|
|
540 <has_n_lines n="0"/>
|
|
|
541 </assert_contents>
|
|
|
542 </output>
|
|
|
543 <output name="out_cm">
|
|
|
544 <assert_contents>
|
|
|
545 <has_n_lines n="0"/>
|
|
|
546 </assert_contents>
|
|
|
547 </output>
|
|
|
548 <output name="out_cl">
|
|
|
549 <assert_contents>
|
|
|
550 <has_n_lines n="0"/>
|
|
|
551 </assert_contents>
|
|
|
552 </output>
|
|
|
553 </test>
|
|
|
554 <!-- #6 chloro, default -->
|
|
|
555 <test expect_num_outputs="8">
|
|
|
556 <conditional name="reads_cond">
|
|
|
557 <param name="reads_sel" value="combined"/>
|
|
|
558 <param name="combined_reads" value="combined.fasta"/>
|
|
|
559 </conditional>
|
|
|
560 <param name="seed_input" value="seed.fasta"/>
|
|
|
561 <conditional name="type_cond">
|
|
|
562 <param name="type_sel" value="chloro"/>
|
|
|
563 </conditional>
|
|
|
564 <!-- Assembly -->
|
|
|
565 <output name="out_c">
|
|
|
566 <assert_contents>
|
|
|
567 <has_n_lines n="110"/>
|
|
|
568 <has_line line=">Contig01+seed"/>
|
|
|
569 </assert_contents>
|
|
|
570 </output>
|
|
|
571 <output name="out_ct">
|
|
|
572 <assert_contents>
|
|
|
573 <has_n_lines n="2"/>
|
|
|
574 <has_line line=">seed"/>
|
|
|
575 </assert_contents>
|
|
|
576 </output>
|
|
|
577 <output name="out_ar1">
|
|
|
578 <assert_contents>
|
|
|
579 <has_n_lines n="10"/>
|
|
|
580 <has_line line=">20/1"/>
|
|
|
581 </assert_contents>
|
|
|
582 </output>
|
|
|
583 <output name="out_ar2">
|
|
|
584 <assert_contents>
|
|
|
585 <has_n_lines n="10"/>
|
|
|
586 <has_line line=">20/2"/>
|
|
|
587 </assert_contents>
|
|
|
588 </output>
|
|
|
589 <output name="out_ua">
|
|
|
590 <assert_contents>
|
|
|
591 <has_n_lines n="0"/>
|
|
|
592 </assert_contents>
|
|
|
593 </output>
|
|
|
594 </test>
|
|
|
595 <!-- #7 mito_plant, default -->
|
|
|
596 <test expect_num_outputs="8">
|
|
|
597 <conditional name="reads_cond">
|
|
|
598 <param name="reads_sel" value="combined"/>
|
|
|
599 <param name="combined_reads" value="combined.fasta"/>
|
|
|
600 </conditional>
|
|
|
601 <param name="seed_input" value="seed.fasta"/>
|
|
|
602 <param name="reference" value="reference.fasta"/>
|
|
|
603 <conditional name="type_cond">
|
|
|
604 <param name="type_sel" value="mito_plant"/>
|
|
|
605 <param name="chloroplast_sequence" value="chloroplast.fasta"/>
|
|
|
606 </conditional>
|
|
|
607 <!-- Assembly -->
|
|
|
608 <output name="out_c">
|
|
|
609 <assert_contents>
|
|
|
610 <has_n_lines n="0"/>
|
|
|
611 </assert_contents>
|
|
|
612 </output>
|
|
|
613 <output name="out_ct">
|
|
|
614 <assert_contents>
|
|
|
615 <has_n_lines n="2"/>
|
|
|
616 </assert_contents>
|
|
|
617 </output>
|
|
|
618 <output name="out_ar1">
|
|
|
619 <assert_contents>
|
|
|
620 <has_n_lines n="10"/>
|
|
|
621 <has_line line=">20/1"/>
|
|
|
622 </assert_contents>
|
|
|
623 </output>
|
|
|
624 <output name="out_ar2">
|
|
|
625 <assert_contents>
|
|
|
626 <has_n_lines n="10"/>
|
|
|
627 <has_line line=">20/2"/>
|
|
|
628 </assert_contents>
|
|
|
629 </output>
|
|
|
630 <output name="out_ua">
|
|
|
631 <assert_contents>
|
|
|
632 <has_text_matching expression=">Contig01.+"/>
|
|
|
633 <has_n_lines n="110"/>
|
|
|
634 </assert_contents>
|
|
|
635 </output>
|
|
|
636 </test>
|
|
|
637 </tests>
|
|
|
638 <help><![CDATA[
|
|
|
639 .. class:: infomark
|
|
|
640
|
|
|
641 **What it does**
|
|
|
642
|
|
|
643 *NOVOPlasty* is a de novo assembler and heteroplasmy/variance caller for short circular genomes.
|
|
|
644
|
|
|
645 **Input**
|
|
|
646
|
|
|
647 - suitable seed: There are different types of seed possible: A single read from the dataset that originates from the organelle genome. A organelle sequence derived from the same or a related species. A complete organelle sequence of a more distant species (recommended when there is no close related sequence available)
|
|
|
648 - input reads: Either two separate files(forward and reverse) or a merged fastq/fasta file. Multiple libraries as input is not yet supported. There is also an Ion Torrent option, but it does not produce the best results.
|
|
|
649 - reference: The assembly will still be de novo, but references of the same genus can be used as a guide to resolve duplicated regions in the plant mitochondria or the inverted repeat in the chloroplast.
|
|
|
650
|
|
|
651 **Output**
|
|
|
652
|
|
|
653 *Assembly*
|
|
|
654
|
|
|
655 - Contigs: This file contains all the assembled contigs.
|
|
|
656 - Circularized assembly: When NOVOPlasty is able to circularize one contig, without any additional contigs being produced, it will just output this circularized fasta file.
|
|
|
657 - Uncircularized assembly
|
|
|
658 - Merged contigs: When there are multiple contigs, NOVOPlasty will try to combine all contigs in to a complete circular genome, all the different possibilities can be found in this file.
|
|
|
659 - Contigs Tmp: If non of the above files are outputted or are empty, you can retrieve some contigs from this file.
|
|
|
660
|
|
|
661 *Variance*
|
|
|
662
|
|
|
663 - Variance: When the variance detection option is selected, an additional VCF file will be outputted.
|
|
|
664
|
|
|
665 *Heteroplasmy*
|
|
|
666
|
|
|
667 - Heteroplasmy: This VCF file contains all the detected heteroplasmy positions.
|
|
|
668 - Heteroplasmy assemblies: NOVOPlasty will assemble around each position of the above VCF output. The resulting assemblies will be outputted in this fasta file.
|
|
|
669 - Linkage table heteroplasmy: This file contains all the polymorphisms that are fully/partially/not linked with each detected heteroplasmic position.
|
|
|
670 - Possible NUMTs / Possible NUMTs assemblies / Linkage table NUMTs: These are the same files as the above three, but then for all the polymorphisms that were identified as NUMT origin.
|
|
|
671 - Circos mutations / Circos links: These files can be used to draw figures with Circos.
|
|
|
672
|
|
|
673 *Others*
|
|
|
674
|
|
|
675 - Log: This is a basic log file with the information that shows up on your terminal.
|
|
|
676 - Extended log: An elaborate log file.
|
|
|
677 - Config file
|
|
|
678
|
|
|
679 .. class:: infomark
|
|
|
680
|
|
|
681 **References**
|
|
|
682
|
|
|
683 More information are available on `GitHub <https://github.com/ndierckx/NOVOPlasty/>`_.
|
|
|
684 ]]></help>
|
|
|
685 <citations>
|
|
|
686 <citation type="doi">10.1093/nar/gkw955</citation>
|
|
|
687 <citation type="doi">10.1093/nargab/lqz011</citation>
|
|
|
688 </citations>
|
|
|
689 </tool>
|
