如何为给定的RNA序列生成具有X个突变数的所有突变体

Question

我想从具有x突变的DNA/RNA序列中找到所有可能的突变子。

我制作了一个代码，用于获得所有具有1个突变的突变体，但我不知道如何将其概括为x个突变数。

def generateAllMutan(baseSeq):
base= ["A", "T", "C", "G"]
mutantList = []
for i in range(len(baseSeq)):
for nuc in base:
if nuc != baseSeq[i]:
seqName = "name" + str(i+1) + "_" + baseSeq[i] + "_to_" + nuc
mutedSeq = baseSeq[0:i] + nuc + baseSeq[i++1:]
mutantList.append(SeqRecord(Seq(mutedSeq), id= seqName, name="", description=""))
return     mutantList

我希望这个函数能处理任意数量的欲望突变如def generateAllMutan(baseSeq, nbMutation):

但我不知道如何有效地做到这一点。

知道吗？

Answer 1

尝试与指定数量的匹配进行比较：n

def function(s, n):
for x in range(len(s)):
for a in b:
if a != s[x] and count <= n:
#code
count +=1