[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