PL field in gvcf file equal to zero
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.
• 1,622 views
•
link
written
by
Researcher •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Cannot use more than 2 sample IDs when demultiplexing using Demuxlet
written by Chilly •
-
My VCF file structure is strange and less compatible to data processing programs.
written by kgwkk2 •
-
Problem with VCF format while using "VariantAnnotation" R package
written by pixie@bioinfo
-
Tassel5v2 vcf file with 'N' in reference and alternate allele
written by Tm
-
Why invariant blocks in GATK consistently have very low quality scores (but not variant sites)
written by yixianshuiesuan •
-
How do genotype likelihoods (PL) work in gVCF files?
written by popantrop •
-
GATK GenotypeGVCFs output is empty gvcf file
written by GK1610
-
Variant calling: GT field disagrees with DP4 field
written by al.bodrug
-
Creation of a large Panel of Normals for Mutect2 analysis
written by filippo.martignano •
-
VCF file PL values for more than 1 alternatives allele
written by bharata1803
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...)
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!?
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