Annotating SNPs from mapping locational information

[JackieBadger]

Hi,

We have a few thousand RAD SNPs which were characterized using a partially annotated genome.


SNP identifies look as follows. gb indicates the genome contig where the SNP is found

gi|354458801|gb|AGKD01000248.1|.0000029886

Can anyone indicate the easiest way to annotate these SNPs?
I imagine this would involve pulling useful annotation information from the genome where possible, and then estimating GO terms by BLASTing unannotated contains?

Thanks in advance