HaplotypeScore not annotated?
Posted in Ask the GATK team | Last updated on 2012-10-23 10:00:08

Comments (4)

Hi there, I'm running into an issue with my Queue pipeline, with variants called with HaplotypeCaller. Once I get the raw VCF file, I use VariantAnnotator before VQSR: however, no HaplotypeScore annotation is being produced, resulting in a failure of the VariantRecalibrator step where 'HaplotypeScore' was indicated as an annotation.

I tried to correct the issue by indicating to VariantAnnotator to use all annotations

  class AnnotationArguments (t: Target) extends VariantAnnotator with UNIVERSAL_GATK_ARGS {
this.reference_sequence = qscript.referenceFile
// Set the memory limit to 7 gigabytes on each command.
    this.memoryLimit = 7
    this.input_file :+= qscript.bamFile
    this.D = qscript.dbSNP_b37

But I still can't get any output in the annotated VCF of that parameter. Here an example of a variant


Any suggestions on what I might be doing wrong?

thanks very much for your help, Francesco

Return to top Comment on this article in the forum