The current GATK version is 3.7-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.
GATK 3.7 is here! Be sure to read the Version Highlights and optionally the full Release Notes.
Register now for the upcoming GATK Best Practices workshop, Feb 20-22 in Leuven, Belgium. Open to all comers! More info and signup at

Haplotype Score Region

chongmchongm Member Posts: 33

Hi Everyone,

I had a few questions about the haplotype score.

In the technical documentation it states that "Higher scores are indicative of regions with bad alignments, often leading to artifactual SNP and indel calls. Note that the Haplotype Score is only calculated for sites with read coverage."

How is the haplotype group for each variant site determined? e.g. Does it take the closest two variants to the query site and then treat the query variant + closest two variants as the haplotype group?

Also, in the case of multiallelic SNPs (>2 SNPs), haplotype score is inappropriate since it only looks at whether a site can be explained by the segregation of two and only two haplotypes, correct? So multiallelic snps will be assigned poor haplotype scores OR will these sites not be annotated at all? If we have a case where there is a truly biallelic SNP and a couple of samples have some reads that are erroneously calling a third allele, this variant site will be assigned a poor haplotype score overall, correct?





Best Answers


  • chongmchongm Member Posts: 33

    Hi Dr. Banks,

    Okay the second part is clear. Regarding the first question I just want to clarify. If there is no connection between previous or successive variants then does the "haplotype" represents 21 bp nucleotide sequences of an individual read? So for an ideal heterozygous site, we would expect two exactly two "haplotypes" or "piles of reads", while in a less than ideal situation some reads will contain SNP artefacts and contribute to additional "haplotypes"? Thus, the probability that a variant site can be explained by 2 segregating piles of reads is diluted?



Sign In or Register to comment.