最近在搞16S，发现了一个实践算法的最佳机会。
见文章：
文章利用了贝叶斯模型，调用了blast和muscle来对OTU进行taxonomy assignment。
可以看一下源代码，非常简单。
Bayesian-based LCA taxonomic classification method
如果你能彻底理解本文的算法，并能看懂其源码，那你基本就打到了生信算法入门的水平。
说不定以后你也能随手发一个算法的文章哦！
原文：http://www.cnblogs.com/leezx/p/7146668.html