Hello! I'm sorry if this is a duplicate question. I've searched the forum and didn't find the answer to what I seek.
I'm working with exome sequencing data and some of the variants are not completely annotated. I have the chromosome positions of each variant, but not the CDS or protein change associated with it.
My question is how to extract the CDS from the chromosome positions of each variant, to translate these variants into Synonymous/Non-Synonymous variations and then get the predicted protein effects of each variant. I thought of using UCSC Genome Browser for that. I there an easy way to extract that information using Python?
Example: I want to get from chr9:34371171 T>C to c.T1773C / p.I591T.