Best software to assemble bacterial genomes

I have sequencing data of five bacteria, which were generated using Illumina MiSeq. Four of them were sequenced using a paired-end 2x300 protocol and one was sequenced using the nextera mate-pair protocol.

My question is: What are the softwares that you recommend me to assemble these genomes (the largest has almost 8 Mbp)?

I have access to a CLC Workbench. It seems quite ease to use, but I dont know if it is the best one. Most of papers that I found that evaluate the performance of assemblers from two-three years ago.

I also have to mention that I have two i7 with 8 threads PCs available for this objective (one with 32 and another with 8 Gb RAM).

Thanks in advance.

illumina assembler assembly bacteria

• 19,377 views • link

updated 4.5 years ago by Ram

• written 11.7 years ago by fhsantanna

Not exactly what you're looking for, but this guide was very useful for me: Beginner’s guide to comparative bacterial genome analysis using next-generation sequence data.

• 0 views • link

written 11.7 years ago by PoGibas

I already read this paper. Very good. They recommended Velvet, but I believe that there are better options. Thank you anyway.

• 0 views • link

written 11.7 years ago by fhsantanna

I like using CLC Workbench. Even if you don't use their assembler you can do the last steps in CLC, it's much more convenient. Try SPAdes, SOAPdenovo and you can compare it to CLC built-in assembler.

• 0 views • link

written 11.7 years ago by marina.v.yurieva

7 answers

You should read in the literature to know which one is the best one for you specific data. Here you have a nice paper comparing some assembly tools, and it is a recent paper (2014): http://genomebiology.com/2014/15/3/R42

In my opinion SOAPdenovo2 and SGA are a good choice. Bambus is quite difficult to install and to understand. SPACE also is nice, but if you want to use the last version you have to pay so...

Hope it helps.

• 0 views • link

written 11.7 years ago by iraun

For bacterial genome, GAGE-B paper (2013) compare 8 genome assemblers:

ABySS v1.3.4
CABOG v7.0
MIRA v3.4.0
MSRCA v1.8.3
SGA v0.9.34
SOAPdenovo2 v2.04 + GapCloser v1.12
SPAdes v2.3.0
Velvet v1.2.08

All GAGE-B data and assembly recipe are available here.

For more recent comparison of genome assemblers, have a look here.

As each bacteria genome size and GC% is different, you need to check these reproducible Benchmarks.

• 0 views • link

updated 4.5 years ago by Ram

• written 11.7 years ago by rtliu

I will second the recommendation of http://nucleotid.es/

It compares numerous assemblers on microbial genomes, with objective metrics as reported by Quast (a tool for evaluating assemblies). And it includes both peak memory usage and CPU-time.

• 0 views • link

updated 4.5 years ago by Ram

• written 11.7 years ago by Brian Bushnell

All the articles mentioned conclude with that there is no single best assembler for bacterial genomes. It depends on the genome and the data. So, you'll have to try a few, then validate them using tools such as FRCBam, REAPR or one of the likelihood methods. If you don't care about all this, use SPAdes. If you want a tool that automates most of this, look at iMetAMOS www.cbcb.umd.edu/software/imetamos

• 0 views • link

written 11.7 years ago by lexnederbragt

Yes, SPAdes performs very well and it's robust: I would reccomend using the --careful option which, according to the nucleotid.es benchmarks reduces the errors while keeping the same N50.

• 0 views • link

written 11.6 years ago by mgalactus

As mentioned SPAdes is great or check out the A5 Assembly pipeline. Following the full GAGE-B paper it has produced better QUAST results than SPAdes for MiSeq data. A lot depends upon your organism, how repetitive, GC content etc.

http://www.plosone.org/article/info%3Adoi%2F10.1371%2Fjournal.pone.0042304

EDIT: Also, check out Abacas for scaffolding if you have a closely related reference genome.

• 1 views • link

updated 6.8 years ago by Ram

• written 11.7 years ago by moorem