溫馨提示×

溫馨提示×

您好,登錄后才能下訂單哦!

密碼登錄×
登錄注冊×
其他方式登錄
點擊 登錄注冊 即表示同意《億速云用戶服務條款》

qiime2分類器如何構建

發布時間:2022-03-21 10:43:34 來源:億速云 閱讀:324 作者: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

亚洲午夜精品一区二区_中文无码日韩欧免_久久香蕉精品视频_欧美主播一区二区三区美女