Annotating variants using latest version fof the genome
2
0
Entering edit mode
4.1 years ago
prasundutta87 ▴ 670

Hi,

I have SNP data from an array which has coordinate information from an older version of the genome (in VCF form). I want to annotate it using GTF/GFF file of the current version of the genome using SnpEff or ANNOVAR or VEP. Do I need to take care of the VCF coordinates somehow (like liftover) so that the variant annotation program can annotate the variants correctly?

Or, the variant annotation program will only annotate the variants with matching coordinates and incorrectly annotate the rest whose coordinates don't match?

I am assuimg its the second, but want other people's opinion on this.

Regards, Prasun

SNP snp VCF • 984 views
ADD COMMENT
3
Entering edit mode
4.1 years ago
GenoMax 147k

You could use crossmap (LINK) to first convert the co-ordinates over to the version of the genome you want to use and then annotate.

ADD COMMENT
0
Entering edit mode

Thanks all..I believe that I have got an answer to my query and I am thinking in the right direction!

ADD REPLY
0
Entering edit mode
4.1 years ago
brunobsouzaa ▴ 830

Just to check if I understood, you have SNP data with hg19 coordinates and you want to annotate with hg38 coordinates? Don't do this, use the same genome version to annotate your files.

ADD COMMENT
0
Entering edit mode

As an hypothetical example, yes, that's exactly what I am asking and I am hoping that the coordinates will all screw up giving all sort of wrong annotations..

ADD REPLY
1
Entering edit mode

If you use different genome versions yes, it will give you wrong annotations!

ADD REPLY

Login before adding your answer.

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