# How can I get the HGVS notation of a mutation from a FASTA?

In an old article, they gave a table of mutations found in a gene and its positions which points to another article containing the FASTA sequence.

This FASTA is the precursor sequence of the gene. How can I know which transcripts those mutations belong to? And how can I obtain their HGVS format?

For example, C-7111-G should have an HGVS notation like c.XXXC>G.

• What's the gene? There are tools you can try, but we'd need to know the gene first. And no matter what you do, this will be a manual process since the old FASTA sequence you have is likely inaccurate anyway.
– terdon
Oct 7 at 9:52
• The gene described in the article is POMC. What I did is the following: first, I copied the same fasta sequence from the old article, then I pinpointed the exact location of the mutations from the table. Second, I copied a couple of unique nucleotides from the left and right of the mutations. Third, I matched those sequences from the actual gene fasta sequence to locate the real nucleotide, and I got its position in the current hg38 build. Then the rest is simple, I searched in the databases what mutations are found in those emplacements. Oct 7 at 15:20
• Yeah, that's probably the best thing you can do. I was going to shamelessly plug the product of the company I work for, varsome.com, which is free and can help you translate variant notations, but if all you have is the position relative to an old fasta sequence and don't know the actual transcript, we can't really help much more. You'll be able to get the proper HGVS notation easily though. Same for variantvalidator.org with which I am not affiliated and which is a great tool for validating HGVS.
– terdon
Oct 7 at 15:36
• Thank you @terdon. I will check those websites. Cheers! Oct 7 at 15:42