This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Some Non-Variant Positions Obtained When Using Bcftools With -V Option. Any Idea Why This Happen?

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 vcftools bcftools

Heads up: This is a 4 year old post.

0 answers

No answers yet.

Log in to answer this question.