-
Notifications
You must be signed in to change notification settings - Fork 4
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
How to improve allotetraploid scaffolding in haphic? #22
Comments
Hi @RezwanCAAS, According to your heatmap, it seems that the homologous chromosomes were incorrectly clustered. How did you assemble the genome and which assembly did you use for scaffolding (e.g., p_utg, p_ctg, or hap*.p_ctg)? Best, |
Hi @zengxiaofei I used hifiasm with following command and got these outputs
output
later, I used yellow_assembly.hic.p_ctg.fasta file for scaffolding with haphic. Please guide some points to improve the scaffolding. |
Hi @RezwanCAAS, You need to concatenate the contigs in the Best, |
@zengxiaofei thank you so much for helping. I will let you know soon after getting the results. |
Hi @zengxiaofei following the above given suggestions. I have this output in form of 44 groups. The groups are shown here in the hic plot. So what do you suggest here? how can I improve it
|
Hi @RezwanCAAS, It seems that the heatmap is clear enough. You can manually adjust it in Juicebox after importing the Best, |
@zengxiaofei thank you for your great help. I will add the final figure here after correction with juicebox as reference for other users as well. |
@RezwanCAAS Thanks for your sharing! |
@zengxiaofei Please check this plot having 44 chromosomes. How does this look? |
Sorry for the delay. I'm quite busy these days. Your contact heatmap shows that there are still many errors in the contig assignment, as well as the ordering and orientation. The Hi-C signals you highlighted with red circles are signals between the homologous chromosomes. They mainly derive from assembly errors (but it's normal for haplotype-phased assemblies). You could check out our manuscript for more information about collapsed contigs, chimeric contigs, and switch errors. You could also have a look at the heatmaps we generated in the tests of real cases, especially those haplotype-phased assemblies. These figures are in the Supplementary Information. I believe it is helpful for you in curating your assembly. |
@zengxiaofei thank you so much for your great support, and shared shared examples will be helpful to improve my genome plot accordingly. |
Hi , I used haphic for allotetraploid genome (2n=4x=44). It makes 22 groups like the following which have huge variations in chromosome sizes as given below and can be seen in juicebox plot. Could you suggest how to improve this
group1 335147534
group2 265712237
group3 259566223
group4 256679759
group5 252275747
group6 223506148
group7 216891361
group8 204149503
group9 178470265
group10 169460411
group11 165539612
group12 151015244
group13 150193252
group14 124757396
group15 106023119
group16 103438627
group17 102100592
group18 94194491
group19 71377479
group20 51412018
group21 41958356
group22 41183465
The text was updated successfully, but these errors were encountered: