I want to count all SNPs from the upstream and downstream of a gene and I have a reference genome sequence and all upstream and downstream sequence of the genes. However, I don't know how to discover the SNPs from these regions.
If you have sequence alignments present then a simple R script would be able to to this. Otherwise you can generate the alignment and then write a script to identify bases that are different to the reference sequence