|
|
@@ -1,9 +1,9 @@ |
|
|
|
task annovar { |
|
|
|
|
|
|
|
File vcf |
|
|
|
String annotated_vcf = basename(vcf,".vcf") |
|
|
|
String basename = basename(vcf,".vcf") |
|
|
|
String hg |
|
|
|
File database |
|
|
|
File annovar_database |
|
|
|
String docker |
|
|
|
String cluster_config |
|
|
|
String disk_size |
|
|
@@ -13,7 +13,14 @@ task annovar { |
|
|
|
set -o pipefail |
|
|
|
set -e |
|
|
|
nt=$(nproc) |
|
|
|
/installations/annovar/table_annovar.pl ${vcf} ${database} -buildver ${hg} -out ${annotated_vcf} -remove -protocol refGene,cytoBand,genomicSuperDups,ljb26_all,dbnsfp35c,intervar_20180118,cosmic70,exac03,gnomad211_exome,clinvar_20200316 -operation g,r,r,f,f,f,f,f,f,f -nastring . -vcfinput -thread $nt |
|
|
|
|
|
|
|
if [ ${hg} == "hg38" ]; then |
|
|
|
/installations/annovar/table_annovar.pl ${vcf} ${annovar_database} -buildver ${hg} -out ${basename} -remove -protocol refGene,cytoBand,genomicSuperDups,ljb26_all,dbnsfp35c,intervar_20180118,cosmic70,exac03,gnomad211_exome,clinvar_20200316 -operation g,r,r,f,f,f,f,f,f,f -nastring . -vcfinput -thread $nt |
|
|
|
fi |
|
|
|
|
|
|
|
if [ ${hg} == "hg19" ]; then |
|
|
|
/installations/annovar/table_annovar.pl ${vcf} ${annovar_database} -buildver ${hg} -out ${basename} -remove -protocol refGene,cytoBand,genomicSuperDups,ljb26_all,snp138,cosmic78,intervar_20170202,popfreq_all_20150413,clinvar_20190305 -operation g,r,r,f,f,f,f,f,f -nastring . -vcfinput -thread $nt |
|
|
|
fi |
|
|
|
>>> |
|
|
|
|
|
|
|
runtime { |
|
|
@@ -24,8 +31,8 @@ task annovar { |
|
|
|
} |
|
|
|
|
|
|
|
output { |
|
|
|
File avinput = "${annotated_vcf}.avinput" |
|
|
|
File multianno_txt = "${annotated_vcf}.${hg}_multianno.txt" |
|
|
|
File multianno = "${annotated_vcf}.${hg}_multianno.vcf" |
|
|
|
File avinput = "${basename}.avinput" |
|
|
|
File multianno_txt = "${basename}.${hg}_multianno.txt" |
|
|
|
File multianno_vcf = "${basename}.${hg}_multianno.vcf" |
|
|
|
} |
|
|
|
} |