###将GATK结果转换成annovar输入文件#####
perl convert2annovar.pl-format vcf4 Hom-972-rawGATK-indel-filter.vcf> Hom-972-indel.avinput
###mm10基因组注释#####
perl annotate_variation.pl-out Het-160-indel -build mm10 Het-160-indel.avinput mousedb
###将GATK结果转换成annovar输入文件#####
perl convert2annovar.pl-format vcf4 Hom-972-rawGATK-indel-filter.vcf> Hom-972-indel.avinput
###mm10基因组注释#####
perl annotate_variation.pl-out Het-160-indel -build mm10 Het-160-indel.avinput mousedb