Effect type prediction for prokaryotes is relatively simple in comparison to eukaryotes because there is no splicing. UCSC genome browser doesn't help you with prokaryote genomes, but you do not need very sophisticated tools anyway.

You need:

• The genome reference sequence for alignment (you have it)
• The genome annotation to get the coordinates of CDS (simply as a .gff file)

You can find the genome annotations for prokaryotes on the NCBI ftp server For the yersinia pestis KIM10 strain this is e.g.: ftp://ftp.ncbi.nih.gov/genbank/genomes/Bacteria/Yersinia_pestis_KIM_10_uid288/AE009952.gff

Make sure the strain is consistent throughout your analysis.

Then:

1. Find the CDS from the gff file for each SNP coordinate if any.
2. extract the reference sequence of all CDS that contain a predicted SNP
3. translate the reference sequence
4. insert the alternative allele at the predicted SNP position
5. translate again, if the AA sequence is different then consequence is non-synonymous, and synonymous otherwise

Hope this helps.

You can try using Annovar, I guess it's a helpful tool where it can annotate variations against UCSC or Refseq gene

I'm not sure whether there is genome sequence for your organism just like Larry mentioned.

You can try annotating your SNPs or indels if your genome sequence is present in the UCSC genome browser to identify whether SNPs cause protein coding changes and the amino acids that are affected. Annovar is open-source:

http://www.openbioinformatics.org/annovar/annovar_faq.html

Just 22 SNPs?

You could do that manually. Get your SNP in context, with maybe 100 bases on either side, and submit it to ncbi's blastx

http://blast.ncbi.nlm.nih.gov/Blast.cgi?PROGRAM=blastx&BLAST_PROGRAMS=blastx&PAGE_TYPE=BlastSearch&SHOW_DEFAULTS=on&LINK_LOC=blasthome

It'll translate your sequence in all 6 frames, and figure out which protein submitted to NCBI best fits. If you have found a SNP whose protein has been submitted, that hit will come up first, so use some caution. But in general, if your SNP causes a non-synonymous change, you'll see it by looking at the top alignment.

Thank you for your answer.

My genome is y. pestis (CO92). I have identified some SNPs and their locations - I want to cross-validate whether the SNPs at these positions are synonymous or non-synonymous. Also, I have a reference sequence of this species as a fasta file. I am relatively new to SNP calling- and still in learning phase. I apologize if these questions are trivial !!!

Thanks.

Hi, listen, I'm pretty new on this blogosphere and Internet thing, so I don't know if there's a sort of "subscription" method that I can use [?]Generic Viagra[?] in order to receive notifications of your new entries...? Thing is I enjoy reading your blog a lot and I'd like to be up to date with your posts!

You can use very simple and useful program Procannot: http://www.platform.genexplain.com/procannot/

hi all of you I want to annotate the SNPs of my genome P.aeruginosa (free software, easy and for windows).