【直播】我的基因组 32:使用annovar注释vcf

ANNOVAR软件用法看我以前的博客: 1.Annovar使用记录 (http://www.bio-info-trainee.com/641.html) 2.用annovar对snp进行注释 (http://www.bio-info-trainee.com/441.html) 3.对感兴趣的基因call variation(http://www.bio-info-trainee.com/2013.html)4.WES(六)用annovar注释(http://www.bio-info-trainee.com/1158.html)

ANNOVAR是一个perl编写的命令行工具,能在安装了perl解释器的多种操作系统上执行。允许多种输入文件格式,包括最常被使用的VCF格式。输出文件也有多种格式,包括注释过的VCF文件、用tab或者逗号分隔的txt文件。ANNOVAR能快速注释遗传变异并预测其功能。类似的variants注释软件还有VEP,snpEff,VAAST,AnnTools等等。

这是annovar的第一讲,照例我们先简单试用一下它的基础功能,因为这样的遗传变异位点的注释软件,功能实在是太多了,但是它最基础的东西是根据变异位点的坐标已经现有研究已知基因结构坐标的信息来将变异进行分类,类别如下:

其中位于外显子上面的突变是尤为重要的,所以会重点分析:

两个步骤即可:

~/biosoft/ANNOVAR/annovar/convert2annovar.pl -format vcf4old realign.vcf >realign.annovar ~/biosoft/ANNOVAR/annovar/annotate_variation.pl -buildver hg19 --geneanno --outfile realign.anno realign.annovar ~/biosoft/ANNOVAR/annovar/humandb/

annovar默认的注释方式是:

Done with 50914 transcripts (including 11516 without coding sequence annotation) for 26271 unique genes

这个是最基础的注释,几分钟就完成了480万个位点的注释。

其中有31741个variation位点在外显子上面,是指的我们重点关注的,其中外显子的variation又可以分成下面几种情况:

一般公司会把上面的数据画成一个好看的饼图,但这不是我们的重点,分析自己的基因组数据是为了探究本质的东西,而不是浮于表面。

文:Jimmy、吃瓜群众

图文编辑:吃瓜群众

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

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

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

发表于

我来说两句

0 条评论
登录 后参与评论

相关文章

来自专栏PHP技术

PEAR简介:用PEAR来写你的下一个php程序

你可能已经是个PHP的老手了,写了很多非常棒的代码。但是,如果你现在要把它们加入到你现在的项目中去,是否有些吃力?你的朋友想使用你的代码作为他的项目中的一个模块...

2607
来自专栏企鹅号快讯

Python的新手指南,教你如何变编程大佬级别的

新的Python?下面是基础知识的细分,包括语言的历史,使用者的语言以及Python 2与3的区别。 Python是一种编程语言,对于编写快速而简单的脚本非常有...

1868
来自专栏瓜大三哥

ROM 单元

在数字系统中,由于ROM掉电后数据不会丢失,因此ROM单元也有着更广泛的应用。 对于容量不大的ROM,在Verilog HDL中可以通过case语句来实现。下面...

1825
来自专栏阿凯的Excel

Pandas实现分列功能(Pandas读书笔记1)

最近发现自己特愿意扯和分享内容本身无关的事情,下述内容纯闲扯: pandas应该怎么分享困扰了我好久,但是看我公众号的朋友更困惑的是这个人为啥要分享panda...

3574
来自专栏浪淘沙

java学习day1--JDK

2018.6.4 今天是实训的第一天,从今天起,我将每天在这里记录每天的学习知识。 1、Java是跨平台的语言。 因为字节码是在虚拟机上运行的,而不...

712
来自专栏FreeBuf

初探Android逆向:通过游戏APP破解引发的安全思考

如今移动互联网已经完全融入到我们的生活中,各类APP也是层出不穷,因此对于安卓APP安全的研究也尤为重要。本文通过对一款安卓APP的破解实例,来引出对于APP安...

643
来自专栏AI科技评论

开发 | 如何在 i5 上实现 20 倍的 Python 运行速度?

Intel Distribution for Python 在今年二月进行了更新——英特尔发布了 Update 2 版本。以“加速”为核心的它,相比原生 Pyt...

3816
来自专栏小L的魔法馆

第13届景驰-埃森哲杯广东工业大学ACM程序设计大赛--K-密码

2646
来自专栏技术博客

VS2013中Python学习笔记[环境搭建]

Python的设计具有很强的可读性,相比其他语言经常使用英文关键字,其他语言的一些标点符号,它具有比其他语言更有特色语法结构。

842
来自专栏ml

catalan---卡特兰数(小结)

(关于卡特兰数的详细介绍)http://baike.baidu.com/view/2499752.htm 下面有练习的题目: 经过测试,_int64/long...

3427

扫码关注云+社区