Mercurial > repos > iuc > hal_halstats

comparison hal_halStats.xml @ 0:25dcde5bf94e draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/main/tools/haltools commit 6244b9d15a5ad97ae20191e2f8fbafe2050c3cac
author iuc
date Fri, 06 Feb 2026 10:39:34 +0000
parents
children
comparison
equal deleted inserted replaced
-1:000000000000 0:25dcde5bf94e
1 <tool id="hal_halstats" name="halStats" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="@PROFILE@">
2 <description>retrieves basic statistics from a HAL file</description>
3 <macros>
4 <import>macros.xml</import>
5 </macros>
6 <expand macro="requirements"/>
7 <expand macro="stdio"/>
8 <command detect_errors="aggressive"><![CDATA[
9 set -o pipefail; ## Sets the pipeline’s exit code to halStats’s on failure.
10 ( ## echo headers for specific numerical data
11 #if $mode.option == '--baseComp':
12 echo -e 'fraction_of_As\tfraction_of_Gs\tfraction_of_Cs\tfraction_of_Ts';
13 #else if $mode.option == '--numSegments':
14 echo -e 'numTopSegments\tnumBottomSegments';
15 #end if
16 halStats
17 #if $mode.option == '--allCoverage':
18 --allCoverage
19 #else if $mode.option == '--branches':
20 --branches
21 #else if $mode.option == '--genomes':
22 --genomes
23 #else if $mode.option == '--metaData':
24 --metaData
25 #else if $mode.option == '--root':
26 --root
27 #else if $mode.option == '--tree':
28 --tree
29 #else if $mode.option == '--baseComp':
30 --baseComp '$mode.baseComp'
31 #else if $mode.option == '--bedSequences':
32 --bedSequences '$mode.bedSequences'
33 #else if $mode.option == '--bottomSegments':
34 --bottomSegments '$mode.bottomSegments'
35 #else if $mode.option == '--branchLength':
36 --branchLength '$mode.branchLength'
37 #else if $mode.option == '--children':
38 --children '$mode.children'
39 #else if $mode.option == '--chromSizes':
40 --chromSizes '$mode.chromSizes'
41 #else if $mode.option == '--coverage':
42 --coverage '$mode.coverage'
43 #else if $mode.option == '--genomeMetaData':
44 --genomeMetaData '$mode.genomeMetaData'
45 #else if $mode.option == '--numSegments':
46 --numSegments '$mode.numSegments'
47 #else if $mode.option == '--parent':
48 --parent '$mode.parent'
49 #else if $mode.option == '--percentID':
50 --percentID '$mode.percentID'
51 #else if $mode.option == '--sequenceStats':
52 --sequenceStats '$mode.sequenceStats'
53 #else if $mode.option == '--sequences':
54 --sequences '$mode.sequences'
55 #else if $mode.option == '--span':
56 --span '$mode.span'
57 #else if $mode.option == '--spanRoot':
58 --spanRoot '$mode.spanRoot'
59 #else if $mode.option == '--topSegments':
60 --topSegments '$mode.topSegments'
61 #end if
62 '$input_hal'
63 ## Pipes specific output to replace commas with tabs. Output is mostly numerical, and Genome names contain no commas, as this would invalidate the HAL Newick tree.
64 #if $mode.option == '--allCoverage' or $mode.option == '--sequenceStats' or $mode.option == '--percentID' or $mode.option == '--coverage':
65 | tr ',' '\t'
66 #else if $mode.option == '--numSegments':
67 | tr ' ' '\t' ## Replace spaces
68 #end if
69 ) > '$out_file'
70 ]]></command>
71 <inputs>
72 <expand macro="input_hal"/>
73 <conditional name="mode">
74 <param name="option" type="select" label="Select the type of statistics you are interested in">
75 <option value="" selected="true">Basic overview</option>
76 <option value="--genomes">List of genomes in alignment (--genomes)</option>
77 <option value="--sequences">List of sequences in a given genome (--sequences)</option>
78 <option value="--bedSequences">List of sequences in a given genome (in BED format) (--bedSequences)</option>
79 <option value="--sequenceStats">Stats for each sequence in a given genome (--sequenceStats)</option>
80 <option value="--tree">Newick tree (--tree)</option>
81 <option value="--branches">List of branches specified by the child genome (--branches)</option>
82 <option value="--span">Branches on path (or spanning tree) between given list of genomes (--span)</option>
83 <option value="--spanRoot">Genomes on path (or spanning tree) with spanning tree root between given list of genomes (--spanRoot)</option>
84 <option value="--children">Names of children of a given genome (--children)</option>
85 <option value="--root">Root genome name (--root)</option>
86 <option value="--parent">Parent name of a given genome (--parent)</option>
87 <option value="--branchLength">Branch length between a given genome and its parent (--branchLength)</option>
88 <option value="--numSegments">Number of top and of bottom segments of a given genome (--numSegments)</option>
89 <option value="--topSegments">Coordinates of all top segments of a given genome (in BED format) (--topSegments)</option>
90 <option value="--bottomSegments">Coordinates of all bottom segments of a given genome (in BED format) (--bottomSegments)</option>
91 <option value="--baseComp">Base composition by sampling every step bases (--baseComp)</option>
92 <option value="--genomeMetaData">Metadata for a given genome (--genomeMetaData)</option>
93 <option value="--metaData">Metadata for the entire alignment (--metaData)</option>
94 <option value="--chromSizes">Name and length of each sequence in a given genome (in format used by wigToBigWig) (--chromSizes)</option>
95 <option value="--percentID">Percent ID of a given genome with all other genomes (--percentID)</option>
96 <option value="--coverage">Histogram of coverage of a given genome with all genomes (--coverage)</option>
97 <option value="--allCoverage">Histogram of coverage from all genomes to all genomes (--allCoverage)</option>
98 </param>
99 <when value=""/>
100 <when value="--allCoverage"/>
101 <when value="--baseComp">
102 <param name="baseComp" type="text" value="" label="Genome and step" help="Parameter value is of the form genome,step. Ex: human,1000">
103 <expand macro="sanitizer_default"/>
104 <validator type="regex" message="Please enter as genome,step without leading or trailing spaces">^[^\s,](?:[^,]*[^\s,])?,[0-9]+$</validator>
105 </param>
106 </when>
107 <when value="--bedSequences">
108 <param name="bedSequences" type="text" value="" label="Genome name">
109 <expand macro="sanitizer_default"/>
110 <expand macro="validator_trim"/>
111 </param>
112 </when>
113 <when value="--topSegments">
114 <param name="topSegments" type="text" value="" label="Genome name">
115 <expand macro="sanitizer_default"/>
116 <expand macro="validator_trim"/>
117 </param>
118 </when>
119 <when value="--bottomSegments">
120 <param name="bottomSegments" type="text" value="" label="Genome name">
121 <expand macro="sanitizer_default"/>
122 <expand macro="validator_trim"/>
123 </param>
124 </when>
125 <when value="--branchLength">
126 <param name="branchLength" type="text" value="" label="Genome name">
127 <expand macro="sanitizer_default"/>
128 <expand macro="validator_trim"/>
129 </param>
130 </when>
131 <when value="--branches"/>
132 <when value="--children">
133 <param name="children" type="text" value="" label="Genome name">
134 <expand macro="sanitizer_default"/>
135 <expand macro="validator_trim"/>
136 </param>
137 </when>
138 <when value="--chromSizes">
139 <param name="chromSizes" type="text" value="" label="Genome name">
140 <expand macro="sanitizer_default"/>
141 <expand macro="validator_trim"/>
142 </param>
143 </when>
144 <when value="--coverage">
145 <param name="coverage" type="text" value="" label="Genome name">
146 <expand macro="sanitizer_default"/>
147 <expand macro="validator_trim"/>
148 </param>
149 </when>
150 <when value="--genomeMetaData">
151 <param name="genomeMetaData" type="text" value="" label="Genome name">
152 <expand macro="sanitizer_default"/>
153 <expand macro="validator_trim"/>
154 </param>
155 </when>
156 <when value="--genomes"/>
157 <when value="--metaData"/>
158 <when value="--numSegments">
159 <param name="numSegments" type="text" value="" label="Genome name">
160 <expand macro="sanitizer_default"/>
161 <expand macro="validator_trim"/>
162 </param>
163 </when>
164 <when value="--parent">
165 <param name="parent" type="text" value="" label="Genome name">
166 <expand macro="sanitizer_default"/>
167 <expand macro="validator_trim"/>
168 </param>
169 </when>
170 <when value="--percentID">
171 <param name="percentID" type="text" value="" label="Genome name" help="Only non-duplicated and unambiguous sites are considered">
172 <expand macro="sanitizer_default"/>
173 <expand macro="validator_trim"/>
174 </param>
175 </when>
176 <when value="--root"/>
177 <when value="--sequenceStats">
178 <param name="sequenceStats" type="text" value="" label="Genome name">
179 <expand macro="sanitizer_default"/>
180 <expand macro="validator_trim"/>
181 </param>
182 </when>
183 <when value="--sequences">
184 <param name="sequences" type="text" value="" label="Genome name">
185 <expand macro="sanitizer_default"/>
186 <expand macro="validator_trim"/>
187 </param>
188 </when>
189 <when value="--span">
190 <param name="span" type="text" value="" label="List of genomes" help="Enter a comma-separated (no spaces) list of genomes">
191 <expand macro="sanitizer_default"/>
192 <expand macro="validator_comma_list"/>
193 </param>
194 </when>
195 <when value="--spanRoot">
196 <param name="spanRoot" type="text" value="" label="List of genomes" help="Enter a comma-separated (no spaces) list of genomes">
197 <expand macro="sanitizer_default"/>
198 <expand macro="validator_comma_list"/>
199 </param>
200 </when>
201 <when value="--tree"/>
202 </conditional>
203 </inputs>
204 <outputs>
205 <data name="out_file" format="txt" label="${tool.name} on ${on_string}: Stats">
206 <change_format>
207 <when input="mode.option" value="--numSegments" format="tabular"/>
208 <when input="mode.option" value="--allCoverage" format="tabular"/>
209 <when input="mode.option" value="--sequenceStats" format="tabular"/>
210 <when input="mode.option" value="--percentID" format="tabular"/>
211 <when input="mode.option" value="--coverage" format="tabular"/>
212 <when input="mode.option" value="--chromSizes" format="tabular"/>
213 <when input="mode.option" value="--baseComp" format="tabular"/>
214 <when input="mode.option" value="--metaData" format="tabular"/>
215 <when input="mode.option" value="--genomeMetaData" format="tabular"/>
216 <when input="mode.option" value="--bedSequences" format="bed"/>
217 <when input="mode.option" value="--topSegments" format="bed"/>
218 <when input="mode.option" value="--bottomSegments" format="bed"/>
219 </change_format>
220 </data>
221 </outputs>
222 <tests>
223 <test expect_num_outputs="1">
224 <param name="input_hal" value="halTest.hal"/>
225 <output name="out_file" ftype="txt">
226 <assert_contents>
227 <has_line line="(Genome_1:1,Genome_2:1,Genome_3:1)Genome_0;"/>
228 <has_line line="GenomeName, NumChildren, Length, NumSequences, NumTopSegments, NumBottomSegments"/>
229 <has_line line="Genome_0, 3, 1758, 1, 0, 8"/>
230 <has_n_lines n="10"/>
231 </assert_contents>
232 </output>
233 </test>
234 <test expect_num_outputs="1">
235 <param name="input_hal" value="halTest.hal"/>
236 <conditional name="mode">
237 <param name="option" value="--allCoverage"/>
238 </conditional>
239 <output name="out_file" ftype="tabular">
240 <assert_contents>
241 <has_line line="FromGenome&#009; ToGenome&#009; sitesCovered1Times&#009; sitesCovered2Times&#009; sitesCovered3Times&#009; sitesCovered4Times&#009; sitesCovered5Times"/>
242 <has_line line="Genome_1&#009; Genome_1&#009; 5472&#009; 4688&#009; 3516&#009; 2637&#009; 1465"/>
243 <has_n_lines n="10"/>
244 </assert_contents>
245 </output>
246 </test>
247 <test expect_num_outputs="1">
248 <param name="input_hal" value="halTest.hal"/>
249 <conditional name="mode">
250 <param name="option" value="--bedSequences"/>
251 <param name="bedSequences" value="Genome_0"/>
252 </conditional>
253 <output name="out_file" ftype="bed">
254 <assert_contents>
255 <has_line line="Genome_0_seq&#009;0&#009;1758"/>
256 <has_n_lines n="1"/>
257 </assert_contents>
258 </output>
259 </test>
260 <test expect_num_outputs="1">
261 <param name="input_hal" value="halTest.hal"/>
262 <conditional name="mode">
263 <param name="option" value="--topSegments"/>
264 <param name="topSegments" value="Genome_1"/>
265 </conditional>
266 <output name="out_file" ftype="bed">
267 <assert_contents>
268 <has_line line="Genome_1_seq&#009;0&#009;293"/>
269 <has_line line="Genome_1_seq&#009;3223&#009;3399"/>
270 <has_line line="Genome_1_seq&#009;5274&#009;5472"/>
271 <has_n_lines n="28"/>
272 </assert_contents>
273 </output>
274 </test>
275 <test expect_num_outputs="1">
276 <param name="input_hal" value="halTest.hal"/>
277 <conditional name="mode">
278 <param name="option" value="--bottomSegments"/>
279 <param name="bottomSegments" value="Genome_0"/>
280 </conditional>
281 <output name="out_file" ftype="bed">
282 <assert_contents>
283 <has_line line="Genome_0_seq&#009;0&#009;293"/>
284 <has_line line="Genome_0_seq&#009;1033&#009;1172"/>
285 <has_line line="Genome_0_seq&#009;1465&#009;1758"/>
286 <has_n_lines n="8"/>
287 </assert_contents>
288 </output>
289 </test>
290 <test expect_num_outputs="1">
291 <param name="input_hal" value="halTest.hal"/>
292 <conditional name="mode">
293 <param name="option" value="--tree"/>
294 </conditional>
295 <output name="out_file" ftype="txt">
296 <assert_contents>
297 <has_line line="(Genome_1:1,Genome_2:1,Genome_3:1)Genome_0;"/>
298 <has_n_lines n="1"/>
299 </assert_contents>
300 </output>
301 </test>
302 <test expect_num_outputs="1">
303 <param name="input_hal" value="halTest.hal"/>
304 <conditional name="mode">
305 <param name="option" value="--spanRoot"/>
306 <param name="spanRoot" value="Genome_0,Genome_1"/>
307 </conditional>
308 <output name="out_file" ftype="txt">
309 <assert_contents>
310 <has_text text="Genome_0 Genome_1"/>
311 <has_n_lines n="1"/>
312 </assert_contents>
313 </output>
314 </test>
315 <test expect_num_outputs="1">
316 <param name="input_hal" value="halTest.hal"/>
317 <conditional name="mode">
318 <param name="option" value="--sequences"/>
319 <param name="sequences" value="Genome_0"/>
320 </conditional>
321 <output name="out_file" ftype="txt">
322 <assert_contents>
323 <has_line line="Genome_0_seq"/>
324 <has_n_lines n="1"/>
325 </assert_contents>
326 </output>
327 </test>
328 <test expect_num_outputs="1">
329 <param name="input_hal" value="halTest.hal"/>
330 <conditional name="mode">
331 <param name="option" value="--sequenceStats"/>
332 <param name="sequenceStats" value="Genome_0"/>
333 </conditional>
334 <output name="out_file" ftype="tabular">
335 <assert_contents>
336 <has_line line="SequenceName&#009; Length&#009; NumTopSegments&#009; NumBottomSegments"/>
337 <has_line line="Genome_0_seq&#009; 1758&#009; 0&#009; 8"/>
338 <has_n_lines n="3"/>
339 </assert_contents>
340 </output>
341 </test>
342 <test expect_num_outputs="1">
343 <param name="input_hal" value="halTest.hal"/>
344 <conditional name="mode">
345 <param name="option" value="--root"/>
346 </conditional>
347 <output name="out_file" ftype="txt">
348 <assert_contents>
349 <has_line line="Genome_0"/>
350 <has_n_lines n="1"/>
351 </assert_contents>
352 </output>
353 </test>
354 <test expect_num_outputs="1">
355 <param name="input_hal" value="halTest.hal"/>
356 <conditional name="mode">
357 <param name="option" value="--parent"/>
358 <param name="parent" value="Genome_1"/>
359 </conditional>
360 <output name="out_file" ftype="txt">
361 <assert_contents>
362 <has_line line="Genome_0"/>
363 <has_n_lines n="1"/>
364 </assert_contents>
365 </output>
366 </test>
367 <test expect_num_outputs="1">
368 <param name="input_hal" value="halTest.hal"/>
369 <conditional name="mode">
370 <param name="option" value="--percentID"/>
371 <param name="percentID" value="Genome_0"/>
372 </conditional>
373 <output name="out_file" ftype="tabular">
374 <assert_contents>
375 <has_line line="Genome&#009; % ID&#009; numID&#009; numSites"/>
376 <has_line line="Genome_0&#009; 1&#009; 1758&#009; 1758"/>
377 <has_n_lines n="5"/>
378 </assert_contents>
379 </output>
380 </test>
381 <test expect_num_outputs="1">
382 <param name="input_hal" value="halTest.hal"/>
383 <conditional name="mode">
384 <param name="option" value="--numSegments"/>
385 <param name="numSegments" value="Genome_1"/>
386 </conditional>
387 <output name="out_file" ftype="tabular">
388 <assert_contents>
389 <has_line line="numTopSegments&#009;numBottomSegments"/>
390 <has_line line="28&#009;0"/>
391 <has_n_lines n="2"/>
392 </assert_contents>
393 </output>
394 </test>
395 <test expect_num_outputs="1">
396 <param name="input_hal" value="halTest.hal"/>
397 <conditional name="mode">
398 <param name="option" value="--genomes"/>
399 </conditional>
400 <output name="out_file" ftype="txt">
401 <assert_contents>
402 <has_line line="Genome_0 Genome_1 Genome_2 Genome_3"/>
403 <has_n_lines n="1"/>
404 </assert_contents>
405 </output>
406 </test>
407 <test expect_num_outputs="1">
408 <param name="input_hal" value="halTest.hal"/>
409 <conditional name="mode">
410 <param name="option" value="--children"/>
411 <param name="children" value="Genome_0"/>
412 </conditional>
413 <output name="out_file" ftype="txt">
414 <assert_contents>
415 <has_line line="Genome_1 Genome_2 Genome_3"/>
416 <has_n_lines n="1"/>
417 </assert_contents>
418 </output>
419 </test>
420 <test expect_num_outputs="1">
421 <param name="input_hal" value="halTest.hal"/>
422 <conditional name="mode">
423 <param name="option" value="--chromSizes"/>
424 <param name="chromSizes" value="Genome_1"/>
425 </conditional>
426 <output name="out_file" ftype="tabular">
427 <assert_contents>
428 <has_line line="Genome_1_seq&#009;5472"/>
429 <has_n_lines n="1"/>
430 </assert_contents>
431 </output>
432 </test>
433 <test expect_num_outputs="1">
434 <param name="input_hal" value="halTest.hal"/>
435 <conditional name="mode">
436 <param name="option" value="--branches"/>
437 </conditional>
438 <output name="out_file" ftype="txt">
439 <assert_contents>
440 <has_line line="Genome_1 Genome_2 Genome_3"/>
441 <has_n_lines n="1"/>
442 </assert_contents>
443 </output>
444 </test>
445 <test expect_num_outputs="1">
446 <param name="input_hal" value="halTest.hal"/>
447 <conditional name="mode">
448 <param name="option" value="--branchLength"/>
449 <param name="branchLength" value="Genome_1"/>
450 </conditional>
451 <output name="out_file" ftype="txt">
452 <assert_contents>
453 <has_line line="1"/>
454 <has_n_lines n="1"/>
455 </assert_contents>
456 </output>
457 </test>
458 <test expect_num_outputs="1">
459 <param name="input_hal" value="halTest.hal"/>
460 <conditional name="mode">
461 <param name="option" value="--baseComp"/>
462 <param name="baseComp" value="Genome_0,1000"/>
463 </conditional>
464 <output name="out_file" ftype="tabular">
465 <assert_contents>
466 <has_line line="fraction_of_As&#009;fraction_of_Gs&#009;fraction_of_Cs&#009;fraction_of_Ts"/>
467 <has_line line="0.5&#009;0&#009;0.5&#009;0"/>
468 <has_n_lines n="2"/>
469 </assert_contents>
470 </output>
471 </test>
472 </tests>
473 <help><![CDATA[
474 halStats prints structural and summary information from a HAL file, which must be provided as input.
475 It can list genomes, sequences, sizes, and relationships, provides sequence level statistics such as coverage histograms and percent ID, and can export sequence or segment information in BED format.
476
477 It is useful for quick inspection of a HAL file and for extracting per genome or per sequence summaries.
478
479 -----
480
481 **Output**
482
483 The tool generates different output formats based on the selected type of statistic:
484
485 - **Tabular** for --coverage, --allCoverage, --sequenceStats, --percentID, --baseComp, --chromSizes, --metaData, --numSegments, or --genomeMetaData
486 - **BED** for --bedSequences, --topSegments, or --bottomSegments
487 - Plain **text** for all other type of statistics
488
489 ]]></help>
490 <expand macro="citation"/>
491 <expand macro="creator"/>
492 </tool>