Python：多个共识序列答案

【问题标题】：Python: Multiple Consensus sequencesPython：多个共识序列
【发布时间】：2016-07-26 10:01:46
【问题描述】：

从 dna 序列列表开始，作为回报，我必须得到所有可能的共识（结果在每个位置具有最高核苷酸频率的序列）序列。如果在某些位置核苷酸有相同的最高频率，我必须获得最高频率的所有可能组合。我还必须返回配置文件矩阵（每个序列的每个核苷酸频率的矩阵）。

这是我目前的代码（但它只返回一个共识序列）：

seqList = ['TTCAAGCT','TGGCAACT','TTGGATCT','TAGCAACC','TTGGAACT','ATGCCATT','ATGGCACT']
n = len(seqList[0])
profile = { 'T':[0]*n,'G':[0]*n ,'C':[0]*n,'A':[0]*n }

for seq in seqList:

    for i, char in enumerate(seq):
        profile[char][i] += 1



consensus = ""
for i in range(n):
    max_count = 0
    max_nt = 'x'
    for nt in "ACGT":
        if profile[nt][i] > max_count:
            max_count = profile[nt][i]
            max_nt = nt
    consensus += max_nt
print(consensus)
for key, value in profile.items():
     print(key,':', " ".join([str(x) for x in value] ))

TTGCAACT
C : 0 0 1 3 2 0 6 1
A : 2 1 0 1 5 5 0 0
G : 0 1 6 3 0 1 0 0
T : 5 5 0 0 0 1 1 6

（如你所见，在位置四，C和G得分相同，这意味着我必须获得两个共识序列）

是否可以修改此代码以获得所有可能的序列，或者你能解释一下逻辑（伪代码）如何获得正确的结果？

非常感谢您！

【问题讨论】：

目前您的代码只打印最后一个序列，而不是共识序列
谢谢你，Ohad，抱歉我的错误，我修复了我的代码，现在应该返回共识。（但只有一个，所以主要问题仍然存在:)。关于您的答案，输出必须显示整个两个序列。
为你解决了这个问题

标签： python bioinformatics rosalind

【解决方案1】：

我确信有更好的方法，但这是一个简单的方法：

bestseqs = [[]]
for i in range(n):
    d = {N:profile[N][i] for N in ['T','G','C','A']}
    m = max(d.values())
    l = [N for N in ['T','G','C','A'] if d[N] == m]
    bestseqs = [ s+[N] for N in l for s in bestseqs ]

for s in bestseqs:
    print(''.join(s))

# output:
ATGGAACT
ATGCAACT

【讨论】：

@3lli0t 欢迎您。你应该接受答案。