Question: How do variant callers get genotype info?
gravatar for afzm
15 months ago by
afzm20 wrote:

How are variant callers able to compute which copy of the chrs in case of a diploid species has the heterozygous detected variants? What information do they use, just paired end reads?

I could find some statistical data in GATK webpage, but I would like to understand if there is other information used, the rationale behind it the accuracy it would have and the factors that affect this process.

Thank you very much

ADD COMMENTlink modified 15 months ago by Jeremy Leipzig18k • written 15 months ago by afzm20
gravatar for Jeremy Leipzig
15 months ago by
Philadelphia, PA
Jeremy Leipzig18k wrote:

Usually the caller has no idea which chromosome homolog a variant is on. It can just see variants that are in the same read or read pair (unlikely for short reads) or it can try to infer which variants are on the same chromosome homolog (phased) using read-backed phasing (as part of the read assembly performed by the haplotype caller).

These in silico methods are spotty at best. Most people who need phasing just use a long-read technology, or they sequence the parents.

ADD COMMENTlink modified 15 months ago • written 15 months ago by Jeremy Leipzig18k
gravatar for Pierre Lindenbaum
15 months ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum122k wrote:

Mathematical Notes on SAMtools Algorithms

.... good luck...

ADD COMMENTlink written 15 months ago by Pierre Lindenbaum122k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 896 users visited in the last hour