Annotating 1000 Genomes data

[ashkot]

Hi all,
I am working with 1KGenomes data abd I have conducted all the analysis on my files of interest.

What I am really after is to annotate my data, I have been using annovar, which works great, however i have come accross an issue.

First I used dbSNP 129 and hg18 for annotating, this time most variants were not annotated. I was then advised to use buildver 37 and hg19 and i used dbSNP 132. This time around most variants where annotated, but my variant was interest was not. Thinking it may not be annotated i looked at the filtered file and even the genomic coordinates for that variant did not appear in that.

Can someone clearly explain with assembly versions, dbSNP builds etc.

Happy Holidays to all !

Ashwin

[laura]

Askhot

Are you looking for functional annotation of your variants?

Do you know what assembly version your variants have their coordinates in?

thanks

[ashkot]

hi, the dataset that i am working on is one of the .bam files that i have downloaded from 1K Genomes.

As for annotation, all i really want is to have dbSNP rs numbers in the VCF files. I did a lot of downstream analysis such as converting the .bam into bcf and then piping it into annovar to generate an annotated vcf file. however, this vcf file does not have any rs id's listed.

the coordinate system that I want to use is GR37. i am just wondering what is the most straight forward way to accomplish this?

thanks,
ashwin

[laura]

You will probably find the dbSNP vcf files useful

ftp://ftp.ncbi.nih.gov/snp/organisms...9606/VCF/v4.0/

[ashkot]

Thanks, i think this might do the trick.