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

Higher QV in a_ctg than p_ctg #57

Open
pcarbone opened this issue Dec 2, 2020 · 1 comment
Open

Higher QV in a_ctg than p_ctg #57

pcarbone opened this issue Dec 2, 2020 · 1 comment

Comments

@pcarbone
Copy link

pcarbone commented Dec 2, 2020

Thank you for the nice and very efficient assembler.

I assembled a highly heterozygous genome and hifiasm did very well separating primary and haplotigs (p_ctg size is only a bit larger than the expected genome size and a_ctg length is 90% of p_ctg size).

I ran Merqury for QC from short-read kmers of the same sample. I always have higher consensus accuracy QV in a_ctg (QV=~49) than p_ctg (QV=~37) independently of the parameters I use in hifiasm.
Is it expected?
Are there some hifiasm options that could improve QV of p_ctg?

Thanks

@pcarbone pcarbone changed the title Higher QV in p_ctg than in a_ctg Higher QV in a_ctg than p_ctg Dec 2, 2020
@chhylp123
Copy link
Owner

chhylp123 commented Dec 2, 2020

The QV of p_ctg should be roughly as same as that of a_ctg, and QV=~37 is too low for HiFi assembly. Could you please evaluate QV with yak qv (https://github.com/lh3/yak)? I'm not sure if it is the problem of hifiasm or Merqury.

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

2 participants