combine hifiasm and hifiasm-ont assemblies.
1
0
Entering edit mode
6 weeks ago
Mark ▴ 60

I am wondering if it's possible to perform genome assemblies by combining PacBio HiFi and ONT data. How I understand it, contigs made by hifiasm from PacBio HiFi reads are more accurate but shorter, and contigs made by hifiasm-ont from ONT reads are more less accurate but more contiguous, I'm wondering if there's some tool that allows you to create de novo assemblies by creating a draft assembly from PacBio contigs, and then scaffolding them using ONT contigs. Does such a pipeline exist yet / is there a better, alternative strategy?

wgs ngs genomics sequencing • 1.1k views
ADD COMMENT
1
Entering edit mode

If you are working with a good model organism then you may be able to use available reference to try and merge the assemblies.

In case you have not seen this hifiasm author hints that a version that can use both PacBio and Nanopore data at the same time may be available in future: https://github.com/chhylp123/hifiasm/issues/849#issuecomment-3019737180

ADD REPLY
3
Entering edit mode
6 weeks ago
shelkmike ★ 1.8k

I suppose that you can just combine Nanopore and PacBio reads in one file and give this file to Hifiasm with the option "--ont". Probably, Hifiasm will treat PacBio reads just as a "more accurate fraction" of Nanopore reads.

Alternatively, you can make two separate assemblies and merge them with a tool like Quickmerge (https://github.com/mahulchak/quickmerge). However, in my experience, merging of assemblies in not a good thing, because if one assembly contains a structural misassembly, the merged assembly will have duplicated sequences. For example, if there is a genomic region A-B, and one of the assemblies contain these sequences in an improper order (B-A), then merging may produce A-B-A or B-A-B. In other words, some sequences will be duplicated. Such improper merges can be detected by an increase in BUSCO's "D".

ADD COMMENT

Login before adding your answer.

Traffic: 3635 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6