Hi GATK Users,

Happy Thanksgiving!
Our staff will be observing the holiday and will be unavailable from 22nd to 25th November. This will cause a delay in reaching out to you and answering your questions immediately. Rest assured we will get back to it on Monday November 26th. We are grateful for your support and patience.
Have a great holiday everyone!!!

Regards
GATK Staff

Why is gatk4 GetPileupSummaries always failed ?

Binfls6321DTUBinfls6321DTU Member
edited October 2017 in Ask the GATK team

I used the following command

/home/projects/cu_10098/apps/src/gatk-4.beta.5/gatk-launch --javaOptions "-Xmx20g" GetPileupSummaries \
-R /home/databases/gatk-legacy-bundles/b37/human_g1k_v37_decoy.fasta \
-I tumor.bam \
-V /home/projects/cu_10098/data/public_data/gnomad.exomes.r2.0.2.sites.vcf.gz \
-O /home/projects/cu_10098/data/ennbin/project/hg37_testWithNewgnomAD/filtration_test/GetPileupSummaries.table

and it always fails, and I cant find where the error is.......
the error message is shown below:

Using GATK jar /home/projects/cu_10098/apps/src/gatk-4.beta.5/gatk-package-4.beta.5-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=1 -Dsnappy.disable=true -Xmx20g -jar /home/projects/cu_10098/apps/src/gatk-4.beta.5/gatk-package-4.beta.5-local.jar GetPileupSummaries -R /home/databases/gatk-legacy-bundles/b37/human_g1k_v37_decoy.fasta -I tumor.bam -V /home/projects/cu_10098/data/public_data/gnomad.exomes.r2.0.2.sites.vcf.gz -O /home/projects/cu_10098/data/ennbin/project/hg37_testWithNewgnomAD/filtration_test/GetPileupSummaries.table
19:19:21.288 INFO NativeLibraryLoader - Loading libgkl_compression.so from jar:file:/home/projects/cu_10098/apps/src/gatk-4.beta.5/gatk-package-4.beta.5-local.jar!/com/intel/gkl/native/libgkl_compression.so
[October 9, 2017 7:19:21 PM CEST] GetPileupSummaries --output /home/projects/cu_10098/data/ennbin/project/hg37_testWithNewgnomAD/filtration_test/GetPileupSummaries.table --variant /home/projects/cu_10098/data/public_data/gnomad.exomes.r2.0.2.sites.vcf.gz --input tumor.bam --reference /home/databases/gatk-legacy-bundles/b37/human_g1k_v37_decoy.fasta --minimumPopulationAlleleFrequency 0.01 --maximumPopulationAlleleFrequency 0.2 --interval_set_rule UNION --interval_padding 0 --interval_exclusion_padding 0 --interval_merging_rule ALL --readValidationStringency SILENT --secondsBetweenProgressUpdates 10.0 --disableSequenceDictionaryValidation false --createOutputBamIndex true --createOutputBamMD5 false --createOutputVariantIndex true --createOutputVariantMD5 false --lenient false --addOutputSAMProgramRecord true --addOutputVCFCommandLine true --cloudPrefetchBuffer 40 --cloudIndexPrefetchBuffer -1 --disableBamIndexCaching false --help false --version false --showHidden false --verbosity INFO --QUIET false --use_jdk_deflater false --use_jdk_inflater false --gcs_max_retries 20 --disableToolDefaultReadFilters false
[October 9, 2017 7:19:21 PM CEST] Executing as [email protected] on Linux 3.10.0-514.10.2.el7.x86_64 amd64; OpenJDK 64-Bit Server VM 1.8.0_144-b01; Version: 4.beta.5
19:19:21.429 INFO GetPileupSummaries - HTSJDK Defaults.COMPRESSION_LEVEL : 1
19:19:21.429 INFO GetPileupSummaries - HTSJDK Defaults.USE_ASYNC_IO_READ_FOR_SAMTOOLS : false
19:19:21.429 INFO GetPileupSummaries - HTSJDK Defaults.USE_ASYNC_IO_WRITE_FOR_SAMTOOLS : true
19:19:21.429 INFO GetPileupSummaries - HTSJDK Defaults.USE_ASYNC_IO_WRITE_FOR_TRIBBLE : false
19:19:21.429 INFO GetPileupSummaries - Deflater: IntelDeflater
19:19:21.429 INFO GetPileupSummaries - Inflater: IntelInflater
19:19:21.429 INFO GetPileupSummaries - GCS max retries/reopens: 20
19:19:21.429 INFO GetPileupSummaries - Using google-cloud-java patch c035098b5e62cb4fe9155eff07ce88449a361f5d from https://github.com/droazen/google-cloud-java/tree/dr_all_nio_fixes
19:19:21.429 INFO GetPileupSummaries - Initializing engine
19:19:22.766 INFO FeatureManager - Using codec VCFCodec to read file file:///home/projects/cu_10098/data/public_data/gnomad.exomes.r2.0.2.sites.vcf.gz
19:19:23.569 INFO FeatureManager - Using codec VCFCodec to read file file:///home/projects/cu_10098/data/public_data/gnomad.exomes.r2.0.2.sites.vcf.gz
19:19:25.109 WARN IndexUtils - Feature file "/home/projects/cu_10098/data/public_data/gnomad.exomes.r2.0.2.sites.vcf.gz" appears to contain no sequence dictionary. Attempting to retrieve a sequence dictionary from the associated index file
19:19:25.184 WARN IndexUtils - Feature file "/home/projects/cu_10098/data/public_data/gnomad.exomes.r2.0.2.sites.vcf.gz" appears to contain no sequence dictionary. Attempting to retrieve a sequence dictionary from the associated index file
19:19:25.219 WARN IndexUtils - Feature file "/home/projects/cu_10098/data/public_data/gnomad.exomes.r2.0.2.sites.vcf.gz" appears to contain no sequence dictionary. Attempting to retrieve a sequence dictionary from the associated index file
19:19:25.254 WARN IndexUtils - Feature file "/home/projects/cu_10098/data/public_data/gnomad.exomes.r2.0.2.sites.vcf.gz" appears to contain no sequence dictionary. Attempting to retrieve a sequence dictionary from the associated index file
19:19:25.292 INFO GetPileupSummaries - Done initializing engine
19:19:25.292 INFO ProgressMeter - Starting traversal
19:19:25.293 INFO ProgressMeter - Current Locus Elapsed Minutes Variants Processed Variants/Minute
19:19:25.310 INFO GetPileupSummaries - Shutting down engine
[October 9, 2017 7:19:25 PM CEST] org.broadinstitute.hellbender.tools.walkers.contamination.GetPileupSummaries done. Elapsed time: 0.07 minutes.
Runtime.totalMemory()=1781006336
java.lang.NumberFormatException: For input string: "."
at sun.misc.FloatingDecimal.readJavaFormatString(FloatingDecimal.java:2043)
at sun.misc.FloatingDecimal.parseDouble(FloatingDecimal.java:110)
at java.lang.Double.parseDouble(Double.java:538)
at java.lang.Double.valueOf(Double.java:502)
at htsjdk.variant.variantcontext.CommonInfo.getAttributeAsDouble(CommonInfo.java:323)
at htsjdk.variant.variantcontext.VariantContext.getAttributeAsDouble(VariantContext.java:732)
at org.broadinstitute.hellbender.tools.walkers.contamination.GetPileupSummaries.alleleFrequencyInRange(GetPileupSummaries.java:111)
at org.broadinstitute.hellbender.tools.walkers.contamination.GetPileupSummaries.apply(GetPileupSummaries.java:94)
at org.broadinstitute.hellbender.engine.VariantWalkerBase.lambda$traverse$0(VariantWalkerBase.java:110)
at java.util.stream.ForEachOps$ForEachOp$OfRef.accept(ForEachOps.java:184)
at java.util.stream.ReferencePipeline$2$1.accept(ReferencePipeline.java:175)
at java.util.Iterator.forEachRemaining(Iterator.java:116)
at java.util.Spliterators$IteratorSpliterator.forEachRemaining(Spliterators.java:1801)
at java.util.stream.AbstractPipeline.copyInto(AbstractPipeline.java:481)
at java.util.stream.AbstractPipeline.wrapAndCopyInto(AbstractPipeline.java:471)
at java.util.stream.ForEachOps$ForEachOp.evaluateSequential(ForEachOps.java:151)
at java.util.stream.ForEachOps$ForEachOp$OfRef.evaluateSequential(ForEachOps.java:174)
at java.util.stream.AbstractPipeline.evaluate(AbstractPipeline.java:234)
at java.util.stream.ReferencePipeline.forEach(ReferencePipeline.java:418)
at org.broadinstitute.hellbender.engine.VariantWalkerBase.traverse(VariantWalkerBase.java:108)
at org.broadinstitute.hellbender.engine.GATKTool.doWork(GATKTool.java:838)
at org.broadinstitute.hellbender.cmdline.CommandLineProgram.runTool(CommandLineProgram.java:119)
at org.broadinstitute.hellbender.cmdline.CommandLineProgram.instanceMainPostParseArgs(CommandLineProgram.j>ava:176)
at org.broadinstitute.hellbender.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:195)
at org.broadinstitute.hellbender.Main.runCommandLineProgram(Main.java:131)
at org.broadinstitute.hellbender.Main.mainEntry(Main.java:152)
at org.broadinstitute.hellbender.Main.main(Main.java:233)

Answers

Sign In or Register to comment.