网创优客建站品牌官网
为成都网站建设公司企业提供高品质网站建设
热线:028-86922220
成都专业网站建设公司

定制建站费用3500元

符合中小企业对网站设计、功能常规化式的企业展示型网站建设

成都品牌网站建设

品牌网站建设费用6000元

本套餐主要针对企业品牌型网站、中高端设计、前端互动体验...

成都商城网站建设

商城网站建设费用8000元

商城网站建设因基本功能的需求不同费用上面也有很大的差别...

成都微信网站建设

手机微信网站建站3000元

手机微信网站开发、微信官网、微信商城网站...

建站知识

当前位置:首页 > 建站知识

如何正确的使用Annovar

如何正确的使用Annovar,很多新手对此不是很清楚,为了帮助大家解决这个难题,下面小编将为大家详细讲解,有这方面需求的人可以来学习下,希望你能有所收获。

成都创新互联是一家朝气蓬勃的网站建设公司。公司专注于为企业提供信息化建设解决方案。从事网站开发,网站制作,网站设计,网站模板,微信公众号开发,软件开发,成都小程序开发,10年建站对茶楼设计等多个领域,拥有丰富的网站设计经验。

  

annovar的下载及安装

Annovar是用perl语言写的,可以在任何安装了perl的系统上运行,且不需要进行安装,直接下载解压就可以使用。但它的下载需要注册,且需要使用教育机构或者科研单位后缀的邮箱。当然,如果你没有注册邮箱也没有关系,后台回复annovar即可得到软件安装包。Annovar主要有三种不同形式的注释方式:

1、Gene-based annotation:根据SNP或者CNV的位置来判断是否会引起蛋白质编码的变化,是否发生了氨基酸的改变。

2、 Region-based annotation: 来鉴定特定基因组区域的突变。

3、Filter-based annotation:用来鉴定特定数据库中的突变。  
 

下载完annovar并且解压之后,主要包括以下文件:

如何正确的使用Annovar

  • example:存放的是示例文件

  • humandb:部分注释数据库的文件,annovar的软件中自带了一部分,根据自己的研究需要也可以自己下载

  • annotate_variation.pl:主程序,用来进行数据库的下载,以及不同形式的注释

  • coding_change.pl:用来推断蛋白质的序列是否发生变化

  • convert2annovar.pl:将其他多种形式转化为annovar可识别的形式(如将vcf文件转化为annovar可识别形式)

  • retrieve_seq_from_fasta.pl:自行建立其它物种的转录本

  • table_annovar.pl:可以一次完成三种不同形式的注释

  • variants_reduction.pl:用来定制过滤注释流程

—  输入文件 —

Annovar的输入文件是一个简单的文本格式文件,其中前五列应分别是染色体号、突变位点在染色体上的起始位置、突变位点的结束位置、该突变位点在参考序列上的碱基以及该位点的突变碱基,其他列的内容可以有也可以没有。

如何正确的使用Annovar

如果输入文件是vcf文件,可以采用annovar的convert2annovar.pl程序将vcf文件转化为annovar可识别的文件形式,具体的命令如下:

   perl convert2annovar.pl -format vcf4 G-001.vcf -outfile G.avinput

输出文件的格式为:

如何正确的使用Annovar

— 数据库下载  —

Annovar的注释主要依赖于数据库,因此在进行分析之前,应将所需的数据库下载到humandb文件夹中,下载的命令如下:

perl annotate_variation.pl -buildver hg19 -downdb -webfrom annovar avsnp147 humandb/

-buildver:对应参考基因组的版本

-downdb –webfrom annovar:从annovar库中下载对应的数据库,如果不知道要下载什么数据库,可以在annovar库中查看对应的数据库以及对应的功能,网址为:(http://annovar.openbioinformatics.org/en/latest/user-guide/download/)

avsnp147:下载的数据库的名称

humandb:下载到humandb文件夹中

— 结果注释—

整理好输入文件格式以及下载好数据库后,就可以进行注释了,下面以table_annovar.pl为例介绍下annovar的注释功能,具体命令如下:

perl table_annovar.pl GCK.avinput annovar/humandb/ -buildver hg19 -out GCK -remove –protocol refGene,1000g2015aug_eas,1000g2015aug_eur,1000g2015aug_sas,1000g2015aug_amr  -operation g,f,f,f ,f -nastring .

table_annovar.pl:输入文件

-buildver:参考序列版本

-out:输出文件

-remove:删掉程d序运行过程中产生的中间文件

–protocol:数据库的名称

-operation:对应顺序的数据库的类型,如千人基因组,dbsnp数据库等(g代表gene-based、r代表region-based、f代表filter-based),与前面数据库一一对应

-nastring .:缺省值用.表示

看完上述内容是否对您有帮助呢?如果还想对相关知识有进一步的了解或阅读更多相关文章,请关注创新互联行业资讯频道,感谢您对创新互联的支持。


新闻标题:如何正确的使用Annovar
文章分享:http://bjjierui.cn/article/ghgeii.html

其他资讯