start_time: 2024-03-29 09:25:59 +0800

SnpEff使用方法

96
LiuYueRR
IP属地: 云南
1.0 2019.10.16 09:49 字数 304

SnpEff使用方法

SnpEff 软件通过基因组结构注释数据(GTF文件),对VCF文件中的SNP/InDel信息进行注释,即主要解释了SNP/InDel是否能够对编码蛋白基因造成影响。

最近在给课题组做一些突变体基因定位的工作(BSA混池测序),得到了最终的VCF文件,然后最终将得到的SNP/InDel注释出来。

1.SnpEff的下载

普通下载(推荐)

#软件下载
wget https://downloads.sourceforge.net/project/snpeff/snpEff_latest_core.zip
#解压缩
unzip snpEff_latest_core.zip -d ~

conda下载

conda activate mutmap
conda install -y snpeff
#conda下载的话,需要自己寻找snpeff软件包的位置

2.SnpEff使用

SnpEff软件的主要程序就是snpEff.jar,该软件需要Java运行程序。SnpEff使用最多的程序就是build和eff,build适用于数据库的构建,eff适用于对SNP/InDel进行注释。

2.1构建SnpEff数据库

SnpEff软件的运行,首先需要基因组fasta序列信息和GTF注释信息,来构建数据库。

配置文件步骤如下:

1.在~/snpEff/目录中,创建一个文件夹:data
2.在~/snpEFF/data目录下,创建两个文件夹
    AT_10/   genomes/
    这两个文件夹中,分别放置了GTF文件和基因组文件
    genes.gtf sequences.fa
3.编辑~/snpEff/snpEff.config文件
    在文件的最后一行添加信息:
    AT_10.genome: AT

构建数据库步骤如下:

java -jar ~/snpEff/snpEff.jar build -c ~/snpEff/snpEff.config -gtf22 -v AT_10

#参数说明
java -jar: Java环境下运行程序
-c snpEff.config配置文件路径
-gtf22 设置输入的基因组注释信息是gtf2.2格式
-gff3 设置输入基因组注释信息是gff3格式
-v 设置在程序运行过程中输出的日志信息
最后的AT_10参数 设置输入的基因组版本信息,和~/snpEff/snpEff.config配置文件中添加的信息一致
2.2使用SnpEff进行注释
java -Xmx10G -jar ~/snpEff/snpEff.jar eff -c ~/snpEff/snpEff.config AT_10 positive.vcf > positive.snp.eff.vcf -csvStats positive.csv -stats positive.html

最终会产生四个文件 positive.snp.eff.vcf positive.html positive.csv positive.genes.txt

可以在positive.snp.eff.vcf文件中,分析自己的后续基因位点了。

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