Qiime安装

参考资料:http://blog.sina.com.cn/s/blog_83f77c940101h2rp.html

Qiime script官方说明http://qiime.org/scripts/index.html

一般ITS数据简要分析过程

reference下载:https://github.com/downloads/qiime/its-reference-otus/its_12_11_otus.tar.gz

其中97_* 和99_*分别是相似度97和99的reference

以下过程仅做参考,具体根据详细项目更改,

1   质控 

对数据进行质量控制,去除低质量数据

2   拼接

针对不同平台数据,采用相应的拼接软件拼接,根据contig长度,方差,N50(主要)等指标使用最佳k-mer值

3   pick_otus.py

官方文档:http://qiime.org/scripts/pick_otus.html

从fasta文件中提取OTUs,提取方法可选,具体有cd-hit,blast,Mothur,usearch等等

具体使用方法:

python pick_otus.py -i seqs.fna -o picked_otus_default

-i参数指定输入序列,-o参数指定输出文件目录,输出文件为seqs_otus.txt和seqs_otus.log,-m可指定提取方法,更多具体参数说明参见官方文档

4   assign_taxonomy.py

官方文档:http://qiime.org/scripts/assign_taxonomy.html

使用下载的reference中的rep_set和 taxonomy文件夹中对应的文件,对序列分类,使用方法:

python assign_taxonomy.py –i seqs.fna -r <otus_dir>/rep_set/97_otus.fasta -t <otus_dir>/taxonomy/97_otu_taxonomy.txt

更多具体参数说明参见官方文档,输出文件示例:tax_assignments.txt

5   make_otu_table.py

官方文档:http://qiime.org/scripts/make_otu_table.html

将提取出的seqs_otus.txt和分类文件tax_assignments.txt制作成out_table,即.biom文件

使用方法:

python make_otu_table.py -i otu_map.txt -t tax_assignments.txt -o otu_table.biom

6   pick_open_reference_otus.py

官方文档:http://qiime.org/scripts/pick_open_reference_otus.html

以上三个步骤可以通过Qiime中自带的pick_open_reference_otus.py流程完成,将需要使用到的参数写到文件params.txt中,并通过-p参数调用,param文件说明文档: http://qiime.org/documentation/qiime_parameters_files.html

7   summarize_taxa_through_plots.py

官方文档:http://qiime.org/scripts/summarize_taxa_through_plots.html

统计各个样品中的OUT中物种丰度信息,并从门纲目科属水平分别绘制物种丰度bar图

python summarize_taxa_through_plots.py -o taxa_summary -i otu_table.biom 在输出文件夹中,双击html文件即可观察bar图

8   alpha_rarefaction.py

官方文档:http://qiime.org/scripts/alpha_rarefaction.html

对拼接后的reads

使用示例:python alpha_rarefaction.py -i otu_table.biom -o arare/ -m map.txt,对reads进行随机抽样,得到的序列数和out数目比较,做出稀释曲线图,得到同样可以通过html文件查看

9   alpha_diversity.py

官方文档:http://qiime.org/scripts/alpha_diversity.html

使用示例:python alpha_diversity.py -i otu_table.biom -m chao1 -o adiv_chao1.txt,单样品多样性分析,绘制Alpha多样性bar图

10 统计.biom文件长度,并选取最短长度作为beta多样性的深度参数

可以借助ipython工具中提供的biom summarize-table方法,选取最短长度n,作为下一步骤中-e 参数(覆盖深度)

11 beta_diversity_through_plots.py

官方文档:http://qiime.org/scripts/beta_diversity_through_plots.html

该流程包含beta_diversity.py principal_coordinates.py make_prefs_file.py三个脚本,能够给出(un)weignted unifrac距离矩阵,并做出2D和3D的PCoA图形

使用示例:python beta_diversity_through_plots.py -i otu_table.biom -o bdiv_even100/ -m map.txt -e 100,使用-p参数指定距离范式

