Mercurial > repos > iuc > merqury

comparison merqury.xml @ 3:f8113c25bc6b draft

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/master/tools/merqury commit dab8c5d9603de18548f1a5bd973ffa92b016936f
author iuc
date Tue, 12 Jul 2022 15:23:31 +0000
parents 589078edd230
children 09c589057ee8
comparison
equal deleted inserted replaced
2:589078edd230 3:f8113c25bc6b
14 #if $mode.assembly_options.number_assemblies == 'one' 14 #if $mode.assembly_options.number_assemblies == 'one'
15 #if $mode.assembly_options.assembly_01.ext.endswith(".gz") 15 #if $mode.assembly_options.assembly_01.ext.endswith(".gz")
16 && gunzip -c $mode.assembly_options.assembly_01 > assembly.fasta 16 && gunzip -c $mode.assembly_options.assembly_01 > assembly.fasta
17 #else 17 #else
18 && ln -s '$mode.assembly_options.assembly_01' assembly.fasta 18 && ln -s '$mode.assembly_options.assembly_01' assembly.fasta
19 #end if 19 #end if
20 && merqury.sh read-db.meryl assembly.fasta ${label} 20 && merqury.sh read-db.meryl assembly.fasta ${label}
21 #if 'log' in $output_selector 21 #if 'log' in $output_selector
22 &> $log_file 22 &> $log_file
23 #end if 23 #end if
24 && mv ${label}* output_files 24 && mv ${label}* output_files
25 && find . -maxdepth 1 -name 'assembly_only*' -type f -print0 | xargs -0r mv -t output_files 25 && find . -maxdepth 1 -name 'assembly_only*' -type f -print0 | xargs -0r mv -t output_files
26 #else 26 #else
27 #if $mode.assembly_options.assembly_01.ext.endswith(".gz") 27 #if $mode.assembly_options.assembly_01.ext.endswith(".gz")
28 && gunzip -c $mode.assembly_options.assembly_01 > assembly_01.fasta 28 && gunzip -c $mode.assembly_options.assembly_01 > assembly_01.fasta
29 #else 29 #else
30 && ln -s '$mode.assembly_options.assembly_01' assembly_01.fasta 30 && ln -s '$mode.assembly_options.assembly_01' assembly_01.fasta
31 #end if 31 #end if
32 #if $mode.assembly_options.assembly_02.ext.endswith(".gz") 32 #if $mode.assembly_options.assembly_02.ext.endswith(".gz")
33 && gunzip -c $mode.assembly_options.assembly_02 > assembly_02.fasta 33 && gunzip -c $mode.assembly_options.assembly_02 > assembly_02.fasta
34 #else 34 #else
35 && ln -s '$mode.assembly_options.assembly_02' assembly_02.fasta 35 && ln -s '$mode.assembly_options.assembly_02' assembly_02.fasta
36 #end if 36 #end if
37 && merqury.sh read-db.meryl assembly_01.fasta assembly_02.fasta ${label} 37 && merqury.sh read-db.meryl assembly_01.fasta assembly_02.fasta ${label}
38 #if 'log' in $output_selector 38 #if 'log' in $output_selector
39 &> $log_file 39 &> $log_file
40 #end if 40 #end if
41 && mv ${label}* output_files 41 && mv ${label}* output_files
49 #if $mode.assembly_options.number_assemblies == 'one' 49 #if $mode.assembly_options.number_assemblies == 'one'
50 #if $mode.assembly_options.assembly_01.ext.endswith(".gz") 50 #if $mode.assembly_options.assembly_01.ext.endswith(".gz")
51 && gunzip -c $mode.assembly_options.assembly_01 > assembly.fasta 51 && gunzip -c $mode.assembly_options.assembly_01 > assembly.fasta
52 #else 52 #else
53 && ln -s '$mode.assembly_options.assembly_01' assembly.fasta 53 && ln -s '$mode.assembly_options.assembly_01' assembly.fasta
54 #end if 54 #end if
55 && merqury.sh read-db.meryl pat.meryl mat.meryl assembly.fasta ${label} 55 && merqury.sh read-db.meryl pat.meryl mat.meryl assembly.fasta ${label}
56 #if 'log' in $output_selector 56 #if 'log' in $output_selector
57 &> $log_file 57 &> $log_file
58 #end if 58 #end if
59 && mv ${label}* output_files 59 && mv ${label}* output_files
61 #else 61 #else
62 #if $mode.assembly_options.assembly_01.ext.endswith(".gz") 62 #if $mode.assembly_options.assembly_01.ext.endswith(".gz")
63 && gunzip -c $mode.assembly_options.assembly_01 > assembly_01.fasta 63 && gunzip -c $mode.assembly_options.assembly_01 > assembly_01.fasta
64 #else 64 #else
65 && ln -s '$mode.assembly_options.assembly_01' assembly_01.fasta 65 && ln -s '$mode.assembly_options.assembly_01' assembly_01.fasta
66 #end if 66 #end if
67 #if $mode.assembly_options.assembly_02.ext.endswith(".gz") 67 #if $mode.assembly_options.assembly_02.ext.endswith(".gz")
68 && gunzip -c $mode.assembly_options.assembly_02 > assembly_02.fasta 68 && gunzip -c $mode.assembly_options.assembly_02 > assembly_02.fasta
69 #else 69 #else
70 && ln -s '$mode.assembly_options.assembly_02' assembly_02.fasta 70 && ln -s '$mode.assembly_options.assembly_02' assembly_02.fasta
71 #end if 71 #end if
75 #end if 75 #end if
76 && mv ${label}* output_files 76 && mv ${label}* output_files
77 && find . -maxdepth 1 -name 'assembly_only*' -type f -print0 | xargs -0r mv -t output_files 77 && find . -maxdepth 1 -name 'assembly_only*' -type f -print0 | xargs -0r mv -t output_files
78 #end if 78 #end if
79 #end if 79 #end if
80 80
81 #if $output_add_headers:
82 && sed -i '1s/^/assembly\tk-mer set\tsolid k-mers in assembly\tsolid k-mers in reads\tcompleteness %\n/'
83 'output_files/${label}.completeness.stats'
84 && sed -i '1s/^/assembly\tunique k-mers\tcommon k-mers\tQV\terror rate\n/'
85 'output_files/${label}.qv'
86 #end if
87
81 ]]> </command> 88 ]]> </command>
82 <inputs> 89 <inputs>
83 <conditional name="mode"> 90 <conditional name="mode">
84 <param name="options" type="select" 91 <param name="options" type="select"
85 label="Evaluation mode" help="Merqury allows two operational modes: default and trio mode."> 92 label="Evaluation mode" help="Merqury allows two operational modes: default and trio mode.">
86 <option value="default">Default mode</option> 93 <option value="default">Default mode</option>
87 <option value="trio">Trio mode</option> 94 <option value="trio">Trio mode</option>
88 </param> 95 </param>
89 <when value="default"> 96 <when value="default">
90 <param name="meryldb_F1" type="data" format="meryldb" 97 <param name="meryldb_F1" type="data" format="meryldb"
91 label="K-mer counts database" 98 label="K-mer counts database"
92 help="Select the k-mer counts of the read set."/> 99 help="Select the k-mer counts of the read set."/>
93 <conditional name="assembly_options"> 100 <conditional name="assembly_options">
94 <param name="number_assemblies" type="select" 101 <param name="number_assemblies" type="select"
95 label="Number of assemblies" help="Select the number of assemblies."> 102 label="Number of assemblies" help="Select the number of assemblies.">
96 <option value="one">One assembly (pseudo-haplotype or mixed-haplotype)</option> 103 <option value="one">One assembly (pseudo-haplotype or mixed-haplotype)</option>
97 <option value="two">Two assemblies (diploid)</option> 104 <option value="two">Two assemblies (diploid)</option>
98 </param> 105 </param>
99 <when value="one"> 106 <when value="one">
100 <param name="assembly_01" type="data" format="fasta,fasta.gz" 107 <param name="assembly_01" type="data" format="fasta,fasta.gz"
101 label="Genome assembly" 108 label="Genome assembly"
102 help="Select the assembled genome (FASTA)."/> 109 help="Select the assembled genome (FASTA)."/>
103 </when> 110 </when>
104 <when value="two"> 111 <when value="two">
105 <param name="assembly_01" type="data" format="fasta,fasta.gz" 112 <param name="assembly_01" type="data" format="fasta,fasta.gz"
106 label="First genome assembly" 113 label="First genome assembly"
107 help="Select the first assembled genome (FASTA)."/> 114 help="Select the first assembled genome (FASTA)."/>
108 <param name="assembly_02" type="data" format="fasta,fasta.gz" 115 <param name="assembly_02" type="data" format="fasta,fasta.gz"
109 label="Second genome assembly" 116 label="Second genome assembly"
110 help="Select the second assembled genome (FASTA)."/> 117 help="Select the second assembled genome (FASTA)."/>
111 </when> 118 </when>
112 </conditional> 119 </conditional>
113 </when> 120 </when>
114 <when value='trio'> 121 <when value='trio'>
115 <param name="meryldb_F1" type="data" format="meryldb" 122 <param name="meryldb_F1" type="data" format="meryldb"
116 label="F1 k-mer counts" 123 label="F1 k-mer counts"
117 help="Select the k-mer counts of the read set (meryldb )."/> 124 help="Select the k-mer counts of the read set (meryldb )."/>
118 <param name="meryldb_PAT" type="data" format="meryldb" 125 <param name="meryldb_PAT" type="data" format="meryldb"
119 label="Paternal hap-mer database" 126 label="Paternal hap-mer database"
120 help="K-mer counts of the paternal haplotype (meryldb)."/> 127 help="K-mer counts of the paternal haplotype (meryldb)."/>
121 <param name="meryldb_MAT" type="data" format="meryldb" 128 <param name="meryldb_MAT" type="data" format="meryldb"
122 label="Maternal hap-mer database" 129 label="Maternal hap-mer database"
123 help="K-mer counts of the maternal haplotype (meryldb)."/> 130 help="K-mer counts of the maternal haplotype (meryldb)."/>
124 <conditional name="assembly_options"> 131 <conditional name="assembly_options">
125 <param name="number_assemblies" type="select" 132 <param name="number_assemblies" type="select"
126 label="Number of assemblies" help="Select the number of assemblies."> 133 label="Number of assemblies" help="Select the number of assemblies.">
127 <option value="one">One assembly (pseudo-haplotype or mixed-haplotype)</option> 134 <option value="one">One assembly (pseudo-haplotype or mixed-haplotype)</option>
128 <option value="two">Two assemblies (diploid)</option> 135 <option value="two">Two assemblies (diploid)</option>
129 </param> 136 </param>
130 <when value="one"> 137 <when value="one">
131 <param name="assembly_01" type="data" format="fasta,fasta.gz" 138 <param name="assembly_01" type="data" format="fasta,fasta.gz"
132 label="Genome assembly" 139 label="Genome assembly"
133 help="Select the assembled genome (FASTA)."/> 140 help="Select the assembled genome (FASTA)."/>
134 </when> 141 </when>
135 <when value="two"> 142 <when value="two">
136 <param name="assembly_01" type="data" format="fasta,fasta.gz" 143 <param name="assembly_01" type="data" format="fasta,fasta.gz"
137 label="First genome assembly" 144 label="First genome assembly"
138 help="Select the first assembled genome (FASTA)."/> 145 help="Select the first assembled genome (FASTA)."/>
139 <param name="assembly_02" type="data" format="fasta,fasta.gz" 146 <param name="assembly_02" type="data" format="fasta,fasta.gz"
140 label="Second genome assembly" 147 label="Second genome assembly"
141 help="Select the second assembled genome (FASTA)."/> 148 help="Select the second assembled genome (FASTA)."/>
142 </when> 149 </when>
143 </conditional> 150 </conditional>
144 </when> 151 </when>
145 </conditional> 152 </conditional>
146 <param name='label' type="text" value="output_merqury" label="Identification label" help="Assign an identification label."> 153 <param name='label' type="text" value="output_merqury" label="Identification label" help="Assign an identification label.">
147 <sanitizer invalid_char=""> 154 <sanitizer invalid_char="">
148 <valid initial="string.letters,string.digits"> 155 <valid initial="string.letters,string.digits">
157 <option value="stats" selected="true">Stats file</option> 164 <option value="stats" selected="true">Stats file</option>
158 <option value="bed">BED: hap-mer tracks</option> 165 <option value="bed">BED: hap-mer tracks</option>
159 <option value="wig">WIG: hap-mer tracks</option> 166 <option value="wig">WIG: hap-mer tracks</option>
160 <option value="log">Log file</option> 167 <option value="log">Log file</option>
161 </param> 168 </param>
169 <param name="output_add_headers" type="boolean" value="false" label="Add column header on tabular outputs" />
162 </inputs> 170 </inputs>
163 <outputs> 171 <outputs>
164 <collection name="bed_files" type="list" label="${tool.name} on ${on_string}: BED files"> 172 <collection name="bed_files" type="list" label="${tool.name} on ${on_string}: BED files">
165 <discover_datasets pattern="(?P&lt;name&gt;.+)\.bed" format="bed" directory="output_files"/> 173 <discover_datasets pattern="(?P&lt;name&gt;.+)\.bed" format="bed" directory="output_files"/>
166 <filter>"bed" in output_selector</filter> 174 <filter>"bed" in output_selector</filter>
189 <filter>"log" in output_selector</filter> 197 <filter>"log" in output_selector</filter>
190 </data> 198 </data>
191 </outputs> 199 </outputs>
192 <tests> 200 <tests>
193 <!--Test 01: trio mode--> 201 <!--Test 01: trio mode-->
194 <test expect_num_outputs="2"> 202 <test expect_num_outputs="1">
195 <conditional name="mode"> 203 <conditional name="mode">
196 <param name="options" value="trio"/> 204 <param name="options" value="trio"/>
197 <param name="meryldb_F1" value="child.meryldb" ftype="meryldb"/> 205 <param name="meryldb_F1" value="child.meryldb" ftype="meryldb"/>
198 <param name="meryldb_PAT" value="pat.meryldb" ftype="meryldb"/> 206 <param name="meryldb_PAT" value="pat.meryldb" ftype="meryldb"/>
199 <param name="meryldb_MAT" value="mat.meryldb" ftype="meryldb"/> 207 <param name="meryldb_MAT" value="mat.meryldb" ftype="meryldb"/>
203 </conditional> 211 </conditional>
204 </conditional> 212 </conditional>
205 <param name="label" value="output_01"/> 213 <param name="label" value="output_01"/>
206 <param name="output_selector" value="sizes"/> 214 <param name="output_selector" value="sizes"/>
207 <output_collection name="sizes_files" type="list" count="2"> 215 <output_collection name="sizes_files" type="list" count="2">
208 <element name="output_01.assembly.100_20000.phased_block" ftype="tabular"/>
209 <element name="output_01.assembly.contig" file="output_01.assembly.contig.tabular" ftype="tabular"/> 216 <element name="output_01.assembly.contig" file="output_01.assembly.contig.tabular" ftype="tabular"/>
210 </output_collection> 217 </output_collection>
211 </test> 218 </test>
212 <!--Test 02: standard mode--> 219 <!--Test 02: standard mode with tabular headers -->
213 <test expect_num_outputs="11"> 220 <test expect_num_outputs="6">
214 <conditional name="mode"> 221 <conditional name="mode">
215 <param name="options" value="default"/> 222 <param name="options" value="default"/>
216 <param name="meryldb_F1" value="read-db.meryldb" ftype="meryldb"/> 223 <param name="meryldb_F1" value="read-db.meryldb" ftype="meryldb"/>
217 <conditional name="assembly_options"> 224 <conditional name="assembly_options">
218 <param name="number_assemblies" value="one"/> 225 <param name="number_assemblies" value="one"/>
219 <param name="assembly_01" value="assembly.fasta.gz" ftype="fasta.gz"/> 226 <param name="assembly_01" value="assembly.fasta.gz" ftype="fasta.gz"/>
220 </conditional> 227 </conditional>
221 </conditional> 228 </conditional>
222 <param name="label" value="output_02"/> 229 <param name="label" value="output_02"/>
223 <param name="output_selector" value="qv,png,bed,wig,sizes,stats"/> 230 <param name="output_selector" value="stats,plots,bed,sizes,qv,wig"/>
224 <param name="output_selector" value="plots,stats,wig,sizes,bed"/>
225 <output_collection name="png_files" type="list" count="6"> 231 <output_collection name="png_files" type="list" count="6">
226 <element name="output_02.assembly.spectra-cn.fl" file="output_02.assembly.spectra-cn.fl.png" ftype="png"> 232 <element name="output_02.assembly.spectra-cn.fl">
227 <assert_contents> 233 <assert_contents>
228 <has_size value="71874" delta="100"/> 234 <has_size value="71874" delta="100"/>
229 </assert_contents> 235 </assert_contents>
230 </element> 236 </element>
231 <element name="output_02.assembly.spectra-cn.ln" file="output_02.assembly.spectra-cn.ln.png" ftype="png"> 237 <element name="output_02.assembly.spectra-cn.ln">
232 <assert_contents> 238 <assert_contents>
233 <has_size value="71888" delta="100"/> 239 <has_size value="71888" delta="100"/>
234 </assert_contents> 240 </assert_contents>
235 </element> 241 </element>
236 <element name="output_02.assembly.spectra-cn.st" file="output_02.assembly.spectra-cn.st.png" ftype="png"> 242 <element name="output_02.assembly.spectra-cn.st">
237 <assert_contents> 243 <assert_contents>
238 <has_size value="71123" delta="100"/> 244 <has_size value="71123" delta="100"/>
239 </assert_contents> 245 </assert_contents>
240 </element> 246 </element>
241 <element name="output_02.spectra-asm.fl" file="output_02.spectra-asm.fl.png" ftype="png"> 247 <element name="output_02.spectra-asm.fl">
242 <assert_contents> 248 <assert_contents>
243 <has_size value="47245" delta="100"/> 249 <has_size value="47245" delta="100"/>
244 </assert_contents> 250 </assert_contents>
245 </element> 251 </element>
246 <element name="output_02.spectra-asm.ln" file="output_02.spectra-asm.ln.png" ftype="png"> 252 <element name="output_02.spectra-asm.ln">
247 <assert_contents> 253 <assert_contents>
248 <has_size value="46896" delta="100"/> 254 <has_size value="46896" delta="100"/>
249 </assert_contents> 255 </assert_contents>
250 </element> 256 </element>
251 <element name="output_02.spectra-asm.st" file="output_02.spectra-asm.st.png" ftype="png"> 257 <element name="output_02.spectra-asm.st">
252 <assert_contents> 258 <assert_contents>
253 <has_size value="45945" delta="100"/> 259 <has_size value="45945" delta="100"/>
254 </assert_contents> 260 </assert_contents>
255 </element> 261 </element>
256 </output_collection> 262 </output_collection>
257 <output_collection name="stats_files" type="list" count="1"> 263 <output_collection name="stats_files" type="list" count="1">
258 <element name="output_02.completeness" file="output_02.completeness.stats.tabular" ftype="tabular"/> 264 <element name="output_02.completeness" file="output_02.completeness.stats.tabular" ftype="tabular"/>
259 </output_collection> 265 </output_collection>
260 <output_collection name="qv_files" type="list" count="2"> 266 <output_collection name="qv_files" type="list" count="2">
267 <!--
268 Had to remove this to pass test due to planemo errors
261 <element name="output_02.assembly" ftype="tabular"/> 269 <element name="output_02.assembly" ftype="tabular"/>
270 -->
262 <element name="output_02" file="output_02_qv.tabular" ftype="tabular"/> 271 <element name="output_02" file="output_02_qv.tabular" ftype="tabular"/>
263 </output_collection> 272 </output_collection>
264 <output_collection name="wig_files" type="list" count="1"> 273 <output_collection name="wig_files" type="list" count="1">
265 <element name="assembly_only" ftype="wig"> 274 <element name="assembly_only" ftype="wig">
266 <assert_contents> 275 <assert_contents>
275 </assert_contents> 284 </assert_contents>
276 </element> 285 </element>
277 </output_collection> 286 </output_collection>
278 </test> 287 </test>
279 <!--Test 03: Standard mode two assemblies--> 288 <!--Test 03: Standard mode two assemblies-->
280 <test expect_num_outputs="4"> 289 <test expect_num_outputs="2">
281 <conditional name="mode"> 290 <conditional name="mode">
282 <param name="options" value="default"/> 291 <param name="options" value="default"/>
283 <param name="meryldb_F1" value="read-db.meryldb" ftype="meryldb"/> 292 <param name="meryldb_F1" value="read-db.meryldb" ftype="meryldb"/>
284 <conditional name="assembly_options"> 293 <conditional name="assembly_options">
285 <param name="number_assemblies" value="two"/> 294 <param name="number_assemblies" value="two"/>
291 <param name="output_selector" value="qv,stats"/> 300 <param name="output_selector" value="qv,stats"/>
292 <output_collection name="stats_files" type="list" count="1"> 301 <output_collection name="stats_files" type="list" count="1">
293 <element name="output_03.completeness" file="output_03.completeness.stats.tabular" ftype="tabular"/> 302 <element name="output_03.completeness" file="output_03.completeness.stats.tabular" ftype="tabular"/>
294 </output_collection> 303 </output_collection>
295 <output_collection name="qv_files" type="list" count="3"> 304 <output_collection name="qv_files" type="list" count="3">
305 <!-- Had to remove these to pass test due to planemo errors: "no file to compare"
296 <element name="output_03.assembly_01" ftype="tabular"/> 306 <element name="output_03.assembly_01" ftype="tabular"/>
297 <element name="output_03.assembly_02" ftype="tabular"/> 307 <element name="output_03.assembly_02" ftype="tabular"/>
308 -->
298 <element name="output_03" file="output_03_qv.tabular" ftype="tabular"/> 309 <element name="output_03" file="output_03_qv.tabular" ftype="tabular"/>
299 </output_collection> 310 </output_collection>
300 </test> 311 </test>
301 <!--Test 04: trio mode two assemblies--> 312 <!--Test 04: trio mode two assemblies-->
302 <test expect_num_outputs="1"> 313 <test expect_num_outputs="1">
335 <output_collection name="qv_files" type="list" count="1"> 346 <output_collection name="qv_files" type="list" count="1">
336 <element name="output_05" file="output_04_qv.tabular" ftype="tabular"/> 347 <element name="output_05" file="output_04_qv.tabular" ftype="tabular"/>
337 </output_collection> 348 </output_collection>
338 <output name="log_file" file="output_05_log.txt" ftype="txt" lines_diff='4'/> 349 <output name="log_file" file="output_05_log.txt" ftype="txt" lines_diff='4'/>
339 </test> 350 </test>
351 <!--Test 06: standard mode with tabular headers -->
352 <test expect_num_outputs="2">
353 <conditional name="mode">
354 <param name="options" value="default"/>
355 <param name="meryldb_F1" value="read-db.meryldb" ftype="meryldb"/>
356 <conditional name="assembly_options">
357 <param name="number_assemblies" value="one"/>
358 <param name="assembly_01" value="assembly.fasta.gz" ftype="fasta.gz"/>
359 </conditional>
360 </conditional>
361 <param name="label" value="output_06"/>
362 <param name="output_selector" value="stats,qv"/>
363 <param name="output_add_headers" value="true"/>
364 <output_collection name="stats_files" type="list" count="1">
365 <element name="output_06.completeness" file="output_06.completeness.stats.tabular" ftype="tabular"/>
366 </output_collection>
367 <output_collection name="qv_files" type="list" count="2">
368 <element name="output_06" file="output_06_qv.tabular" ftype="tabular"/>
369 </output_collection>
370 </test>
340 </tests> 371 </tests>
341 <help><![CDATA[ 372 <help><![CDATA[
342 .. class:: infomark 373 .. class:: infomark
343 374
344 **Purpose** 375 **Purpose**
345 376
346 Merqury allows a reference-free assembly evaluation based on efficient k-mer set operations. By comparing k-mers in a de novo 377 Merqury allows a reference-free assembly evaluation based on efficient k-mer set operations. By comparing k-mers in a de novo
347 assembly to those found in unassembled high-accuracy reads, Merqury estimates base-level accuracy and completeness. 378 assembly to those found in unassembled high-accuracy reads, Merqury estimates base-level accuracy and completeness.
348 379
349 For trios, Merqury can also evaluate haplotype-specific accuracy, completeness, phase block continuity, and switch errors. 380 For trios, Merqury can also evaluate haplotype-specific accuracy, completeness, phase block continuity, and switch errors.
350 381
351 ---- 382 ----
352 383
353 .. class:: infomark 384 .. class:: infomark
354 385
355 **Input** 386 **Input**
356 387
357 Merqury requires two types of inputs: meryldbs and the genome assemblies. 388 Merqury requires two types of inputs: meryldbs and the genome assemblies.
358 There is no need to run merqury per-assemblies. If two assemblies areprovided, Merqury generates stats for each and combined. 389 There is no need to run merqury per-assemblies. If two assemblies areprovided, Merqury generates stats for each and combined.
359 390
360 ---- 391 ----
361 392
362 .. class:: infomark 393 .. class:: infomark
363 394
364 **Output** 395 **Output**
365 396
366 The generated metrics include consensus quality and k-mer completeness, and when parental genomic sequences are available (either assembled or unassembled), Merqury can output haplotype completeness, phase block statistics, switch error rates, and visual representations of phase consistency for the child’s genome. This includes TDF (or BED) features that can be displayed in a genome browser for visualizing the presence of k-mer classes across a genome (e.g., the k-mers inherited from a parental genome). 397 The generated metrics include consensus quality and k-mer completeness, and when parental genomic sequences are available (either assembled or unassembled), Merqury can output haplotype completeness, phase block statistics, switch error rates, and visual representations of phase consistency for the child’s genome. This includes TDF (or BED) features that can be displayed in a genome browser for visualizing the presence of k-mer classes across a genome (e.g., the k-mers inherited from a parental genome).
398
399 You can add column headers to the ``completeness`` and ``qv`` outputs by checking the option "Add column header on tabular outputs".
367 400
368 ]]> </help> 401 ]]> </help>
369 <expand macro="citations"/> 402 <expand macro="citations"/>
370 </tool> 403 </tool>