Hi all,
Can anybody suggest how much pacbio coverage do we need for a genome of 1.7 giga bases (de novo assembly)? We already have about 47X illumina reads (paired end, mate pairs 9 and 11 Kb) and the assembly that we generated is not so great (scaffold N50 of 36Kb, but 50% Ns). Most importantly, the genome has ~50% repeat content. We have a shared budget that we need to split for both pacbio DNA sequencing and Iso-seq (for future genome annotation: this spp doesn't have any transcriptome data). What do you guys think is the ideal coverage we need aim for to get a quality assembly?
Any suggestions on this will be greatly appreciated!
1 answer
It depends a lot on what you are planning to do with your PacBio data. For assembly and error correction, there are different strategies:
- Ideal coverage required for PacBio error correction using HGAP + PacBio assembly (HGAP, FALCON, etc..) would requires >50X coverage.
- What tools you use or know for PacBio Long Read error correction? + some form of hybrid or low coverage PacBio assembly (FALCON). This works with 10-20X coverage + >50X Illumina. However correction can be computationally intensive for a 1.5GB genome and I haven't seen a adequate hybrid assembler for this size of genome.
- Scaffolding + Gapclosing (pbJelly2, sspace-longread, finisherSC) can be done with raw or error corrected PacBio reads and 10X coverage should be sufficient
Also, I've seen a lot of issues with large/repetitive plant genomes and library prep/size selection for PacBio. Make sure to test your protocols. Good subread length is very important if you do low coverage PacBio. With short subreads you may get a theoretical coverage of 10X, but since all subreads from one read stack onto each other, you'll have much less real overlapping fragments than you need for assembly/scaffolding
Log in to answer this question.
-
How to compare the quality of assemblies
written by kirillkirilenko •
-
Why does low-coverage library give more contiguous assembly than high-coverage library?
written by Timotheus •
-
BBmap statistic values
written by Dũng •
-
About 10x chromium Gene assembly
written by r00628112 •
-
Assembly statistics: ABySS & BBMap
written by el97004 •
-
How calculate CG content for each contig in an transcript assembly?
written by luzglongoria •
-
Why we need 100X coverage to get a high-quality assembly?
written by lamz138138 •
-
a question about -unitigs.fa and -scaffolds.fa
written by xiaofeiwang198266 •
-
i got two different trinity.fasta by using two version of trinity
written by Kurban
-
Larger-Than-Expected Total Scaffolds Size Using Nextera Mate Pairs
written by Rayan Chikhi
just as a gut feeling, to complement your Illumina-Seq I'd say 10-15x coverage. But I might be way off here!