VCF_format, a point in REF/ALT columns
2
0
Entering edit mode
9.8 years ago
iraun 6.2k

Hi all,

Maybe this is a basic question but I am a bit confused about vcf-format... My question is, what's the meaning of a point ".", in the REF or ALT columns? Ex.:

#CHROM    POS    ID    REF    ALT    QUAL    FILTER    INFO
20     1230237 .         T      .       47   PASS   NS=3;DP=13;AA=T

Thanks in advance.

vcf-format • 6.4k views
ADD COMMENT
0
Entering edit mode

I can see a '.' in ID and FILTER, not in REF and ALT.

ADD REPLY
0
Entering edit mode

Oh sorry, my mistake, see my new edit.

ADD REPLY
5
Entering edit mode
9.8 years ago

from the SPEC: http://samtools.github.io/hts-specs/VCFv4.2.pdf

If there are no alternative alleles,then the missing value (.) should be used

(no variation was called here)

ADD COMMENT
0
Entering edit mode

So, there is a line in vcf file suggesting that a variant has been found in that position, but actually there is no variant?

ADD REPLY
3
Entering edit mode

Depending on how you called your variants, you may have output all positions to the VCF. This is the default (I believe) in GATK UnifiedGenotyper for instance, you have to explicitly tell it to only output variant positions.

ADD REPLY
0
Entering edit mode

This is correct and to add, this also happens in multi sample VCFs where a locus is successfully genotyped in one sample and lacks sufficient coverage in another sample.

ADD REPLY
0
Entering edit mode

Single base deletion?

ADD REPLY
0
Entering edit mode

deletion would be .eg: TA/T

ADD REPLY
0
Entering edit mode

I see, Thanks all for the explanations, now it is a little more clear. :)

ADD REPLY

Login before adding your answer.

Traffic: 2771 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