arithmetic exception
Posted in Ask the team | Last updated on 2012-10-18 01:16:20


Comments (2)

Dear GATK team,

I'm trying to call variants on some 'haploid' human data (Illumina reads from a mix of clones). I did exactly the following:

crd8% cd /wga/dev/jaffe
crd8% mkdir GATK2
crd8% cd GATK2
crd8% (got GenomeAnalysisTK-2.1-12.tar.bz2 from the GATK site on the web)
crd8% bunzip2 GenomeAnalysisTK-2.1-12.tar.bz2
crd8% cat GenomeAnalysisTK-2.1-12.tar | tar xf -

crd8% cd /local/scratch/jaffe/BroadCRD/fixed

crd8% java -jar /wga/dev/jaffe/GATK2/GenomeAnalysisTK-2.1-12-ga99c19d/GenomeAnalysisTK.jar -R /wga/scr2/bigrefs/human19/genome.fasta -T UnifiedGenotyper -I frag.list -o raw2.vcf -U -baq CALCULATE_AS_NECESSARY -nt 48

crd8% java -jar /wga/dev/jaffe/GATK2/GenomeAnalysisTK-2.1-12-ga99c19d/GenomeAnalysisTK.jar -R /wga/scr2/bigrefs/human19/genome.fasta -T VariantFiltration -U -V raw2.vcf -o var.vcf --filterExpression "QD<5.0||AC<2||DP<6" --filterName junk

The last command failed, with output

... 
INFO  10:38:07,429 TraversalEngine -        Location processed.sites  runtime per.1M.sites completed total.runtime remaining 
WARN  10:38:08,453 Interpreter - ![12,18]: 'QD < 5.0 || AC < 2 || DP < 6;' < error 
java.lang.ArithmeticException: Long coercion: java.util.ArrayList:([1, 1])
    at org.apache.commons.jexl2.JexlArithmetic.toLong(JexlArithmetic.java:914)
    at org.apache.commons.jexl2.JexlArithmetic.compare(JexlArithmetic.java:718)
    at org.apache.commons.jexl2.JexlArithmetic.lessThan(JexlArithmetic.java:774)
    at org.apache.commons.jexl2.Interpreter.visit(Interpreter.java:967)
    at org.apache.commons.jexl2.parser.ASTLTNode.jjtAccept(ASTLTNode.java:18)
...

Can you please suggest a solution? You should be able to access my data, if that would help. Well actually you would have to login to crd8, but you could copy files from there.

Thank you very much.

David


Return to top Comment on this article in the forum