Recreate p_ctg from p_utg
diego-rt opened this issue · comments
Hello,
I'm trying to curate the assembly produced by hifiasm by reassembling the reads at the assembly tips with different parameters and assemblers. This would be bit more efficient if I could edit the assembly graph (manually add/remove or update nodes and edges), and recreate the contigs as hifiasm(UL) does. Do you by any chance have any code that you could share that could take the p_utg from the UL pipeline as input and output the primary p_ctg/hap1/hap2 graphs?
Thank you so much in advance!
@chhylp123 @lh3 Or alternatively, if hifiasm could output the order, orientation and coordinates of the unitigs used to generate each contig/haplotig, it would be super helpful for curation. In other words, similar to the 'A' lines in the current GFA but with unitigs rather than reads.
Thanks a lot in advance and sorry for the shameless feature request...
@diego-rt Sorry for the late reply since I was quite busy this month. I will add this very soon.
Amazing, thank you so much @chhylp123 !!