Fast and accurate chromosome-scale phased assembly about dipasm HOT 7 CLOSED

Dear Shilpagarg,

Thanks for the instructions. I am following the pstools commands.

Is there a way to estimate the amount of RAM pstools resolve_haplotypes will need?

My "hic_name_connection.output" is around 70G when raw HiC reads were provided. "pstools resolve_haplotypes" ran out of RAM with a 1 TB server.

Should the analysis starts with the set of HiC reads that proven to be "valid tags"? Or the validation of HiC tags were included in the workflow?

Many thanks in advance. Have a nice weekend,

Weihong

from dipasm.

<150 Gb for humans with ~30x Hi-C.

from dipasm.

Dear Shilpagarg,

Many thanks for your prompt response.

For the new analysis, only valid HiC tags after filtering were used for mapping, and the resulting "hic_name_connection.output" is only 7G. But "pstools resolve_haplotypes" still failed when 3 TB RAM was allocated.

The genome is highly heterozygous. Usually the primary assembly needs to be purged aggressively at the end. But each of the haploid genome size is less than 1 Gbp.

The stderr of "pstools resolve_haplotypes":
start main
start get bubbles
rm: missing operand
Try 'rm --help' for more information.
finish get bubble chain
utg...
... (here total 687 lines of utg ids)
start get bubbles
/var/lib/slurm/slurmd/job3498159/slurm_script: line 8: 1537222 Killed pstools resolve_haplotypes -t32 hic_name_connection.output s60.asm.r_utg.gfa ./

The tmp output "pstools.clean_graph.temp.0.out" is about 1.6 G.

If you have any suggestions on how to further proceed, they will be highly appreciated. Many thanks in advance.

Kind regards,

Weihong

from dipasm.

Thank you. If you could send me properties of the genome as well as datasets used at [email protected], I will be happy to help.

from dipasm.

@tcb72: Please download binaries and use them as mentioned in the issue. We have tested on several non-human genomes and it seems to work fine. I look forward to helping you with your experiments using pstools (You don't need to install docker etc.).

from dipasm.

Thanks for the information here, I just wondered would DipAsm work for a triploid species? What's its main advantage comparing with hifiasm?

from dipasm.

is the link https://pstools.s3.us-east-2.amazonaws.com/pstools_1 still working?

from dipasm.