QIIME2学习笔记（二）
1. 序列质控及Feature表构建
导入数据之后便是质控了。
这次不用测试数据了，用实际数据跑一下，所以同样重复之前的步骤，把fastq文件压缩下，然后，生成样本数据列表（ps.不知道fastq文件不压缩可不可以用，有空试下）。依然参考自上篇的统计咨询公众号的文章代码如下：
gzip * #gzip压缩文件夹里所有的数据
source activate qiime2-2018.2 #激活qiime2环境
##1## import data
qiime tools import \
--type 'SampleData[PairedEndSequencesWithQuality]' \
--input-path 171213_16s-manifest \
--output-path 171213_16s.qza \
--source-format PairedEndFastqManifestPhred33
##2## quality control
#visualization
qiime demux summarize \
--i-data 171213_16s.qza\
--o-visualization 171213_16s.qzv
##3##filter 聚类
qiime dada2 denoise-paired \
--i-demultiplexed-seqs 171213_16s.qza \
--p-trunc-len-f 0 \
--p-trunc-len-r 0 \
--o-representative-sequences rep-seqs-dada2.qza \ #输出feature表代表序列
--o-table table-dada2.qza \ #输出表
--p-n-threads 36 #36线程，服务器是40线程的
#qzv 可视化
qiime feature-table summarize \
--i-table table-dada2.qza \
--o-visualization table-dada2.qzv主要作用是去除低质量序列、嵌合体；再生成OTU表，现在叫Feature表，因为不再使用聚类方法，相当于QIIME时代100%相似度的OTU表。–QIIME2中文帮助文档(https://forum.qiime2.org/t/qiime2-chinese-manual/838)
2.导出数据
这里顺便学习一下导出数据，毕竟qiime2全部是专有格式的文件，不导出没办法在外部查看。到这里，没有中文的相关参考了，直接去学习qiime2官方文档。命令比较简单，就是qiime tools export 后面加上文件等参数就行了。
qiime tools export \
feature-table.qza \
--output-dir exported-feature-table运行后导出了一个biom格式的文件。还是比较亲切的格式。
3.metadata 实验设计
文档介绍说类似于qiime1的mappingfile，只不过是tsv格式的，还兼容mapping文件的。文件长得这样：
这个文件对于我来说好你暂时没有用到，所以暂时没有创建，不知道后边会不会卡脖子，因为我用到的是分好的样品数据，没有进行分数据的步骤。
本文参与 腾讯云自媒体同步曝光计划，分享自微信公众号。
原始发表：2018-04-24，如有侵权请联系 cloudcommunity@tencent.com 删除
评论
登录后参与评论
推荐阅读
目录
相关产品与服务
文件存储
文件存储（Cloud File Storage，CFS）为您提供安全可靠、可扩展的共享文件存储服务。文件存储可与腾讯云服务器、容器服务、批量计算等服务搭配使用，为多个计算节点提供容量和性能可弹性扩展的高性能共享存储。腾讯云文件存储的管理界面简单、易使用，可实现对现有应用的无缝集成；按实际用量付费，为您节约成本，简化 IT 运维工作。