The current GATK version is 3.7-0
Examples: Monday, today, last week, Mar 26, 3/26/04

#### Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

You can opt in to receive email notifications, for example when your questions get answered or when there are new announcements, by following the instructions given here.

#### ☞ Did you remember to?

1. Search using the upper-right search box, e.g. using the error message.
3. Include tool and Java versions.
4. Tell us whether you are following GATK Best Practices.
5. Include relevant details, e.g. platform, DNA- or RNA-Seq, WES (+capture kit) or WGS (PCR-free or PCR+), paired- or single-end, read length, expected average coverage, somatic data, etc.
6. For tool errors, include the error stacktrace as well as the exact command.
7. For format issues, include the result of running ValidateSamFile for BAMs or ValidateVariants for VCFs.
8. For weird results, include an illustrative example, e.g. attach IGV screenshots according to Article#5484.
9. For a seeming variant that is uncalled, include results of following Article#1235.

#### ☞ Formatting tip!

Wrap blocks of code, error messages and BAM/VCF snippets--especially content with hashes (#)--with lines with three backticks ( ` ) each to make a code block as demonstrated here.

GATK 3.7 is here! Be sure to read the Version Highlights and optionally the full Release Notes.

# JexlExpression with VariantEval

Member Posts: 5

Hello,

I'm currently running variantEval to count up variants per individual stratified by a variety of annotations.

My GATK call looks like:

java -jar /humgen/gsa-hpprojects/GATK/bin/current/GenomeAnalysisTK.jar \
-T VariantEval \
-R Homo_sapiens_assembly19.fasta \
-o output.txt \
-L input.vcf \
-eval input.vcf \
-ST Sample -noST \
-noEV -EV CountVariants \
-ST JexlExpression --select_names "nonsynon" --select_exps "resource.VAT_CDS == 'nonsynonymous' && resource.FOUNDERS_FRQ > 0.05" \
-ST JexlExpression --select_names "synon" --select_exps "resource.VAT_CDS == 'synonymous' && resource.FOUNDERS_FRQ > 0.05" ...

where the VAT_CDS section of the INFO field in the VCF has a functional annotation or is set to "na" if an annotation is unavailable. I'm getting the following error:

##### ERROR ------------------------------------------------------------------------------------------

but weirdly the error is not consistent (my data is stratified by chromosome and most chromosomes will run without throwing the error while one or two chromosomes do exhibit the error. Do you have any ideas what's causing this behavior?

Thanks!

Tagged:

That is a little weird. Can you confirm that your expression is working on a small subsample of your file? If so, can you narrow down the records that are failing?

Geraldine Van der Auwera, PhD

• Member Posts: 5

Hi Geraldine,

I can confirm that the expression was working on a subset of the data which would seem to indicate that it was indeed a specific record causing the error.

I tried narrowing down which variant in particular was causing the problem so I reran the script to do this and, lo and behold, did not get an error this time around. Sorry I can't be more helpful here (though I'm glad everything is error free at the moment).

The only difference between the run that caused the error and the one that did not was that I updated the INFO field in the VCF using the variantAnnotator the first time around but the second time wrote my own script to do it. Not sure why that would have changed things.

If I stumble across the error again, I'll try digging into it a bit more.

Sara