K-mer应用实例分析

We started with 600 million 100mer reads from a genomic library and plotted its K-mer distribution. It looked like the following plot.

K-mer应用实例分析

The peak around 52 suggests that non-repetitive regions of the genome are covered ~52 times by the reads. For details about the basic philosophy of what is being plotted and why, please read our earlier commentary on Maximizing Utility of Available RAMs in K-mer World.

When we added another 800 million 51mer reads into the pot, the K-mer distribution changed to the following image. The peak of two combined libraries moved to 58.

K-mer应用实例分析

Does that make sense? The new library has about 70% as many sequences as the old one. Shouldn’t the peak shift from 52 to 88?

Article From:

http://www.homolog.us/blogs/2011/10/27/sharing-a-k-mer-story/

  • 本文由 整理发表
  • 网站部分文章源自互联网,若未正确标注来源,请联系管理员更新。文章转载,请务必保留本文链接

您必须才能发表评论!