在QIIME2中,DADA2的聚类过程和质量控制过程是紧密结合的,无法直接分离。如果你希望仅进行聚类过程,可以通过将质量控制参数设置为宽松的方式来实现。

以下是一个示例代码,其中我们将trunc_len设置为0(表示不截断序列),max_ee设置为999(表示接受所有序列),trunc_q设置为0(表示不剪切序列):

qiime dada2 denoise-paired \
  --i-demultiplexed-seqs demux-paired-end.qza \
  --p-trim-left-f 0 \
  --p-trim-left-r 0 \
  --p-trunc-len-f 0 \
  --p-trunc-len-r 0 \
  --p-max-ee-f 999 \
  --p-max-ee-r 999 \
  --p-trunc-q 0 \
  --o-table table.qza \
  --o-representative-sequences rep-seqs.qza \
  --o-denoising-stats denoising-stats.qza

请注意,这种方式可能会导致较多的错误序列被保留,因此需要根据实际情况进行评估和分析。


原文地址: https://www.cveoy.top/t/topic/bR3C 著作权归作者所有。请勿转载和采集!

免费AI点我,无需注册和登录