How does one calculate the information content of DNA sequence like ATCGGCT where mutation rate of G's is 10% and the most common mutation product binds with C and A with equal frequency.
I know that the individual information of a sequence is the dot product of the sequence and the weights matrix for each base. Essentially I = sum (base1xweight + base2xweight + ...)
What I can't figure out is how to determine the the weight of each base so as to incorporate the 10% probability of mutation of G's.. Can anyone help me figure that out
EDITED since wording may be confusing: Original problem says: " Spontaneous deamination of exocyclic amine is dG (deoxyguanosine) can occur occasionally in DNA strand. Given probability of this occurring during the time it takes DNA polymerase to copy DNA is 20% per dG in a particular DNA. Given that most common deamination product pairs with dC and dA with equal frequency calculate information content during DNA replication in the given sequence. " As an example I choose ATCGGCT and changed it to 10% deamination of dG.
PROPOSED SOLUTION: For sequence ATCGGCT, I am only consider the initial strand bases A,T,C are 2 bits a base, given by I=-1/ln(2)ln(n_before/n_after) ,where n_before = 1, and n_after = 4 In the case of mutated G (Gm) , I say, I = -1/ln(2)ln(1/3)= 1.58 bits , n_after=3 since A,C are equivalent to Gm
Perfect copy would give us 7 bases x 2 bits each = 14 bits content of strand. Since there are 2 G's, with 10% change of mutating --> 0.02 of G in strand will mutate.
Thus I content during replication, Information content I(ideal) - I(inc. mutation) = 14-0.03 = 13.97 bits during time of replication