如何用Python生成指定碱基量的随机DNA序列、查找CpG岛并打乱序列？

一、修正现有代码：打乱列表中每个DNA序列的碱基顺序 #

你的代码问题在于random.sample()返回了打乱后的字符列表，但没有保存结果，且直接遍历列表元素无法修改原列表。以下是两种可行的修正方案：

方案1：生成新的打乱后序列列表 #

import random

finalList = ['AGTCTTTGTGT', 'ACGCGTAAACG', 'AGGGTCTC']
shuffled_list = []

for seq in finalList:
    # 将打乱后的字符列表转回字符串，加入新列表
    shuffled_seq = ''.join(random.sample(seq, k=len(seq)))
    shuffled_list.append(shuffled_seq)

print(shuffled_list)

方案2：直接修改原列表（通过索引遍历） #

import random

finalList = ['AGTCTTTGTGT', 'ACGCGTAAACG', 'AGGGTCTC']

for i in range(len(finalList)):
    seq = finalList[i]
    # 替换原列表中的对应元素
    finalList[i] = ''.join(random.sample(seq, k=len(seq)))

print(finalList)

注：random.sample(seq, len(seq))会随机抽取序列的所有元素实现打乱，返回字符列表，''.join()用于将列表转回字符串。

二、生成包含指定数量A/G/T/C的随机DNA序列 #

如果需要生成指定碱基数量的DNA序列（比如每个序列含3个A、2个G、1个T、2个C），可以先构造对应碱基列表再打乱：

import random

def generate_custom_dna(a_count, g_count, t_count, c_count):
    # 构造包含指定数量碱基的列表
    bases = ['A']*a_count + ['G']*g_count + ['T']*t_count + ['C']*c_count
    # 原地打乱列表顺序
    random.shuffle(bases)
    # 转成字符串返回
    return ''.join(bases)

# 示例：生成5个含3A、2G、1T、2C的DNA序列
generated_sequences = [generate_custom_dna(3,2,1,2) for _ in range(5)]
print(generated_sequences)

注：random.shuffle()直接原地打乱列表，比random.sample()更高效。

内容的提问来源于stack exchange，提问作者user19278046