基于snpeff进行snp注释

#总体思路:基于conda下载snpeff,会配置好java环境问题,但是conda不好修改snpEff.config,在此环境配置好的基础上,单独下载snpeff,开始运行
#snpeff下载
wget -c https://snpeff.blob.core.windows.net/versions/snpEff_latest_core.zip
#解压
unzip snpEff_latest_core.zip
#vim 编辑器修改snpEff.config文件1.修改data路径 2.加入目标物种
vim snpEff.config
#下两行粘贴到snpEff.config最后
#Arabidopsis,version Arabidopsis
Arabidopsis.genome : Arabidopsis
#SnpEff下面新建data文件夹,并在里面新建你添加的物种的文件夹
#data下分别新建文件夹Arabidopsis(放入注释文命名为genes.gff)和genomes(放入Arabidopsis.fa)
#添加执行权限
chmod +x snpEff.jar
chmod +x SnpSift.jar
#建库
java -jar snpEff.jar build -gff3 -c snpEff.config -v Arabidopsis -d -noCheckCds -noCheckProtein
#注释
java -Xmx8g -jar ../snpEff.jar -o gatk -v Arabidopsis ./merged_snp.vcf > merged_snp.ann.vcf

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值