SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Inconsistent genotype in VCF file rlegendre Bioinformatics 1 04-16-2012 06:17 AM
vcf snp genotype M&M Bioinformatics 1 09-02-2011 10:52 AM
using vcf tools to extract genotype information rna_seeker Bioinformatics 3 07-10-2011 06:25 PM
From 1000 genome project, how can I get genotype information of European population? wwmm933 Bioinformatics 0 06-08-2011 08:27 AM
vcf: difference between position DP and genotype DP melo Bioinformatics 0 05-10-2011 08:50 AM

Reply
 
Thread Tools
Old 08-20-2012, 06:40 AM   #1
Karen Chait Berman
Junior Member
 
Location: Israel

Join Date: May 2012
Posts: 4
Default Incorrect Genotype information in VCF

Hi,

I am working with samtools and I get variants with genotypes that do not make sense. The DP4 shows high number of reads that match the reference G (101) and only 25 reads that match the Alt A (16,9) but the genotype is homozygote AA acording to the PL.

SL2.40ch05 14606139 . G A 151 . DP=129;VDB=0.0365;AF1=1;AC1=2;DP4=101,0,16,9;MQ=30;FQ=-28;PV4=1.2e-07,1,1,1 GT:PL:GQ 1/1:180,1,0:4

another example:
SL2.40ch05 14606550 . T A 125 . DP=186;VDB=0.0399;AF1=1;AC1=2;DP4=1,146,15,22;MQ=29;FQ=-53;PV4=3.3e-11,1,1,1 GT:PL:GQ 1/1:158,26,0:49

does anybody know how do get the correct genotype?
Karen Chait Berman is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 03:42 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO