#来源于https://mp.weixin.qq.com/s/j901X2DYQhyFXGRNEFaoaQ
#eggnog-mapper软件来对基因做功能注释,由于需要注释的基因有10万条,在线版单次只支持5000
#条序列因此使用本地版来进行注释
conda create -n eggnog-mapper -c bioconda eggnog-mapper
#下载数据库解压缩
axel -an50 http://eggnog5.embl.de/download/emapperdb-5.0.2/eggnog.db.gz
axel -an50 http://eggnog5.embl.de/download/emapperdb-5.0.2/eggnog_proteins.dmnd.gz
#使用
python /home/abc/miniconda3/bin/emapper.py -i /home/abc/Desktop/pep.fa
--output /home/abc/Desktop/At -m diamond --cpu 35
--seed_ortholog_evalue 1e-5
--dmnd_db /data/public/database/eggnog/eggnog_proteins.dmnd
--data_dir /data/public/database/eggnog/eggnog.db
#注意需要写对软件的路径及数据库相对应得路径,本次使用35个线程来进行注释还是相当快的,此##外该程序对内存得占用不是很高
暂无评论