This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Convert from Protein position to ENSEMBL Transcript's coding sequence position

Hi, Firstly, I am a Computer Science Engineer and a newbie in Bioinformatics.

I have Human mutation data in the following format: UniProt ID, Position, Original Amino Acid, New Amino Acid

I am trying to match this with ExAC data which contains: ENSEMBL Transcript ID, ..., CDS_Start, CDS_End, number of coding base pairs, ... etc

I have the UniProtKB ID to ENSEMBL Transcript ID mapping from one of the sites. But I am unable to map the positions.

How do I convert Protein Position to Transcript's CDS_start and CDS_end (beginning and end of Transcript's coding sequence positions?

I request your help.

Thank you, Rashmi

protein exac uniprot cds ensembl
• 3,735 views • link
updated by Emily Emily • written by rashmi.bangalore123 • rashmi.bangalore123

2 answers

The string 'original amino acid, position, new amino acid' looks like HGVS notation to me. You can use the Ensembl VEP for your variants and find out what their cDNA and CDS positions are. In addition the VEP will tell you whether the variants are observed in ExAC, 1000Genomes and ESP, plus clinical significance and plenty more.

• 1 views • link
written by Denise CS Denise CS

You can do this using the Transcript object from the Ensembl API. It has coding_region_start and coding_region_end methods. Or if that doesn't work for you, you can also use the TranscriptMapper object.

• 1 views • link
written by Jean-Karim Heriche Jean-Karim Heriche

Log in to answer this question.

More posts like this