脚本1:
perl /mnt/sdb/share/bin/find_identical.pl
数据路径:/mnt/sdb/share/bin/linux-demo/program/table
查看两个table表之间基因共有ID,各自特有ID,如果共有分别输出各自共有的结果和合并的共有结果,数据1.DEG.xls,2.DEG.xls
查找pep.faa文件和1.DEG.xls第一列共有的编号,分别输出共有,特有的结果
3.1. DEG.xls第一列基因编号加上功能描述, 没有功能的用 - 替代; 基因功能文件:all.gene.func.txt
4.找出snp.anno.xls文件内targetScanS ,esp6500si_all , GnomAD_EAS_AF, CADD表头分别位于第几列,并输出
5去除snp.anno.xls文件内Ori_REF,Ori_ALT, shared_hom, shared_het, OMIM这些列