This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Retrieve mutation position and ID for a mutation in hgvs format

How can I find a mutation's chr pos and id represented in HGVS format?

Eg:

Gene: TMEM231    cdna_Change: NM_001077418.1:c.582+3A>G    protein_change: p.?

The mutations are represented in HGVS format. How and where can I find the rs#, chr and pos for this particular mutation.

I have a set of 10000 mutations and would like to annotate then with their chr, pos and rs#

mutation hgvs annotation dbsnp
• 10,242 views • link
updated by Ram Ram • written by vigprasud vigprasud

What programming languages do you know? This could be done in R (and presumably biopython/bioperl) relatively easily.

• 0 views • link
written by Devon Ryan Devon Ryan

I know python and R

• 0 views • link
written by vigprasud vigprasud

If VEP doesn't work for you, then you can do this in R. The general steps would be to:

  1. Load this file as a dataframe and parse the cdna information to split the ID from the position information.
  2. Load a txdb that contains these IDs (they don't all).
  3. You can then just apply a function to each transcript to calculate the cDNA position of each exon (you'd just use the 5' or 3' most coordinate).
  4. Now you have numbers you can compare, so you'll need to apply a function to extract the appropriate transcript and then just determine (A) which exon it would be in (or intron following an exon as in your example) and then (B) increment/decrement the genomic position of said exon by the appropriate offset.
• 0 views • link
written by Devon Ryan Devon Ryan

3 answers

I would try VEP

• 0 views • link
updated by Ram Ram • written by Jeremy Leipzig Jeremy Leipzig

I second VEP, it works quite well for this. I forget what the limit is for the number of variants through the online web interface but you can either do it that way in batches or do it through the command-line version. You just have to switch from the default and you can put in HGVS mutations using RefSeq sequences

• 0 views • link
written by DG DG

Thank you. That helps.

• 0 views • link
written by vigprasud vigprasud

Also consider the Python hgvs package. [Disclosure: I'm one of the authors.]

• 0 views • link
updated by Ram Ram • written by Reece Reece

Well as a non-python user I didn't like this tool initially but given it's well-written documentation, I was able to follow. Posting an example just in case if someone like me is struggling.

I had a NC IG as follows: "NC_000002.11:g.113890610C>T" and I wanted NP IDs for the same

## Initializing hgvs shell
hgvs-shell
 var_g = parse("NC_000002.11:g.113890610C>T") 
transcripts = am37.relevant_transcripts(var_g)

In [27]: for ac in sorted(transcripts):
    ...:     var_t = g_to_t(var_g, ac)
    ...:     var_p = t_to_p(var_t)
    ...:     print("-> " + str(var_t) + " (" + str(var_p) + ") ")

This returned all the NP IDS

• 0 views • link
written by rohitsatyam102 rohitsatyam102

The lesser mentioned yet totally awesome webpage

• 0 views • link
updated by Ram Ram • written by Vivek Vivek

I see that mutalyzer works for converting rsIDs to HGvs but not the other way around.

• 0 views • link
written by Ritika • Ritika

Log in to answer this question.

More posts like this