ALLHiC icon indicating copy to clipboard operation
ALLHiC copied to clipboard

AllHiC on Hifiasm assembly

Open sandipmkale opened this issue 3 years ago • 2 comments

Hello,

Thanks for great tool. I have question regarding development of phased assembly. I used hifiasm to generate haplotype level contain assemblies. As the phases are interchanged in the hifiasm assembly, can I merge the *hap1.p.ctg.fa and *hap2.p.ctg.fa and run the AllHiC to generate phased assemblies?

Thanks and regards

Sandip

sandipmkale avatar Nov 26 '21 13:11 sandipmkale

Yes, It should be totally OK to merge the hap1.p_ctg.fa and hap2.p_ctg.fa and run ALLHiC for haplotype phasing. Alternatively, you can also use p_utg.fa for ALLHiC phasing if you do not have trio data.

tangerzhang avatar Nov 30 '21 06:11 tangerzhang

Thanks alot for reply. The haploid genome size is 2.5 Gb, but the *p_utg.fa file is 4.5 Gb so around 500 Mb data is merged in the *p_utg.fa file. I run the AllHiC but didn't get satisfactory results. I calculated coverage for each units using the pipeline mentioned in the tea genome (Nature Genetics) paper but the maximum depth in *read_depth.txt file is 1. so decided to merge hap1.p_ctg.fa and hap2.p_ctg.fa.

sandipmkale avatar Nov 30 '21 08:11 sandipmkale