GATK v4.0.1.1 HaplotypeCaller

I am using GATK v4.0.1.1 HaplotypeCaller for variant analysis. (paired-end DNA sequenced data mapped to the reference using BWA mem).

The command I used;
“gatk HaplotypeCaller –R Reference.fna –I input.bam –O output.vcf”

It runs for a while (couple of seconds) but does not produce an output. No error message was given. Am I doing this right? Any help is appreciated.


  • SheilaSheila Broad InstituteMember, Broadie, Moderator


    What do you mean "does not produce an output"? Can you post the entire log output (all the lines after the tool starts to run)?


  • SKWSKW Member

    Hi Sheila,
    After running the command there's no output file to be found. Any advice would be appreciated.
    Thank you!

    Here is the entire log output;

    [[email protected] project]$ /local/cluster/bin/gatk HaplotypeCaller -R Reference_genome.fna -I GH-50-a_aln_pairs_mapped_sorted.bam -O GH-50-a_GATK_output.vcf
    Using GATK jar /local/cluster/gatk-
    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 -jar /local/cluster/gatk- HaplotypeCaller -R Reference_genome.fna -I GH-50-a_aln_pairs_mapped_sorted.bam -O GH-50-a_GATK_output.vcf
    18:50:41.895 INFO NativeLibraryLoader - Loading from jar:file:/local/cluster/gatk-!/com/intel/gkl/native/
    18:50:42.173 INFO HaplotypeCaller - ------------------------------------------------------------
    18:50:42.174 INFO HaplotypeCaller - The Genome Analysis Toolkit (GATK) v4.0.1.1
    18:50:42.174 INFO HaplotypeCaller - For support and documentation go to
    18:50:42.175 INFO HaplotypeCaller - Executing as [email protected] on Linux v3.10.0-327.el7.x86_64 amd64
    18:50:42.175 INFO HaplotypeCaller - Java runtime: Java HotSpot(TM) 64-Bit Server VM v1.8.0_71-b15
    18:50:42.176 INFO HaplotypeCaller - Start Date/Time: February 19, 2018 6:50:41 PM PST
    18:50:42.176 INFO HaplotypeCaller - ------------------------------------------------------------
    18:50:42.176 INFO HaplotypeCaller - ------------------------------------------------------------
    18:50:42.177 INFO HaplotypeCaller - HTSJDK Version: 2.14.1
    18:50:42.177 INFO HaplotypeCaller - Picard Version: 2.17.2
    18:50:42.177 INFO HaplotypeCaller - HTSJDK Defaults.COMPRESSION_LEVEL : 1
    18:50:42.177 INFO HaplotypeCaller - HTSJDK Defaults.USE_ASYNC_IO_READ_FOR_SAMTOOLS : false
    18:50:42.178 INFO HaplotypeCaller - HTSJDK Defaults.USE_ASYNC_IO_WRITE_FOR_SAMTOOLS : true
    18:50:42.178 INFO HaplotypeCaller - HTSJDK Defaults.USE_ASYNC_IO_WRITE_FOR_TRIBBLE : false
    18:50:42.178 INFO HaplotypeCaller - Deflater: IntelDeflater
    18:50:42.178 INFO HaplotypeCaller - Inflater: IntelInflater
    18:50:42.178 INFO HaplotypeCaller - GCS max retries/reopens: 20
    18:50:42.178 INFO HaplotypeCaller - Using google-cloud-java patch 6d11bef1c81f885c26b2b56c8616b7a705171e4f from
    18:50:42.179 INFO HaplotypeCaller - Initializing engine
    18:50:43.288 INFO HaplotypeCaller - Done initializing engine
    18:50:43.393 INFO HaplotypeCallerEngine - Disabling physical phasing, which is supported only for reference-model confidence output
    18:50:44.112 INFO HaplotypeCaller - Shutting down engine
    [February 19, 2018 6:50:44 PM PST] done. Elapsed time: 0.04 minutes.
    java.lang.IllegalArgumentException: samples cannot be empty
    at org.broadinstitute.hellbender.utils.Utils.validateArg(
    at org.broadinstitute.hellbender.engine.GATKTool.doWork(
    at org.broadinstitute.hellbender.cmdline.CommandLineProgram.runTool(
    at org.broadinstitute.hellbender.cmdline.CommandLineProgram.instanceMainPostParseArgs(
    at org.broadinstitute.hellbender.cmdline.CommandLineProgram.instanceMain(
    at org.broadinstitute.hellbender.Main.runCommandLineProgram(
    at org.broadinstitute.hellbender.Main.mainEntry(
    at org.broadinstitute.hellbender.Main.main(
    14.625u 1.297s 0:05.17 307.7% 0+0k 0+1184io 0pf+0w

  • SkyWarriorSkyWarrior TurkeyMember
    edited February 20

    There is no output because there is an error.

    java.lang.IllegalArgumentException: samples cannot be empty

    This is an error statement.

    There seems to be a problem with the command or with your input files.

  • SheilaSheila Broad InstituteMember, Broadie, Moderator
    edited February 23


    Indeed @SkyWarrior is correct. Can you run ValidateSamFile on your input BAM file.


Sign In or Register to comment.