ITS简要分析流程(using Qiime)的更多相关文章

  1. RxJava && Agera 从源码简要分析基本调用流程(2)

    版权声明:本文由晋中望原创文章,转载请注明出处: 文章原文链接:https://www.qcloud.com/community/article/124 来源:腾云阁 https://www.qclo ...

  2. Activity源码简要分析总结

    Activity源码简要分析总结 摘自参考书籍,只列一下结论: 1. Activity的顶层View是DecorView,而我们在onCreate()方法中通过setContentView()设置的V ...

  3. Android Hal层简要分析

    Android Hal层简要分析 Android Hal层(即 Hardware Abstraction Layer)是Google开发的Android系统里上层应用对底层硬件操作屏蔽的一个软件层次, ...

  4. 转:InnoDB多版本(MVCC)实现简要分析

    InnoDB多版本(MVCC)实现简要分析 基本知识 假设对于多版本(MVCC)的基础知识,有所了解.InnoDB为了实现多版本的一致读,采用的是基于回滚段的协议. 行结构 InnoDB表数据的组织方 ...

  5. Android 5.1 Settings源代码简要分析

    转载请注明出处,谢谢~http://blog.csdn.net/u011974987/article/details/51004854. 概述: 先声明:本人工作快两年了,仍是菜鸟级别的.羞愧啊!曾经 ...

  6. Google发布SSLv3漏洞简要分析报告

    今天上午,Google发布了一份关于SSLv3漏洞的简要分析报告.根据Google的说法,该漏洞贯穿于所有的SSLv3版本中,利用该漏洞,黑客可以通过中间人攻击等类似的方式(只要劫持到的数据加密两端均 ...

  7. 构建ASP.NET MVC4+EF5+EasyUI+Unity2.x注入的后台管理系统(34)-文章发布系统①-简要分析

    原文:构建ASP.NET MVC4+EF5+EasyUI+Unity2.x注入的后台管理系统(34)-文章发布系统①-简要分析 系列目录 最新比较闲,为了学习下Android的开发构建ASP.NET ...

  8. CVE-2015-5122 简要分析(2016.4)

    CVE-2015-5122 简要分析 背景 最近在学习Flash漏洞的分析,其与IE漏洞的分析还是有诸多的不同(不便)之处,折腾了一阵子终于克服了没有符号表.Flash的超时定时器等问题.所以找到了去 ...

  9. Java7中的ForkJoin并发框架初探(中)——JDK中实现简要分析

    原文发表于 2013 年 8 月 28 日 由 三石 根据前文描述的Doug Lea的理论基础,在JDK1.7中已经给出了Fork Join的实现.在Java SE 7的API中,多了ForkJoin ...

随机推荐

  1. jmeter接口系列:时间戳、加密

    JMeter安装配置 从官网下载JMeter的软件包apache-jmeter-x.x.zip,下载完成之后解压打开jmeter.bat即可. 说明 这里使用的jmeter版本是3.0,jdk版本是j ...

  2. navigator.userAgent浏览器检测(前端基础系列)

    对于前端来说,浏览器检测已经不陌生了,在做一些页面是,需要针对不同的浏览器进行处理不同的逻辑,最简单的就是区分pc和移动端的浏览器,或是android 和ios下的浏览器. 一.浏览器检测的由来?  ...

  3. java学习总结篇二--3 种简单排序

    本篇文章,先从数据结构开始,一边总结,一边反思,寻求最优解. 本文简单温习下最基础的三类算法:选择,冒泡,插入.先定义一个交换数组作为备用: /** * 交换数组元素 * @param arr * @ ...

  4. [亲测]ASP.NET Core 2.0怎么发布/部署到Ubuntu Linux服务器并配置Nginx反向代理实现域名访问

    前言 ASP.NET Core 2.0 怎么发布到Ubuntu服务器?又如何在服务器上配置使用ASP.NET Core网站绑定到指定的域名,让外网用户可以访问呢? 步骤 第1步:准备工作 一台Liun ...

  5. tomcat线程初探

    博主:handsomecui,希望路过的各位大佬留下你们宝贵的意见,在这里祝大家冬至快乐. 缘由: 初探缘由,在业务层想要通过(当前线程的栈)来获取到控制层的类名,然后打日志,可是发现并不能通过当前线 ...

  6. 使用CEF的JSON解析功能

    Cef提供了JSON解析功能,在cef_parser.h文件内有三个JSON相关的方法: CefParseJSON CefParseJSONAndReturnError CefWriteJSON 以最 ...

  7. 技术债务管理以及Firefox/Chromium的债务评价

    如今的软件开发是在遍地敏捷,人人讲唯快不破的时代,哪有人有时间思考代码质量,设计的质量? 哪个又不是从一堆代码中杀出血路来实现还有一个功能?一个产品都存活不了几年,何必考虑什么可维护性? 我们追求进度 ...

  8. 【转】JAVA处理线程超时

    在实际业务中,由其是多线程并开业务中,经常会遇到某个线程执行超时.而程序如果不捕获这类情况,就会导致程序一直处于等待状态,从而影响后续线程的运行.比如说网络通迅.单任务下的复杂数据库查询等,通常处理这 ...

  9. Java并发编程实践读书笔记(1)线程安全性和对象的共享

    2.线程的安全性 2.1什么是线程安全 在多个线程访问的时候,程序还能"正确",那就是线程安全的. 无状态(可以理解为没有字段的类)的对象一定是线程安全的. 2.2 原子性 典型的 ...

  10. Ignite与Spark集成时,ClassNotFoundException问题解决

    参考文章:https://apacheignite-fs.readme.io/docs/installation-deployment Spark application deployment mod ...