Question: How to use pacbio reads to join contigs
0
gravatar for 934963534
2.3 years ago by
93496353420
93496353420 wrote:

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!

pacbio long assembly tgs • 1.2k views
ADD COMMENTlink modified 2.3 years ago by h.mon30k • written 2.3 years ago by 93496353420

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

ADD REPLYlink written 2.3 years ago by Rohit1.4k
2
gravatar for h.mon
2.3 years ago by
h.mon30k
Brazil
h.mon30k wrote:

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.

ADD COMMENTlink written 2.3 years ago by h.mon30k
Please log in to add an answer.

Help
Access

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