This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Protein differences in HGVS format

I have a bunch of protein sequences (~3.000) that I wish to compare with reference protein sequences, and I want the result in HGVS format. Example:

   ARNDCEQGHILKMFPSTWYV <=> ARNMCEQGHILKMFPSTYV => p.[Asp4Met; Trp18del]

I have written and rewritten routines for this, but with multiple deletions/insertions and/or frameshifts it gets complicated... What are my options?

Thanks.

protein alignment protein comparison hgvs
• 1,847 views • link
updated by Ram Ram • written by M. Möller • M. Möller

1 answer

Oh boy. This is ONE TOUGH task, I'm afraid - it's beyond what most tools are capable of. You might want to check the hgvs python module - if your protein sequence is a known reference protein, you might luck out with their VariantMapper or some such and get the appropriate notation for your mutant protein.

As far as I can see, the module does not support custom alignments, but you might be able to use this module to build your own tool.

• 0 views • link
written by Ram Ram

Log in to answer this question.

More posts like this