VEP chromosome name query
1
0
Entering edit mode
4 months ago
prasundutta87 ▴ 580

Hi,

I have been annotating my VCFs which has chromosome names such as chr1, chr2 etc. using a reference clinvar structural variation VCF which has chromosome names as 1, 2, etc. I am using the structural variation plugin of VEP.

Theoritically, VEP should throw an error due to wrong chromosome names, but it does not. Does VEP take care of the chromosome name by itself?

ensembl VCF annotation VEP • 486 views
ADD COMMENT
0
Entering edit mode

Are you using VEP command line? If yes, what is your command?

Side note: It's ensEMBL, not ensemble. There's no e at the end.

ADD REPLY
0
Entering edit mode

Oh! my bad..thanks for correcting me..

My command is:

vep -i input.vcf -o output.vcf --offline --dir_cache VEP_genomes --vcf --force_overwrite --plugin StructuralVariantOverlap,file=nstd102.GRCh38.variant_call.vcf.gz
ADD REPLY
1
Entering edit mode

It looks like the plugin is designed to check both with and without the chr prefix: https://github.com/Ensembl/VEP_plugins/commit/7c4c055bdb60bc02af460763cb83d48c62246947

If the output looks fine (you should manually test a few entries), you can rest assured that the prefix doesn't really matter.

ADD REPLY
0
Entering edit mode

Thanks a lot for sharing this! Although, I did rename the chromosomes just to be sure. But its good to know this. Thanks a lot again!

ADD REPLY
2
Entering edit mode
4 months ago

VEP is very clever, you're fine.

ADD COMMENT
0
Entering edit mode

Haha! I am sure of that, Emily! :-)

ADD REPLY

Login before adding your answer.

Traffic: 2562 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