【Python】求根据给定GC比率生成的DNA字符串与给定字符串完全匹配的概率

题目要求

解决方法

def random_genome(dna, gc_content):
    log_probabilities = []
    for gc_pro in gc_content:
        GC_frequency = gc_pro/2
        AT_frequency = (1 - gc_pro)/2
        probabilities = 1
        for ch in dna:
            if ch == 'A' or ch == 'T':
                probabilities *= AT_frequency
            else:
                probabilities *= GC_frequency

        log_probabilities.append(round(math.log10(probabilities), 3))

    return log_probabilities

#test
print(random_genome("ACGATACAA",[0.129, 0.287, 0.423, 0.476, 0.641, 0.742, 0.783]))

 

输出结果

 

成功！！！