To celebrate the release of GATK 4.0, we are giving away free credits for running the GATK4 Best Practices pipelines in FireCloud, our secure online analysis portal. It’s first come first serve, so sign up now to claim your free credits worth $250. Sponsored by Google Cloud. Learn more at

Haplotype caller in GVCF mode still taking a very long time. Can I possibly speed up the process?

apolyakapolyak State College, PA, USAMember

I have 13 whole exome sequencing samples, and unfortunately, I'm having a hard time getting HaplotypeCaller to complete within the time frame the cluster I use allows (150 hours). I use 10 nodes at a time with 10gb ram with 8 cores per node. Is there any way to speed up this rate? I tried using HaplotypeCaller in GVCF mode with the following command:

java -d64 -Xmx8g -jar $GATKDIR/GenomeAnalysisTK.jar \
-T HaplotypeCaller \
-R $REF --dbsnp $DBSNP \
-I 7-27_realigned.bam \
-o 7-27_hg19.vcf \
-gt_mode DISCOVERY \
-mbq 20 \
-stand_emit_conf 20 -G Standard -A AlleleBalance -nct 16 \
--emitRefConfidence GVCF --variant_index_type LINEAR --variant_index_parameter 128000

Am I doing something incorrectly? Is there anything I can tweak to minimize the runtime? What is the expected runtime for WES on a standard setup (a few cores and some ram)?

Best Answer


Sign In or Register to comment.