Rosalind第28题——ros_bio28_PROB

如果第一次阅读，请查看写在前面

import math

#读取文件
with open('../examples/ros_bio28_PROB.txt') as f:
    file = f.readlines()

dna = file[0].rstrip()
probability = file[1]
probability = probability.split(' ')

#根据提供的各个概率计算序列相加的对数
all_lo = []
for i in probability:
    i = float(i)
    gpro = i/2
    apro = (1-i)/2
    loga = 0
    for j in dna:
        if j == 'A' or j == 'T':
            loga += math.log(apro, 10)
        else:
            loga += math.log(gpro, 10)
    all_lo.append(loga)

#格式化为三位小数
result = ''
for line in all_lo:
    line = format(line, '.3f')
    result += line + ' '
print(result)