Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Hi-C mode kat spectra-cn validation fails #98

Open
chklopp opened this issue Apr 16, 2021 · 5 comments
Open

Hi-C mode kat spectra-cn validation fails #98

chklopp opened this issue Apr 16, 2021 · 5 comments

Comments

@chklopp
Copy link

chklopp commented Apr 16, 2021

I've assembled a 1.3 Gb genome with hifiasm (version 0.14 e6e6dbf), merged both haplotypes and compared the kmer content to the read kmers with kat spectra-cn.
Because of a small heterozygocity of 0.5% I was expecting to see the major gaussian (homozygous kmers, centered at 25X) in purple and the smaller gaussian to the left (heterozygous kmers, centered at 12X) in red.
This in not what the graph shows :
http://genoweb.toulouse.inra.fr/~klopp/tmp/hifiasm_0_14_Mrh_both_hap_new-main.mx.spectra-cn.png

Their are unexpected 3X and 4X kmer sets attesting the some parts of the genome are over represented. The central red curve is also showing that a large set of homozygous kmers are missing a copy in the merged assembly.

@lh3
Copy link
Collaborator

lh3 commented Apr 16, 2021

merged both haplotypes

Don't merge p_ctg and a_ctg except for dup purging.

@chklopp
Copy link
Author

chklopp commented Apr 16, 2021

I did not. I merged hap1 and hap2.

@lh3
Copy link
Collaborator

lh3 commented Apr 16, 2021

Do you have Hi-C or parental data? If Hi-C, you may wait for the next release.

@chklopp
Copy link
Author

chklopp commented Apr 16, 2021

It is Hi-C...I will wait...

@chhylp123
Copy link
Owner

chhylp123 commented Apr 17, 2021

Please have a try with the new version: https://github.com/chhylp123/hifiasm/releases/tag/0.15. Both Hi-C and non-Hi-C modules should be able to output two balanced haplotypes.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants