-
Notifications
You must be signed in to change notification settings - Fork 0
/
annovar_db_download.sh
executable file
·94 lines (72 loc) · 4.64 KB
/
annovar_db_download.sh
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
# Date: 201908
# Annovar database download
# Reference genome version: GRCh37/hg19 | GRCh38/hg38
# URL: http://annovar.openbioinformatics.org/en/latest/user-guide/download/
# Hang-on database path: /home/workshop/ngs_tools/annovar/humandb/
#---------------------------
# Creat humandb folder
if [ -d ~/humandb ]; then
echo 'humandb folder exist!!'
else mkdir -p ~/humandb;
echo 'Create humandb folder done!'
fi
export DB_PATH=$HOME/humandb/
##### Gene-based #######
# NCBI RefSeq v.201706
annotate_variation.pl -buildver hg19 -downdb -webfrom annovar refGene $DB_PATH
annotate_variation.pl -buildver hg38 -downdb -webfrom annovar refGene $DB_PATH
# UCSC KnownGene v.201706
annotate_variation.pl -buildver hg19 -downdb -webfrom annovar knownGene $DB_PATH
annotate_variation.pl -buildver hg38 -downdb -webfrom annovar knownGene $DB_PATH
# Ensembl Gene v.201706
annotate_variation.pl -buildver hg19 -downdb -webfrom annovar ensGene $DB_PATH
annotate_variation.pl -buildver hg38 -downdb -webfrom annovar ensGene $DB_PATH # based on Gencode v26 Basic collection v.20170912
##### Region-based #####
# UCSC cytoBand
annotate_variation.pl -buildver hg19 -downdb cytoBand $DB_PATH
# UCSC cytoBand hg38
cd $DB_PATH && wget -c http://hgdownload.cse.ucsc.edu/goldenpath/hg38/database/cytoBand.txt.gz && gunzip -d cytoBand.txt.gz && mv cytoBand.txt hg38_cytoBand.txt
# UCSC genomicSuperDups
annotate_variation.pl -buildver hg19 -downdb genomicSuperDups $DB_PATH
# UCSC genomicSuperDups hg38
cd $DB_PATH && wget -c http://hgdownload.cse.ucsc.edu/goldenpath/hg38/database/genomicSuperDups.txt.gz && gunzip -d genomicSuperDups.txt.gz && mv genomicSuperDups.txt hg38_genomicSuperDups.txt
# UCSC phastConsElements100way
annotate_variation.pl -buildver hg19 -downdb phastConsElements100way $DB_PATH
cd $DB_PATH && wget -c http://hgdownload.cse.ucsc.edu/goldenpath/hg38/database/phastConsElements100way.txt.gz && gunzip -d phastConsElements100way.txt.gz && mv phastConsElements100way.txt hg38_phastConsElements100way.txt
##### Filter-based #####
# 1000 Genome v.20150824
annotate_variation.pl -buildver hg19 -downdb -webfrom annovar 1000g2015aug $DB_PATH
annotate_variation.pl -buildver hg38 -downdb -webfrom annovar 1000g2015aug $DB_PATH
# NCBI dbSNP150 v.20170929
annotate_variation.pl -buildver hg19 -downdb -webfrom annovar avsnp150 $DB_PATH
annotate_variation.pl -buildver hg38 -downdb -webfrom annovar avsnp150 $DB_PATH
# NCI60 v.20150428
annotate_variation.pl -buildver hg19 -downdb -webfrom annovar nci60 $DB_PATH
annotate_variation.pl -buildver hg38 -downdb -webfrom annovar nci60 $DB_PATH
# NCBI ClinVar v.20190305
annotate_variation.pl -buildver hg19 -downdb -webfrom annovar clinvar_20190305 $DB_PATH
annotate_variation.pl -buildver hg38 -downdb -webfrom annovar clinvar_20190305 $DB_PATH
# gnomAD genome collection
annotate_variation.pl -buildver hg19 -downdb -webfrom annovar gnomad211_genome $DB_PATH
annotate_variation.pl -buildver hg38 -downdb -webfrom annovar gnomad211_genome $DB_PATH
# gnomAD exome collection
annotate_variation.pl -buildver hg19 -downdb -webfrom annovar gnomad211_exome $DB_PATH
annotate_variation.pl -buildver hg38 -downdb -webfrom annovar gnomad211_exome $DB_PATH
# ExAC 65000 exome allele frequency data for ALL, AFR (African), AMR (Admixed American), EAS (East Asian), FIN (Finnish), NFE (Non-finnish European), OTH (other), SAS (South Asian)). version 0.3. Left normalization done.
annotate_variation.pl -buildver hg19 -downdb -webfrom annovar exac03 $DB_PATH
annotate_variation.pl -buildver hg38 -downdb -webfrom annovar exac03 $DB_PATH
# ESP6500 v.20141222
annotate_variation.pl -buildver hg19 -downdb -webfrom annovar esp6500siv2_all $DB_PATH
annotate_variation.pl -buildver hg38 -downdb -webfrom annovar esp6500siv2_all $DB_PATH
# InterVar: clinical interpretation of missense variants (indels not supported) 20180325
annotate_variation.pl -buildver hg19 -downdb -webfrom annovar intervar_20180118 $DB_PATH
annotate_variation.pl -buildver hg38 -downdb -webfrom annovar intervar_20180118 $DB_PATH
# whole-exome SIFT, PolyPhen2 HDIV, PolyPhen2 HVAR, LRT, MutationTaster, MutationAssessor, FATHMM, MetaSVM, MetaLR, VEST, CADD, GERP++, DANN, fitCons, PhyloP and SiPhy scores from dbNSFP (protein domain for variants)
annotate_variation.pl -buildver hg19 -downdb -webfrom annovar dbnsfp31a_interpro $DB_PATH
annotate_variation.pl -buildver hg38 -downdb -webfrom annovar dbnsfp31a_interpro $DB_PATH
# gwasCatalog
annotate_variation.pl -buildver hg19 -downdb gwasCatalog $DB_PATH
annotate_variation.pl -buildver hg38 -downdb gwasCatalog $DB_PATH
# LoF tools
wget -c http://www.openbioinformatics.org/annovar/download/LoFtool_scores.txt.gz -P $DB_PATH
gzip -d $DB_PATH/LoFtool_scores.txt.gz