Question: calculating genotype quality from genotype likelihood
gravatar for Ana
3.3 years ago by
Ana170 wrote:

Hi all, I have faced with a question and I would appreciate if anyone has any comment or suggestion. I have got a vcf file made by FreeBayes, genotype quality is missing. I want to filter out those below the cutoff value. Therefore I want to know how can I get (calculate) genotype quality from genotype likelihood? This is my genotype field:

GT:DP:DPR:RO:QR:AO:QA:GL    ./0:11:4,0:7:255:0:0:-7.21754,-9.32475,-22.6408


ADD COMMENTlink modified 3.3 years ago by youcai50 • written 3.3 years ago by Ana170
gravatar for youcai
3.3 years ago by
youcai50 wrote:

bcftools may be help

bcftools view

Consensus/variant calling options:

       -c        SNP calling (force -e)
       -d FLOAT  skip loci where less than FLOAT fraction of samples covered [0]
       -e        likelihood based analyses
       -g        call genotypes at variant sites (force -c)
       -i FLOAT  indel-to-substitution ratio [-1]
       -I        skip indels
       -m FLOAT  alternative model for multiallelic and rare-variant calling, include if P(chi^2)>=FLOAT
       -p FLOAT  variant if P(ref|D)<FLOAT [0.5]
       -P STR    type of prior: full, cond2, flat [full]
       -t FLOAT  scaled substitution mutation rate [0.001]
       -T STR    constrained calling; STR can be: pair, trioauto, trioxd and trioxs (see manual) [null]
       -v        output potential variant sites only (force -c)
ADD COMMENTlink written 3.3 years ago by youcai50

I actually did find some ideas from the post by Philip How to interpret genotype likelihoods? They suggested to calculate genotype likelihood ratio. Does any one have any idea about good threshold value for this?

ADD REPLYlink written 3.3 years ago by Ana170
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: 1203 users visited in the last hour