std::string _get_mutated_aa(Tx & tx, std::string base, std::string codon, int intra_codon) { /** find the amino acid resulting from a base change to a codon @tx transcript object for a gene @base alternate base (e.g. 'G') to introduce @codon DNA sequence of a single codon @intra_codon position within the codon to be altered (0-based) @returns single character amino acid code translated from the altered codon. */ // figure out what the mutated codon is codon.replace(intra_codon, 1, base); return tx.translate(codon); }