Shasta long read assembler
This is the documentation for the
Shasta long read assembler.
If you are seeing this documentation on
it applies to the latest version of the Shasta assembler
on GitHub (not necessarily the same as the latest release).
Documentation for any version of Shasta is available
in the source code in the
and in any build under the
The goal of the Shasta long read assembler is to rapidly produce accurate assembled sequence using as input DNA reads generated by Oxford Nanopore flow cells.
Computational methods used by the Shasta assembler include:
- Using a run-length representation of the read sequence. This makes the assembly process more resilient to errors in homopolymer repeat counts, which are the most common type of errors in Oxford Nanopore reads.
- Using in some phases of the computation a representation of the read sequence based on markers, a fixed subset of short k-mers (k ≈ 10).
An initial implementation of the Shasta assembler is complete and functional,
but significant improvements in several areas are possible.
As currently implemented, it can run an assembly of a human genome at coverage around 60x
in about 5 hours using a single, large machine (AWS instance type
x1.32xlarge, with 128 virtual processors and 1952 GB of memory).
The compute cost of such an assembly is around $20 at AWS spot market or reserved prices.
The accuracy of assembled sequence is being analyzed. Early indications are that Shasta is similar or better in assembly quality when compared to other long read assemblers.
Command line options
- Run an assembly
- Run an assembly in Docker
- Maximize assembly performance
- Build the code from source
- Explore assembly results
- Contribute to Shasta
- Report problems or ask questions