溫馨提示×
您好，登錄后才能下訂單哦！
點擊 登錄注冊 即表示同意《億速云用戶服務條款》
您好，登錄后才能下訂單哦！
本篇內容介紹了“qiime2怎么使用”的有關知識，在實際案例的操作過程中，不少人都會遇到這樣的困境，接下來就讓小編帶領大家學習一下如何處理這些情況吧！希望大家仔細閱讀，能夠學有所成！
qiime2 vsearch OTU pip 整理
### 方法1. DADA2 feature table 雙端合并，去除嵌合體，截去接頭序列降噪生成feature table 一步完成 cd $workdir mkdir dada2 echo dada2 start date cd dada2 qiime dada2 denoise-paired \ --i-demultiplexed-seqs $workdir/demux.qza \ --p-n-threads 1 \ --p-trim-left-f 29 --p-trim-left-r 18 \ --p-trunc-len-f 0 --p-trunc-len-r 0 \ --o-table dada2-table_biom.qza \ --o-representative-sequences dada2-repset-seqs.qza \ --o-denoising-stats denoising-stats.qza echo dada2 end date
### 方法2. 外部導入特征表和代表序列(常用) cd $workdir mkdir import_table cd import_table # 上傳我們生成的OTU表otu_table.txt和代表序列rep_set.fa # 轉換文本為Biom1.0，注意biom --version 2.1.5/8可以，2.1.7報錯 qiime tools import --input-path $workdir/otu_table.biom \ --type 'FeatureTable[Frequency]' --input-format BIOMV100Format \ --output-path table_biom.qza qiime tools import --input-path $workdir/otus.fa \ --type 'FeatureData[Sequence]' \ --output-path repset-seqs.qza
### 方法3. deblur feature table deblur只能輸入合并好的序列， 速度比dada2快 #因此先做序列處理，去除引物，合并，過濾低質量等等 # Trim amplicon primers 去除引物序列 # 341f # 806r qiime cutadapt trim-paired \ --i-demultiplexed-sequences $workdir/1.import_data/demux.qza \ --p-cores 4 \ --p-no-indels \ --p-front-f CCTAYGGGRBGCASCAG \ --p-front-r GGACTACNNGGGTATCTAAT \ --o-trimmed-sequences primer-trimmed-demux.qza # Summarise the reads 查看結果 qiime demux summarize \ --i-data primer-trimmed-demux.qza \ --o-visualization primer-trimmed-demux.qzv \ #雙端合并 qiime vsearch join-pairs \ --i-demultiplexed-seqs primer-trimmed-demux.qza \ --p-threads 4 \ --o-joined-sequences demux-joined.qza #查看合并結果 qiime demux summarize \ --i-data demux-joined.qza \ --o-visualization demux-joined.qzv #質控過濾數據 qiime quality-filter q-score \ --i-demux demux-joined.qza \ --o-filtered-sequences demux-joined-filtered.qza \ --o-filter-stats demux-joined-filter-stats.qza qiime metadata tabulate --m-input-file demux-joined-filter-stats.qza \ --o-visualization demux-joined-filter-stats.qzv #生產feature table 自帶去除嵌合體 #deblur在denoising時需要輸入整齊一樣長度的序列，所以需要trim成相同的長度。這里需要用到前面join summary里我們記下來的join起來的質量有保障的序列大概有多長的數值。D#eblur的開發者們建議設置一個質量分數開始迅速下降的長度（recommend setting this value to a length where the median quality score begins to drop too low）。于#是本例中--p-trim-length為240。 qiime deblur denoise-16S \ --i-demultiplexed-seqs demux-joined-filtered.qza \ --p-trim-length 240 \ --p-sample-stats \ --p-jobs-to-start 2 \ --p-min-reads 1 \ --o-representative-sequences rep-seqs.qza \ --o-table feature-table.qza \ --o-stats deblur-stats.qza #結果可視化 qiime feature-table summarize \ --i-table feature-table.qza \ --o-visualization feature-table.qzv qiime deblur visualize-stats \ --i-deblur-stats deblur-stats.qza \ --o-visualization deblur-stats.qzv
### 方法4. q2 vsearch otu echo vsearch start date cd $workdir mkdir vsearch cd vsearch qiime vsearch dereplicate-sequences \ --i-sequences $workdir/deblur/demux-joined-filtered.qza \ --o-dereplicated-table table.qza \ --o-dereplicated-sequences rep-seqs.qza echo vsearch denovo start date #denovo pick otu qiime vsearch cluster-features-de-novo \ --i-table table.qza \ --i-sequences rep-seqs.qza \ --p-perc-identity 0.99 \ --o-clustered-table table-dn-99.qza \ --o-clustered-sequences rep-seqs-dn-99.qza echo vsearch close start date #close reference qiime vsearch cluster-features-closed-reference \ --i-table table.qza \ --i-sequences rep-seqs.qza \ --i-reference-sequences 85_otus.qza \ --p-perc-identity 0.85 \ --o-clustered-table table-cr-85.qza \ --o-clustered-sequences rep-seqs-cr-85.qza \ --o-unmatched-sequences unmatched-cr-85.qza echo vsearch open start date #open referenced qiime vsearch cluster-features-open-reference \ --i-table table.qza \ --i-sequences rep-seqs.qza \ --i-reference-sequences 85_otus.qza \ --p-perc-identity 0.85 \ --o-clustered-table table-or-85.qza \ --o-clustered-sequences rep-seqs-or-85.qza \ --o-new-reference-sequences new-ref-seqs-or-85.qza echo vsearch open end date #denovo 去除嵌合體 qiime vsearch uchime-denovo \ --i-table atacama-table.qza \ --i-sequences atacama-rep-seqs.qza \ --output-dir uchime-dn-out ### 特征表和代表序列統計 qiime feature-table summarize \ --i-table table_biom.qza \ --o-visualization table.qzv \ --m-sample-metadata-file $fastmap qiime feature-table tabulate-seqs \ --i-data repset-seqs.qza \ --o-visualization rep-seqs.qzv # 下載qzv在線查看，dada2只有1千多個ASV
“qiime2怎么使用”的內容就介紹到這里了，感謝大家的閱讀。如果想了解更多行業相關的知識可以關注億速云網站，小編將為大家輸出更多高質量的實用文章！
免責聲明：本站發布的內容（圖片、視頻和文字）以原創、轉載和分享為主，文章觀點不代表本網站立場，如果涉及侵權請聯系站長郵箱：is@yisu.com進行舉報，并提供相關證據，一經查實，將立刻刪除涉嫌侵權內容。
