How can I revert back the Ensembl vep annotated positions to vcf positions?

Sample VCF positions:

CHROM   POS     ID  REF  ALT
chr1    939570   .  T   TCCCTGGAGGACC
chr1    3411794  .  T   C
chr1    1043223  .  CCT C
chr1    22511093 .  CTT CT
chr1    24454031 .  CT  CTT

VEP annotted positions

Positions            Allele
1:939570-939570      CCCTGGAGGACC
1:3411794-3411794        C
1:1043223-1043225        -
1:22511093-22511095      -
1:22511093-22511095      T
1:24454031-24454032      -
1:24454031-24454032      TT

I wanted to revert back VEP positions to VCF positions to remap certain parameters in the VCF file.

Any help appreciated!!!

ensembl vep vcf

• 1,812 views • link

updated 5.3 years ago by ATpoint

• written 5.3 years ago by DareDevil

Hello,

can you please describe what your goal is? The Position in the vcf file describes the starting position of the variant, the VEP Position describes the start and end position of the variant. If you transfer this range into your vcf, this is no valid vcf anymore.

fin swimmer

• 0 views • link

written 5.3 years ago by finswimmer

1 answer

Just get your output in VCF. If using the offline tool, use --vcf. If using the online tool, use the Export options. The VEP annotation will appear as CSQ in the INFO field. If offline, it will also keep all the other stuff from your original VCF.

• 0 views • link

written 5.3 years ago by Emily

More posts like this

BCFtools isec output vs BCFtools query
written 2.9 years ago by gtho123

I have 5 samples and I have performed variant calling using the `Strelka2` software. The output of this software is a `VCF`file for each sample …
Filtering of tricky overlapping sites in VCF
written 3.4 years ago by dr_cfoster

I'm working on a project where I'm working with variants ranging from low frequency (AF~=0.03) up to high frequency (AF=1) after variant calling with `lofreq`. …
how to seperate VEP INFO column into seperate columns
written 3.5 years ago by minoo •

I have a vcf files like below: ``` #CHROM POS ID REF ALT QUAL FILTER INFO FORMAT treatmentSample chr1 857100 . C T 1756.06 PASS …
bcftools filter by protein prediction
written 5.2 years ago by storm1907 •

Hello, I have VEP annotated vcf files with following content: #CHROM POS ID REF ALT QUAL FILTER INFO FORMAT file chr1 183937 . G A …
htsjdk.tribble.TribbleException: The provided VCF file is malformed at approximately line number 92…
written 5.4 years ago by DareDevil

I run variantsToTable module of gatk as follows: gatk VariantsToTable -V 001.vcf -F CHROM -F POS -F ID -F REF -F ALT -F QUAL -F …
VEP output has no gene names
written 8.0 years ago by Gene_MMP8

I am trying to annotate a variant file(generated using strelka) from mice WGS data. This is the command I used: ./vep -i /path/to/somatic.snvs.vcf \ --cache …
NCBI Genome Remapping Service- clinical remap
written 8.1 years ago by Sudhir Jadhao

Dear All, I am converting LRG coding sequence co-ordinate to hg19 genomic. In "NCBI Genome Remapping Service- clinical remap " gives two types of output …
Grep Position List
written 9.2 years ago by BAGeno

Hi, I have VCF file and list of position file. I want to grep those records from VCF which have position corresponding to position present …
Extracting Specific Reads From Vcf File
written 14.1 years ago by Bnfoguy •

Hello, I am trying to extract indels with coverage ~1 and they are in a VCF file. How do I write a python script or …
Joining Snpeff And Vcf
written 14.5 years ago by Pierre Lindenbaum

<p>I've extracted the distinct mutations from a set of <strong>VCF</strong> files:</p> <pre><code>#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT SAMPLE chr1 66480 . AT …