106 qiime tools import \
107 --type 'FeatureData[Sequence]' \
108 --input-path gg_13_8_otus/rep_set/99_otus.fasta \
109 --output-path 99_otus.qza
110 导入物种分类信息
111 qiime tools import \
112 --type 'FeatureData[Taxonomy]' \
113 --input-format HeaderlessTSVTaxonomyFormat \
114 --input-path gg_13_8_otus/taxonomy/99_otu_taxonomy.txt \
115 --output-path ref-taxonomy.qza
116 5. 训练分类器—全长 (通用) ,耗时约半小时
117 time qiime feature-classifier fit-classifier-naive-bayes \
118 --i-reference-reads 99_otus.qza \
119 --i-reference-taxonomy ref-taxonomy.qza \
120 --o-classifier classifier_gg_13_8_99.qza
121 如果提示版本错误,详见常见问题 3。
122 6. (可选) 训练分类器—指定 V 区分类器
123 如果扩增了指定的 16S 区域,还可以构建特异区域的分类器,可进一步提高分类精
124 度。常用 GreenGenes 13_8 按 99%聚类操作分类单元 (Operational taxonomic un
125 its, OTUs) 序列中的 V4 区域 (341F CCTACGGGNGGCWGCAG/805R GACTA
126 CHVGGGTATCTAATCC) 构建分类器。确定此处使用的引物与扩增引物保持一致。
127 本次使用与测试数据对应的 V5 (799F) - V7 (1193R) 引物为例进行提取序列,
128 耗时约 9 分钟。
129 time qiime feature-classifier extract-reads \
130 --i-sequences 99_otus.qza \
131 --p-f-primer AACMGGATTAGATACCCKG \
132 --p-r-primer ACGTCATCCCCACCTTCC \
133 --o-reads ref-seqs.qza
134 基于筛选的指定区段,生成实验特异的分类器,耗时约 8 分钟。
135 time qiime feature-classifier fit-classifier-naive-bayes \
136 --i-reference-reads ref-seqs.qza \
评论0