We've moved!
This site is now read-only. You can find our new documentation site and support forum for posting questions here.
Be sure to read our welcome blog!

AD field

Hi,
Part of my GATK result is below. I don't understand why those called as 2/2 were not called as 1/1. I thought the second position in AD field is the depth for the first alternate allele (C in this case) and the third position in the AD field is the depth for the second alternate allele (G in this case). Do I misunderstand something?

1 745049 . A C,G 39258.66 PASS AC=29,240;AF=0.094,0.779;AN=308;BaseQRankSum=-3.240e-01;ClippingRankSum=0.509;DP=1690;FS=0.000;GQ_MEAN=33.06;GQ_STDDEV=22.58;InbreedingCoeff=0.8847;MLEAC=28,246;MLEAF=0.091,0.799;MQ=60.00;MQ0=0;MQRankSum=0.417;NCC=0;QD=26.96;ReadPosRankSum=0.509 GT:AD:DP:GQ:PL 1/2:0,2,2:4:57:135,57,78,84,0,78 0/0:13,0,0:13:7:0,7,317,7,317,317 2/2:0,0,16:16:48:464,513,765,48,51,0 2/2:0,17,0:17:51:495,495,495,51,51,0 1/2:0,1,8:9:18:311,281,333,42,0,18 2/2:0,6,0:6:18:170,170,170,18,18,0 2/2:0,14,0:14:42:398,398,398,42,42,0 2/2:0,12,0:12:36:327,327,327,36,36,0 2/2:0,8,0:8:24:205,205,205,24,24,0 2/2:0,17,0:17:51:480,480,480,51,51,0 2/2:0,10,0:10:30:287,287,287,30,30,0 0/0:11,0,0:11:20:0,20,348,20,348,3481/2:0,1,5:6:27:161,133,207,42,0,27 2/2:0,14,0:14:42:396,396,396,42,42,0 2/2:0,9,0:9:27:257,257,257,27,27,0 2/2:0,11,0:11:33:317,317,317,33,33,0 2/2:0,7,0:7:21:200,200,200,21,21,0 2/2:0,16,0:16:48:454,454,454,48,48,0 1/2:0,4,14:18:99:527,403,618,165,0,123 2/2:0,18,0:18:54:514,514,514,54,54,0 2/2:0,26,0:26:78:721,721,721,78,78,0 2/2:0,11,0:11:33:304,304,304,33,33,0 2/2:0,11,0:11:33:318,318,318,33,33,0
Thank you,
Ginger

Best Answer

Answers

Sign In or Register to comment.