by using Python
*According to junkes cantor evolution model
For this question, generate a random transition rate matrix ? to model the rates of nucleic acid substitutions in such a way that the off-diagonal entries vary between 0 and 1 and the diagonal entries are set so that each row sums to 0.
*Calculate the transition probability matrix ?(?) corresponding to the transition rate matrix ? above and plot the corresponding transition probabilities between each nucleotide pair in a 4 × 4 table for ? ∈ [0,5].
(Hint: You need to figure out how to calculate matrix exponentials.)