Question: Some Non-Variant Positions Obtained When Using Bcftools With -V Option. Any Idea Why This Happen?
gravatar for dapregi
5.2 years ago by
dapregi50 wrote:

Hello, I am wondering why, if I have used the -v option in bcftools (output only variant sites), I have in the output some reference homozygotes. Could it be because the genotype likelihood is too low that could be also a het? Does the -v option look at the most probable genotype (the one written down) or at the genotype likelihoods?

bcftools view -Svc myfile.vcf.gz

#CHROM  POS     ID      REF     ALT     QUAL    FILTER  INFO    FORMAT  MET_03912_4_1_TAGCTT
1       15906   .       A       .       3.55    .       DP=3;QS=0.000000,0.906977,0.093023,0.000000;VDB=4.034274e-02;AF1=1;AC1=2;DP4=0,0,0,3;MQ=16;FQ=-33;AF1=1;AC1=2;FQ=-30    PL:DP   31,3,0:3
vcf bcftools vcftools • 1.6k views
ADD COMMENTlink modified 8 months ago by Biostar ♦♦ 20 • written 5.2 years ago by dapregi50

Heads up: This is a 4 year old post.

ADD REPLYlink written 8 months ago by RamRS21k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1914 users visited in the last hour