I have a list of SNPs that consists of their rsids or SNP names, like rs6699871, for example. From that list, I would like to retrieve a list of gene variations that result from that SNP. I would expect something like this as output:
rs6699871 var1 CTACATCGATTTGCAGCACCCAGCTGA[C]CCAGAAATCGACAAGTCGACCCGTGCTAGATCGACGA <- (this is the COMPLETE gene or ORF sequence, not the X bp flanking sequence) rs6699871 var2 CTACATCGATTTGCAGCACCCAGCTGA[T]CCAGAAATCGACAAGTCGACCCGTGCTAGATCGACGA rs6699871 var3 CTACATCGATTTGCAGCACCCAGCTGA[G]CCAGAAATCGACAAGTCGACCCGTGCTAGATCGACGA
where the [X] refers to the location of the SNP and it's different forms (C, T, or G). Each line represents a variation in the sequence of a given gene resulting from the SNP rs6699871. This variation could be either in Introns or Exons.
So my question is: what database would have this SNP versus gene variations map that I'm looking for, or what database could have enough information on SNPs that could allow me to retrieve this output easily?
I have experience in linux terminal, perl, and c++, so if I know where to look I can make a script to retrieve this data automatically for a large list of SNPs.