When researchers discover a new protein, they would like to infer the strand of mRNA from which this protein could have been translated, thus allowing them to locate genes associated with this protein on the genome.
Unfortunately, although any RNA string can be translated into a unique protein string, reversing the process yields a huge number of possible RNA strings from a single protein string because most amino acids correspond to multiple RNA codons (see the RNA codon table).
Because of memory considerations, most data formats that are built into languages have upper bounds on how large an integer can be: in some versions of Python, a variable of type int may be required to be no larger than $$2^{31} - 1 = 2,147,483,647$$. As a result, to deal with very large numbers in, we need to devise a system that allows us to manipulate large numbers without actually having to store large numbers.
For positive integers $$a$$ and $$n$$, $$a$$ modulo $$n$$ (written $$a \mod n$$ in shorthand) is the remainder when $$a$$ is divided by $$n$$. For example, $$29\mod 11 = 7$$ because $$29 = 11 \times 2 + 7$$.
Modular arithmetic is the study of addition, subtraction, multiplication, and division with respect to the modulo operation. We say that $$a$$ and $$b$$ are congruent modulo $$n$$ if $$a\mod n = b \mod n$$. In this case, we use the notation $$a \equiv b\mod n$$.
Two useful facts in modular arithmetic are that if $$a \equiv b\mod n$$
and $$c \equiv d\mod n$$ then $$a + c \equiv b + d\mod n$$ and $$a \times
c \equiv b \times d\mod n$$. To check your understanding of these rules,
you may wish to verify these relationships for $$a = 29$$, $$b = 73$$, $$c
= 10$$, $$d = 32$$ and $$n = 11$$.
As you will see in this exercise, some assignments will ask for a (very large) integer solution modulo a smaller number to avoid the computational pitfalls that arise with storing such large numbers.
Write a function reverseTranslations that takes a protein string. The function must return total number of different RNA strings from which the protein could have been translated, modulo 1,000,000. Don't neglect the importance of the stop codon in protein translation.
>>> reverseTranslations('MA') 12 >>> from Bio import SeqIO >>> reverseTranslations(*SeqIO.parse('data.fna', 'fasta')) 115264