Header Leaderboard Ad

Collapse

Reading a vcf - genotypes

Collapse

Announcement

Collapse

SEQanswers June Challenge Has Begun!

The competition has begun! We're giving away a $50 Amazon gift card to the member who answers the most questions on our site during the month. We want to encourage our community members to share their knowledge and help each other out by answering questions related to sequencing technologies, genomics, and bioinformatics. The competition is open to all members of the site, and the winner will be announced at the beginning of July. Best of luck!

For a list of the official rules, visit (https://www.seqanswers.com/forum/sit...wledge-and-win)
See more
See less
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Reading a vcf - genotypes

    Hi, I am trying to understand the genotype calls in my vcf, a trio.
    Looking at one variant of interest, the calls are:
    Dad Mum Kid
    0/2 2/1 1/1

    How can the kid be 1/1?
    The full information is:

    FORMAT
    GT:ADP:GQ:PL

    Dad
    0/2:9,3,0:16:19.63:134,0,98,27,20,47

    Mum
    145 2/1:8,4,0:20:15.38:197,15,52,44,0,41

    Child
    1/1:7,2,0:11:1.37:113,31,61,9,1,0


    Thanks for any help.

  • #2
    Note the GQ value, which is quite low for the child. You might consider filtering by quality.

    You will run into this sort of thing on occasion. The general causes are (1) sequencing error (2) genotype call error (3) sample swap (this should be pretty obvious though) (4) non-paternity. The last one is more common than you would think (1-3% of cases where there's otherwise no dispute of paternity), though you'll again end up with a lot of sites like this.

    Comment

    Latest Articles

    Collapse

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 06-01-2023, 08:56 PM
    0 responses
    9 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 06-01-2023, 07:33 AM
    0 responses
    9 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 05-31-2023, 07:50 AM
    0 responses
    4 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 05-26-2023, 09:22 AM
    0 responses
    11 views
    0 likes
    Last Post seqadmin  
    Working...
    X