[khmer] zero abundance kmers
C. Titus Brown
ctb at msu.edu
Thu Jan 30 06:19:10 PST 2014
On Wed, Jan 29, 2014 at 09:15:24PM -0800, Aaron Liston wrote:
> I am using load-into-counting.py and abundance-dist.py to plot the kmer
> abundance distribution of read data at various stages in the digital
> normalization process. Some of the plots have a small fraction (<0.001)
> of kmers with zero abundance. Is this an error, or does it mean
> something?
Hi Aaron,
that should only happen if you're either missing some reads (e.g. if you
load data set A into a counting hash, and then count all the k-mers in data
set B) _or_ if you have Ns in your reads somewhere.
Gang, we should verify this and put it in a FAQ... thanks for pointing it out!
cheers,
--titus
--
C. Titus Brown, ctb at msu.edu
More information about the khmer
mailing list