Kevin2015的个人博客分享 http://blog.sciencenet.cn/u/Kevin2015

博文

全外显子组生信分析流程-13-table_annovar

已有 3036 次阅读 2019-3-25 21:43 |个人分类:全外显子项目|系统分类:科研笔记

#purpose: re-annotate VCF with table_annovar.pl 
perl annotate_variation.pl -buildver hg19 -downdb -webfrom annovar refGene humandb/
perl annotate_variation.pl -buildver hg19 -downdb cytoBand humandb/
perl annotate_variation.pl -buildver hg19 -downdb -webfrom annovar exac03 humandb/ 
perl annotate_variation.pl -buildver hg19 -downdb -webfrom annovar avsnp147 humandb/ 
perl annotate_variation.pl -buildver hg19 -downdb -webfrom annovar dbnsfp30a humandb/


#set env
annovar_adr="/home/zhanghl/supporting_softwares/annovar/annovar"
infolder="/path/mutect2_result"
outfolder="/path/annovar_result"
annovardb=$annovar_adr/humandb

for sample in `cat /path/sample_list`;do 
/home/zhanghl/local/perl_install/bin/perl $annovar_adr/table_annovar.pl  $infolder/${sample}.mutect2.filter.vcf \
$annovardb/ \
-buildver hg19 \
-out $outfolder/$sample \
-remove \
-protocol refGene,cytoBand,exac03,avsnp147,dbnsfp30a \
-operation gx,r,f,f,f \
-nastring . \
-csvout \
-polish \
-xref /home/zhanghl/supporting_softwares/annovar/annovar/example/gene_fullxref.txt
done




https://blog.sciencenet.cn/blog-2609994-1169595.html

上一篇:全外显子组生信分析流程-12-Control-Freec
下一篇:全外显子组生信分析流程-14-annovar_to_maftools
收藏 IP: 159.226.149.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-2-29 10:55

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部