Mercurial > repos > bgruening > hifiasm

comparison hifiasm.xml @ 19:284a1cb82b61 draft

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/hifiasm commit 47637031b6d38391762d49ee5c8317810aa8d766
author iuc
date Thu, 11 Jul 2024 14:44:22 +0000
parents df25d4fb79b7
children 5161f204ba63
comparison
equal deleted inserted replaced
18:df25d4fb79b7 19:284a1cb82b61
1 <tool id="hifiasm" name="Hifiasm" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@"> 1 <tool id="hifiasm" name="Hifiasm" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@">
2 <description>haplotype-resolved de novo assembler for PacBio Hifi reads</description> 2 <description>haplotype-resolved de novo assembler for PacBio Hifi reads</description>
3 <macros> 3 <macros>
4 <token name="@TOOL_VERSION@">0.19.8</token> 4 <token name="@TOOL_VERSION@">0.19.9</token>
5 <token name="@VERSION_SUFFIX@">1</token> 5 <token name="@VERSION_SUFFIX@">0</token>
6 <token name="@FORMATS@">fasta,fasta.gz,fastq,fastq.gz</token> 6 <token name="@FORMATS@">fasta,fasta.gz,fastq,fastq.gz</token>
7 <xml name="reads"> 7 <xml name="reads">
8 <param name="reads" type="data" format="@FORMATS@" multiple="true" label="Input reads" /> 8 <param name="reads" type="data" format="@FORMATS@" multiple="true" label="Input reads"/>
9 </xml> 9 </xml>
10 </macros> 10 </macros>
11 <requirements> 11 <requirements>
12 <requirement type="package" version="@TOOL_VERSION@">hifiasm</requirement> 12 <requirement type="package" version="@TOOL_VERSION@">hifiasm</requirement>
13 <requirement type="package" version="0.1">yak</requirement> 13 <requirement type="package" version="0.1">yak</requirement>
182 <param name="mode_selector" type="select" label="Assembly mode"> 182 <param name="mode_selector" type="select" label="Assembly mode">
183 <option value="standard">Standard</option> 183 <option value="standard">Standard</option>
184 <option value="trio">Trio mode</option> 184 <option value="trio">Trio mode</option>
185 </param> 185 </param>
186 <when value="standard"> 186 <when value="standard">
187 <expand macro="reads" /> 187 <expand macro="reads"/>
188 </when> 188 </when>
189 <when value="trio"> 189 <when value="trio">
190 <expand macro="reads" /> 190 <expand macro="reads"/>
191 <conditional name="trioinput"> 191 <conditional name="trioinput">
192 <param name="trio_input_selector" type="select" label="What parental information are you using?"> 192 <param name="trio_input_selector" type="select" label="What parental information are you using?">
193 <option value="reads">Parental reads (FASTQ files, gzipped or otherwise)</option> 193 <option value="reads">Parental reads (FASTQ files, gzipped or otherwise)</option>
194 <option value="lists">Lists of reads assigned by parent (text files)</option> 194 <option value="lists">Lists of reads assigned by parent (text files)</option>
195 </param> 195 </param>
196 <when value="reads"> 196 <when value="reads">
197 <param name="hap1_reads" type="data" format="fastq,fastq.gz" multiple="true" label="Haplotype 1 reads" /> 197 <param name="hap1_reads" type="data" format="fastq,fastq.gz" multiple="true" label="Haplotype 1 reads"/>
198 <param name="hap2_reads" type="data" format="fastq,fastq.gz" multiple="true" label="Haplotype 2 reads" /> 198 <param name="hap2_reads" type="data" format="fastq,fastq.gz" multiple="true" label="Haplotype 2 reads"/>
199 </when> 199 </when>
200 <when value="lists"> 200 <when value="lists">
201 <param name="hap1_list" type="data" format="txt,tabular" label="Haplotype 1 read list" /> 201 <param name="hap1_list" type="data" format="txt,tabular" label="Haplotype 1 read list"/>
202 <param name="hap2_list" type="data" format="txt,tabular" label="Haplotype 2 read list" /> 202 <param name="hap2_list" type="data" format="txt,tabular" label="Haplotype 2 read list"/>
203 </when> 203 </when>
204 </conditional> 204 </conditional>
205 <param name="max_kmers" argument="-c" type="integer" value="2" label="Lower bound of the binned k-mer's frequency" /> 205 <param name="max_kmers" argument="-c" type="integer" value="2" label="Lower bound of the binned k-mer's frequency"/>
206 <param name="min_kmers" argument="-d" type="integer" value="5" label="Upper bound of the binned k-mer's frequency" /> 206 <param name="min_kmers" argument="-d" type="integer" value="5" label="Upper bound of the binned k-mer's frequency"/>
207 <param name="yak_kmer_length" type="integer" min="0" max="64" value="31" label="Yak counter k-mer length" /> 207 <param name="yak_kmer_length" type="integer" min="0" max="64" value="31" label="Yak counter k-mer length"/>
208 <param name="trio_dual" argument="--trio-dual" type="boolean" truevalue="--trio-dual" falsevalue="" label="Utilize homology information to correct trio-phasing errors" /> 208 <param argument="--trio-dual" type="boolean" truevalue="--trio-dual" falsevalue="" label="Utilize homology information to correct trio-phasing errors"/>
209 </when> 209 </when>
210 </conditional> 210 </conditional>
211 <param name="filter_bits" argument="-f" type="integer" min="0" value="37" label="Bits for bloom filter" help="A value of 0 disables the bloom filter" /> 211 <param name="filter_bits" argument="-f" type="integer" min="0" value="37" label="Bits for bloom filter" help="A value of 0 disables the bloom filter"/>
212 <conditional name="assembly_options"> 212 <conditional name="assembly_options">
213 <param name="assembly_selector" type="select" label="Assembly options"> 213 <param name="assembly_selector" type="select" label="Assembly options">
214 <option value="blank">Leave default</option> 214 <option value="blank">Leave default</option>
215 <option value="set">Specify</option> 215 <option value="set">Specify</option>
216 </param> 216 </param>
217 <when value="blank" /> 217 <when value="blank"/>
218 <when value="set"> 218 <when value="set">
219 <param name="cleaning_rounds" argument="-a" type="integer" value="4" label="Cleaning rounds" /> 219 <param name="cleaning_rounds" argument="-a" type="integer" value="4" label="Cleaning rounds"/>
220 <param name="adapter_length" argument="-z" type="integer" min="0" value="0" label="Length of adapters to be removed" /> 220 <param name="adapter_length" argument="-z" type="integer" min="0" value="0" label="Length of adapters to be removed"/>
221 <param name="pop_contigs" argument="-m" type="integer" value="10000000" label="Minimum contig bubble size" help="Pop contig graph bubbles smaller than this value" /> 221 <param name="pop_contigs" argument="-m" type="integer" value="10000000" label="Minimum contig bubble size" help="Pop contig graph bubbles smaller than this value"/>
222 <param name="pop_unitigs" argument="-p" type="integer" value="100000" label="Minimum unitig bubble size" help="Pop unitig graph bubbles smaller than this value" /> 222 <param name="pop_unitigs" argument="-p" type="integer" value="100000" label="Minimum unitig bubble size" help="Pop unitig graph bubbles smaller than this value"/>
223 <param name="remove_tips" argument="-n" type="integer" value="3" label="Tip unitigs" help="Keep only tip unitigs with a number of reads greater than or equal to this value" /> 223 <param name="remove_tips" argument="-n" type="integer" value="3" label="Tip unitigs" help="Keep only tip unitigs with a number of reads greater than or equal to this value"/>
224 <param name="max_overlap" argument="-x" type="float" min="0" max="1" value="0.8" label="Maximum overlap drop ratio" help="This option is used with -r. Given a node N in the assembly graph, let max(N) be the length of the largest overlap of N. Hifiasm iteratively drops overlaps of N if their length/max(N) are below a threshold controlled by -x. Hifiasm applies -r rounds of short overlap removal with an increasing threshold between -x and -y"/> 224 <param name="max_overlap" argument="-x" type="float" min="0" max="1" value="0.8" label="Maximum overlap drop ratio" help="This option is used with -r. Given a node N in the assembly graph, let max(N) be the length of the largest overlap of N. Hifiasm iteratively drops overlaps of N if their length/max(N) are below a threshold controlled by -x. Hifiasm applies -r rounds of short overlap removal with an increasing threshold between -x and -y"/>
225 <param name="min_overlap" argument="-y" type="float" min="0" max="1" value="0.2" label="Minimum overlap drop ratio" help="This option is used with -r. Given a node N in the assembly graph, let max(N) be the length of the largest overlap of N. Hifiasm iteratively drops overlaps of N if their length/max(N) are over a threshold controlled by -y. Hifiasm applies -r rounds of short overlap removal with an increasing threshold between -x and -y"/> 225 <param name="min_overlap" argument="-y" type="float" min="0" max="1" value="0.2" label="Minimum overlap drop ratio" help="This option is used with -r. Given a node N in the assembly graph, let max(N) be the length of the largest overlap of N. Hifiasm iteratively drops overlaps of N if their length/max(N) are over a threshold controlled by -y. Hifiasm applies -r rounds of short overlap removal with an increasing threshold between -x and -y"/>
226 <param name="disable_post_join" argument="-u" type="boolean" truevalue="-u" falsevalue="" label="Skip post join contigs step" help="May improve N50" /> 226 <param name="disable_post_join" argument="-u" type="boolean" truevalue="-u" falsevalue="" label="Skip post join contigs step" help="May improve N50"/>
227 <param name="ignore_error_corrected" argument="-i" type="boolean" truevalue="-i" falsevalue="" value="False" label="Ignore error corrected reads and overlaps" help="Ignore error corrected reads and overlaps saved in prefix.*.bin files. Apart from assembly graphs, hifiasm also outputs three binary files that save alloverlap information during assembly step. With these files, hifiasm can avoid the time-consuming all-to-all overlap calculation step, and do the assembly directly and quickly. This might be helpful when users want to get an optimized assembly by multiple rounds of experiments with different parameters." /> 227 <param name="ignore_error_corrected" argument="-i" type="boolean" truevalue="-i" falsevalue="" value="False" label="Ignore error corrected reads and overlaps" help="Ignore error corrected reads and overlaps saved in prefix.*.bin files. Apart from assembly graphs, hifiasm also outputs three binary files that save alloverlap information during assembly step. With these files, hifiasm can avoid the time-consuming all-to-all overlap calculation step, and do the assembly directly and quickly. This might be helpful when users want to get an optimized assembly by multiple rounds of experiments with different parameters."/>
228 <param argument="--hom-cov" type="integer" optional="True" value="" label="Homozygous read coverage" /> 228 <param argument="--hom-cov" type="integer" optional="True" value="" label="Homozygous read coverage"/>
229 </when> 229 </when>
230 </conditional> 230 </conditional>
231 <conditional name="purge_options"> 231 <conditional name="purge_options">
232 <param name="purge_selector" type="select" label="Options for purging duplicates"> 232 <param name="purge_selector" type="select" label="Options for purging duplicates">
233 <option value="blank">Leave default</option> 233 <option value="blank">Leave default</option>
234 <option value="set">Specify</option> 234 <option value="set">Specify</option>
235 </param> 235 </param>
236 <when value="blank" /> 236 <when value="blank"/>
237 <when value="set"> 237 <when value="set">
238 <param name="purge_level" argument="-l" type="select" label="Purge level"> 238 <param name="purge_level" argument="-l" type="select" label="Purge level">
239 <option value="0" selected="true">None (0)</option> 239 <option value="0" selected="true">None (0)</option>
240 <option value="1">Light (1)</option> 240 <option value="1">Light (1)</option>
241 <option value="2">Aggressive (2)</option> 241 <option value="2">Aggressive (2)</option>
242 <option value="3">Aggressive - high heterozygosity rate (3)</option> 242 <option value="3">Aggressive - high heterozygosity rate (3)</option>
243 </param> 243 </param>
244 <param name="similarity_threshold" argument="-s" type="float" min="0" max="1" value="0.75" label="Similarity threshold for duplicate haplotigs" /> 244 <param name="similarity_threshold" argument="-s" type="float" min="0" max="1" value="0.75" label="Similarity threshold for duplicate haplotigs"/>
245 <param name="minimum_overlap" argument="-O" type="integer" value="1" label="Minimum overlapped reads for duplicate haplotigs" /> 245 <param name="minimum_overlap" argument="-O" type="integer" value="1" label="Minimum overlapped reads for duplicate haplotigs"/>
246 <param argument="--purge-max" type="integer" optional="true" label="Coverage upper bound" help="If not set, this will be determined automatically" /> 246 <param argument="--purge-max" type="integer" optional="true" label="Coverage upper bound" help="If not set, this will be determined automatically"/>
247 <param argument="--n-hap" type="integer" min="0" value="" optional="true" label="Assumtion of haplotype number" help="A haplotype is defined as the combination of alleles for different polymorphisms that occur on the same chromosome." /> 247 <param argument="--n-hap" type="integer" min="0" value="" optional="true" label="Assumtion of haplotype number" help="A haplotype is defined as the combination of alleles for different polymorphisms that occur on the same chromosome."/>
248 </when> 248 </when>
249 </conditional> 249 </conditional>
250 <conditional name="hic_partition"> 250 <conditional name="hic_partition">
251 <param name="hic_partition_selector" type="select" label="Options for Hi-C-partition"> 251 <param name="hic_partition_selector" type="select" label="Options for Hi-C-partition">
252 <option value="blank">Leave default</option> 252 <option value="blank">Leave default</option>
253 <option value="set">Specify</option> 253 <option value="set">Specify</option>
254 </param> 254 </param>
255 <when value="blank" /> 255 <when value="blank"/>
256 <when value="set"> 256 <when value="set">
257 <param argument="--h1" type="data" format="fastq,fastq.gz" multiple="true" label="Hi-C R1 reads" /> 257 <param argument="--h1" type="data" format="fastq,fastq.gz" multiple="true" label="Hi-C R1 reads"/>
258 <param argument="--h2" type="data" format="fastq,fastq.gz" multiple="true" label="Hi-C R2 reads" /> 258 <param argument="--h2" type="data" format="fastq,fastq.gz" multiple="true" label="Hi-C R2 reads"/>
259 <param argument="--seed" type="integer" min="1" value="" optional="true" label="RNG seed" /> 259 <param argument="--seed" type="integer" min="1" value="" optional="true" label="RNG seed"/>
260 <param argument="--n-weight" type="integer" min="1" value="" optional="true" label="Rounds of reweighting Hi-C links. Increasing this may improves phasing results but takes longer time" /> 260 <param argument="--n-weight" type="integer" min="1" value="" optional="true" label="Rounds of reweighting Hi-C links. Increasing this may improves phasing results but takes longer time"/>
261 <param argument="--n-perturb" type="integer" min="1" value="" optional="true" label="Rounds of perturbation. Increasing this may improves phasing results but takes longer time" /> 261 <param argument="--n-perturb" type="integer" min="1" value="" optional="true" label="Rounds of perturbation. Increasing this may improves phasing results but takes longer time"/>
262 <param argument="--f-perturb" type="float" min="0" max="1" value="" optional="true" label="Fraction to flip for perturbation. Increasing this may improves phasing results but takes longer time" /> 262 <param argument="--f-perturb" type="float" min="0" max="1" value="" optional="true" label="Fraction to flip for perturbation. Increasing this may improves phasing results but takes longer time"/>
263 <param argument="--l-msjoin" type="integer" min="0" value="500000" label="Detect misjoined unitigs of greater than or equal to specified size" help="A value of 0 disables this filter"/> 263 <param argument="--l-msjoin" type="integer" min="0" value="500000" label="Detect misjoined unitigs of greater than or equal to specified size" help="A value of 0 disables this filter"/>
264 </when> 264 </when>
265 </conditional> 265 </conditional>
266 <conditional name="ont_integration"> 266 <conditional name="ont_integration">
267 <param name="ont_integration_selector" type="select" label="Options for ONT integration"> 267 <param name="ont_integration_selector" type="select" label="Options for ONT integration">
268 <option value="blank">Leave default</option> 268 <option value="blank">Leave default</option>
269 <option value="set">Specify</option> 269 <option value="set">Specify</option>
270 </param> 270 </param>
271 <when value="blank" /> 271 <when value="blank"/>
272 <when value="set"> 272 <when value="set">
273 <param argument="--ul" type="data" format="fastq,fastq.gz,fasta,fasta.gz" multiple="true" label="Ultra-long (>100kb) reads" /> 273 <param argument="--ul" type="data" format="fastq,fastq.gz,fasta,fasta.gz" multiple="true" label="Ultra-long (>100kb) reads"/>
274 <param argument="--ul-rate" type="float" value="0.2" optional="true" label="error rate of ultra-long reads (default 0.2)" /> 274 <param argument="--ul-rate" type="float" value="0.2" optional="true" label="error rate of ultra-long reads (default 0.2)"/>
275 <param argument="--ul-tip" type="integer" value="6" optional="true" label="remove tip unitigs composed of less than (INT) reads for the UL assembly (default 6)" /> 275 <param argument="--ul-tip" type="integer" value="6" optional="true" label="remove tip unitigs composed of less than (INT) reads for the UL assembly (default 6)"/>
276 </when> 276 </when>
277 </conditional> 277 </conditional>
278 <conditional name="advanced_options"> 278 <conditional name="advanced_options">
279 <param name="advanced_selector" type="select" label="Advanced options"> 279 <param name="advanced_selector" type="select" label="Advanced options">
280 <option value="blank">Leave default</option> 280 <option value="blank">Leave default</option>
281 <option value="set">Specify</option> 281 <option value="set">Specify</option>
282 </param> 282 </param>
283 <when value="blank" /> 283 <when value="blank"/>
284 <when value="set"> 284 <when value="set">
285 <param name="hifiasm_kmer_length" argument="-k" type="integer" min="0" max="64" value="51" label="Hifiasm k-mer length" /> 285 <param name="hifiasm_kmer_length" argument="-k" type="integer" min="0" max="64" value="51" label="Hifiasm k-mer length"/>
286 <param name="window_size" argument="-w" type="integer" min="0" value="51" label="Minimizer window size" /> 286 <param name="window_size" argument="-w" type="integer" min="0" value="51" label="Minimizer window size"/>
287 <param name="drop_kmers" argument="-D" type="float" value="5.0" label="Drop k-mers" help="K-mers that occur more than this value multiplied by the coverage will be discarded" /> 287 <param name="drop_kmers" argument="-D" type="float" value="5.0" label="Drop k-mers" help="K-mers that occur more than this value multiplied by the coverage will be discarded"/>
288 <param name="max_overlaps" argument="-N" type="integer" value="100" label="Maximum overlaps to consider" help="The software selects the larger of this value and the k-mer count multiplied by coverage" /> 288 <param name="max_overlaps" argument="-N" type="integer" value="100" label="Maximum overlaps to consider" help="The software selects the larger of this value and the k-mer count multiplied by coverage"/>
289 <param name="correction_rounds" argument="-r" type="integer" value="3" label="Correction rounds" /> 289 <param name="correction_rounds" argument="-r" type="integer" value="3" label="Correction rounds"/>
290 <param argument="--min-hist-cnt" type="integer" min="0" value="" optional="true" label="Minimum count threshold" help="When analyzing the k-mer spectrum, ignore counts below this value" /> 290 <param argument="--min-hist-cnt" type="integer" min="0" value="" optional="true" label="Minimum count threshold" help="When analyzing the k-mer spectrum, ignore counts below this value"/>
291 <param argument="--max-kocc" type="integer" min="0" value="20000" label="Maximum k-mer ocurrence" help="Employ k-mers occurring less than INT times to rescue repetitive overlaps" /> 291 <param argument="--max-kocc" type="integer" min="0" value="20000" label="Maximum k-mer ocurrence" help="Employ k-mers occurring less than INT times to rescue repetitive overlaps"/>
292 <param argument="--hg-size" type="text" value="" optional="true" label="Estimated haploid genome size" 292 <param argument="--hg-size" type="text" value="" optional="true" label="Estimated haploid genome size"
293 help="Estimated haploid genome size used for inferring read coverage. If not provided, this parameter will be infered by hifism. Common suffices are required, for example, 100m or 3g"> 293 help="Estimated haploid genome size used for inferring read coverage. If not provided, this parameter will be infered by hifism. Common suffices are required, for example, 100m or 3g">
294 <sanitizer invalid_char=""> 294 <sanitizer invalid_char="">
295 <valid initial="string.digits"> 295 <valid initial="string.digits">
296 <add value="k" /> 296 <add value="k"/>
297 <add value="K" /> 297 <add value="K"/>
298 <add value="m" /> 298 <add value="m"/>
299 <add value="M" /> 299 <add value="M"/>
300 <add value="G" /> 300 <add value="G"/>
301 <add value="g" /> 301 <add value="g"/>
302 </valid> 302 </valid>
303 </sanitizer> 303 </sanitizer>
304 <validator type="regex">[0-9kKmMGg]+</validator> 304 <validator type="regex">[0-9kKmMGg]+</validator>
305 </param> 305 </param>
306 </when> 306 </when>
308 <conditional name="scaffolding_options"> 308 <conditional name="scaffolding_options">
309 <param name="scaffold_selector" type="select" label="Options for hifiasm scaffolding"> 309 <param name="scaffold_selector" type="select" label="Options for hifiasm scaffolding">
310 <option value="blank">No scaffolding (default)</option> 310 <option value="blank">No scaffolding (default)</option>
311 <option value="set">Turn on hifiasm scaffolding</option> 311 <option value="set">Turn on hifiasm scaffolding</option>
312 </param> 312 </param>
313 <when value="blank" /> 313 <when value="blank"/>
314 <when value="set"> 314 <when value="set">
315 <param argument="--scaf-gap" type="integer" min="1" value="3000000" optional="true" label="Max scaffolding gap size" help="Set the max gap size of scaffolds (default is 3,000,000)" /> 315 <param argument="--scaf-gap" type="integer" min="1" value="3000000" optional="true" label="Max scaffolding gap size" help="Set the max gap size of scaffolds (default is 3,000,000)"/>
316 </when> 316 </when>
317 </conditional> 317 </conditional>
318 <param name="log_out" type="boolean" label="Output log file?" truevalue="yes" falsevalue="no" /> 318 <param name="log_out" type="boolean" label="Output log file?" truevalue="yes" falsevalue="no"/>
319 <param name="bins_out" type="boolean" label="Output .bin files (used for development and debugging)?" truevalue="yes" falsevalue="no" /> 319 <param name="bins_out" type="boolean" label="Output .bin files (used for development and debugging)?" truevalue="yes" falsevalue="no"/>
320 </inputs> 320 </inputs>
321 <outputs> 321 <outputs>
322 <!--Standard mode--> 322 <!--Standard mode-->
323 <data name="raw_unitigs" format="gfa1" from_work_dir="output.r_utg.gfa" label="${tool.name} on ${on_string}: haplotype-resolved raw unitig graph for pseudohaplotype assembly"> 323 <data name="raw_unitigs" format="gfa1" from_work_dir="output.r_utg.gfa" label="${tool.name} on ${on_string}: haplotype-resolved raw unitig graph for pseudohaplotype assembly">
324 <filter>mode['mode_selector'] == 'standard' and hic_partition['hic_partition_selector'] == 'blank'</filter> 324 <filter>mode['mode_selector'] == 'standard' and hic_partition['hic_partition_selector'] == 'blank'</filter>
371 <!--Log, noseq, and bin output--> 371 <!--Log, noseq, and bin output-->
372 <data name="log_file" format="txt" from_work_dir="output.log" label="${tool.name} ${on_string}: log file"> 372 <data name="log_file" format="txt" from_work_dir="output.log" label="${tool.name} ${on_string}: log file">
373 <filter>log_out</filter> 373 <filter>log_out</filter>
374 </data> 374 </data>
375 <collection name="noseq_files" type="list" label="${tool.name} on ${on_string}: noseq files"> 375 <collection name="noseq_files" type="list" label="${tool.name} on ${on_string}: noseq files">
376 <discover_datasets pattern="__name_and_ext__" format="gfa1" directory="noseq_files" /> 376 <discover_datasets pattern="__name_and_ext__" format="gfa1" directory="noseq_files"/>
377 </collection> 377 </collection>
378 <collection name="bin_files" type="list" label="${tool.name} on ${on_string}: bin files"> 378 <collection name="bin_files" type="list" label="${tool.name} on ${on_string}: bin files">
379 <filter>bins_out</filter> 379 <filter>bins_out</filter>
380 <discover_datasets pattern="__name_and_ext__" format="gfa1" directory="bin_files" /> 380 <discover_datasets pattern="__name_and_ext__" format="gfa1" directory="bin_files"/>
381 </collection> 381 </collection>
382 </outputs> 382 </outputs>
383 <tests> 383 <tests>
384 <!-- TEST 1 --> 384 <!-- TEST 1 -->
385 <test expect_num_outputs="5"> 385 <test expect_num_outputs="5">
386 <param name="reads" value="hifiasm-in1.fa.gz" ftype="fasta.gz" /> 386 <param name="reads" value="hifiasm-in1.fa.gz" ftype="fasta.gz"/>
387 <param name="filter_bits" value="0" /> 387 <param name="filter_bits" value="0"/>
388 <param name="mode_selector" value="standard" /> 388 <param name="mode_selector" value="standard"/>
389 <output name="raw_unitigs" file="hifiasm-out1-raw.gfa" ftype="gfa1" /> 389 <output name="raw_unitigs" file="hifiasm-out1-raw.gfa" ftype="gfa1"/>
390 <output name="processed_unitigs" file="hifiasm-out1-processed.gfa" ftype="gfa1" /> 390 <output name="processed_unitigs" file="hifiasm-out1-processed.gfa" ftype="gfa1"/>
391 <output name="primary_contig_graph" file="hifiasm-out1-primary.gfa" ftype="gfa1" /> 391 <output name="primary_contig_graph" file="hifiasm-out1-primary.gfa" ftype="gfa1"/>
392 </test> 392 </test>
393 <!-- TEST 2 --> 393 <!-- TEST 2 -->
394 <test expect_num_outputs="5"> 394 <test expect_num_outputs="5">
395 <param name="reads" value="hifiasm-in2-0.fa.gz,hifiasm-in2-1.fa.gz,hifiasm-in2-2.fa.gz,hifiasm-in2-3.fa.gz,hifiasm-in2-4.fa.gz" ftype="fasta.gz" /> 395 <param name="reads" value="hifiasm-in2-0.fa.gz,hifiasm-in2-1.fa.gz,hifiasm-in2-2.fa.gz,hifiasm-in2-3.fa.gz,hifiasm-in2-4.fa.gz" ftype="fasta.gz"/>
396 <param name="filter_bits" value="0" /> 396 <param name="filter_bits" value="0"/>
397 <param name="mode_selector" value="standard" /> 397 <param name="mode_selector" value="standard"/>
398 <output name="raw_unitigs" file="hifiasm-out2-raw.gfa" ftype="gfa1" /> 398 <output name="raw_unitigs" file="hifiasm-out2-raw.gfa" ftype="gfa1"/>
399 <output name="processed_unitigs" file="hifiasm-out2-processed.gfa" ftype="gfa1" /> 399 <output name="processed_unitigs" file="hifiasm-out2-processed.gfa" ftype="gfa1"/>
400 <output name="primary_contig_graph" file="hifiasm-out2-primary.gfa" ftype="gfa1" /> 400 <output name="primary_contig_graph" file="hifiasm-out2-primary.gfa" ftype="gfa1"/>
401 </test> 401 </test>
402 <!-- TEST 3: Test logfile out--> 402 <!-- TEST 3: Test logfile out-->
403 <test expect_num_outputs="6"> 403 <test expect_num_outputs="6">
404 <param name="reads" value="hifiasm-in1.fa.gz" ftype="fasta.gz" /> 404 <param name="reads" value="hifiasm-in1.fa.gz" ftype="fasta.gz"/>
405 <param name="filter_bits" value="0" /> 405 <param name="filter_bits" value="0"/>
406 <param name="mode_selector" value="standard" /> 406 <param name="mode_selector" value="standard"/>
407 <param name="log_out" value="yes"/> 407 <param name="log_out" value="yes"/>
408 <output name="raw_unitigs" file="hifiasm-out1-raw.gfa" ftype="gfa1" /> 408 <output name="raw_unitigs" file="hifiasm-out1-raw.gfa" ftype="gfa1"/>
409 <output name="processed_unitigs" file="hifiasm-out1-processed.gfa" ftype="gfa1" /> 409 <output name="processed_unitigs" file="hifiasm-out1-processed.gfa" ftype="gfa1"/>
410 <output name="primary_contig_graph" file="hifiasm-out1-primary.gfa" ftype="gfa1" /> 410 <output name="primary_contig_graph" file="hifiasm-out1-primary.gfa" ftype="gfa1"/>
411 <output name="log_file" ftype="txt"> 411 <output name="log_file" ftype="txt">
412 <assert_contents> 412 <assert_contents>
413 <has_text text="-o output -f 0 --primary input_0.fasta.gz"/> 413 <has_text text="-o output -f 0 --primary input_0.fasta.gz"/>
414 </assert_contents> 414 </assert_contents>
415 </output> 415 </output>
416 </test> 416 </test>
417 <!--TEST 4: Test Hi-C reads--> 417 <!--TEST 4: Test Hi-C reads-->
418 <test expect_num_outputs="6"> 418 <test expect_num_outputs="6">
419 <param name="reads" value="hifiasm-in1.fa.gz" ftype="fasta.gz" /> 419 <param name="reads" value="hifiasm-in1.fa.gz" ftype="fasta.gz"/>
420 <param name="filter_bits" value="0" /> 420 <param name="filter_bits" value="0"/>
421 <param name="mode_selector" value="standard" /> 421 <param name="mode_selector" value="standard"/>
422 <conditional name="hic_partition"> 422 <conditional name="hic_partition">
423 <param name="hic_partition_selector" value="set"/> 423 <param name="hic_partition_selector" value="set"/>
424 <param name="h1" value="hic_1.fastq.gz"/> 424 <param name="h1" value="hic_1.fastq.gz"/>
425 <param name="h2" value="hic_2.fastq.gz"/> 425 <param name="h2" value="hic_2.fastq.gz"/>
426 <param name="n_weight" value="1"/> 426 <param name="n_weight" value="1"/>
427 <param name="n_perturb" value="1"/> 427 <param name="n_perturb" value="1"/>
428 <param name="l_perturb" value="0"/> 428 <param name="l_perturb" value="0"/>
429 <param name="l_msjoin" value="0"/> 429 <param name="l_msjoin" value="0"/>
430 </conditional> 430 </conditional>
431 <output name="hic_pcontig_graph" file="hifiasm-out-hifi-p.gfa" ftype="gfa1" /> 431 <output name="hic_pcontig_graph" file="hifiasm-out-hifi-p.gfa" ftype="gfa1"/>
432 <output name="hic_acontig_graph" file="hifiasm-out-hifi-a.gfa" ftype="gfa1" /> 432 <output name="hic_acontig_graph" file="hifiasm-out-hifi-a.gfa" ftype="gfa1"/>
433 <output name="hic_balanced_contig_hap1_graph" ftype="gfa1"> 433 <output name="hic_balanced_contig_hap1_graph" ftype="gfa1">
434 <assert_contents> 434 <assert_contents>
435 <has_text_matching expression="^S" /> 435 <has_text_matching expression="^S"/>
436 <has_size value="83914" delta="500"/> 436 <has_size value="83914" delta="500"/>
437 </assert_contents> 437 </assert_contents>
438 </output> 438 </output>
439 <output name="hic_balanced_contig_hap1_graph" ftype="gfa1" > 439 <output name="hic_balanced_contig_hap1_graph" ftype="gfa1" >
440 <assert_contents> 440 <assert_contents>
441 <has_text_matching expression="^S" /> 441 <has_text_matching expression="^S"/>
442 <has_size value="83914" delta="500"/> 442 <has_size value="83914" delta="500"/>
443 </assert_contents> 443 </assert_contents>
444 </output> 444 </output>
445 <output name="hic_raw_initig" ftype="gfa1" > 445 <output name="hic_raw_initig" ftype="gfa1" >
446 <assert_contents> 446 <assert_contents>
447 <has_text_matching expression="^S" /> 447 <has_text_matching expression="^S"/>
448 <has_size value="83904" delta="500"/> 448 <has_size value="83904" delta="500"/>
449 </assert_contents> 449 </assert_contents>
450 </output> 450 </output>
451 </test> 451 </test>
452 <!-- TEST 5: Test trio mode --> 452 <!-- TEST 5: Test trio mode -->
462 </conditional> 462 </conditional>
463 <param name="max_kmers" value="2"/> 463 <param name="max_kmers" value="2"/>
464 <param name="min_kmers" value="5"/> 464 <param name="min_kmers" value="5"/>
465 </conditional> 465 </conditional>
466 <assert_command> 466 <assert_command>
467 <has_text text="-1 hap1.yak" /> 467 <has_text text="-1 hap1.yak"/>
468 <has_text text="-2 hap2.yak" /> 468 <has_text text="-2 hap2.yak"/>
469 <has_text text="--primary" /> 469 <has_text text="--primary"/>
470 </assert_command> 470 </assert_command>
471 </test> 471 </test>
472 <!-- TEST 6: Test ignore-error-corrected option --> 472 <!-- TEST 6: Test ignore-error-corrected option -->
473 <test expect_num_outputs="5"> 473 <test expect_num_outputs="5">
474 <param name="reads" value="hifiasm-in1.fa.gz" ftype="fasta.gz" /> 474 <param name="reads" value="hifiasm-in1.fa.gz" ftype="fasta.gz"/>
475 <param name="filter_bits" value="0" /> 475 <param name="filter_bits" value="0"/>
476 <param name="mode_selector" value="standard" /> 476 <param name="mode_selector" value="standard"/>
477 <conditional name="assembly_options"> 477 <conditional name="assembly_options">
478 <param name="assembly_selector" value="set"/> 478 <param name="assembly_selector" value="set"/>
479 <param name="ignore_error_corrected" value="True"/> 479 <param name="ignore_error_corrected" value="True"/>
480 </conditional> 480 </conditional>
481 <output name="raw_unitigs" file="hifiasm-out3-raw.gfa" ftype="gfa1" /> 481 <output name="raw_unitigs" file="hifiasm-out3-raw.gfa" ftype="gfa1"/>
482 <output name="processed_unitigs" file="hifiasm-out3-processed.gfa" ftype="gfa1" /> 482 <output name="processed_unitigs" file="hifiasm-out3-processed.gfa" ftype="gfa1"/>
483 <output name="primary_contig_graph" file="hifiasm-out3-primary.gfa" ftype="gfa1" /> 483 <output name="primary_contig_graph" file="hifiasm-out3-primary.gfa" ftype="gfa1"/>
484 </test> 484 </test>
485 <!-- TEST 7: Test expected haplotype number --> 485 <!-- TEST 7: Test expected haplotype number -->
486 <test expect_num_outputs="5"> 486 <test expect_num_outputs="5">
487 <param name="reads" value="hifiasm-in1.fa.gz" ftype="fasta.gz" /> 487 <param name="reads" value="hifiasm-in1.fa.gz" ftype="fasta.gz"/>
488 <param name="filter_bits" value="0" /> 488 <param name="filter_bits" value="0"/>
489 <param name="mode_selector" value="standard" /> 489 <param name="mode_selector" value="standard"/>
490 <conditional name="purge_options"> 490 <conditional name="purge_options">
491 <param name="purge_selector" value="set"/> 491 <param name="purge_selector" value="set"/>
492 <param name="n_hap" value="1"/> 492 <param name="n_hap" value="1"/>
493 </conditional> 493 </conditional>
494 <output name="raw_unitigs" file="hifiasm-out4-raw.gfa" ftype="gfa1" /> 494 <output name="raw_unitigs" file="hifiasm-out4-raw.gfa" ftype="gfa1"/>
495 <output name="processed_unitigs" file="hifiasm-out4-processed.gfa" ftype="gfa1" /> 495 <output name="processed_unitigs" file="hifiasm-out4-processed.gfa" ftype="gfa1"/>
496 <output name="primary_contig_graph" file="hifiasm-out4-primary.gfa" ftype="gfa1" /> 496 <output name="primary_contig_graph" file="hifiasm-out4-primary.gfa" ftype="gfa1"/>
497 </test> 497 </test>
498 <!-- TEST 8: Test min_hist_cnt option --> 498 <!-- TEST 8: Test min_hist_cnt option -->
499 <test expect_num_outputs="5"> 499 <test expect_num_outputs="5">
500 <param name="reads" value="hifiasm-in1.fa.gz" ftype="fasta.gz" /> 500 <param name="reads" value="hifiasm-in1.fa.gz" ftype="fasta.gz"/>
501 <param name="filter_bits" value="0" /> 501 <param name="filter_bits" value="0"/>
502 <param name="mode_selector" value="standard" /> 502 <param name="mode_selector" value="standard"/>
503 <conditional name="advanced_options"> 503 <conditional name="advanced_options">
504 <param name="advanced_selector" value="set"/> 504 <param name="advanced_selector" value="set"/>
505 <param name="min_hist_cnt" value="1"/> 505 <param name="min_hist_cnt" value="1"/>
506 </conditional> 506 </conditional>
507 <output name="raw_unitigs" file="hifiasm-out5-raw.gfa" ftype="gfa1" /> 507 <output name="raw_unitigs" file="hifiasm-out5-raw.gfa" ftype="gfa1"/>
508 <output name="processed_unitigs" file="hifiasm-out5-processed.gfa" ftype="gfa1" /> 508 <output name="processed_unitigs" file="hifiasm-out5-processed.gfa" ftype="gfa1"/>
509 <output name="primary_contig_graph" file="hifiasm-out5-primary.gfa" ftype="gfa1" /> 509 <output name="primary_contig_graph" file="hifiasm-out5-primary.gfa" ftype="gfa1"/>
510 </test> 510 </test>
511 <!-- TEST 9: Test max_kooc option --> 511 <!-- TEST 9: Test max_kooc option -->
512 <test expect_num_outputs="5"> 512 <test expect_num_outputs="5">
513 <param name="reads" value="hifiasm-in1.fa.gz" ftype="fasta.gz" /> 513 <param name="reads" value="hifiasm-in1.fa.gz" ftype="fasta.gz"/>
514 <param name="filter_bits" value="0" /> 514 <param name="filter_bits" value="0"/>
515 <param name="mode_selector" value="standard" /> 515 <param name="mode_selector" value="standard"/>
516 <conditional name="advanced_options"> 516 <conditional name="advanced_options">
517 <param name="advanced_selector" value="set"/> 517 <param name="advanced_selector" value="set"/>
518 <param name="max_kooc" value="21000"/> 518 <param name="max_kooc" value="21000"/>
519 </conditional> 519 </conditional>
520 <output name="raw_unitigs" file="hifiasm-out6-raw.gfa" ftype="gfa1" /> 520 <output name="raw_unitigs" file="hifiasm-out6-raw.gfa" ftype="gfa1"/>
521 <output name="processed_unitigs" file="hifiasm-out6-processed.gfa" ftype="gfa1" /> 521 <output name="processed_unitigs" file="hifiasm-out6-processed.gfa" ftype="gfa1"/>
522 <output name="primary_contig_graph" file="hifiasm-out6-primary.gfa" ftype="gfa1" /> 522 <output name="primary_contig_graph" file="hifiasm-out6-primary.gfa" ftype="gfa1"/>
523 </test> 523 </test>
524 <!-- TEST 10: Test hg-size option --> 524 <!-- TEST 10: Test hg-size option -->
525 <test expect_num_outputs="5"> 525 <test expect_num_outputs="5">
526 <param name="reads" value="hifiasm-in1.fa.gz" ftype="fasta.gz" /> 526 <param name="reads" value="hifiasm-in1.fa.gz" ftype="fasta.gz"/>
527 <param name="filter_bits" value="0" /> 527 <param name="filter_bits" value="0"/>
528 <param name="mode_selector" value="standard" /> 528 <param name="mode_selector" value="standard"/>
529 <conditional name="advanced_options"> 529 <conditional name="advanced_options">
530 <param name="advanced_selector" value="set"/> 530 <param name="advanced_selector" value="set"/>
531 <param name="hg_size" value="1k"/> 531 <param name="hg_size" value="1k"/>
532 </conditional> 532 </conditional>
533 <output name="raw_unitigs" file="hifiasm-out7-raw.gfa" ftype="gfa1" /> 533 <output name="raw_unitigs" file="hifiasm-out7-raw.gfa" ftype="gfa1"/>
534 <output name="processed_unitigs" file="hifiasm-out7-processed.gfa" ftype="gfa1" /> 534 <output name="processed_unitigs" file="hifiasm-out7-processed.gfa" ftype="gfa1"/>
535 <output name="primary_contig_graph" file="hifiasm-out7-primary.gfa" ftype="gfa1" /> 535 <output name="primary_contig_graph" file="hifiasm-out7-primary.gfa" ftype="gfa1"/>
536 </test> 536 </test>
537 <!-- TEST 11: Test ignore-error-corrected option --> 537 <!-- TEST 11: Test ignore-error-corrected option -->
538 <test expect_num_outputs="5"> 538 <test expect_num_outputs="5">
539 <param name="reads" value="hifiasm-in1.fa.gz" ftype="fasta.gz" /> 539 <param name="reads" value="hifiasm-in1.fa.gz" ftype="fasta.gz"/>
540 <param name="filter_bits" value="0" /> 540 <param name="filter_bits" value="0"/>
541 <param name="mode_selector" value="standard" /> 541 <param name="mode_selector" value="standard"/>
542 <conditional name="assembly_options"> 542 <conditional name="assembly_options">
543 <param name="assembly_selector" value="set"/> 543 <param name="assembly_selector" value="set"/>
544 <param name="hom-cov" value="1000"/> 544 <param name="hom-cov" value="1000"/>
545 </conditional> 545 </conditional>
546 <output name="raw_unitigs" file="hifiasm-out8-raw.gfa" ftype="gfa1" /> 546 <output name="raw_unitigs" file="hifiasm-out8-raw.gfa" ftype="gfa1"/>
547 <output name="processed_unitigs" file="hifiasm-out8-processed.gfa" ftype="gfa1" /> 547 <output name="processed_unitigs" file="hifiasm-out8-processed.gfa" ftype="gfa1"/>
548 <output name="primary_contig_graph" file="hifiasm-out8-primary.gfa" ftype="gfa1" /> 548 <output name="primary_contig_graph" file="hifiasm-out8-primary.gfa" ftype="gfa1"/>
549 </test> 549 </test>
550 <!-- TEST 12: test nanopore input --> 550 <!-- TEST 12: test nanopore input -->
551 <test expect_num_outputs="5"> 551 <test expect_num_outputs="5">
552 <param name="reads" value="hifiasm-in1.fa.gz" ftype="fasta.gz" /> 552 <param name="reads" value="hifiasm-in1.fa.gz" ftype="fasta.gz"/>
553 <param name="mode_selector" value="standard" /> 553 <param name="mode_selector" value="standard"/>
554 <param name="filter_bits" value="0" /> 554 <param name="filter_bits" value="0"/>
555 <conditional name="ont_integration"> 555 <conditional name="ont_integration">
556 <param name="ont_integration_selector" value="set" /> 556 <param name="ont_integration_selector" value="set"/>
557 <param name="ul" value="nanopore.fasta.gz" /> 557 <param name="ul" value="nanopore.fasta.gz"/>
558 <param name="ul_tip" value="1" /> 558 <param name="ul_tip" value="1"/>
559 </conditional> 559 </conditional>
560 <output name="primary_contig_graph" file="hifiasm-out11-primary.gfa" ftype="gfa1" /> 560 <output name="primary_contig_graph" file="hifiasm-out11-primary.gfa" ftype="gfa1"/>
561 <output name="raw_unitigs" file="hifiasm-out11-raw.gfa" ftype="gfa1" /> 561 <output name="raw_unitigs" file="hifiasm-out11-raw.gfa" ftype="gfa1"/>
562 </test> 562 </test>
563 <!-- TEST 13: test multi-file nanopore input --> 563 <!-- TEST 13: test multi-file nanopore input -->
564 <test expect_num_outputs="6"> 564 <test expect_num_outputs="6">
565 <param name="reads" value="hifiasm-in1.fa.gz" ftype="fasta.gz" /> 565 <param name="reads" value="hifiasm-in1.fa.gz" ftype="fasta.gz"/>
566 <param name="mode_selector" value="standard" /> 566 <param name="mode_selector" value="standard"/>
567 <param name="filter_bits" value="0" /> 567 <param name="filter_bits" value="0"/>
568 <param name="log_out" value="yes" /> 568 <param name="log_out" value="yes"/>
569 <conditional name="ont_integration"> 569 <conditional name="ont_integration">
570 <param name="ont_integration_selector" value="set" /> 570 <param name="ont_integration_selector" value="set"/>
571 <param name="ul" value="nanopore.fasta.gz,nanopore.fasta.gz" /> 571 <param name="ul" value="nanopore.fasta.gz,nanopore.fasta.gz"/>
572 <param name="ul_tip" value="1" /> 572 <param name="ul_tip" value="1"/>
573 </conditional> 573 </conditional>
574 <output name="log_file" ftype="txt"> 574 <output name="log_file" ftype="txt">
575 <assert_contents> 575 <assert_contents>
576 <has_text text="--ul ./ultralong/input_0.fasta.gz,./ultralong/input_1.fasta.gz"/> 576 <has_text text="--ul ./ultralong/input_0.fasta.gz,./ultralong/input_1.fasta.gz"/>
577 </assert_contents> 577 </assert_contents>
578 </output> 578 </output>
579 </test> 579 </test>
580 <!-- TEST 14: test bin files --> 580 <!-- TEST 14: test bin files -->
581 <test expect_num_outputs="6"> 581 <test expect_num_outputs="6">
582 <param name="reads" value="hifiasm-in1.fa.gz" ftype="fasta.gz" /> 582 <param name="reads" value="hifiasm-in1.fa.gz" ftype="fasta.gz"/>
583 <param name="filter_bits" value="0" /> 583 <param name="filter_bits" value="0"/>
584 <param name="mode_selector" value="standard" /> 584 <param name="mode_selector" value="standard"/>
585 <param name="bins_out" value="yes" /> 585 <param name="bins_out" value="yes"/>
586 <output_collection name="bin_files" type="list" count="3" /> 586 <output_collection name="bin_files" type="list" count="3"/>
587 </test> 587 </test>
588 <!-- TEST 15: Test trio LIST mode --> 588 <!-- TEST 15: Test trio LIST mode -->
589 <test expect_num_outputs="6"> 589 <test expect_num_outputs="6">
590 <param name="filter_bits" value="0"/> 590 <param name="filter_bits" value="0"/>
591 <param name="log_out" value="yes"/> 591 <param name="log_out" value="yes"/>