码迷,mamicode.com
首页 > 其他好文 > 详细

如何使用SnpEff 对SNP结果进行分析

时间:2014-09-16 23:26:11      阅读:1462      评论:0      收藏:0      [点我收藏+]

标签:http   io   os   使用   java   ar   for   文件   数据   

SnpEff is a variant annotation and effect prediction tool. It annotates and predicts the effects of variants on genes 

详细的说明请阅读:

http://snpeff.sourceforge.net/SnpEff_manual.html

 

一, 安装:

首先在家目录下, 下载安装包

 wget http://sourceforge.net/projects/snpeff/files/snpEff_latest_core.zip

然后进行解压

unzip snpEff_latest_core.zip

会产生一个snpEff目录 所有的程序都在这里面

二, 配置自己的基因组和注释文件, 官方的数据库中有大量的参考基因组,一般都不需要配置。如果在官方的database中没有找到就需要自己配置

假如我现在有一个参考基因组: Osativa_204.fa

有个这个基因组的注释文件: Osativa_204_gene.gff3

首先编辑配置文件,加入新基因组的entry

配置文件在snpEff目录下, 配置文件名为snpEff.config

用vi进行编辑 加入如下两行

# Rice genome, version Osativa_204
Osativa_204.genome : Rice

然后保存退出

还是在snpEff文件下, 创建目录data

mkdir data

cd data

创建Osativa_204 和 genomes目录

mkdir Osativa_204

mkdir genoems

将你的gff3注释文件放在Osativa_204目录下

将你的参考序列文件放在genomes目录下

注意,要将注释文件重新命名为genes.gff

完成后回到 snpEff 目录, 执行命令:

java -jar snpEff.jar build -gff3 -v Osativa_204

 

执行:

先将vcf文件copy到data目录下

然后在snpEff目录下执行命令:

java -Xmx8g -jar snpEff.jar Osativa_204 data/testgroup.filtered.ordered.vcf > test.eff.vcf

命令执行完后在snpEff目录下会产生三个文件

snpEff_genes.txt

snpEff_summary.html

test.eff.vcf

然后将snpEff_summary.html用浏览器打开就可以看到结果的汇总情况

 

 

 

 

by freemao

FAFU.

free_mao@qq.com

 

如何使用SnpEff 对SNP结果进行分析

标签:http   io   os   使用   java   ar   for   文件   数据   

原文地址:http://www.cnblogs.com/freemao/p/3975927.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!