Question: HGAP assembly algorithm(help)
gravatar for midox
2.6 years ago by
midox180 wrote:

In the assembler HGAP, they used long and short reads or reads only long?

because they wrote in their article "we developed a consensus algorithm that preassembles long and highly accurate overlapping sequences by correcting errors on the longest reads using shorter reads from the same library."

"using go short reads from the Saami library" long or short reads reads from the Saami library?
thank you for your reply

ADD COMMENTlink modified 2.6 years ago • written 2.6 years ago by midox180
gravatar for Medhat
2.6 years ago by
Medhat6.9k wrote:

only long reads 

We present a hierarchical genome-assembly process (HGAP) for high-quality de novo microbial genome assemblies using only a single, long-insert shotgun DNA library in conjunction with Single Molecule, Real-Time (SMRT) DNA sequencing. Our method uses the longest reads as seeds to recruit all other reads for construction of highly accurate preassembled reads through a directed acyclic graph–based consensus procedure, which we follow with assembly using off-the-shelf long-read assemblers. In contrast to hybrid approaches, HGAP does not require highly accurate raw reads for error correction. 


The work flow as follow 

A workflow to first preassemble reads, assemble the preassembled reads using Celera® Assembler, then polish using Quiver


The longest reads are selected as 'seed' reads, to which all other reads are mapped.

as you can see in the figure the longest read in blue first used to map  short once to it that is called the preassembles then they used this preassembled to form assembly  as you can see in the second part of the figure , but in your consideration that The best de novo assembly results are generated by HGAP if at least 60X-100X coverage generated 



I hope that make it clearer  


ADD COMMENTlink modified 2.6 years ago • written 2.6 years ago by Medhat6.9k

thank you for your reply.
 can you explain the term "seed reads"?
 thank you

ADD REPLYlink written 2.6 years ago by midox180
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1372 users visited in the last hour