Variant Effect Predictor (VEP) warning on insertions: coordinates are not (start = end + 1)
0
0
Entering edit mode
6 days ago
beketamas • 0

I'm trying to run Variant Effect Predictor (VEP) on TCGA VCF(-like) files that were converted from MAFs. I get the following warnings on insertions: WARNING: Alleles look like an insertion (-/GCAAT) but coordinates are not start = end + 1 (START=56620027, END=56620027) for variant (.) on line 19.

Therefore, the insertions are not present in the output files...

My VCF file looks like the following:

VCF-like file

Columns are Chromosome, Start position, ID, Ref. seq, Alt. seq, VAF.

The problem is with the Start position column but I couldn't figure out how to modify it to work. If I'm correct, I should specify the End position as well, but I got no luck in finding the correct format that VEP accepts. Can someone please help me with how to resolve it? Thanks!

VCF VEP • 156 views
ADD COMMENT
0
Entering edit mode

why a screenshot when you can just copy+paste the data ?

My VCF file looks like the following

'-' is not valid: this is not how a indel should be defined in a vcf. https://samtools.github.io/hts-specs/VCFv4.2.pdf

For simple insertions and deletions in which either the REF or one of the ALT alleles would otherwise be null/empty, the REF and ALT Strings must include the base before the event (which must be reflected in the POS field)

ADD REPLY

Login before adding your answer.

Traffic: 684 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6