使用 for 循环替换 DNA 序列中的坏核苷酸

Question

我有一个序列列表（为简单起见，如下所示）

seqList=["ACCTGCCSSSTTTCCT","ACCTGCCFFFTTTCCT"]

我想使用 for 循环将除 ["A","C","G","T"] 以外的每个核苷酸实例替换为 "N"

到目前为止我的代码

seqList=["ACCTGCCSSSTTTCCT","ACCTGCCFFFTTTCCT"]
for x in range(len(seqList)):
    for i in range(len(seqList[x])):
        if seqList[x][i] not in ["A","C","G","T"]:
            seqList[x][i].replace(seqList[x][i],"N")
            print(seqList)

问题是，核苷酸没有被替换，原始序列没有任何变化，我不知道原因！！！

Answer 1

python 中的字符串是不可变的。 你可以让ot像这样工作

seqList=  ["ACCTGCCSSSTTTCCT","ACCTGCCFFFTTTCCT"]
for x in range(len(seqList)):
    stringl=list(seqList[x])
    for i in range(len(seqList[x])):
        if seqList[x][i] not in ["A","C","G","T"]:
            stringl[i].m="N"
    seqList[x]="".join(stringl)

Answer 2

不循环所有字母的方法将替换所有不是ACGT的字母

def replace_bad(seq):
    unique = [
        letter
        for letter in set(seq)
        if letter not in "ACGT"
    ]

    for each in unique:
        seq = seq.replace(each, "N")

    return seq


if __name__ == '__main__':
    for seq in seqList:
        print(replace_bad(seq))