Wednesday, November 21, 2012

Basic Installation of Quake - Package to correct substitution sequencing errors in experiments with deep coverage

What is Quake?
(Taken from Quake Site)

Quake is a package to correct substitution sequencing errors in experiments with deep coverage (e.g. >15X), specifically intended for Illumina sequencing reads. Quake adopts the k-mer error correction framework, first introduced by the EULER genome assembly package. Unlike EULER and similar progams, Quake utilizes a robust mixture model of erroneous and genuine k-mer distributions to determine where errors are located. Then Quake uses read quality values and learns the nucleotide to nucleotide error rates to determine what types of errors are most likely. This leads to more corrections and greater accuracy, especially with respect to avoiding mis-corrections, which create false sequence unsimilar to anything in the original genome sequence from which the read was taken.

Setting up is quite straight-forward, just untar in an appropriate directory.
# tar -zxvf quake-0.3.4.tar.gz
# cd Quake\src

Edit the Makefile if you are using Linux (Link CFLAGS to Boost Directory). Boot Software can be downloaded at Boost C++ Libraries
CC=g++
CFLAGS=-O3 -fopenmp -I/usr/local/boost/include/boost -I.
LDFLAGS=-L. -lgzstream -lz
.....
.....

To complete the installation, do a make at the src
 # make

You should see executable in the src

No comments: