可以使用qiime2的命令行工具qiime metadata edit来进行批量重命名ASV号。 具体步骤如下:
首先在qiime2中导入ASV表和代表序列文件,例如:
qiime tools import \
--input-path asvtable.qza \
--type 'FeatureTable[Frequency]' \
--input-format QIIME2 \
--output-path asvtable.qza qiime tools import \
--input-path repseqs.qza \
--type 'FeatureData[Sequence]' \
--input-format QIIME2 \
--output-path repseqs.qza
使用qiime metadata edit命令修改ASV号,例如:
qiime metadata edit \
--type FeatureTable[Frequency] \
--input-path asvtable.qza \
--output-path asvtable_renamed.qza \
--mapping 'ASV1,ASV2\n1,2\n2,3\n3,4\n4,5'
这里的mapping就是用来指定新旧ASV号的对应关系,格式为CSV文件。例如,上面的mapping表示将ASV1改为ASV2,将ASV2改为ASV3,以此类推。
使用qiime metadata edit命令修改代表序列中的ID,例如:
qiime metadata edit \
--type FeatureData[Sequence] \
--input-path repseqs.qza \
--output-path repseqs_renamed.qza \
--mapping 'ASV1,ASV2\nASV2,ASV3\nASV3,ASV4\nASV4,ASV5'
这里的mapping也是用来指定新旧ID的对应关系,格式同上。
最后再使用修改后的ASV表和代表序列进行后续分析,例如:
qiime feature-table summarize \
--i-table asvtablerenamed.qza \
--o-visualization asvtablesummary.qzv qiime feature-classifier classify-sklearn \
--i-classifier classifier.qza \
--i-reads repseqsrenamed.qza \
--o-classification taxonomy.qza
这样就可以方便快捷地批量重命名qiime2做出的ASV号了。注意,修改后的文件不能再与原始文件混淆,建议另存为新的文件名。
这家伙很懒,还没有设置简介