【直播】我的基因组 31:vcf文件标记dbSNP的rsID号

vcf文件标记dbSNP的rsID号的这个问题非常多的人问过,大部分的variation calling软件给出的vcf文件里面第3列都是一个纯粹的dot占位符,如下:

但是,我们往往想知道,这个染色体的这个坐标发现的这个变异,是否在dbSNP数据库里面出现呢?这里,我们可以用snpEFF软件套装里面的SnpSift工具,具体安装教程见前面第5讲。http://snpeff.sourceforge.net/SnpSift.html#annotate

首先,需要下载dbSNP数据库的文件,并且要保证基因组版本是一致的!

java -jar ~/biosoft/SnpEff/snpEff/SnpSift.jar annotate ~/annotation/variation/human/dbSNP/All_20160601.vcf.gz realign.vcf >realign.dbsnp.vcf

当然,仅仅是注释一个dbSNP 的ID意义不大,我们得理解为什么这样注释,没有基础知识背景的小伙伴还需要理解上面的dbSNP的含义。我们还可以分析480万variation里面哪一些是可以注释到dbSNP的,哪一些不能,它们这两种在测序深度或者杂合纯合上面有显著的区别吗?有染色体分布的偏好性吗?

我们都会在下一讲慢慢说明~~~

文:Jimmy、阿尔的太阳

图文编辑:吃瓜群众

原文发布于微信公众号 - 生信技能树(biotrainee)

原文发表时间:2016-12-31

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

发表于

我来说两句

0 条评论
登录 后参与评论

相关文章

来自专栏陈仁松博客

ASP.NET Core 'Microsoft.Win32.Registry' 错误修复

今天在发布Asp.net Core应用到Azure的时候出现错误InvalidOperationException: Cannot find compilati...

4878
来自专栏pangguoming

Spring Boot集成JasperReports生成PDF文档

由于工作需要,要实现后端根据模板动态填充数据生成PDF文档,通过技术选型,使用Ireport5.6来设计模板,结合JasperReports5.6工具库来调用渲...

1.3K7
来自专栏杨龙飞前端

scrollto 到指定位置

2564
来自专栏我和未来有约会

Silverlight第三方控件专题

这里我收集整理了目前网上silverlight第三方控件的专题,若果有所遗漏请告知我一下。 名称 简介 截图 telerik 商 RadC...

4105
来自专栏C#

DotNet加密方式解析--非对称加密

    新年新气象,也希望新年可以挣大钱。不管今年年底会不会跟去年一样,满怀抱负却又壮志未酬。(不过没事,我已为各位卜上一卦,卦象显示各位都能挣钱...)...

5008
来自专栏张善友的专栏

Silverlight + Model-View-ViewModel (MVVM)

     早在2005年,John Gossman写了一篇关于Model-View-ViewModel模式的博文,这种模式被他所在的微软的项目组用来创建Expr...

3038
来自专栏一个爱瞎折腾的程序猿

sqlserver使用存储过程跟踪SQL

USE [master] GO /****** Object: StoredProcedure [dbo].[sp_perfworkload_trace_s...

2210
来自专栏闻道于事

js登录滑动验证,不滑动无法登陆

js的判断这里是根据滑块的位置进行判断,应该是用一个flag判断 <%@ page language="java" contentType="text/html...

7288
来自专栏Golang语言社区

【Golang语言社区】GO1.9 map并发安全测试

var m sync.Map //全局 func maintest() { // 第一个 YongHuomap := make(map[st...

4918
来自专栏落花落雨不落叶

canvas画简单电路图

67011

扫码关注云+社区