SNP genotyping algorithm consraints on +/- read proportion

gene

Hi,

I am observing the following scenario at one particular SNP (C/G) using two different enrichment technologies: (I am using IGV syntax: ALLELE|number of reads w/ allele|%of total reads|+strand reads|- strand reads) - technology1:
C: 15 47% 15+,0- G: 17 53% 17+,0- - technology2:
C: 17 37% 13+,4- G: 29 63% 26+,3- As you can see both technologies have good coverage of the SNP and also good representation of each allele. SNP(C/G) does not get called in technology1.

My questions are: 1- Does the GATK algorithm have some sort of constraint on the proportion of reads coming from only one strand (as with technology1) in order to try to predict or discard duplicates? 2- I know that the base call of a particular base is bounded by the mapping quality of its read. If my --stand_call_conf is 30 and one of the bases at this SNP position has MQ<30 does this avoid this position getting called? Or is it more like the avg(MQ) has to be >30 (meaning more than one read at this position is taken into account)?

Thanks for any clarification, Gene

ebanks

1) No, we don't mark duplicates (that's done in a previous data processing step as per our best practices). 2) No, the calling confidence threshold is not the same thing as a base quality cutoff; it is the value of the QUAL in the VCF record (see the documentation for how they are related).

gene

Hi,

Thanks for your answer.

I looked in the paper "A framework for variation discovery and genotyping using next-generation DNA sequencing data" but could not find the relationship between QUAL and the base calls for a specific SNP position. Would it be possible to point me to the documentation of this relationship?

Thx, Gene

Geraldine_VdAuwera

Hi Gene,

The math is detailed in this document. Note that it covers single-allele cases only; for multi-allelic calls, see this slide.

newbie16

Hi I also wanted to check out the documents you mentioned above. However the links don't work. Could you post new links? Thanks

Geraldine_VdAuwera

I fixed the links in the above comment, thanks for reporting them.

newbie16

Thanks !