======简介====== meta分析流程中,会遇到客户需求按分组绘制CARD双圈图。 此处整理补充了meta CARD分组双圈图分析思路。 ======功能====== CARD分组双圈图分析。 ======数据准备====== **1. unigene丰度文件准备:** 双圈图的作图文件为Unigenes.absolute.total.tax.xls和Unigenes.protein.rgi.txt,整体性思路是按照分组将2个表格处理成分组表格, 首先处理Unigenes.absolute.total.tax.xls将同一个分组的样本列提取出来,然后使用脚本处理该文件的行,脚本如下: perl /TJPROJ5/META_ASS/16s/chenjiawei/script/META/CARD/fenzu-circos/shaixuan.pl Unigenes.absolute.total.tax.xls 此处挑选出分组unigene id形成一个list,命名为shaixuan.list。 **2. rgi预测文件准备:** 根据挑选出的unigene id list,即shaixuan.list对rgi文件进行筛选,选出该分组的rgi预测结果。 perl shaixuan2.pl Unigenes.protein.rgi.txt shaixuan.list Unigenes.protein.rgi-new.txt ======数据分析====== **1. 生成分组双圈图各层级绘图文件:** perl /PUBLIC/software/MICRO/share/MetaGenome_pipeline/MetaGenome_pipeline_V5.1/lib/05.Function/lib/CARD/lib/mechanism.tax.pl --mech /PUBLIC/software/MICRO/share/MetaGenome_pipeline/MetaGenome_pipeline_V5.1/lib/05.Function/lib/CARD/lib/../database/ARO_Mechnism.list --anno Unigenes.protein.rgi.txt --table Unigenes.absolute.total.tax.xls --outdir Mech_tax --prefix meta.tax **2. 绘制双圈图:** perl /PUBLIC/software/MICRO/share/MetaGenome_pipeline/MetaGenome_pipeline_V5.1/lib/05.Function/lib/CARD/lib/circos_overview4.0.pl --table Mech_tax/meta.tax.p.xls 完整脚本路径: /TJPROJ5/META_ASS/16s/chenjiawei/script/META/CARD/fenzu-circos