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

Resolving switching error (?) #650

Open
s-kyungyong opened this issue May 14, 2024 · 0 comments
Open

Resolving switching error (?) #650

s-kyungyong opened this issue May 14, 2024 · 0 comments

Comments

@s-kyungyong
Copy link

Hi!

I am assembling a hetrozygous diploid plant genome with 110X HiFi reads and Hi-C data. The hetrozygosity is about 1%, and the genome size is 1 Gb per haplotype. After scaffolding and comparison to the reference genome, I found some sudden changes in homology between two corresponding chromosomes. In the graph below, the x-axis is my genome, and the y-axis is the reference, the green lines indicate each contig. The red box indicates such a sudden shift in homology that happened within a single contig.

image

After mapping the HiFi reads and investigating the regions around such breaks, I always identify a region in which haplotypes don't seem to be fully resolved as below. What would be the best way to resolve this? Some scaffolds are already good as they are now. But some others have too many to fix manually.
image

I used v0.19.9-r616 with the following command:
hifiasm --telo-m GGGTAAA --dual-scaf -l3 -t 56 --h1 .. --h2 .. ...

I am playing with -u and --b-cov, but in my first attempt, N50 dropped from 15 Mb to 3.2 Mb. I am not sure if this would be helpful.

Thank you!

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

1 participant