hifiasm icon indicating copy to clipboard operation
hifiasm copied to clipboard

Two hap are larger than survey and flow cytometry result

Open zhang144999 opened this issue 6 months ago • 1 comments

Dear Haoyu

I have a species that survey and flow cytometry are 2.4G.

This is PacBio data survey. WechatIMG2163

This is next generation sequence data survey. WechatIMG2152

But the hifiasm default parameter results, hap1 is 5.5G, hap2 is 3.4G, p.ctg is 4.7G, p.utg is 8.7G.

hap1 BUSCO is

C:97.2%[S:25.2%,D:72.0%],F:0.9%,M:1.9%,n:425	   
413	Complete BUSCOs (C)			   
107	Complete and single-copy BUSCOs (S)	   
306	Complete and duplicated BUSCOs (D)	   
4	Fragmented BUSCOs (F)			   
8	Missing BUSCOs (M)			   
425	Total BUSCO groups searched		

hap2 BUSCO is

C:87.0%[S:53.6%,D:33.4%],F:5.6%,M:7.4%,n:425	   
370	Complete BUSCOs (C)			   
228	Complete and single-copy BUSCOs (S)	   
142	Complete and duplicated BUSCOs (D)	   
24	Fragmented BUSCOs (F)			   
31	Missing BUSCOs (M)			   
425	Total BUSCO groups searched	

Here is log.

12716149.txt

Next, I add '--hom-cov 31' this parameter. I'm still not satisfied with the result. Here is new log 12723223.txt

hap1 is 4.8G, hap2 is 3.5G, p.ctg is 3.8G, p.utg is 8.4G. I haven't done a busco evaluation yet, but I don't think it will make a difference.

I'm confused. Can you give me some advice? I am looking forward to your response and will greatly appreciate any guidance or advice you provide in this situation.

Best Zhang

zhang144999 avatar Dec 16 '23 01:12 zhang144999