This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Wrong ZYGOSITY by ensembl vep on DeepVariant vcf file

Hello,
I'm annotating a deepvariant vcf file using vep ensembl version 110.1. however although the allele depth shows obvious heterozygous variant the ZYG output of vep when using --individual all shows HOM.
Here is the vcf (here I just put one variant) :

##fileformat=VCFv4.2
##FILTER=<ID=PASS,Description="All filters passed">
##contig=<ID=chr1,length=248956422,assembly=38>
##contig=<ID=chr2,length=242193529,assembly=38>
##contig=<ID=chr3,length=198295559,assembly=38>
##contig=<ID=chr4,length=190214555,assembly=38>
##contig=<ID=chr5,length=181538259,assembly=38>
##contig=<ID=chr6,length=170805979,assembly=38>
##contig=<ID=chr7,length=159345973,assembly=38>
##contig=<ID=chr8,length=145138636,assembly=38>
##contig=<ID=chr9,length=138394717,assembly=38>
##contig=<ID=chr10,length=133797422,assembly=38>
##contig=<ID=chr11,length=135086622,assembly=38>
##contig=<ID=chr12,length=133275309,assembly=38>
##contig=<ID=chr13,length=114364328,assembly=38>
##contig=<ID=chr14,length=107043718,assembly=38>
##contig=<ID=chr15,length=101991189,assembly=38>
##contig=<ID=chr16,length=90338345,assembly=38>
##contig=<ID=chr17,length=83257441,assembly=38>
##contig=<ID=chr18,length=80373285,assembly=38>
##contig=<ID=chr19,length=58617616,assembly=38>
##contig=<ID=chr20,length=64444167,assembly=38>
##contig=<ID=chr21,length=46709983,assembly=38>
##contig=<ID=chr22,length=50818468,assembly=38>
##contig=<ID=chrX,length=156040895,assembly=38>
##contig=<ID=chrY,length=57227415,assembly=38>
##contig=<ID=chrM,length=16569,assembly=38>
#CHROM  POS     ID      REF     ALT     QUAL    FILTER  INFO    FORMAT  SAMPLE
chr3    195784506       .       A       AAG     18.6    PASS    .       GT:GQ:DP:AD:VAF:PL      1/1:12:45:14,27:0.6:18,13,0

The vep command I used:

 vep \
    --verbose --cache --offline --species homo_sapiens --assembly GRCh38 \
    --use_given_ref \
    --tab \
    --force_overwrite \
    --dir /opt/vep/.vep \
    --input_file /opt/vep/files/${inputVcf_file} \
    --output_file /opt/vep/files/${output_file} \
    --fasta /opt/vep/.vep/custom/references/Homo_sapiens_assembly38.fasta \
    --allele_number \
    --individual all \
    --per_gene

Why is this happening?
Thanks for any help in Advance

annotation vcf vep zygosity deepvariant

I found out that vep is only checking the GT field therefore assigning HOM to the Zygosity. Therefore the problem should be from deepvariant side

0 answers

No answers yet.

Log in to answer this question.