annovar注释SNP

1.软件安装

1.1使用conda安装gtfToGenePred软件

给conda添加下载channels

conda config --add channels defaults
conda config --add channels bioconda
conda config --add channels conda-forge

下载gtfToGenePred

conda install ucsc-gtftogenepred
conda update ucsc-gtftogenepred

usage

gtfToGenePred -genePredExt Gbarbadense_gene_model.gtf  Gbarbadense_gene_model.refGene.txt

annovar 注释

  1. 将参考基因组文件转换格式

  2. --format指定要转换的格式

  3. --seqfile后面接参考基因组序列文件

  4. --outfile输出文件名

Gbarbadense_gene_model.Pred文件为gtgtfToGenePred软件生成的文件

module load annovar
retrieve_seq_from_fasta.pl --format refGene --seqfile Gbarbadense_genome_HAU_v2.0.fasta Gbarbadense_gene_model.refGene.txt --outfile Gbarbadense_refGEneMrna.fa

  1. 将vcf文件转换为annovar格式

    6G的vcf文件大概跑了

convert2annovar.pl --includeinfo --allsample  --withfreq --format vcf4 ./../Gbarbadense_genome.snp.filter.recode.vcf >Gbarbadence.avinput

  1. table_annovar.pl进行注释

    gtf转换后的文件和基因序列转换后的文件都要放在Gbarbadense/目录下

    • --protocol指定数据库类型

    • --operation注释类型 g、r、f分别只按照基因、region、filter进行注释，对应的数据库--protocol参数也有指明

    • --thread线程数

    • --maxgenethread当线程数超过6时，需要声明，不然最多就是6个线程在跑

    • --outfile输出文件前缀

    • Gbarbadense/文件夹中包含Gbarbadense_gene_model.refGene.txt文件

   table_annovar.pl --maxgenethread 10  --thread 10  Gbarbadense.avinput  Gbarbadense/ -buildver Gbarbadense --outfile Gbarbadense_annovar --protocol refGene,refGene,refGene --operation g,r,f

只对基因区域进行SNP的注释

   table_annovar.pl --maxgenethread 10  --thread 10  Gbarbadense.avinput  Gbarbadense/ -buildver Gbarbadense --outfile Gbarbadense_annovar --protocol refGene --operation g

  1. 最终生成文件

    由于--protocol参数我用的都是refGene数据库类型，所以region、fileter模式的注释应该都有问题；没放出来

    ├── Gbarbadense_annovar.refGene.exonic_variant_function
├── Gbarbadense_annovar.refGene.invalid_input
├── Gbarbadense_annovar.refGene.log
├── Gbarbadense_annovar.refGene.variant_function

参考

  1. 非人类https://blog.csdn.net/u013816205/article/details/51262289

  2. 非人类https://blog.csdn.net/g863402758/article/details/75304391

  3. https://anjingwd.github.io/AnJingwd.github.io/2018/01/20/ANNOVAR%E8%BF%9B%E8%A1%8C%E7%AA%81%E5%8F%98%E6%B3%A8%E9%87%8A/

  4. gtfTogGenePred安装 https://bioconda.github.io/recipes/ucsc-gtftogenepred/README.html

PreviousbedtoolsNextHiC-Pro安装

Last updated