Variant Filtration error "Number format exception"
Can you please provide
a) GATK version used : 4.1.5.0
b) Exact GATK commands used
$GATK VariantFiltration -R $REF -V ${TRIO}.snps.hardfilt.fordenov.hg19_multianno.vcf -O ${TRIO}.snps.hardfilt.fordenov.varfilt.vcf --filter-name "rare_ExAC" --filter-expression 'vc.getAttribute("ExAC_ALL") >0.0005' --filter-name "rare_1000G" --filter-expression 'vc.getAttribute("1000g2015aug_all") > 0.0005' --filter-name "rare_esp" --filter-expression 'vc.getAttribute("esp6500siv2_all") > 0.0005'
c) The entire error log if applicable.
Using GATK jar /usr/bin/gatk-4.1.5.0/gatk-package-4.1.5.0-local.jar
Running:
java -Dsamjdk.use_async_io_read_samtools=false -Dsamjdk.use_async_io_write_samtools=true -Dsamjdk.use_async_io_write_tribble=false -Dsamjdk.compression_level=2 -jar /usr/bin/gatk-4.1.5.0/gatk-package-4.1.5.0-local.jar VariantFiltration -R /media/sheetalshetty/NewDrive/Pipeline_New//Data//Inputs//Reference_Genome/hg19.p13.plusMT.no_alt_analysis_set/hg19.p13.plusMT.full_analysis_set.fa -V /media/sheetalshetty/NewDrive/Pipeline_New//Data//Outputs//M_F750-003-A_Trio/M_F750-003-A.snps.hardfilt.fordenov.hg19_multianno.vcf -O /media/sheetalshetty/NewDrive/Pipeline_New//Data//Outputs//M_F750-003-A_Trio/M_F750-003-A.snps.hardfilt.fordenov.varfilt.vcf --filter-name rare_ExAC --filter-expression vc.getAttribute("ExAC_ALL") >0.0005 --filter-name rare_1000G --filter-expression vc.getAttribute("1000g2015aug_all") > 0.0005 --filter-name rare_esp --filter-expression vc.getAttribute("esp6500siv2_all") > 0.0005
16:29:14.956 INFO NativeLibraryLoader - Loading libgkl_compression.so from jar:file:/usr/bin/gatk-4.1.5.0/gatk-package-4.1.5.0-local.jar!/com/intel/gkl/native/libgkl_compression.so
Aug 04, 2020 4:29:15 PM shaded.cloud_nio.com.google.auth.oauth2.ComputeEngineCredentials runningOnComputeEngine
INFO: Failed to detect whether we are running on Google Compute Engine.
16:29:15.109 INFO VariantFiltration - ------------------------------------------------------------
16:29:15.109 INFO VariantFiltration - The Genome Analysis Toolkit (GATK) v4.1.5.0
16:29:15.109 INFO VariantFiltration - For support and documentation go to https://software.broadinstitute.org/gatk/
16:29:15.109 INFO VariantFiltration - Executing as sheetalshetty@KruerLab-Neurology on Linux v4.15.0-106-generic amd64
16:29:15.109 INFO VariantFiltration - Java runtime: OpenJDK 64-Bit Server VM v11.0.8+10-post-Ubuntu-0ubuntu118.04.1
16:29:15.110 INFO VariantFiltration - Start Date/Time: August 4, 2020 at 4:29:14 PM MST
16:29:15.110 INFO VariantFiltration - ------------------------------------------------------------
16:29:15.110 INFO VariantFiltration - ------------------------------------------------------------
16:29:15.111 INFO VariantFiltration - HTSJDK Version: 2.21.2
16:29:15.111 INFO VariantFiltration - Picard Version: 2.21.9
16:29:15.111 INFO VariantFiltration - HTSJDK Defaults.COMPRESSION_LEVEL : 2
16:29:15.111 INFO VariantFiltration - HTSJDK Defaults.USE_ASYNC_IO_READ_FOR_SAMTOOLS : false
16:29:15.111 INFO VariantFiltration - HTSJDK Defaults.USE_ASYNC_IO_WRITE_FOR_SAMTOOLS : true
16:29:15.111 INFO VariantFiltration - HTSJDK Defaults.USE_ASYNC_IO_WRITE_FOR_TRIBBLE : false
16:29:15.111 INFO VariantFiltration - Deflater: IntelDeflater
16:29:15.111 INFO VariantFiltration - Inflater: IntelInflater
16:29:15.111 INFO VariantFiltration - GCS max retries/reopens: 20
16:29:15.111 INFO VariantFiltration - Requester pays: disabled
16:29:15.111 INFO VariantFiltration - Initializing engine
16:29:15.411 INFO FeatureManager - Using codec VCFCodec to read file file:///media/sheetalshetty/NewDrive/Pipeline_New/Data/Outputs/M_F750-003-A_Trio/M_F750-003-A.snps.hardfilt.fordenov.hg19_multianno.vcf
16:29:15.461 INFO VariantFiltration - Done initializing engine
16:29:15.579 INFO ProgressMeter - Starting traversal
16:29:15.579 INFO ProgressMeter - Current Locus Elapsed Minutes Variants Processed Variants/Minute
16:29:15.610 INFO VariantFiltration - Shutting down engine
[August 4, 2020 at 4:29:15 PM MST] org.broadinstitute.hellbender.tools.walkers.filters.VariantFiltration done. Elapsed time: 0.01 minutes.
Runtime.totalMemory()=2155872256
java.lang.NumberFormatException: For input string: "."
at java.base/jdk.internal.math.FloatingDecimal.readJavaFormatString(FloatingDecimal.java:2054)
at java.base/jdk.internal.math.FloatingDecimal.parseDouble(FloatingDecimal.java:110)
at java.base/java.lang.Double.parseDouble(Double.java:543)
at org.apache.commons.jexl2.JexlArithmetic.toDouble(JexlArithmetic.java:1016)
at org.apache.commons.jexl2.JexlArithmetic.compare(JexlArithmetic.java:699)
at org.apache.commons.jexl2.JexlArithmetic.greaterThan(JexlArithmetic.java:790)
at org.apache.commons.jexl2.Interpreter.visit(Interpreter.java:796)
at org.apache.commons.jexl2.parser.ASTGTNode.jjtAccept(ASTGTNode.java:18)
at org.apache.commons.jexl2.Interpreter.interpret(Interpreter.java:232)
at org.apache.commons.jexl2.ExpressionImpl.evaluate(ExpressionImpl.java:65)
at htsjdk.variant.variantcontext.JEXLMap.evaluateExpression(JEXLMap.java:186)
at htsjdk.variant.variantcontext.JEXLMap.get(JEXLMap.java:95)
at htsjdk.variant.variantcontext.JEXLMap.get(JEXLMap.java:15)
at htsjdk.variant.variantcontext.VariantContextUtils.match(VariantContextUtils.java:338)
at org.broadinstitute.hellbender.tools.walkers.filters.VariantFiltration.matchesFilter(VariantFiltration.java:427)
at org.broadinstitute.hellbender.tools.walkers.filters.VariantFiltration.filter(VariantFiltration.java:381)
at org.broadinstitute.hellbender.tools.walkers.filters.VariantFiltration.apply(VariantFiltration.java:340)
at org.broadinstitute.hellbender.engine.VariantWalker.lambda$traverse$0(VariantWalker.java:104)
at java.base/java.util.stream.ForEachOps$ForEachOp$OfRef.accept(ForEachOps.java:183)
at java.base/java.util.stream.ReferencePipeline$3$1.accept(ReferencePipeline.java:195)
at java.base/java.util.stream.ReferencePipeline$2$1.accept(ReferencePipeline.java:177)
at java.base/java.util.stream.ReferencePipeline$3$1.accept(ReferencePipeline.java:195)
at java.base/java.util.Iterator.forEachRemaining(Iterator.java:133)
at java.base/java.util.Spliterators$IteratorSpliterator.forEachRemaining(Spliterators.java:1801)
at java.base/java.util.stream.AbstractPipeline.copyInto(AbstractPipeline.java:484)
at java.base/java.util.stream.AbstractPipeline.wrapAndCopyInto(AbstractPipeline.java:474)
at java.base/java.util.stream.ForEachOps$ForEachOp.evaluateSequential(ForEachOps.java:150)
at java.base/java.util.stream.ForEachOps$ForEachOp$OfRef.evaluateSequential(ForEachOps.java:173)
at java.base/java.util.stream.AbstractPipeline.evaluate(AbstractPipeline.java:234)
at java.base/java.util.stream.ReferencePipeline.forEach(ReferencePipeline.java:497)
at org.broadinstitute.hellbender.engine.VariantWalker.traverse(VariantWalker.java:102)
at org.broadinstitute.hellbender.engine.GATKTool.doWork(GATKTool.java:1048)
at org.broadinstitute.hellbender.cmdline.CommandLineProgram.runTool(CommandLineProgram.java:139)
at org.broadinstitute.hellbender.cmdline.CommandLineProgram.instanceMainPostParseArgs(CommandLineProgram.java:191)
at org.broadinstitute.hellbender.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:210)
at org.broadinstitute.hellbender.Main.runCommandLineProgram(Main.java:163)
at org.broadinstitute.hellbender.Main.mainEntry(Main.java:206)
at org.broadinstitute.hellbender.Main.main(Main.java:292)
-
Hi Sheetal Shetty, it looks like the tool is expecting a number but there is a period. (java.lang.NumberFormatException: For input string: "."). To troubleshoot this issue you can run ValidateVariants or look through your filters to make sure they contain the correct values to be evaluated.
Please sign in to leave a comment.
1 comment