Banana Slug Genomics

This page is read only. You can view the source, but not change it. Ask your administrator if you think this is wrong.
====== Quake ======
[[http://www.cbcb.umd.edu/software/quake/index.html|Quake]] is a sequencing error correction program that uses kmer count distributions to make the distinction between trusted (suspected true) and untrusted (suspected erroneous) reads.[(cite:quake>Kelley, David R; Schatz, Michael C; and Salzberg, Steven L. Quake: quality-aware detection and correction of sequencing errors. Genome Biology 2010, 11:R116 doi:[[http://dx.doi.org/10.1186/gb-2010-11-11-r116|10.1186/gb-2010-11-11-r116]])]

===== Running Quake =====
[[http://www.cbcb.umd.edu/software/quake/manual.html|Quake manual]]

To run Quake with the Jellyfish counts, it seems we will want to first use stats to dump the kmers as tab-separated and column (-t and -c respectively):

   jellyfish stats -t -c -o [counts file]

Then, build the covariance model to find the cutoff, default is 200:1 likelihood ratio of untrusted to trusted reads (same as plotting the histogram from Jellyfish):

   cov_model.py [counts file]

Finally, correct the reads:

   correct -f [fastq list file] -k [k-mer size] -c [cutoff] -m [counts file] -p [number of cores]

[Kevin] I think that we want to select the k-mer size manually, rather than relying on quake.  Their default cutoff is very conservative, and we'll probably do better over-correcting than under-correcting.
===== Methods =====
The paper gives an example of the distribution typically seen in kmer counting below.  They fit a Gamma distribution to the untrusted reads, and a Gaussian + Zeta (for the high frequency repeats) mixture for trusted reads.  The distribution of the trusted reads is actually expected to be Poisson, but the variance is significantly larger than the mean due to sequencing biases.  Their example is a high-coverage run, and they note that the distributions will overlap--as we've seen--with lower coverages.  This of course makes finding the cutoff trickier.  They chose a point where there was a high likelihood ratio of untrusted to trusted read, around 1000:1.

{{:bioinformatic_tools:quake_kmer_distribution.jpg|}}

They also recommend a kmer probability of 0.01 in a random sequence that is as long as the genome. That is, 2*G/4^k ~ 0.01, where G is the size of the sequenced genome and k is the size of the kmer. Simplified, k ~ log4(200*G), which is about 19 for our prediction of the banana slug genome size, {{https://banana-slug.soe.ucsc.edu/bioinformatic_tools:jellyfish|2.042e+09}}.
Banana Slug Genomics

User Tools

Site Tools

Page Tools
