<p>Try <a href="http://www.fishbrowser.org/software/LR_Gapcloser/index.php/Home/Index/Downloads.html" rel="nofollow noreferrer">LR_Gapcloser</a>.</p>

<p>I've used L_RNA_Scaffolder for trying to scaffold a genome (which turned out to be more complex than I had expected). It looks like LR_Gapcloser (written by the same people) is similar, but designed specifically for scaffolding using long-reads.</p>

<p>That page also suggests <a href="https://sourceforge.net/projects/pb-jelly/" rel="nofollow noreferrer">PBJelly</a> and <a href="https://sourceforge.net/projects/gmcloser/" rel="nofollow noreferrer">GMCloser</a> as competing tools.</p>

CC BY-SA 3.0
2017-11-27T18:08:37.770Z
2017-11-27T18:08:37.770Z
gringer
<p>I am assembling a ~500MB genome, and have ~150x long reads and ~200x 150bp PE short reads, with a ~400bp insert size.</p>

<p>I've done a lot of work with minimap+miniasm, and have what I think is a good set of unitigs and a .gfa file from miniasm. Depending on the parameters, I have a few hundred to ~1K unitigs. I can easily polish these unitigs with pilon, racon, and nanopolish, so that they are as error free as possible.</p>

<p>The question I have is: given access to</p>

<ol>
<li>Polished (or not) unitigs</li>
<li>150x long reads (long enough to span almost all repeats in the genome)</li>
<li>200x 150bp PE short reads (~400bp insert size)</li>
</ol>

<p>What would you do to scaffold a genome? And what are the pros and cons of your suggested approach?</p>

<p>What I need here is a draft genome. I would prefer a conservative approach that is more fragmented and has fewer misassemblies. So of course one option is just to keep the unitigs and make no attempt at scaffolding. But it does strike me there may be suitably conservative approaches I haven't considered. E.g. perhaps I could use ABYSS to scaffold my polished unitigs using my short-read data, thereby leveraging whatever information I can from the PE short reads?</p>

CC BY-SA 3.0
2017-11-27T16:11:13.663Z
You don't specify how long the long reads are, but if they are longer than 700 bp, then the PE reads won't be much help for scaffolding as they will contain less long-range info than the long reads. The short reads will probably be most useful for polishing the long reads.
heathobrien
CC BY-SA 3.0
2017-11-27T19:48:27.227Z
Good point. The long reads are almost all >10KB. I guess I had a vague notion that there might be some useful information in the PE reads.
roblanf
<assembly><long-reads><illumina><scaffold>
Scaffolding a genome with hybrid data
# Scaffolding a genome with hybrid data

(Asked by: username_0 on 2017-11-27)

I am assembling a ~500MB genome, and have ~150x long reads and ~200x 150bp PE short reads, with a ~400bp insert size.

I've done a lot of work with minimap+miniasm, and have what I think is a good set of unitigs and a .gfa file from miniasm. Depending on the parameters, I have a few hundred to ~1K unitigs. I can easily polish these unitigs with pilon, racon, and nanopolish, so that they are as error free as possible.

The question I have is: given access to

1. Polished (or not) unitigs
2. 150x long reads (long enough to span almost all repeats in the genome)
3. 200x 150bp PE short reads (~400bp insert size)

What would you do to scaffold a genome? And what are the pros and cons of your suggested approach?

What I need here is a draft genome. I would prefer a conservative approach that is more fragmented and has fewer misassemblies. So of course one option is just to keep the unitigs and make no attempt at scaffolding. But it does strike me there may be suitably conservative approaches I haven't considered. E.g. perhaps I could use ABYSS to scaffold my polished unitigs using my short-read data, thereby leveraging whatever information I can from the PE short reads?

**Tags:** assembly, long-reads, illumina, scaffold
> **Comments:**
> **username_1**: You don't specify how long the long reads are, but if they are longer than 700 bp, then the PE reads won't be much help for scaffolding as they will contain less long-range info than the long reads. The short reads will probably be most useful for polishing the long reads.  
> **username_0**: Good point. The long reads are almost all >10KB. I guess I had a vague notion that there might be some useful information in the PE reads.  

---

## 1 answers

### Answer by username_2 on 2017-11-27. Score: 2 | Accepted answer

Try [LR\_Gapcloser](http://www.fishbrowser.org/software/LR_Gapcloser/index.php/Home/Index/Downloads.html).

I've used L\_RNA\_Scaffolder for trying to scaffold a genome (which turned out to be more complex than I had expected). It looks like LR\_Gapcloser (written by the same people) is similar, but designed specifically for scaffolding using long-reads.

That page also suggests [PBJelly](https://sourceforge.net/projects/pb-jelly/) and [GMCloser](https://sourceforge.net/projects/gmcloser/) as competing tools.


---

