2020-02-11 eggnog-mapper的使用及结果

序自《2020-02-08 wget(含参数技巧)下载eggNOG数据库、安装以及碰到的问题》

一、使用

1、数据库内容

在下载好eggnog数据库之后,

1)og2level.tsv.gz;;2)eggnog.db.gz;;3)OG_fasta.tar.gz;;4)eggnog_proteins.dmnd.gz;;

5)HMM database(s): bact(存在了HMM_dblevels里)

包括:index.html;;bact_50.hmm.h3f;;bact_50.hmm.h3i;;bact_50.hmm.h3m;;bact_50.hmm.h3p;;bact_50.hmm.idmap;;bact_50.info;;bact_50.pkl

2、安装好后测试

在软件安装目录下

emapper.py -h        #查看使用帮助

python emapper.py -i test/polb.fa --output polb_bact -d bact --data_dir ~/data/eggnog.db        #按照别人的操作,我的报的报错了,后来发现是不能到eggnog,db下,因为他不是一个目录,这是一个文件(也看是跟我下载的我有关系)

python emapper.py -i test/polb.fa --output polb_bact -d bact --data_dir ~/data        #这样就可以了


运行的比对过程
这是给出的引用信息

出现这些就证明可以正常使用了!

二、使用

注意:

python emapper.py -i /share/data1//lx_genome/genepredict/SZDP190917115/115.pep --resume --output 115_bact -d bact --data_dir /share/data1//software/eggnog-mapper-1.0.3/data

参数的选择以及解释:

参考:https://www.jianshu.com/p/e646c0fa6443

eggnog-mapper默认是以HMMER进行序列搜索,尽管这可以通过-m diamond更改成DIAMOND,但结果中就会缺少一些信息列。

-d/--database表示搜索数据库的类型,既可以是"euk, bact, arch" 三大类的其中一种,也可以是"eggnog-mapper/data/OG_fasta"中的小类,例如"maNOG"表示的就是哺乳动物的NOG。此外也可以是自定义HMM数据库,-d /path/to/pfam.hmm

--usemem和--cpu 线程数提高运行速度(如果服务器内存比较大,线程比较多的话)

--output表示输出文件的前缀,默认输出在当前文件夹下,

--output_dir可以更改为其他文件路径;

--resume表示任务重启后可以跳过之前已经完成的部分, 

--override则表示覆盖原先的输出结果。

其他几个基因的分析命令同上

至于如何同时使用这个在后台同时运行多条,还需要学习,以后会补上。

三、结果

注意:生成的文件默认存在当前的目录下


参数参考:

1、使用eggNOG数据库进行基因功能注释 http://www.chenlianfu.com/?p=2804

2、序列功能注释神器:eggNOG-mapper,KEGG/COG/KOG/GO/BiGG 一网打尽http://www.biostack.org/?p=698

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 199,478评论 5 467
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 83,825评论 2 376
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 146,482评论 0 330
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 53,726评论 1 271
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 62,633评论 5 359
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,018评论 1 275
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,513评论 3 390
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,168评论 0 254
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,320评论 1 294
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,264评论 2 317
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,288评论 1 328
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 32,995评论 3 315
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,587评论 3 303
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,667评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 30,909评论 1 255
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,284评论 2 345
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 41,862评论 2 339

推荐阅读更多精彩内容