My objective is to find the frequency of in-frame codons within 50 bp of a specified location. I have the ensembl transcript ID and the genomic coordinates,e.g.:
chromosome position strand ensemblTID
chr2 219130603 - ENST00000538028
I have been trying to use the getBM and getSequence functions, but I'm not even close to getting what I want.
Any help to point me in the correct direction would be appreciated!