Excess values in VCF INFO/AD field
I have a multisample VCF4.2 file in which variants are called in two different samples against a reference genome. For variant site, the INFO fields DP;ADF;ADR;AD;AC1;DP4 . are returned as below:
CHROM POS ID REF ALT QUAL FILTER INFO
chr1A 100 . G C 30 . DP=5;ADF=1,1,0;ADR=0,3,0;AD=1,4,0;AC1=3;DP4=1,0,1,3
I am however getting 3 values in the ADF, ADR, and AD fields instead of two values i.e ref and alternate allele depth as expected for biallelic sites. And this is the same for all my ~20,000 biallelic sites. Can anyone kindly explain why this might be.
Many thanks
Olu
• 3,295 views
•
link
written
by
oluwaseyi.shorinola •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Bcftools how to add SNP type to the INFO field (DEL,INS,INV,DUP..)
written by Michal Nevo •
-
How to correctly calculate Alternative Allele Frequency (AAF) from Samtools mpileup
written by Daniel •
-
Could bcftools output GQ for homozygous REF call?
written by wang.yiguan •
-
Mpileup Multiple Samples - Individual DP Values
written by Ks-seq •
-
Filtering for homozygous reference alleles
written by AGE •
-
bcftools call error "--gvcf output mode requires FORMAT/DP tag, which is not present in the input h…
written by AGE •
-
Estimating cross contamination in a set of BAMS
written by Pierre Lindenbaum
-
samtools v1.3: AD, ADF and ADR fields.
written by iraun
-
Some Non-Variant Positions Obtained When Using Bcftools With -V Option. Any Idea Why This Happen?
written by dapregi •
-
Vcf File From Mpileup ... What Next
written by Dataminer
Splitting mulit-allelic sites into multiple vcf records will not clean up excess allele fields.
did you discover why this happens? I have the same problem