Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • genotype likelihood inconsist with genotype field in mpileup output

    hi, all
    I call SNPs using mpileup from samtools, I got strange results like:

    Gm01 2923425 . A G,C 999 . DP=308;VDB=0.0224;AF1=0.6;G3=0.4,1.742e-06,0.6;HWE=0.00948;AC1=8;DP4=46,45,111,93;MQ=55;FQ=999;PV4=0.61,1,4.6e-0
    9,1 GT:PLP:SP:GQ 1/1:255,157,0,255,157,255:52:0:99 0/1:255,255,255,169,169,0:56:0:3 0/0:200,255,255,0,255,255:59:2:48 1/1:255,102,0,255,102,25
    5:34:0:99 0/0:0,154,255,154,255,255:51:0:99 0/1:255,255,255,90,90,0:30:0:3 1/1:255,39,0,255,39,255:13:0:37

    According VCF format specification:
    1000genomes.org is your first and best source for all of the information you’re looking for. From general topics to more of what you would expect to find here, 1000genomes.org has it all. We hope you find what you are searching for!


    the ordering of genotypes for the likelihoods is given by: F(j/k) = (k*(k+1)/2)+j. In other words, for biallelic sites the ordering is: AA,AB,BB; for triallelic sites the ordering is: AA,AB,BB,AC,BC,CC, etc.

    However, in my case, for example,
    GT:PLP:SP:GQ 0/1:255,255,255,169,169,0:56:0:3
    the six field from PL string has the lowest value-"0", thus the genotype might be CC, so the GT field should be 1/1, but I get 0/1 in this case.

    another example:
    0/0:200,255,255,0,255,255:59:2:48
    the forth field was "0", thus genotype might be A/C, and GT field must be 0/1, I get 0/0 in this case.

    I am confused about this.

Latest Articles

Collapse

  • seqadmin
    Understanding Genetic Influence on Infectious Disease
    by seqadmin




    During the COVID-19 pandemic, scientists observed that while some individuals experienced severe illness when infected with SARS-CoV-2, others were barely affected. These disparities left researchers and clinicians wondering what causes the wide variations in response to viral infections and what role genetics plays.

    Jean-Laurent Casanova, M.D., Ph.D., Professor at Rockefeller University, is a leading expert in this crossover between genetics and infectious...
    09-09-2024, 10:59 AM
  • seqadmin
    Addressing Off-Target Effects in CRISPR Technologies
    by seqadmin






    The first FDA-approved CRISPR-based therapy marked the transition of therapeutic gene editing from a dream to reality1. CRISPR technologies have streamlined gene editing, and CRISPR screens have become an important approach for identifying genes involved in disease processes2. This technique introduces targeted mutations across numerous genes, enabling large-scale identification of gene functions, interactions, and pathways3. Identifying the full range...
    08-27-2024, 04:44 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Today, 06:25 AM
0 responses
13 views
0 likes
Last Post seqadmin  
Started by seqadmin, Yesterday, 01:02 PM
0 responses
12 views
0 likes
Last Post seqadmin  
Started by seqadmin, 09-18-2024, 06:39 AM
0 responses
14 views
0 likes
Last Post seqadmin  
Started by seqadmin, 09-11-2024, 02:44 PM
0 responses
14 views
0 likes
Last Post seqadmin  
Working...
X