This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to use pacbio reads to join contigs

Hello everyone!

I have some pre-assembled contigs and want to join and extend the contigs using pacbio reads. I tried DBG2OLC on the contigs but the result shows that the contigs are even shorter.

Does anybody know reasons or are there any other tools to join contigs with pacbio reads?

Thanks for your answer!

assembly tgs pacbio
• 2,948 views • link
updated by Wyclif • Wyclif • written by 934963534 • 934963534

Do you have enough coverage i.e. 10X for pacbio. Since DBG2OLC requires multiple alignment of the pacbio reads, if there is lower coverage, It might be easier to clarify if you provide the actual command used, since it doesn't help you have a higher kmer_coverage threshold and low sequencing coverage. If you have pacbio corrected data, then PBJelly scaffolds

• 1 views • link
written by Rohit Rohit

Hi, did you find the solution for this? I am new to bioinformatics here and I would like to use DBG2OLC to combine the contigs into scaffolds using PaBio RSII sequences. iny installation procedures or pipeline will be very helpful. Thanks.

• 0 views • link
written by Wyclif • Wyclif

1 answer

There are indeed some tools, PacBio site points to at least two, FGAP and GMcloser:

https://www.pacb.com/publications/fgap-an-automated-gap-closing-tool/

https://www.pacb.com/publications/gmcloser-closing-gaps-in-assemblies-accurately-with-a-likelihood-based-selection-of-contig-or-long-read-alignments/

What is the quality of the assembly? Are your PacBio reads error-corrected? Did you align the pre- and post-correction contigs? Although I don't know DBG2OLC, one possible reason for the "even shorter" contigs is that DBG2OLC could be breaking chimeric contigs, as per PacBio reads indication - but again, I don't even know if DBG2OLC performs such correction.

• 1 views • link
written by h.mon h.mon

Log in to answer this question.

More posts like this