This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Calculate allelic frequency from VEP output vcf file

I am calculating allelic frequency for somatic variants. VEP gives a VCF file, for example that has:

9761707 G>A and outputs

DP FDP SDP SUBDP  AU  CU    GU    TU
99   0  0   0    0,0  0,0   99,99  0,0

or

121467604 T > C

DP FDP SDP SUBDP  AU    CU   GU    TU
153  0  0   0    0,0  77,79  0,0  76,76

How can I calculate the allele frequency from this information? I don't quite understand why in the first case there is nothing on AU, but the second case has something on both CU and TU. Based on my readings (recommendation here) it should be:

AU1/(AU1+GU1) for the first one and CU1/(CU1+TU1) for the second one. Yet, on the first case this would mean an allelic frequency of zero, but then why is the variant being reported at all??

Why is this?

vep variant
• 1,035 views • link
written by ramiro.barrantes • ramiro.barrantes

1 answer

Found the answer here!!

https://github.com/Illumina/strelka/blob/master/docs/userGuide/README.md#somatic

Although this leads to a strange situation in which some variants might have zero for allele frequency, which has been commented on before: https://github.com/Illumina/strelka/issues/65

• 0 views • link
written by ramiro.barrantes • ramiro.barrantes

Log in to answer this question.

More posts like this
  • Manta: What does all this mean
    written by josh.mannheimer • josh.mannheimer
    Hi, I am trying to understand these Manta somatic variant VCFs that I have annotated with VEP. chr10 11974892 MantaBND:6118:0:1:0:0:0:0 G ]chr1:41850074]G . PASS SVTYPE=BND;MATEID=MantaBND:6118:0:1:0:0:0:1;CIPOS=0,2;HOMLEN=2;HOMSEQ=AG;SOMATIC;SOMATICSCORE=54;BND_DEPTH=62;MATE_BND_DEPTH=0;CHR2=chr1;ENDPOSSV=41850074;SVLEN=0;RO=-+;TRA;LENSV=-1;RGENUPS=ENSCAFG00000045430___None,-9947,+;RGENISEC=ENSCAFG00000000435___BEST3,0,-;RGENDNS=ENSCAFG00000000436___RAB3IP,59491,+;LGENUPS=ENSCAFG00000000405___PLEKHG1,-146330,+;LGENISEC=ENSCAFG00000000413___MTHFD1L,0,+;LGENDNS=ENSCAFG00000000419___AKAP12,173336,+;CSQ=]chr1:41850074]G|feature_truncation|HIGH|MTHFD1L|MTHFD1L_3|Transcript|XM_038526323.1|mRNA||||||||||||1||||Canis_familiaris.UU_Cfam_GSD_1.0.gtf.gz|,]chr1:41850074]G|feature_truncation&intron_variant|HIGH|BEST3|BEST3_2|Transcript|XM_038549568.1|mRNA||4/5||||||||||-1||||Canis_familiaris.UU_Cfam_GSD_1.0.gtf.gz| …
  • BSSNPer vs. BisSNP
    written by aroso491 • aroso491
    Hello everyone, I am trying to conduct a methylation analysis on some bisulfite-seq data. I have to follow (as much as I can) the procedure …
  • Problem with VCF format while using "VariantAnnotation" R package
    written by pixie@bioinfo pixie@bioinfo
    Hello, I am using VariantAnnotation_1.40.0 to read in a VCF file using vcf <- readVcf("220106_MN01111_0064_v2_A000H3MT5V.tomato.filtered.txt") However, I get this error: Error in h(simpleError(msg, call)) : …
  • What is the correct way of setting the genotype after splitting multi-allelic sites in a VCF file?
    written by terdon terdon
    Consider this minimal VCF file which contains a multi-allelic call made by a somatic variant caller on a human tumor sample: ``` ##fileformat=VCFv4.3 ##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype"> ##FORMAT=<ID=AD,Number=R,Type=Integer,Description="Allelic …
  • Help understanding bcftools call genotyping decision
    written by Fran Rodriguez-Algarra • Fran Rodriguez-Algarra
    Hi everyone! I'm new here, so don't hesitate to let me know if I need to improve my post in anyway. Although I know that …
  • Tassel5v2 vcf file with 'N' in reference and alternate allele
    written by Tm Tm
    Hello All, Can anyone please help me understand why VCF file obtained using tassel5v2 pipeline "DiscooverySNPcallerV2" is giving 'N' in reference as well as alternate …
  • Variant calling on bacterial population using Freebayes gave wrong allele frequency and extremely l…
    written by xuye.july • xuye.july
    I am working on genomic data of 8 evolving bacterial populations. Each one of them is a community of millions non-clonal individuals and was pooled …
  • Strelka output: read count of REF and ALT, SNVs
    written by noeD noeD
    Hello! I have used Strelka as variant caller for SNV. I'm having trouble with the output... I would like to obtained the read counts of …
  • Good filtering practices for VCF (Bcftools 1.3)
    written by umn_bist • umn_bist
    I have a tumor sample without matches. This sample is most likely impure without a certain ploidy, copy #. I am wondering for such samples …
  • Need help understanding read depth for somaticsniper vcf's
    written by Jordan Jordan
    Hi, I'm having trouble understanding read depth in vcf's called by `somatic-sniper`. Here is a sample vcf file: ``` ##fileformat=VCFv4.1 ##fileDate=20140428 ##phasing=none ##reference=human_hg19.fa ##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype"> ##FORMAT=<ID=IGT,Number=1,Type=String,Description="Genotype …