温馨提示×

温馨提示×

您好,登录后才能下订单哦!

密码登录×
登录注册×
其他方式登录
点击 登录注册 即表示同意《亿速云用户服务条款》

qiime2分类器如何构建

发布时间:2022-03-21 10:43:34 来源:亿速云 阅读:277 作者:iii 栏目:开发技术

本文小编为大家详细介绍“qiime2分类器如何构建”,内容详细,步骤清晰,细节处理妥当,希望这篇“qiime2分类器如何构建”文章能帮助大家解决疑惑,下面跟着小编的思路慢慢深入,一起来学习新知识吧。

qiime2 18S 分类器构建

##wget -c https://www.arb-silva.de/fileadmin/silva_databases/qiime/SILVA_132_QIIME_release.zip
####unzip SILVA_132_QIIME_release.zip
qiime tools import \
 --type 'FeatureData[Sequence]' \
 --input-path SILVA_132_QIIME_release/rep_set/rep_set_18S_only/99/silva_132_99_18S.fna \
 --output-path 18S_SILVA132_99_otus.qza
qiime tools import \
 --type 'FeatureData[Taxonomy]' \
 --input-format HeaderlessTSVTaxonomyFormat \
 --input-path SILVA_132_QIIME_release/taxonomy/18S_only/99/majority_taxonomy_all_levels.txt \
 --output-path 18S_SILVA132_99_taxonomy.qza

#长度过滤
qiime rescript filter-seqs-length-by-taxon \
    --i-sequences  18S_SILVA132_99_otus.qza \
    --i-taxonomy 18S_SILVA132_99_taxonomy.qza \
    --p-labels  Eukaryota \
    --p-min-lens  1400 \
    --o-filtered-seqs 18S_SILVA132_nr99-seqs-filt.qza \
    --o-discarded-seqs 18S_SILVA132_nr99-seqs-discard.qza
#重复序列合并
qiime rescript dereplicate \
    --i-sequences 18S_SILVA132_nr99-seqs-filt.qza  \
    --i-taxa 18S_SILVA132_99_taxonomy.qza \
    --p-rank-handles 'silva' \
    --p-mode 'uniq' \
    --o-dereplicated-sequences 18S_SILVA132-nr99-seqs-derep-uniq.qza \
    --o-dereplicated-taxa 18S_SILVA132-tax-derep-uniq.qza
#全长分类器构建
qiime feature-classifier fit-classifier-naive-bayes \
  --i-reference-reads  18S_SILVA132-nr99-seqs-derep-uniq.qza \
  --i-reference-taxonomy 18S_SILVA132-nr99-tax-derep-uniq.qza \
  --o-classifier 18S_SILVA132-nr99-classifier.qza

###http://www.earthmicrobiome.org/protocols-and-standards/18s/ 39
qiime feature-classifier extract-reads  \
   --i-sequences 18S_SILVA132_99_otus.qza \
 --p-f-primer GTACACACCGCCCGTC \
 --p-r-primer TGATCCTTCTGCAGGTTCACCTAC \
 --p-max-length 150 \
 --o-reads ref-seqs_18S_99_SILVA132_RL150.qza
qiime feature-classifier fit-classifier-naive-bayes \
  --i-reference-reads ref-seqs_18S_99_SILVA132_RL150.qza \
  --i-reference-taxonomy 18S_SILVA132_99_taxonomy.qza \
  --o-classifier 18S_EMB_SILVA132_classifier.qza
qiime feature-classifier fit-classifier-naive-bayes \
 --i-reference-reads 18S_SILVA132_99_otus.qza \
 --i-reference-taxonomy 18S_SILVA132_99_taxonomy.qza \
 --o-classifier 18S_full_SILVA132_classifier.qza

读到这里,这篇“qiime2分类器如何构建”文章已经介绍完毕,想要掌握这篇文章的知识点还需要大家自己动手实践使用过才能领会,如果想了解更多相关内容的文章,欢迎关注亿速云行业资讯频道。

向AI问一下细节

免责声明:本站发布的内容(图片、视频和文字)以原创、转载和分享为主,文章观点不代表本网站立场,如果涉及侵权请联系站长邮箱:is@yisu.com进行举报,并提供相关证据,一经查实,将立刻删除涉嫌侵权内容。

AI