We've moved!
This site is now read-only. You can find our new documentation site and support forum for posting questions here.
Be sure to read our welcome blog!

HaplotypeCaller: Record size mismatch! Something went wrong in the merging of alleles


I've trying to use HaplotypeCaller with queue. The qscript is attached. I managed to run it on 2 samples without problem. It also succeed with 44 samples for calling variants on chromosome 1. But when I run the script with 44 samples for the whole exome calling, the following error occurred:

ERROR ------------------------------------------------------------------------------------------
ERROR stack trace

org.broadinstitute.sting.utils.exceptions.ReviewedStingException: Record size mismatch! Something went wrong in the merging of alleles.
at org.broadinstitute.sting.gatk.walkers.haplotypecaller.GenotypingEngine.assignGenotypeLikelihoods(GenotypingEngine.java:217)
at org.broadinstitute.sting.gatk.walkers.haplotypecaller.HaplotypeCaller.map(HaplotypeCaller.java:500)
at org.broadinstitute.sting.gatk.walkers.haplotypecaller.HaplotypeCaller.map(HaplotypeCaller.java:132)
at org.broadinstitute.sting.gatk.traversals.TraverseActiveRegions.processActiveRegion(TraverseActiveRegions.java:552)
at org.broadinstitute.sting.gatk.traversals.TraverseActiveRegions.processActiveRegions(TraverseActiveRegions.java:512)
at org.broadinstitute.sting.gatk.traversals.TraverseActiveRegions.traverse(TraverseActiveRegions.java:244)
at org.broadinstitute.sting.gatk.traversals.TraverseActiveRegions.traverse(TraverseActiveRegions.java:69)
at org.broadinstitute.sting.gatk.executive.LinearMicroScheduler.execute(LinearMicroScheduler.java:100)
at org.broadinstitute.sting.gatk.GenomeAnalysisEngine.execute(GenomeAnalysisEngine.java:283)
at org.broadinstitute.sting.gatk.CommandLineExecutable.execute(CommandLineExecutable.java:113)
at org.broadinstitute.sting.commandline.CommandLineProgram.start(CommandLineProgram.java:245)
at org.broadinstitute.sting.commandline.CommandLineProgram.start(CommandLineProgram.java:152)
at org.broadinstitute.sting.gatk.CommandLineGATK.main(CommandLineGATK.java:91)

ERROR ------------------------------------------------------------------------------------------
ERROR A GATK RUNTIME ERROR has occurred (version 2.4-9-g532efad):
ERROR Please visit the wiki to see if this is a known problem
ERROR If not, please post the error, with stack trace, to the GATK forum
ERROR Visit our website and forum for extensive documentation and answers to
ERROR commonly asked questions http://www.broadinstitute.org/gatk
ERROR MESSAGE: Record size mismatch! Something went wrong in the merging of alleles.
ERROR ------------------------------------------------------------------------------------------

The command line I use:

java -Djava.io.tmpdir=$TEMP_DIR -Xmx12g -Xms12g -jar Queue-2.4-9/Queue.jar -S HaplotypeCaller.scala -I $BAM_LIST_FILE -R $REF_FILE -sg $SCATTER_NUMBER -jobRunner GridEngine -V_out $OUTPUT_FILE -retry 2 -run

BAM_LIST_FILE contains a list of paths of 44 recalibrated bam files ( not reduced reads bams, I think HaplotypeCaller does support it).

Could anyone tell me if the error is bug or I did anything wrong? Many thanks,


Best Answer


  • byb121byb121 UKMember

    I will update you once I have more findings on this. Thanks.

  • byb121byb121 UKMember
    edited May 2013

    Update: The cause of the problem could not be found, but the lasted version of Queue (2.5.2) has fixed it.

  • CarneiroCarneiro Charlestown, MAMember admin
Sign In or Register to comment.