# VEP command abruptly stops running

**URL:** <https://discuss.hail.is/t/vep-command-abruptly-stops-running/1453>\
**Category:** Hail Query & hailctl\
**Created:** [June 11, 2020, 11:31am UTC](https://discuss.hail.is/t/vep-command-abruptly-stops-running/1453 "2020-06-11T11:31:09Z")\
**Posts on this page:** 10\
**Page:** 1

<div class="post-metadata">

**Author:** ![wonu](https://avatars.discourse-cdn.com/v4/letter/w/f17d59/32.png) [@wonu](https://discuss.hail.is/u/wonu)\
**Post date:** [June 11, 2020, 11:31am UTC](https://discuss.hail.is/t/vep-command-abruptly-stops-running/1453/1 "2020-06-11T11:31:09Z")

</div>

Hi,

I am having some trouble annotating variants using VEP. I am running hail locally and using the following commands:

MT = ‘Analyses/scz\_maf\_filtered\_rare.mt’  
ht = hl.read\_matrix\_table(MT).rows()  
ht\_vep = hl.vep(ht, “vep\_configuration\_wonu.json”)  
ht\_vep.write(“annotations/vep\_annotate\_wonu.ht”, overwrite=True)

My json configuration file below is based on the example here: [https://hail.is/docs/0.2/methods/genetics.html?highlight=vep#hail.methods.vep](https://hail.is/docs/0.2/methods/genetics.html?highlight=vep#hail.methods.vep) with a different plugin and additional variables from the plugin in the json schema.

{  
“command”: [  
“/media/veracrypt10/Analyses/ensembl-vep/vep”,  
“–format”, “vcf”,  
“ **OUTPUT\_FORMAT\_FLAG** ”,  
“–everything”,  
“–allele\_number”,  
“–no\_stats”,  
“–cache”, “–offline”,  
“–port”, “3337”,  
“–minimal”,  
“–assembly”, “GRCh37”,  
“–plugin”, “dbNSFP,/media/veracrypt10/Analyses/ensembl-vep/dbNSFPv3.5a/dbNSFP\_hg19.gz,SIFT\_score,SIFT\_pred,Polyphen2\_HDIV\_score,Polyphen2\_HDIV\_pred,Polyphen2\_HVAR\_score,Polyphen2\_HVAR\_pred,LRT\_score,LRT\_pred,MutationTaster\_score,MutationTaster\_pred,MutationAssessor\_score,MutationAssessor\_pred,PROVEAN\_score,PROVEAN\_pred”,  
“–dir\_plugins”, “/home/wonu/.vep/Plugins/”,  
“–output\_file”, “STDOUT”  
],  
“env”: {  
“PERL5LIB”: “/home/wonu/.vep/Plugins/dbNSFP.pm”  
},  
“vep\_json\_schema”: “Struct{assembly\_name:String,allele\_string:String,ancestral:String,ensp:String,colocated\_variants:Array[Struct{af:Float64,afr\_af:Float64,amr\_af:Float64,eas\_af:Float64,eur\_af:Float64,sas\_af:Float64,aa\_af:Float64,ea\_af:Float64,gnomAD\_AF:Float64,gnomAD\_AFR\_AF:Float64,gnomAD\_AMR\_AF:Float64,gnomAD\_ASJ\_AF:Float64,gnomAD\_EAS\_AF:Float64,gnomAD\_FIN\_AF:Float64,gnomAD\_NFE\_AF:Float64,gnomAD\_OTH\_AF:Float64,gnomAD\_SAS\_AF:Float64,MAX\_AF:Float64,MAX\_AF\_POPS:Float64,aa\_allele:String,aa\_maf:Float64,afr\_allele:String,afr\_maf:Float64,allele\_string:String,amr\_allele:String,amr\_maf:Float64,clin\_sig:Array[String],end:Int32,eas\_allele:String,eas\_maf:Float64,ea\_allele:String,ea\_maf:Float64,eur\_allele:String,eur\_maf:Float64,exac\_adj\_allele:String,exac\_adj\_maf:Float64,exac\_allele:String,exac\_afr\_allele:String,exac\_afr\_maf:Float64,exac\_amr\_allele:String,exac\_amr\_maf:Float64,exac\_eas\_allele:String,exac\_eas\_maf:Float64,exac\_fin\_allele:String,exac\_fin\_maf:Float64,exac\_maf:Float64,exac\_nfe\_allele:String,exac\_nfe\_maf:Float64,exac\_oth\_allele:String,exac\_oth\_maf:Float64,exac\_sas\_allele:String,exac\_sas\_maf:Float64,id:String,minor\_allele:String,minor\_allele\_freq:Float64,phenotype\_or\_disease:Int32,pubmed:Array[Int32],sas\_allele:String,sas\_maf:Float64,somatic:Int32,start:Int32,strand:Int32}],context:String,end:Int32,id:String,input:String,intergenic\_consequences:Array[Struct{allele\_num:Int32,consequence\_terms:Array[String],impact:String,minimised:Int32,variant\_allele:String}],most\_severe\_consequence:String,motif\_feature\_consequences:Array[Struct{allele\_num:Int32,consequence\_terms:Array[String],high\_inf\_pos:String,impact:String,minimised:Int32,motif\_feature\_id:String,motif\_name:String,motif\_pos:Int32,motif\_score\_change:Float64,strand:Int32,variant\_allele:String}],regulatory\_feature\_consequences:Array[Struct{allele\_num:Int32,biotype:String,consequence\_terms:Array[String],impact:String,minimised:Int32,regulatory\_feature\_id:String,variant\_allele:String}],seq\_region\_name:String,start:Int32,strand:Int32,transcript\_consequences:Array[Struct{allele\_num:Int32,amino\_acids:String,biotype:String,canonical:Int32,ccds:String,cdna\_start:Int32,cdna\_end:Int32,cds\_end:Int32,cds\_start:Int32,codons:String,consequence\_terms:Array[String],distance:Int32,domains:Array[Struct{db:String,name:String}],exon:String,gene\_id:String,gene\_pheno:Int32,symbol:String,symbol\_source:String,hgnc\_id:String,hgvsc:String,hgvsp:String,hgvs\_offset:Int32,impact:String,intron:String,lof:String,lof\_flags:String,lof\_filter:String,lof\_info:String,minimised:Int32,polyphen\_prediction:String,polyphen\_score:Float64,protein\_end:Int32,protein\_start:Int32,protein\_id:String,sift\_prediction:String,sift\_score:Float64,strand:Int32,swissprot:String,transcript\_id:String,trembl:String,uniparc:String,variant\_allele:String,SIFT\_score:Float64,SIFT\_pred:String,Polyphen2\_HDIV\_score:Float64,Polyphen2\_HDIV\_pred:String,Polyphen2\_HVAR\_score:Float64,Polyphen2\_HVAR\_pred:String,LRT\_score:Float64,LRT\_pred:String,MutationTaster\_score:Float64,MutationTaster\_pred:String,MutationAssessor\_score:Float64,MutationAssessor\_pred:String,PROVEAN\_score:Float64,PROVEAN\_pred:String}],variant\_class:String}”  
}

The vep command seems to run for about 10 mins and then abruptly stops with no error message, and the resulting ht file is incomplete and unreadable as no metadata.json.gz file is produced. It also does not produce any output file. I am not sure what the problem is or how to proceed. I’m running hail version 0.2.44-6cfa355a1954

Thanks!

---

<div class="post-metadata">

**Author:** ![danking](https://yyz2.discourse-cdn.com/flex036/user_avatar/discuss.hail.is/danking/32/43_2.png) [@danking](https://discuss.hail.is/u/danking)\
**Post date:** [June 11, 2020, 2:05pm UTC](https://discuss.hail.is/t/vep-command-abruptly-stops-running/1453/2 "2020-06-11T14:05:23Z")

</div>

Hi @wonu, sorry you’re having trouble. Can you share the hail log file? That will have information necessary to diagnose the issue.

---

<div class="post-metadata">

**Author:** ![wonu](https://avatars.discourse-cdn.com/v4/letter/w/f17d59/32.png) [@wonu](https://discuss.hail.is/u/wonu)\
**Post date:** [June 11, 2020, 2:13pm UTC](https://discuss.hail.is/t/vep-command-abruptly-stops-running/1453/3 "2020-06-11T14:13:44Z")

</div>

Hi, what’s the best way to do this? I’ve tried just uploading the file but it doesn’t seem to work.

---

<div class="post-metadata">

**Author:** ![johnc1231](https://yyz2.discourse-cdn.com/flex036/user_avatar/discuss.hail.is/johnc1231/32/286_2.png) [@johnc1231](https://discuss.hail.is/u/johnc1231)\
**Post date:** [June 11, 2020, 2:24pm UTC](https://discuss.hail.is/t/vep-command-abruptly-stops-running/1453/4 "2020-06-11T14:24:50Z")

</div>

If you can’t upload you can email to [hail@broadinstitute.org](mailto:hail@broadinstitute.org)

---

<div class="post-metadata">

**Author:** ![wonu](https://avatars.discourse-cdn.com/v4/letter/w/f17d59/32.png) [@wonu](https://discuss.hail.is/u/wonu)\
**Post date:** [June 11, 2020, 3:41pm UTC](https://discuss.hail.is/t/vep-command-abruptly-stops-running/1453/5 "2020-06-11T15:41:57Z")

</div>

email sent. I think the file was too large to upload

---

<div class="post-metadata">

**Author:** ![johnc1231](https://yyz2.discourse-cdn.com/flex036/user_avatar/discuss.hail.is/johnc1231/32/286_2.png) [@johnc1231](https://discuss.hail.is/u/johnc1231)\
**Post date:** [June 11, 2020, 4:24pm UTC](https://discuss.hail.is/t/vep-command-abruptly-stops-running/1453/6 "2020-06-11T16:24:10Z")

</div>

So I’m pretty sure that what’s happening here is that your json schema isn’t consistent with what the plugin is doing. Admittedly, it’s not really easy to add new VEP plugins in our current config file model. I am pasting a few lines of the log below:

```auto
2020-06-09 16:47:18 Hail: WARN: struct{allele_num: int32, amino_acids: str, biotype: str, canonical: int32, ccds: str, cdna_start: int32, cdna_end: int32, cds_end: int32, cds_start: int32, codons: str, consequence_terms: array<str>, distance: int32, domains: array<struct{db: str, name: str}>, exon: str, gene_id: str, gene_pheno: int32, symbol: str, symbol_source: str, hgnc_id: str, hgvsc: str, hgvsp: str, hgvs_offset: int32, impact: str, intron: str, lof: str, lof_flags: str, lof_filter: str, lof_info: str, minimised: int32, polyphen_prediction: str, polyphen_score: float64, protein_end: int32, protein_start: int32, protein_id: str, sift_prediction: str, sift_score: float64, strand: int32, swissprot: str, transcript_id: str, trembl: str, uniparc: str, variant_allele: str, SIFT_score: float64, SIFT_pred: str, Polyphen2_HDIV_score: float64, Polyphen2_HDIV_pred: str, Polyphen2_HVAR_score: float64, Polyphen2_HVAR_pred: str, LRT_score: float64, LRT_pred: str, MutationTaster_score: float64, MutationTaster_pred: str, MutationAssessor_score: float64, MutationAssessor_pred: str, PROVEAN_score: float64, PROVEAN_pred: str} has no field flags at <root>.transcript_consequences[element] for value JArray(List(JString(cds_start_NF)))
2020-06-09 16:47:18 Hail: WARN: struct{allele_num: int32, amino_acids: str, biotype: str, canonical: int32, ccds: str, cdna_start: int32, cdna_end: int32, cds_end: int32, cds_start: int32, codons: str, consequence_terms: array<str>, distance: int32, domains: array<struct{db: str, name: str}>, exon: str, gene_id: str, gene_pheno: int32, symbol: str, symbol_source: str, hgnc_id: str, hgvsc: str, hgvsp: str, hgvs_offset: int32, impact: str, intron: str, lof: str, lof_flags: str, lof_filter: str, lof_info: str, minimised: int32, polyphen_prediction: str, polyphen_score: float64, protein_end: int32, protein_start: int32, protein_id: str, sift_prediction: str, sift_score: float64, strand: int32, swissprot: str, transcript_id: str, trembl: str, uniparc: str, variant_allele: str, SIFT_score: float64, SIFT_pred: str, Polyphen2_HDIV_score: float64, Polyphen2_HDIV_pred: str, Polyphen2_HVAR_score: float64, Polyphen2_HVAR_pred: str, LRT_score: float64, LRT_pred: str, MutationTaster_score: float64, MutationTaster_pred: str, MutationAssessor_score: float64, MutationAssessor_pred: str, PROVEAN_score: float64, PROVEAN_pred: str} has no field mutationassessor_pred at <root>.transcript_consequences[element] for value JString(M)
2020-06-09 16:47:18 root: WARN: struct{allele_num: int32, amino_acids: str, biotype: str, canonical: int32, ccds: str, cdna_start: int32, cdna_end: int32, cds_end: int32, cds_start: int32, codons: str, consequence_terms: array<str>, distance: int32, domains: array<struct{db: str, name: str}>, exon: str, gene_id: str, gene_pheno: int32, symbol: str, symbol_source: str, hgnc_id: str, hgvsc: str, hgvsp: str, hgvs_offset: int32, impact: str, intron: str, lof: str, lof_flags: str, lof_filter: str, lof_info: str, minimised: int32, polyphen_prediction: str, polyphen_score: float64, protein_end: int32, protein_start: int32, protein_id: str, sift_prediction: str, sift_score: float64, strand: int32, swissprot: str, transcript_id: str, trembl: str, uniparc: str, variant_allele: str, SIFT_score: float64, SIFT_pred: str, Polyphen2_HDIV_score: float64, Polyphen2_HDIV_pred: str, Polyphen2_HVAR_score: float64, Polyphen2_HVAR_pred: str, LRT_score: float64, LRT_pred: str, MutationTaster_score: float64, MutationTaster_pred: str, MutationAssessor_score: float64, MutationAssessor_pred: str, PROVEAN_score: float64, PROVEAN_pred: str} has no field mutationassessor_pred at <root>.transcript_consequences[element] for value JString(M)
2020-06-09 16:47:18 Hail: WARN: struct{allele_num: int32, amino_acids: str, biotype: str, canonical: int32, ccds: str, cdna_start: int32, cdna_end: int32, cds_end: int32, cds_start: int32, codons: str, consequence_terms: array<str>, distance: int32, domains: array<struct{db: str, name: str}>, exon: str, gene_id: str, gene_pheno: int32, symbol: str, symbol_source: str, hgnc_id: str, hgvsc: str, hgvsp: str, hgvs_offset: int32, impact: str, intron: str, lof: str, lof_flags: str, lof_filter: str, lof_info: str, minimised: int32, polyphen_prediction: str, polyphen_score: float64, protein_end: int32, protein_start: int32, protein_id: str, sift_prediction: str, sift_score: float64, strand: int32, swissprot: str, transcript_id: str, trembl: str, uniparc: str, variant_allele: str, SIFT_score: float64, SIFT_pred: str, Polyphen2_HDIV_score: float64, Polyphen2_HDIV_pred: str, Polyphen2_HVAR_score: float64, Polyphen2_HVAR_pred: str, LRT_score: float64, LRT_pred: str, MutationTaster_score: float64, MutationTaster_pred: str, MutationAssessor_score: float64, MutationAssessor_pred: str, PROVEAN_score: float64, PROVEAN_pred: str} has no field lrt_pred at <root>.transcript_consequences[element] for value JString(U)
2020-06-09 16:47:18 Hail: WARN: struct{allele_num: int32, amino_acids: str, biotype: str, canonical: int32, ccds: str, cdna_start: int32, cdna_end: int32, cds_end: int32, cds_start: int32, codons: str, consequence_terms: array<str>, distance: int32, domains: array<struct{db: str, name: str}>, exon: str, gene_id: str, gene_pheno: int32, symbol: str, symbol_source: str, hgnc_id: str, hgvsc: str, hgvsp: str, hgvs_offset: int32, impact: str, intron: str, lof: str, lof_flags: str, lof_filter: str, lof_info: str, minimised: int32, polyphen_prediction: str, polyphen_score: float64, protein_end: int32, protein_start: int32, protein_id: str, sift_prediction: str, sift_score: float64, strand: int32, swissprot: str, transcript_id: str, trembl: str, uniparc: str, variant_allele: str, SIFT_score: float64, SIFT_pred: str, Polyphen2_HDIV_score: float64, Polyphen2_HDIV_pred: str, Polyphen2_HVAR_score: float64, Polyphen2_HVAR_pred: str, LRT_score: float64, LRT_pred: str, MutationTaster_score: float64, MutationTaster_pred: str, MutationAssessor_score: float64, MutationAssessor_pred: str, PROVEAN_score: float64, PROVEAN_pred: str} has no field mutationassessor_score at <root>.transcript_consequences[element] for value JDouble(2.215)
2020-06-09 16:47:18 root: WARN: struct{allele_num: int32, amino_acids: str, biotype: str, canonical: int32, ccds: str, cdna_start: int32, cdna_end: int32, cds_end: int32, cds_start: int32, codons: str, consequence_terms: array<str>, distance: int32, domains: array<struct{db: str, name: str}>, exon: str, gene_id: str, gene_pheno: int32, symbol: str, symbol_source: str, hgnc_id: str, hgvsc: str, hgvsp: str, hgvs_offset: int32, impact: str, intron: str, lof: str, lof_flags: str, lof_filter: str, lof_info: str, minimised: int32, polyphen_prediction: str, polyphen_score: float64, protein_end: int32, protein_start: int32, protein_id: str, sift_prediction: str, sift_score: float64, strand: int32, swissprot: str, transcript_id: str, trembl: str, uniparc: str, variant_allele: str, SIFT_score: float64, SIFT_pred: str, Polyphen2_HDIV_score: float64, Polyphen2_HDIV_pred: str, Polyphen2_HVAR_score: float64, Polyphen2_HVAR_pred: str, LRT_score: float64, LRT_pred: str, MutationTaster_score: float64, MutationTaster_pred: str, MutationAssessor_score: float64, MutationAssessor_pred: str, PROVEAN_score: float64, PROVEAN_pred: str} has no field mutationassessor_score at <root>.transcript_consequences[element] for value JDouble(2.215)
2020-06-09 16:47:18 Hail: WARN: struct{allele_num: int32, amino_acids: str, biotype: str, canonical: int32, ccds: str, cdna_start: int32, cdna_end: int32, cds_end: int32, cds_start: int32, codons: str, consequence_terms: array<str>, distance: int32, domains: array<struct{db: str, name: str}>, exon: str, gene_id: str, gene_pheno: int32, symbol: str, symbol_source: str, hgnc_id: str, hgvsc: str, hgvsp: str, hgvs_offset: int32, impact: str, intron: str, lof: str, lof_flags: str, lof_filter: str, lof_info: str, minimised: int32, polyphen_prediction: str, polyphen_score: float64, protein_end: int32, protein_start: int32, protein_id: str, sift_prediction: str, sift_score: float64, strand: int32, swissprot: str, transcript_id: str, trembl: str, uniparc: str, variant_allele: str, SIFT_score: float64, SIFT_pred: str, Polyphen2_HDIV_score: float64, Polyphen2_HDIV_pred: str, Polyphen2_HVAR_score: float64, Polyphen2_HVAR_pred: str, LRT_score: float64, LRT_pred: str, MutationTaster_score: float64, MutationTaster_pred: str, MutationAssessor_score: float64, MutationAssessor_pred: str, PROVEAN_score: float64, PROVEAN_pred: str} has no field polyphen2_hdiv_pred at <root>.transcript_consequences[element] for value JString(D)
2020-06-09 16:47:18 root: WARN: struct{allele_num: int32, amino_acids: str, biotype: str, canonical: int32, ccds: str, cdna_start: int32, cdna_end: int32, cds_end: int32, cds_start: int32, codons: str, consequence_terms: array<str>, distance: int32, domains: array<struct{db: str, name: str}>, exon: str, gene_id: str, gene_pheno: int32, symbol: str, symbol_source: str, hgnc_id: str, hgvsc: str, hgvsp: str, hgvs_offset: int32, impact: str, intron: str, lof: str, lof_flags: str, lof_filter: str, lof_info: str, minimised: int32, polyphen_prediction: str, polyphen_score: float64, protein_end: int32, protein_start: int32, protein_id: str, sift_prediction: str, sift_score: float64, strand: int32, swissprot: str, transcript_id: str, trembl: str, uniparc: str, variant_allele: str, SIFT_score: float64, SIFT_pred: str, Polyphen2_HDIV_score: float64, Polyphen2_HDIV_pred: str, Polyphen2_HVAR_score: float64, Polyphen2_HVAR_pred: str, LRT_score: float64, LRT_pred: str, MutationTaster_score: float64, MutationTaster_pred: str, MutationAssessor_score: float64, MutationAssessor_pred: str, PROVEAN_score: float64, PROVEAN_pred: str} has no field polyphen2_hdiv_pred at <root>.transcript_consequences[element] for value JString(D)
2020-06-09 16:47:18 Hail: WARN: struct{allele_num: int32, amino_acids: str, biotype: str, canonical: int32, ccds: str, cdna_start: int32, cdna_end: int32, cds_end: int32, cds_start: int32, codons: str, consequence_terms: array<str>, distance: int32, domains: array<struct{db: str, name: str}>, exon: str, gene_id: str, gene_pheno: int32, symbol: str, symbol_source: str, hgnc_id: str, hgvsc: str, hgvsp: str, hgvs_offset: int32, impact: str, intron: str, lof: str, lof_flags: str, lof_filter: str, lof_info: str, minimised: int32, polyphen_prediction: str, polyphen_score: float64, protein_end: int32, protein_start: int32, protein_id: str, sift_prediction: str, sift_score: float64, strand: int32, swissprot: str, transcript_id: str, trembl: str, uniparc: str, variant_allele: str, SIFT_score: float64, SIFT_pred: str, Polyphen2_HDIV_score: float64, Polyphen2_HDIV_pred: str, Polyphen2_HVAR_score: float64, Polyphen2_HVAR_pred: str, LRT_score: float64, LRT_pred: str, MutationTaster_score: float64, MutationTaster_pred: str, MutationAssessor_score: float64, MutationAssessor_pred: str, PROVEAN_score: float64, PROVEAN_pred: str} has no field mutationtaster_score at <root>.transcript_consequences[element] for value JDouble(0.999949)

```

if you look at the end of the lines, it’s complaining that, for example, the json schema has no field called `polyphen2_hvar_pred`. That’s because you called it " `Polyphen2_HVAR_pred`. I think if you keep changing names in your schema to squash these warnings in the logs, that’ll fix your problem.

---

<div class="post-metadata">

**Author:** ![wonu](https://avatars.discourse-cdn.com/v4/letter/w/f17d59/32.png) [@wonu](https://discuss.hail.is/u/wonu)\
**Post date:** [June 11, 2020, 5:04pm UTC](https://discuss.hail.is/t/vep-command-abruptly-stops-running/1453/7 "2020-06-11T17:04:03Z")

</div>

Thank you! I’ll try that.

---

<div class="post-metadata">

**Author:** ![wonu](https://avatars.discourse-cdn.com/v4/letter/w/f17d59/32.png) [@wonu](https://discuss.hail.is/u/wonu)\
**Post date:** [June 25, 2020, 1:29pm UTC](https://discuss.hail.is/t/vep-command-abruptly-stops-running/1453/8 "2020-06-25T13:29:55Z")

</div>

Hi again,

I’ve been able to edit the more straight forward warnings but one of them is as follows:

2020-06-24 15:26:46 Hail: WARN: Can’t convert JSON value JObject(List((C,JObject(List((gnomad\_fin,JInt(0)), (gnomad\_nfe,JDouble(1.453E-4)), (gnomad\_asj,JDouble(2.167E-4)), (gnomad\_oth,JInt(0)), (gnomad,JDouble(1.24E-4)), (gnomad\_eas,JInt(0)), (gnomad\_sas,JDouble(3.702E-4)), (gnomad\_afr,JInt(0)), (gnomad\_amr,JInt(0))))))) to type array\<struct{gnomad: float64, gnomad\_afr: float64, gnomad\_amr: float64, gnomad\_asj: float64, gnomad\_eas: float64, gnomad\_fin: float64, gnomad\_nfe: float64, gnomad\_oth: float64, gnomad\_sas: float64}\> at .colocated\_variants[element].frequencies.

I’m not familiar with JSON formatting, so I was wondering if you could assist me with the correct structure for this?

The complete configuration is below:  
{  
“command”: [  
“/media/veracrypt10/Analyses/ensembl-vep/vep”,  
“–format”, “vcf”,  
“ **OUTPUT\_FORMAT\_FLAG** ”,  
“–everything”,  
“–allele\_number”,  
“–no\_stats”,  
“–cache”, “–offline”,  
“–port”, “3337”,  
“–minimal”,  
“–assembly”, “GRCh37”,  
“–plugin”, “dbNSFP,/media/veracrypt10/Analyses/ensembl-vep/dbNSFPv3.5a/dbNSFP\_hg19.gz,SIFT\_score,SIFT\_pred,Polyphen2\_HDIV\_score,Polyphen2\_HDIV\_pred,Polyphen2\_HVAR\_score,Polyphen2\_HVAR\_pred,LRT\_score,LRT\_pred,MutationTaster\_score,MutationTaster\_pred,MutationAssessor\_score,MutationAssessor\_pred,PROVEAN\_score,PROVEAN\_pred”,  
“–dir\_plugins”, “/home/wonu/.vep/Plugins/”,  
“–output\_file”, “STDOUT”  
],  
“env”: {  
“PERL5LIB”: “/home/wonu/.vep/Plugins/dbNSFP.pm”  
},  
“vep\_json\_schema”: “Struct{assembly\_name:String,allele\_string:String,ancestral:String,ensp:String,minimised:Int32,colocated\_variants:Array[Struct{af:Float64,afr\_af:Float64,amr\_af:Float64,eas\_af:Float64,eur\_af:Float64,sas\_af:Float64,aa\_af:Float64,ea\_af:Float64,gnomAD\_AF:Float64,gnomAD\_AFR\_AF:Float64,gnomAD\_AMR\_AF:Float64,gnomAD\_ASJ\_AF:Float64,gnomAD\_EAS\_AF:Float64,gnomAD\_FIN\_AF:Float64,gnomAD\_NFE\_AF:Float64,gnomAD\_OTH\_AF:Float64,gnomAD\_SAS\_AF:Float64,MAX\_AF:Float64,MAX\_AF\_POPS:Float64,aa\_allele:String,aa\_maf:Float64,afr\_allele:String,afr\_maf:Float64,allele\_string:String,amr\_allele:String,amr\_maf:Float64,frequencies:Array[Struct{gnomad:Float64,gnomad\_afr:Float64,gnomad\_amr:Float64,gnomad\_asj:Float64,gnomad\_eas:Float64,gnomad\_fin:Float64,gnomad\_nfe:Float64,gnomad\_oth:Float64,gnomad\_sas:Float64}],clin\_sig:Array[String],end:Int32,eas\_allele:String,eas\_maf:Float64,ea\_allele:String,ea\_maf:Float64,eur\_allele:String,eur\_maf:Float64,exac\_adj\_allele:String,exac\_adj\_maf:Float64,exac\_allele:String,exac\_afr\_allele:String,exac\_afr\_maf:Float64,exac\_amr\_allele:String,exac\_amr\_maf:Float64,exac\_eas\_allele:String,exac\_eas\_maf:Float64,exac\_fin\_allele:String,exac\_fin\_maf:Float64,exac\_maf:Float64,exac\_nfe\_allele:String,exac\_nfe\_maf:Float64,exac\_oth\_allele:String,exac\_oth\_maf:Float64,exac\_sas\_allele:String,exac\_sas\_maf:Float64,id:String,minor\_allele:String,minor\_allele\_freq:Float64,phenotype\_or\_disease:Int32,pubmed:Array[Int32],sas\_allele:String,sas\_maf:Float64,somatic:Int32,start:Int32,strand:Int32,seq\_region\_name:String}],context:String,end:Int32,id:String,input:String,intergenic\_consequences:Array[Struct{allele\_num:Int32,consequence\_terms:Array[String],impact:String,minimised:Int32,variant\_allele:String}],most\_severe\_consequence:String,motif\_feature\_consequences:Array[Struct{allele\_num:Int32,consequence\_terms:Array[String],high\_inf\_pos:String,impact:String,minimised:Int32,motif\_feature\_id:String,motif\_name:String,motif\_pos:Int32,motif\_score\_change:Float64,strand:Int32,variant\_allele:String}],regulatory\_feature\_consequences:Array[Struct{allele\_num:Int32,biotype:String,consequence\_terms:Array[String],impact:String,minimised:Int32,regulatory\_feature\_id:String,variant\_allele:String}],seq\_region\_name:String,start:Int32,strand:Int32,transcript\_consequences:Array[Struct{allele\_num:Int32,amino\_acids:String,biotype:String,canonical:Int32,ccds:String,cdna\_start:Int32,cdna\_end:Int32,cds\_end:Int32,cds\_start:Int32,codons:String,consequence\_terms:Array[String],distance:Int32,domains:Array[Struct{db:String,name:String}],exon:String,gene\_id:String,gene\_pheno:Int32,gene\_symbol:String,gene\_symbol\_source:String,hgnc\_id:String,hgvsc:String,hgvsp:String,hgvs\_offset:Int32,impact:String,intron:String,lof:String,lof\_flags:String,lof\_filter:String,lof\_info:String,minimised:Int32,polyphen\_prediction:String,polyphen\_score:Float64,protein\_end:Int32,protein\_start:Int32,protein\_id:String,sift\_prediction:String,strand:Int32,swissprot:String,transcript\_id:String,trembl:String,uniparc:String,variant\_allele:String,sift\_score:Float64,sift\_pred:String,polyphen2\_hdiv\_score:Float64,polyphen2\_hdiv\_pred:String,polyphen2\_hvar\_score:Float64,polyphen2\_hvar\_pred:String,lrt\_score:Float64,lrt\_pred:String,mutationtaster\_score:Float64,mutationtaster\_pred:String,mutationassessor\_score:Float64,mutationassessor\_pred:String,provean\_score:Float64,provean\_pred:String}],variant\_class:String}”  
}

---

<div class="post-metadata">

**Author:** ![johnc1231](https://yyz2.discourse-cdn.com/flex036/user_avatar/discuss.hail.is/johnc1231/32/286_2.png) [@johnc1231](https://discuss.hail.is/u/johnc1231)\
**Post date:** [June 25, 2020, 1:51pm UTC](https://discuss.hail.is/t/vep-command-abruptly-stops-running/1453/9 "2020-06-25T13:51:54Z")

</div>

Sorry this is so hard to use. If I’m reading this correctly, the JSON that’s coming out of VEP looks like:

```auto
JObject(
    List(
        (C, JObject(
                List(
                    (gnomad_fin,JInt(0)), 
                    (gnomad_nfe,JDouble(1.453E-4)), 
                    (gnomad_asj,JDouble(2.167E-4)), 
                    (gnomad_oth,JInt(0)), 
                    (gnomad,JDouble(1.24E-4)), 
                    (gnomad_eas,JInt(0)), 
                    (gnomad_sas,JDouble(3.702E-4)), 
                    (gnomad_afr,JInt(0)), 
                    (gnomad_amr,JInt(0))
                )
            )
        )
    )
) 

```

So to translate, `JObject` is equivalent to `Struct`. So this suggests that the type of `colocated_variants.frequencies` is

```auto
Struct{C:Struct{gnomad:Float64,gnomad_afr:Float64,gnomad_amr:Float64,gnomad_asj:Float64,gnomad_eas:Float64,gnomad_fin:Float64,gnomad_nfe:Float64,gnomad_oth:Float64,gnomad_sas:Float64}}

```

You currently have that type as an `Array[Struct{.......}]`, and I think it should be a `Struct{C: Struct{...}}`

EDIT: Got it wrong the first time, changed it around.

---

<div class="post-metadata">

**Author:** ![wonu](https://avatars.discourse-cdn.com/v4/letter/w/f17d59/32.png) [@wonu](https://discuss.hail.is/u/wonu)\
**Post date:** [June 25, 2020, 4:32pm UTC](https://discuss.hail.is/t/vep-command-abruptly-stops-running/1453/10 "2020-06-25T16:32:00Z")

</div>

No problem! Your suggestion seems to have worked. Thank you!
