kmer估计基因大小

网上实践教程很多，这里只讲理论部分。
base_coverage_depth = Total_base_num/genome_size=(read_num * read_length)/genome_size

Kmer_coverage_depth= Total_kmer_num/genome_size = read_num * (read_length - kmer_size + 1)/genome_size

Kmer_coverage_depth = base_coverage_depth * (read_length - Kmer_size + 1)/read_length
base_coverage_depth = Kmer_coverage_depth * (read_length - Kmer_size + 1) /read_length
由kmer频率分布直方图近似得到Kmer_coverage_depth，再计算出base_coverage_depth，最后得到genome_size。
https://www.jianshu.com/p/f03789bd4d87