Hi, what is everyone's experience in automated gap closure?
I have found a few papers/software that show methods of automated gap closure using paired end reads. What is everyone's experience? I've read about IMAGE, CloG, and the Reconciliator. What other software is available to try out?
In my opinion, two really competitive stand-alone gap closers are:
Your question needs to be reframed though (it appears that you want to do scaffolding). Gaps, in the assembly context, refer to sequences of undetermined nucleotides (NNN...NNN) inside scaffolds.
If you only have contigs (which do not have gaps by the previous definition) and paired-end reads, what you want is to perform scaffolding (another type of software does that). Gap-closing is a later step which consists in replacing these undetermined regions by the actual sequence.
For scaffolding, there are many options, I would recommend SSPACE or a recent one that I haven't tested, SCARPA.
FinIS seems to require a license to a library by Mosek, so I am waiting on that.
GapCloser... are there any examples out there? It requires a scaffold file, but it seems like a catch-22. I would like to join contigs of a de novo assembly that do not have a scaffold yet.
Thank you for introducing these to me!
Yes that is right, you need to get a free academic license for Mosek.. Indeed, that isn't really user friendly. I just know that FinIS has a really nice algorithm.
Concerning GapCloser, the web page (http://soap.genomics.org.cn/about.html) mention that there is a manual inside the package. I successfully used a prior version, it needs the same configuration file that is used for Soapdenovo.
Your question may need to be reframed though. See my updated answer.
What is generally referred to as gap-closing software is a program that takes scaffolds as input, not contigs. Using PE reads, it fills the gaps present in scaffolds. Indeed, scaffolds are produced by some scaffolding software (not by user's guess), and they contain contigs in correct order and orientation.
Yes. Except that IMAGE is actually an hybrid method which does scaffolding and gap-filling at the same time. You did not mention the type of data you have. But if you have Illumina, I am unsure whether IMAGE will do a good job compared to newer methods. Thus, I'd recommend using the gapclosers mentioned in the responses.
Log in to answer this question.
Does anyone have access to CloG and is willing to share? Only if the license permits though, please.