The current GATK version is 3.6-0
Examples: Monday, today, last week, Mar 26, 3/26/04

Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

Powered by Vanilla. Made with Bootstrap.
Last chance to register for the GATK workshop next week in Basel, Switzerland! http://www.sib.swiss/training/upcoming-training-events/training/gatk-workshop-lecture

Expected file size - Haplotype Caller

bvecchiobvecchio Posts: 7Member
edited October 2012 in Ask the GATK team

Hi All,
I've been attempting to use the haplotype caller on my 50x coverage exome data. The bam being parsed is about 12G. Each time, the caller runs for many hours and then the output is only the header of the VCF - no errors seen. I'm wondering if this is due to limited space on my drives or if the expected file size is much larger than I am anticipating.

Command:

GenomeAnalysisTK.jar -T HaplotypeCaller -R  Homo_sapiens_assembly19.fasta -I input.bam --dbsnp dbsnp_132.b37.nochr.vcf  -stand_call_conf 30 -stand_emit_conf 10 -o output.Haplotypes.vcf

Answers

  • Geraldine_VdAuweraGeraldine_VdAuwera Posts: 10,469Administrator, Dev admin

    Hi there,

    Can you clarify -- is the run completing "successfully" (the program finishes without error message) or are you saying you're not seeing anything when you look in the file after some hours but before the run is complete?

    Geraldine Van der Auwera, PhD

  • bvecchiobvecchio Posts: 7Member

    The run completes successfully, without errors. I view the file after completion, and there is no data other than the header.

  • Geraldine_VdAuweraGeraldine_VdAuwera Posts: 10,469Administrator, Dev admin

    I see. Have you tried running on a subset of your file? Since you're working with an exome, I assume you have target intervals? You can pass those with -L; if your problem is due to file size that may work around it.

    Geraldine Van der Auwera, PhD

Sign In or Register to comment.