我们知道,不管是16S等扩增子测序,还是宏基因组,最后最重要的结果,就是物种的丰度情况了,qiime2给出的16S丰度结果是一个计数,对于许多软件来说这是可用的,那么如果我们想获得一个直接的百分比数据应该怎样做呢?
这几天随便搜索snp2hla软件的参考数据集的时候发现一个韩国科学家写了一个数据集合并脚本,在使用韩国人样本测试时准确性较分别只用两个未合并的数据集准确性有所提高,于是,就找到了论文提供的脚本合并了一下。
之前做了一个简单的网页,使用23andme格式+snp2hla软件获得hla分型数据,当然准确性不咋的,也就玩玩,上线后为大约100+人提供了服务,这是伯值得骄傲的事,因为第一次能为大家提供服务。代码我是放在了gitHub的,数据是脚本处理完后自动删除。看网页是不是有点眼熟,这是谷歌中国网页框架,直接搬来的。