Test-drive the GATK tools and Best Practices pipelines on Terra


Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.

Stage 7 of CNVDiscoveryPipeline - all commands fail with java.lang.NullPointerException

d_a_parryd_a_parry EdinburghMember

I am attempting to run GenomeSTRiP's CNVDiscoveryPipeline. I have run SVPreprocess successfully and am now running CNVDiscoveryPipeline on a single chromosome. Each command at stage 7 fails with the following stack trace:

Exception in thread "main" java.lang.RuntimeException
    at org.broadinstitute.sv.commandline.CommandLineProgram.execute(CommandLineProgram.java:65)
    at org.broadinstitute.gatk.utils.commandline.CommandLineProgram.start(CommandLineProgram.java:248)
    at org.broadinstitute.gatk.utils.commandline.CommandLineProgram.start(CommandLineProgram.java:155)
    at org.broadinstitute.sv.commandline.CommandLineProgram.runAndReturnResult(CommandLineProgram.java:29)
    at org.broadinstitute.sv.commandline.CommandLineProgram.run(CommandLineProgram.java:25)
    at org.broadinstitute.sv.genotyping.RefineCNVBoundaries.main(RefineCNVBoundaries.java:133)
Caused by: java.lang.NullPointerException
    at org.broadinstitute.sv.genotyping.BoundaryRefinementAlgorithm.computeNextCoordinate(BoundaryRefinementAlgorithm.java:886)
    at org.broadinstitute.sv.genotyping.BoundaryRefinementAlgorithm.getRefinedInterval(BoundaryRefinementAlgorithm.java:704)
    at org.broadinstitute.sv.genotyping.BoundaryRefinementAlgorithm.refineOneBoundary(BoundaryRefinementAlgorithm.java:612)
    at org.broadinstitute.sv.genotyping.BoundaryRefinementAlgorithm.refineBoundaryStep(BoundaryRefinementAlgorithm.java:548)
    at org.broadinstitute.sv.genotyping.BoundaryRefinementAlgorithm.refineBoundaries(BoundaryRefinementAlgorithm.java:536)
    at org.broadinstitute.sv.genotyping.BoundaryRefinementAlgorithm.processVariant(BoundaryRefinementAlgorithm.java:232)
    at org.broadinstitute.sv.genotyping.RefineCNVBoundaries.run(RefineCNVBoundaries.java:204)
    at org.broadinstitute.sv.commandline.CommandLineProgram.execute(CommandLineProgram.java:54)
    ... 5 more

An example command for this stage is:

'java'  '-Xmx2048m'  '-XX:+UseParallelOldGC'  '-XX:ParallelGCThreads=4'  '-XX:GCTimeLimit=50'  '-XX:GCHeapFreeLimit=10' \
-Djava.io.tmpdir=/exports/eddie/scratch/dparry/bams/.queue/tmp'  \
'-cp' '/exports/igmm/eddie/aitman-lab/svtoolkit/lib/SVToolkit.jar:/exports/igmm/eddie/aitman-lab/svtoolkit/lib/gatk/GenomeAnalysisTK.jar:/exports/igmm/eddie/aitman-lab/svtoolkit/lib/gatk/Queue.jar' \
'-cp' '/exports/igmm/eddie/aitman-lab/svtoolkit/lib/SVToolkit.jar:/exports/igmm/eddie/aitman-lab/svtoolkit/lib/gatk/GenomeAnalysisTK.jar:/exports/igmm/eddie/aitman-lab/svtoolkit/lib/gatk/Queue.jar' \
'org.broadinstitute.sv.genotyping.RefineCNVBoundaries' \
'-I' '/exports/eddie/scratch/dparry/sids_bams/per_chrom/cnv_discovery_chr22/cnv_stage6/seq_chr22/seq_chr22.merged_headers.bam' \
'-O' '/exports/eddie/scratch/dparry/sids_bams/per_chrom/cnv_discovery_chr22/cnv_stage7/seq_chr22/P0354/seq_chr22.merged.brig.vcf' \
'-R' '/exports/igmm/software/pkg/el7/apps/bcbio/share2/genomes/Hsapiens/hg38-noalt/seq/hg38-noalt.fa' \
'-md' 'sv_meta_data'  \
'-configFile' '/exports/igmm/eddie/aitman-lab/svtoolkit/conf/genstrip_parameters.txt' \
'-P' 'depth.readCountCacheIgnoreGenomeMask:true'  \
'-genderMapFile' 'sv_meta_data/sample_gender.report.txt'  \
'-ploidyMapFile' '/exports/igmm/eddie/aitman-lab/svtoolkit/Homo_sapiens_assembly38/Homo_sapiens_assembly38.ploidymap.txt'  \
'-vcf' 'per_chrom/cnv_discovery_chr22/cnv_stage4/seq_chr22/seq_chr22.merged.genotypes.vcf.gz' \
'-site' 'per_chrom/cnv_discovery_chr22/cnv_stage7/seq_chr22/P0354.sites.list'  \
'-boundaryPrecision' '100' \
'-minimumRefinedLength' '500'  \
'-maximumReferenceGapLength' '1000'         

I have tested this on a different chromosome and get the same error. Any help figuring out what is going wrong would be much appreciated!

Thanks

Dave

Answers

  • nmkeurnmkeur NetherlandsMember
    edited December 2016
    
    Exception in thread "main" java.lang.RuntimeException
        at org.broadinstitute.sv.commandline.CommandLineProgram.execute(CommandLineProgram.java:65)
        at org.broadinstitute.gatk.utils.commandline.CommandLineProgram.start(CommandLineProgram.java:248)
        at org.broadinstitute.gatk.utils.commandline.CommandLineProgram.start(CommandLineProgram.java:155)
        at org.broadinstitute.sv.commandline.CommandLineProgram.runAndReturnResult(CommandLineProgram.java:29)
        at org.broadinstitute.sv.commandline.CommandLineProgram.run(CommandLineProgram.java:25)
        at org.broadinstitute.sv.genotyping.RefineCNVBoundaries.main(RefineCNVBoundaries.java:133)
    Caused by: java.lang.NullPointerException
        at org.broadinstitute.sv.genotyping.BoundaryRefinementAlgorithm.computeNextCoordinate(BoundaryRefinementAlgorithm.java:886)
        at org.broadinstitute.sv.genotyping.BoundaryRefinementAlgorithm.getRefinedInterval(BoundaryRefinementAlgorithm.java:704)
        at org.broadinstitute.sv.genotyping.BoundaryRefinementAlgorithm.refineOneBoundary(BoundaryRefinementAlgorithm.java:612)
        at org.broadinstitute.sv.genotyping.BoundaryRefinementAlgorithm.refineBoundaryStep(BoundaryRefinementAlgorithm.java:548)
        at org.broadinstitute.sv.genotyping.BoundaryRefinementAlgorithm.refineBoundaries(BoundaryRefinementAlgorithm.java:536)
        at org.broadinstitute.sv.genotyping.BoundaryRefinementAlgorithm.processVariant(BoundaryRefinementAlgorithm.java:232)
        at org.broadinstitute.sv.genotyping.RefineCNVBoundaries.run(RefineCNVBoundaries.java:204)
        at org.broadinstitute.sv.commandline.CommandLineProgram.execute(CommandLineProgram.java:54)
        ... 5 more 
    

    I am having the exact same problem as above. The pipeline seems to fail on stage7.
    I really appreciate any help you can provide.

    Thanks

    Nick

Sign In or Register to comment.