hifiasm icon indicating copy to clipboard operation
hifiasm copied to clipboard

How to classify hifi reads into haplotype genomes? minimap2 did not work well

Open DayTimeMouse opened this issue 8 months ago • 3 comments

Hi,

I used HiFi + ONT + HiC reads to assemble haplotype genomes, then I want to classfiy hifi reads into haplotype genome separately.

I have tried to use minimap2, due to the similar haplotype genomes, it did not work well.

I cat hap1 and hap2 into hap1_hap2.fa, use "minimap2(2.27-r1193) -c --secondary=no hap1_hap2.fa hifi.fastq.gz >aln.paf", then use reads ID from paf file to calculate the precision, only ~0.56. The HiFi reads were simulated, so I know the reads ID real belong to which haplotype.

So, can hifiasm do this? Or, do you have other method for this?

Best wishes.

DayTimeMouse avatar Jun 13 '24 06:06 DayTimeMouse