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.

HaplotypeCaller producing output without any vcf

Hello,
I am using GATK 4.0.8.1in Ubuntu 18.04 LTS and javac -version => javac 1.8.0_181.

I am basically trying to call variants from Targetted DNA sequencing data using HaplotypeCaller. Though HaplotypeCaller is executing but it is not calling any variants. Now seeing some other posts I did ensure that i have java 8 and validate bam file using picard ValidateSam, and it gave no error. I'm confused. Below is the output of the haplotypecaller run.

Using GATK jar /opt/gatk-4.0.8.1/gatk-package-4.0.8.1-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 /opt/gatk-4.0.8.1/gatk-package-4.0.8.1-local.jar HaplotypeCaller --input BV-15-186_mapto_genom_sortedRG.bam --output BV-15-186_bwa_vars.vcf --reference hg19.fa 17:31:19.630 INFO NativeLibraryLoader - Loading libgkl_compression.so from jar:file:/opt/gatk-4.0.8.1/gatk-package-4.0.8.1-local.jar!/com/intel/gkl/native/libgkl_compression.so 17:31:29.729 INFO HaplotypeCaller - ------------------------------------------------------------ 17:31:29.730 INFO HaplotypeCaller - The Genome Analysis Toolkit (GATK) v4.0.8.1 17:31:29.730 INFO HaplotypeCaller - For support and documentation go to https://software.broadinstitute.org/gatk/ 17:31:29.731 INFO HaplotypeCaller - Executing as [email protected] on Linux v4.15.0-20-generic amd64 17:31:29.731 INFO HaplotypeCaller - Java runtime: OpenJDK 64-Bit Server VM v10.0.2+13-Ubuntu-1ubuntu0.18.04.2 17:31:29.732 INFO HaplotypeCaller - Start Date/Time: 15 September 2018 at 5:31:19 PM IST 17:31:29.732 INFO HaplotypeCaller - ------------------------------------------------------------ 17:31:29.732 INFO HaplotypeCaller - ------------------------------------------------------------ 17:31:29.733 INFO HaplotypeCaller - HTSJDK Version: 2.16.0 17:31:29.734 INFO HaplotypeCaller - Picard Version: 2.18.7 17:31:29.734 INFO HaplotypeCaller - HTSJDK Defaults.COMPRESSION_LEVEL : 2 17:31:29.734 INFO HaplotypeCaller - HTSJDK Defaults.USE_ASYNC_IO_READ_FOR_SAMTOOLS : false 17:31:29.734 INFO HaplotypeCaller - HTSJDK Defaults.USE_ASYNC_IO_WRITE_FOR_SAMTOOLS : true 17:31:29.734 INFO HaplotypeCaller - HTSJDK Defaults.USE_ASYNC_IO_WRITE_FOR_TRIBBLE : false 17:31:29.734 INFO HaplotypeCaller - Deflater: IntelDeflater 17:31:29.735 INFO HaplotypeCaller - Inflater: IntelInflater 17:31:29.735 INFO HaplotypeCaller - GCS max retries/reopens: 20 17:31:29.735 INFO HaplotypeCaller - Using google-cloud-java fork https://github.com/broadinstitute/google-cloud-java/releases/tag/0.20.5-alpha-GCS-RETRY-FIX 17:31:29.735 INFO HaplotypeCaller - Initializing engine 17:31:29.884 INFO HaplotypeCaller - Done initializing engine 17:31:29.891 INFO HaplotypeCallerEngine - Disabling physical phasing, which is supported only for reference-model confidence output 17:31:29.899 INFO NativeLibraryLoader - Loading libgkl_utils.so from jar:file:/opt/gatk-4.0.8.1/gatk-package-4.0.8.1-local.jar!/com/intel/gkl/native/libgkl_utils.so 17:31:29.901 INFO NativeLibraryLoader - Loading libgkl_pairhmm_omp.so from jar:file:/opt/gatk-4.0.8.1/gatk-package-4.0.8.1-local.jar!/com/intel/gkl/native/libgkl_pairhmm_omp.so 17:31:29.939 WARN IntelPairHmm - Flush-to-zero (FTZ) is enabled when running PairHMM 17:31:29.940 INFO IntelPairHmm - Available threads: 4 17:31:29.940 INFO IntelPairHmm - Requested threads: 4 17:31:29.940 INFO PairHMM - Using the OpenMP multi-threaded AVX-accelerated native PairHMM implementation 17:31:29.965 INFO ProgressMeter - Starting traversal 17:31:29.965 INFO ProgressMeter - Current Locus Elapsed Minutes Regions Processed Regions/Minute 17:31:29.971 INFO VectorLoglessPairHMM - Time spent in setup for JNI call : 0.0 17:31:29.971 INFO PairHMM - Total compute time in PairHMM computeLogLikelihoods() : 0.0 17:31:29.972 INFO SmithWatermanAligner - Total compute time in java Smith-Waterman : 0.00 sec 17:31:29.972 INFO HaplotypeCaller - Shutting down engine [15 September 2018 at 5:31:29 PM IST] org.broadinstitute.hellbender.tools.walkers.haplotypecaller.HaplotypeCaller done. Elapsed time: 0.17 minutes. Runtime.totalMemory()=314572800 Exception in thread "main" java.lang.IncompatibleClassChangeError: Inconsistent constant pool data in classfile for class org/broadinstitute/hellbender/transformers/ReadTransformer. Method lambda$identity$d67512bf$1(Lorg/broadinstitute/hellbender/utils/read/GATKRead;)Lorg/broadinstitute/hellbender/utils/read/GATKRead; at index 65 is CONSTANT_MethodRef and should be CONSTANT_InterfaceMethodRef at org.broadinstitute.hellbender.transformers.ReadTransformer.identity(ReadTransformer.java:30) at org.broadinstitute.hellbender.engine.GATKTool.makePreReadFilterTransformer(GATKTool.java:290) at org.broadinstitute.hellbender.engine.AssemblyRegionWalker.traverse(AssemblyRegionWalker.java:262) at org.broadinstitute.hellbender.engine.GATKTool.doWork(GATKTool.java:979) at org.broadinstitute.hellbender.cmdline.CommandLineProgram.runTool(CommandLineProgram.java:137) at org.broadinstitute.hellbender.cmdline.CommandLineProgram.instanceMainPostParseArgs(CommandLineProgram.java:182) at org.broadinstitute.hellbender.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:201) at org.broadinstitute.hellbender.Main.runCommandLineProgram(Main.java:160) at org.broadinstitute.hellbender.Main.mainEntry(Main.java:203) at org.broadinstitute.hellbender.Main.main(Main.java:289)

Answers

Sign In or Register to comment.