Genome Analysis Toolkit

Variant Discovery in High-Throughput Sequencing Data

GATK process banner

Need Help?

Search our documentation

Community Forum

Hi, How can we help?

Developed in the Data Sciences Platform at the Broad Institute, the toolkit offers a wide variety of tools with a primary focus on variant discovery and genotyping. Its powerful processing engine and high-performance computing features make it capable of taking on projects of any size. Learn more

HaplotypeCaller Follow

5 comments

  • 0
    Comment actions Permalink
  • Avatar
    Devendra Waikul

    I am helping setting up GATK4 pipeline. I am getting Java fatal error for individual chromosome intervals on genomicsDBimport. I am guessing I must have set the wrong output file extension for haplotypecaller output. 
    I have set the output to be vcf_haplotype_output = output_haplotype.vcf with flag set to -ERC GVCF

    Can anyone suggest what should be an output extension for this command?

     

    0
    Comment actions Permalink
  • Avatar
    Monete Rajão Gomes

    Hi Devendra,

     

    I think when you use -ERC GVCF option, your output should be "output_haplotype.g.vcf" or (better) "output_haplotype.g.vcf.gz"

    Like described above on session: Single-sample GVCF calling (outputs intermediate GVCF)

    Hope this helps.

    0
    Comment actions Permalink
  • Avatar
    Devendra Waikul

    Hello Monete Rajão Gomes,

    Thank you for providing the quick response. 

    I have one more question. Did you face java errors after the genomicsDBimport command?

    15:33:00.356 INFO ProgressMeter - Starting traversal
    15:33:00.356 INFO ProgressMeter - Current Locus Elapsed Minutes Batches Processed Batches/Minute
    15:33:00.495 INFO GenomicsDBImport - Starting batch input file preload
    15:33:48.560 INFO GenomicsDBImport - Finished batch preload
    15:33:48.563 INFO GenomicsDBImport - Importing batch 1 with 31 samples
    03:08:43.526 INFO ProgressMeter - chr22:1 695.7 1 0.0
    03:08:43.531 INFO GenomicsDBImport - Done importing batch 1/1
    03:08:43.533 INFO ProgressMeter - chr22:1 695.7 1 0.0
    03:08:43.533 INFO ProgressMeter - Traversal complete. Processed 1 total batches in 695.7 minutes.
    03:08:43.533 INFO GenomicsDBImport - Import completed!
    03:08:43.533 INFO GenomicsDBImport - Shutting down engine
    [June 18, 2020 3:08:43 AM EDT] org.broadinstitute.hellbender.tools.genomicsdb.GenomicsDBImport done. Elapsed time: 695.79 minutes.
    Runtime.totalMemory()=28427943936
    Tool returned:
    true
    #
    # A fatal error has been detected by the Java Runtime Environment:
    #
    # SIGBUS (0x7) at pc=0x00007fca5fb9d330, pid=388337, tid=0x00007fda500b6740
    #
    # JRE version: OpenJDK Runtime Environment (8.0_181-b13) (build 1.8.0_181-b13)
    # Java VM: OpenJDK 64-Bit Server VM (25.181-b13 mixed mode linux-amd64 )
    # Problematic frame:
    #

    0
    Comment actions Permalink
  • Avatar
    Monete Rajão Gomes

    Hi Devendra Waikul

    Unfortunately, I've never used genomicsDBImport before. I've been using combinegvcf, instead.

    Maybe you should look for this error in other sites or in read-only GATK v3 previous site.

    0
    Comment actions Permalink

Please sign in to leave a comment.

Powered by Zendesk