View Single Post
Old 04-11-2018, 07:35 AM   #2
Senior Member
Location: San Francisco

Join Date: Aug 2012
Posts: 322

It's possible that ploidy could account for some of the difference, but it is unlikely that it is that simple. Using default parameters I wouldn't expect Canu to not collapse haplotypes and Falcon to fully collapse, although with the correct parameters and enough divergence between haplotypes both assemblers theoretically could.
Some things to try:
1. Align assemblies to one another to see what the difference actually is, for a quick visulization,
2. Look at the alternate contigs in the HGAP.4 assembly, <HGAP.4 job dir>/tasks/falcon_ns.tasks.task_falcon2_run_asm-0/a_ctg.fa. HGAP.4 is based on Falcon so does generate contigs for possible alternate haplotypes
Some more resources on diploid assembly:
rhall is offline   Reply With Quote