Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Inconsistent genotype in VCF file

    Hi everybody !

    I try to analyse my SNP calling results given by samtools and bcftools for my Illumina paired-end data.
    My first problem was identify the order in list of Phred-scaled genotype likelihoods (PL). In a old post, they say :

    "(...) A T,G (...) DP=35;AF1=1;CI95=0.5,1;DP4=0,0,1,8;MQ=60 PL:GT:GQ 96,47,70,35,0,70:1/1:72

    ref = A, alt = T,G right?

    So, genotype may be AA,AT,AG,TT,TG,GG"

    http://seqanswers.com/forums/showthread.php?t=9345&highlight=genotype+vcf

    but when I look my results in a genome browser, I find AA,AT,TT,AG,TG,GG for this example instead of AA,AT,AG,TT,TG,GG, unless I'm very much mistaken...
    So, if it's right, I don't understand how samtools/bcftools calculate individual genotype when there are two alternates alleles.. GT and PL fields are in conflict when they allocate a genotype. I give one example :


    Code:
    chrX	3506519		T	G,C 	DP=114;VDB=0.0324;AF1=0.6132;G3=0.4281,3.194e-06,0.5719;HWE=0.00569;AC1=12;DP4=9,15,43,17;MQ=60;FQ=999;PV4=0.0058,1,1,0.39 1/1:173,18,0,173,18,173:6:17    0/1:81,81,81,9,9,0:3:3  0/0:0,24,182,24,182,182:9:19    0/1:0,3,26,3,26,26:1:4  0/1:124,124,124,18,18,0:6:3

    So genotype may be TT, TG, GG, TC, GC, CC

    sample 1 = 1/1:173,18,0,173,18,173:6:17 => this sample is GG...
    sample 2 = 0/1:81,81,81,9,9,0:3:3 => this sample appears CC with the PL, and TG or TC with the GT... and it is GC in my viewer !!
    sample 3 = 0/0:0,24,182,24,182,182:9:19 => this sample is TT
    sample 4 = 0/1:0,3,26,3,26,26:1:4 => this sample appears TT in the PL and the genome browser, so homozygous for the reference allele, but the GT indicate a heterozygous genotype !
    sample 5 = 0/1:124,124,124,18,18,0:6:3 => this sample is CC


    I dont understand.. if somebody can help me, it would be fabulous !

    Thanks,

    Rachel

  • #2
    I also realized this behaviour of samtools / bcftools and have come across this post:

    http://www.biostars.org/post/show/12...uality-scores/.

    However, I still don't know how to get my genotypes instead :/

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Strategies for Sequencing Challenging Samples
      by seqadmin


      Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
      03-22-2024, 06:39 AM
    • seqadmin
      Techniques and Challenges in Conservation Genomics
      by seqadmin



      The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

      Avian Conservation
      Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
      03-08-2024, 10:41 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Yesterday, 06:37 PM
    0 responses
    11 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, Yesterday, 06:07 PM
    0 responses
    10 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-22-2024, 10:03 AM
    0 responses
    51 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-21-2024, 07:32 AM
    0 responses
    68 views
    0 likes
    Last Post seqadmin  
    Working...
    X