Good afternoon,
I have a list of 100 genes for which (in fact, for one of their transcripts in particular) I would like to get "synonymous coding" and "non-synonymous coding" SNPs that are observed in 1000G data (n=629).
Moreover, it would be fantastic to somehow extract the heterozygosity status for those SNPs.
I tried the ENSEMBLE 1000G browser, however, there are inconsistencies, that is, some SNPs that appear in the VCF file do not show up in the browser view. In addition, I do not want to mess with the dbSNP but am only interested in the SNPs observed in 1000G.
Any help would be much appreciated.
•
written
by
Sarah Tyrell •
1 answer
Sarah,
A pretty direct pipeline should allow you to do so.
-
Get the chromosomal coordinates of your genes
-
Extract the SNPs contained in the regions found in 1 (make sure that you use the same genome annotation version)
- Use a variant annotation software to annotate the SNPs (IE : ANNOVAR)
•
written
by
Simon P •
Log in to answer this question.
-
SNP count per gene, without disease mutations and observation bias
written by Suicyte •
-
What could account for these strange SNP discrepancies?
written by nchuang
-
Visualize Gene Descriptions On The Maps Of Ensembl Genome Browser
written by Anima Mundi
-
Dn/Ds Ratio - Paml : How Is It Possible To Obtain Negative Values For Dn Or Ds ?
written by Francois Olivier Hébert •
-
Conservation At Synonymous Sites
written by Sankalp Sharma •
-
How To Map Snps To Genes Based On Recombination Hotspots And Ld
written by Tafelplankje
-
Missing Snps For Hba1-Hba2 In 1000 Genome Data??
written by Tarbem •
-
Basic Queries Using Ucsc Table Browser
written by User 6659 •
-
How To Find The List Of All Coding Snps / Exomic Variants In A Given Gene ?
written by Khader Shameer
-
1000G And Dbsnp Build 132 In Ucsc Genome Browser?
written by Ryan D
Do you have the VCF file describing the 1000G variants that you want to use?
The inconsistencies you see are probably caused by the fact that there are different 1000genomes releases. In particular, they have published a new one in October 2011, including almost 2000 individuals (http://www.1000genomes.org/announcements/october-2011-integrated-variant-set-release-ichg2011-2011-10-12). Which release are you interested to?