shilpagarg / pstools

MIT License
33 stars 6 forks source link

pred_haplotypes.fa includes only 2/3 of the genome assembly? #11

Open Alteroldis opened 10 months ago

Alteroldis commented 10 months ago

Dear @shilpagarg, I have 3GB genome assembly graph and fasta, but after I ran commands below, I got 2GB file pred_haplotypes.fa. What principles be the basis of filtration edges? I used version from github. Initial genome assembly was created with Flye from ONT reads. Genome of these species is highly heterozygous, repetitive and with enormous count of recombinations (?), because of that I can`t assemble mosaic haploid genome. Commands: pstools hic_mapping -t48 -o map-git.out assembly_graph.fa <(zcat $reads1) <(zcat $reads2) pstools resolve_haplotypes -t48 map-git.out assembly_graph.gfa git-woi

P.S. Command from README with "-i true" argument is crushed with SEGFault.