<p>My understanding is that <a href="https://sourceforge.net/p/pb-jelly/wiki/Home/" rel="nofollow noreferrer">PBJelly</a> is the tool intended to solve this problem.</p>
<p><a href="https://eagle.fish.washington.edu/trilobite/notebooks/sam/onsnetwork.org/kubu4/tag/pbjelly/index.html" rel="nofollow noreferrer">This person</a> seems to have used it starting with Illumina scaffolds. I have not personally used it but have heard that it works pretty well.</p>
<p>Edit: note that the PacBio ecosystem is not super easy to use from scratch. However, they have made pretty good efforts to facilitate setup using conda, as for instance documented <a href="https://github.com/PacificBiosciences/pbbioconda" rel="nofollow noreferrer">here</a>. <a href="https://github.com/cgjosephlee/PBJelly" rel="nofollow noreferrer">This person</a> shows a conda recipe for pbjelly.</p>

CC BY-SA 4.0
2020-08-16T23:21:20.277Z
2020-08-16T23:21:20.277Z
Maximilian Press
<p>I would like to gap-fill (or correct) a mammalian-sized assembly for which we have BACs that have been sequenced with PacBio CLR data. I have seen there are pbm22/gcpp tools available for correcting assemblies, but looking at how they've been used, they seem to correct existing assemblies built with FALCON or similar.</p>
<p>Can pbmm2 and gcpp be used with a reference constructed without pacbio data? E.g. given a mammalian-sized reference and a .bam file of PacBio reads for a BAC, what's the best way to gap-fill the current reference? Thx</p>

<assembly><reference-genome><pacbio>
Gap-fill assembly with PacBio reads
# Gap-fill assembly with PacBio reads

(Asked by: username_0 on 2020-08-10)

I would like to gap-fill (or correct) a mammalian-sized assembly for which we have BACs that have been sequenced with PacBio CLR data. I have seen there are pbm22/gcpp tools available for correcting assemblies, but looking at how they've been used, they seem to correct existing assemblies built with FALCON or similar.

Can pbmm2 and gcpp be used with a reference constructed without pacbio data? E.g. given a mammalian-sized reference and a .bam file of PacBio reads for a BAC, what's the best way to gap-fill the current reference? Thx

**Tags:** assembly, reference-genome, pacbio

---

## 1 answers

### Answer by username_1 on 2020-08-16. Score: 1

My understanding is that [PBJelly](https://sourceforge.net/p/pb-jelly/wiki/Home/) is the tool intended to solve this problem.

[This person](https://eagle.fish.washington.edu/trilobite/notebooks/sam/onsnetwork.org/kubu4/tag/pbjelly/index.html) seems to have used it starting with Illumina scaffolds. I have not personally used it but have heard that it works pretty well.

Edit: note that the PacBio ecosystem is not super easy to use from scratch. However, they have made pretty good efforts to facilitate setup using conda, as for instance documented [here](https://github.com/PacificBiosciences/pbbioconda). [This person](https://github.com/cgjosephlee/PBJelly) shows a conda recipe for pbjelly.


---

