PL field in gvcf file equal to zero
0
0
Entering edit mode
19 months ago
Negin ▴ 10

I have a genomics data from Bacteria (E. Coli.) and used GATK 4 for variant calling. I have noticed that for some samples and some loci, PL is 0,0. Can anyone please explain this?

eg. NC_000913 4473651 . C <non_ref> . . . GT:AD:DP:GQ:PL 0:2,0:2:0:0,0

How can we interpret that? becuase in this case, the GT in later steps will be "." as it will be driven from PL, when we combine and genotype the samples.

thanks in advance.

GATK BP_RESOLUTION GVCF VCF • 423 views
ADD COMMENT
1
Entering edit mode

there is no called genotype here: DP is 2, GQ is 0. this variant could be a node between two statistically different regions, or a variant with an ugly environment (clipping...)

ADD REPLY
1
Entering edit mode

Hi Pierre, as there is no called GT at this loci, in only some of the samples, could that means absence of a gene in those samples? Regards

ADD REPLY
0
Entering edit mode

Thanks for your answer, but could you please explain your answer more as I am quite new in this field and didn't understand you totally!?

ADD REPLY

Login before adding your answer.

Traffic: 1589 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6