Complete this survey about your research needs and be entered to win an Amazon gift card or FireCloud credit.
Read more about it here!
Download the latest Picard release at https://github.com/broadinstitute/picard/releases.
GATK version 4.beta.6 is out. See the GATK4 beta page for download and details.

Haplotype Score Region

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?

Thanks,

MC

Best Answers

Answers

  • 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.