最近在搞16S，发现了一个实践算法的最佳机会。

见文章：

A Bayesian taxonomic classification method for 16S rRNA gene sequences with improved species-level accuracy.

文章利用了贝叶斯模型，调用了blast和muscle来对OTU进行taxonomy assignment。

可以看一下源代码，非常简单。

Bayesian-based LCA taxonomic classification method

如果你能彻底理解本文的算法，并能看懂其源码，那你基本就打到了生信算法入门的水平。

说不定以后你也能随手发一个算法的文章哦！

生信算法实践