SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Pacific Biosciences
Similar Threads
Thread Thread Starter Forum Replies Last Post
Genome Res De novo bacterial genome sequencing: millions of very short reads assembly b_seite Literature Watch 1 10-04-2017 11:26 PM
DeNovo assembly using pacBio data krittika.sasmal Pacific Biosciences 50 06-05-2013 09:56 AM
De novo assembly of PacBio with short Illumina data boetsie Pacific Biosciences 1 10-06-2012 12:35 PM
hybrid assembly for PacBio and Illumina data laelaps Bioinformatics 1 05-01-2012 05:35 AM
Hybrid assembly of PacBio and Illumina allo Bioinformatics 3 05-01-2012 05:27 AM

Reply
 
Thread Tools
Old 10-22-2012, 04:12 AM   #1
lemur2
Junior Member
 
Location: China

Join Date: Aug 2010
Posts: 3
Default Have a genome assembly, what should I do with 15x Pacbio reads?

Hi everybody,
I am working to improve the assembly of a 400 Mb or so eukaryotic genome. We've done a fair amount of Illumina sequencing already, something like 60X coverage, and we've assembled the genome to an N50 of around 100,000 nt.

So, now we've got another 15X coverage of Pacbio reads, and I'm wondering what to do with them. Should we use them to pull the already assembled contigs together? If so, what software should we use?

Or, should we start from zero and do another de novo assembly, this time a hybrid assembly using both Illumina and Pacbio reads?

If it matters at all, we had a fair amount of bacterial contamination in the DNA sample, but I think we've managed to get rid of most of it using depth of coverage (the bacterial reads were of much higher coverage so were fairly easy to spot).

Any advice welcome. Thanks!
lemur2 is offline   Reply With Quote
Old 10-23-2012, 07:12 AM   #2
jbingham
Member
 
Location: Silicon Valley

Join Date: Jul 2011
Posts: 24
Default

At this time, there are really only two good options given the genome size and combination of data: pacBioToCA with Celera Assembler, or pb-jelly to fill gaps and join contigs.

The pacBioToCA approach first aligns the Illumina reads against the longest PacBio reads, then makes a consensus. After that, it assembles the raw Illumina data plus consensus PacBio data together.

pb-jelly uses your existing Illumina assembly and fills in the gaps with PacBio.
jbingham is offline   Reply With Quote
Old 10-25-2012, 09:06 PM   #3
lemur2
Junior Member
 
Location: China

Join Date: Aug 2010
Posts: 3
Default

Thanks jbingham - I think pb-jelly seems most appealing at the moment, I'll give that a try. Cheers!
lemur2 is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off



All times are GMT -8. The time now is 06:37 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2022, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO