Test-drive the GATK tools and Best Practices pipelines on Terra

Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.

Seems Not The Right Genotype of My variants

Hello this is my GT:AD:DP:GQ:PL information of my variants : 0/1:4,42:46:11:1429,0,11 . It seems that "1/1" might the best genotype of my variant, can anybody tell me the reason.
This is my command: genotyping_mode DISCOVERY --emitRefConfidence GVCF -stand_call_conf 10



  • SheilaSheila Broad InstituteMember, Broadie, Moderator admin
    edited March 2018


    If you look at the PLs, you see the tool is not super confident whether the genotype is het or hom var. In this case, there was a bit more confidence in the het call, probably due to high base qualities for the ref allele. If you can use the Genotype Refinement Workflow, that may help with cases like this.


    P.S. Have a look in the Methods and Algorithms section for an explanation of how genotypes are called.

Sign In or Register to comment.