snpEff安装及建库

snpEff安装及使用:

1.安装(推荐这种方法而不是conda)

参考: http://pcingola.github.io/SnpEff/download/

# Download latest version
wget https://snpeff.blob.core.windows.net/versions/snpEff_latest_core.zip

# Unzip file
unzip snpEff_latest_core.zip

2. 配置自己要注释的基因组及注释数据库

参考:http://pcingola.github.io/SnpEff/se_faq/(可以通过 GTF, GFF, RefSeq or GenBank files建库,官方推荐GTF文件,方法最简单)

下载方法:
从下面网址找自己要的基因组:https://ftp.ncbi.nlm.nih.gov/(因为这个下面的文件夹总是变来变去,所以从这个总网站进去自己翻最方便)找refseq再找到物种的fna.gz文件以及gtf.gz文件

我要下载的是肺炎支原体(https://ftp.ncbi.nlm.nih.gov/genomes/refseq/bacteria/Mycoplasma_pneumoniae/latest_assembly_versions/GCF_900660465.1_50648_A01-3/):

  • 建库
    文件夹结构:
    snpEff
    ├── SnpSift.jar
    ├── data
    │ ├── mp2019
    │ │ ├── genes.gff
    │ │ └── snpEffectPredictor.bin #建库成功后会显示
    │ ├── genomes
    │ │ └── mp2019.fa
    ├── examples/
    ├── galaxy/
    ├── scripts/
    ├── snpEff.config
    └── snpEff.jar
# 1.解压改名
mkdir data
cd data
mkdir genomes  #fa文件都放这个文件夹
gzip -d GCF_900660465.1_50648_A01-3_genomic.fna.gz
mv GCF_900660465.1_50648_A01-3_genomic.fna mp2019.fa
cd ..

mkdir mp2019
cd mp2019
gzip -d GCF_900660465.1_50648_A01-3_genomic.gtf.gz
mv GCF_900660465.1_50648_A01-3_genomic.gtf  genes.gtf #必须要这样改!

# 2. 将genome添加到配置文件中
cd  snpEff
vi snpEff.config 
#(shift+g到最后一行进行编辑)
#将以下部分写进去
#---------------------------------------------------------------
# my own databases
# Mycolasma_pneumoniae genome,version mp2019
mp2019.genome : Mycolasma_pneumoniae
# mp2019.chromosomes: NZ_LR214945.1 这一步是为了指定密码子版本(常见的物种比如人不需要加这步)
mp2019.NZ_LR214945.1.codonTable : Mycoplasma

# 3. 建库
java -jar ./snpEff.jar build -gtf22 -v mp2019

不同物种密码子版本请参考
https://www.ncbi.nlm.nih.gov/Taxonomy/taxonomyhome.html/index.cgi?chapter=tgencodes#SG11

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

推荐阅读更多精彩内容