20171103-用single_rarefaction.py为每个样本截取相同的序列数

代码示例:

1. single_rarefaction.py -i OTU_table.biom -o OTU_table_even10000.biom -d 10000

参数: -i 后接输入文件,通常为biom格式

-o 后接输出文件,即为稀释后的文件

-d 表示要稀释为的深度,即序列数,这里表示将每个样本稀释到10000个序列数

2. 将biom格式转为txt文本

/home/Hemaozhang/anaconda3/envs/qiime1/bin/biom convert -i OTU_table_even10000.biom -o OTU_table_even10000.txt --to-tsv --header-key taxonomy

参数:  -i, -o 分别接输入输出文件

--to-tsv  指定了输入格式,这个参数一定要有,否则会报错

--header-ley taxonomy 表示在输出文件的最后一列加上微生物分类信息

最终生成如下文件:

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容

  • 《ijs》速成开发手册3.0 官方用户交流:iApp开发交流(1) 239547050iApp开发交流(2) 10...
    叶染柒丶阅读 5,330评论 0 7
  • Spring Cloud为开发人员提供了快速构建分布式系统中一些常见模式的工具(例如配置管理,服务发现,断路器,智...
    卡卡罗2017阅读 134,973评论 19 139
  • 1. Java基础部分 基础部分的顺序:基本语法,类相关的语法,内部类的语法,继承相关的语法,异常的语法,线程的语...
    子非鱼_t_阅读 31,767评论 18 399
  • 《ilua》速成开发手册3.0 官方用户交流:iApp开发交流(1) 239547050iApp开发交流(2) 1...
    叶染柒丶阅读 10,997评论 0 11
  • linux资料总章2.1 1.0写的不好抱歉 但是2.0已经改了很多 但是错误还是无法避免 以后资料会慢慢更新 大...
    数据革命阅读 12,239评论 2 33