所以我有一个8个字符长的序列TTTCGCGC。然后我有一个这样的字典:
d = {[A:prob1,prob2,... prob8],[T:prob1 ... prob8],...} 这给出了序列中每个可能位置的每个碱基的概率。
Ex:d ['C'] [4] =某个数字,表示C在序列中第5位的概率。
我想知道整个序列的概率。
这就是我所拥有的:
i=0
a=[]
while i<8:
if seq1[i]=='A':
a.append(d['A'][i])
i+8
print(a)
一旦我可以使用它,我将其扩展到所有四个核苷酸,但是我不确定这是怎么回事...请帮助:)