view test-data/output_05_log.txt @ 4:09c589057ee8 draft

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/merqury commit 2c0e1d33bb5a462cd31bb4ad1ec4821520caccd4
author iuc
date Thu, 10 Aug 2023 22:57:52 +0000
parents f8113c25bc6b
children
line wrap: on
line source

read: read-db.meryl

Haplotype dbs provided.
Running Merqury in trio mode...

hap1: pat.meryl
hap2: mat.meryl
asm1: assembly_01.fasta
asm2: assembly_02.fasta
out : output_05


Get spectra-cn plots and QV stats

Get blob plots

Get haplotype specfic spectra-cn plots

Get phase blocks


Get block N plots
No modules available..
# Generate assembly_01.fasta.fai

*** # Found assembly_01.gaps.bed ***

# No gaps found. This is a contig set.

# Generate assembly_02.fasta.fai

*** # Found assembly_02.gaps.bed ***

# No gaps found. This is a contig set.

# Convert output_05.assembly_01.100_20000.phased_block.bed to sizes
 Result saved as output_05.assembly_01.100_20000.phased_block.sizes

# Plot output_05.assembly_01.100_20000.phased_block.bed
Rscript /usr/local/share/merqury/plot/plot_block_N.R -b output_05.assembly_01.100_20000.phased_block.sizes -c output_05.assembly_01.contig.sizes  -o output_05.assembly_01 
Loading required package: argparse
Loading required package: ggplot2
Loading required package: scales
Error in read.table(dat, header = F) : no lines available in input
Calls: block_n -> attach_n -> read.table
Execution halted

# Convert output_05.assembly_02.100_20000.phased_block.bed to sizes
 Result saved as output_05.assembly_02.100_20000.phased_block.sizes

# Plot output_05.assembly_02.100_20000.phased_block.bed
Rscript /usr/local/share/merqury/plot/plot_block_N.R -b output_05.assembly_02.100_20000.phased_block.sizes -c output_05.assembly_02.contig.sizes  -o output_05.assembly_02 
Loading required package: argparse
Loading required package: ggplot2
Loading required package: scales
Error in read.table(dat, header = F) : no lines available in input
Calls: block_n -> attach_n -> read.table
Execution halted