You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
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.
The text was updated successfully, but these errors were encountered:
I noticed that both hap1.gfa and hap2.gfa in the output of hifiasm contain HG:A:m and HG:A:p, why do they contain two kinds of tags? Another problem is that the haplotype tag in gfa only has some reads ids, so how do I get all the haplotype partition reads ids?
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.
The text was updated successfully, but these errors were encountered: