在列表之间随意交换项目

问题描述 投票:0回答:2

我正在编写一个小的脚本,涉及将一个密码子(我列表中的一个项目)更改为另一个密码子。首先,将DNA序列转换为密码子列表,然后我希望脚本选择一个随机的密码子,然后从密码子列表中将其更改为另一个密码子。问题在于,该项目只会在列表上第一次失效。 (即使选择了最后一个项目,它也会用列表中与该项目相同的第一个项目替换它)。我不知道我是否足够清楚,我还在学习。这是代码:

import random 


codones = ['ATA', 'ATC', 'ATT', 'ATG', 
           'ACA', 'ACC', 'ACG', 'ACT', 
           'AAC', 'AAT', 'AAA', 'AAG', 
           'AGC', 'AGT', 'AGA', 'AGG',                  
           'CTA', 'CTC', 'CTG', 'CTT', 
           'CCA', 'CCC', 'CCG', 'CCT', 
           'CAC', 'CAT', 'CAA', 'CAG', 
           'CGA', 'CGC', 'CGG', 'CGT', 
           'GTA', 'GTC', 'GTG', 'GTT', 
           'GCA', 'GCC', 'GCG', 'GCT', 
           "GAC", 'GAT', 'GAA', 'GAG', 
           'GGA', 'GGC', 'GGG', 'GGT', 
           'TCA', 'TCC', 'TCG', 'TCT', 
           'TTC', 'TTT', 'TTA', 'TTG', 
           'TAC', 'TAT', 'TGC', 'TGT', 
           'TGG']


dna = "atgaaaagcatgaaaagc" 
DNA = dna.upper()

print(DNA)

print("El gen tiene " + str(int(len(dna)/3)) + " Aa:")

def codones_dna(seq):
    codon = []
    for i in range(0, len(seq), 3):
        codon.append(seq[i:i +3])
    return codon
lista_codones = codones_dna(DNA)
print(lista_codones)

def translate(seq):
    table = { 
        'ATA':'I', 'ATC':'I', 'ATT':'I', 'ATG':'M', 
        'ACA':'T', 'ACC':'T', 'ACG':'T', 'ACT':'T', 
        'AAC':'N', 'AAT':'N', 'AAA':'K', 'AAG':'K', 
        'AGC':'S', 'AGT':'S', 'AGA':'R', 'AGG':'R',                  
        'CTA':'L', 'CTC':'L', 'CTG':'L', 'CTT':'L', 
        'CCA':'P', 'CCC':'P', 'CCG':'P', 'CCT':'P', 
        'CAC':'H', 'CAT':'H', 'CAA':'Q', 'CAG':'Q', 
        'CGA':'R', 'CGC':'R', 'CGG':'R', 'CGT':'R', 
        'GTA':'V', 'GTC':'V', 'GTG':'V', 'GTT':'V', 
        'GCA':'A', 'GCC':'A', 'GCG':'A', 'GCT':'A', 
        'GAC':'D', 'GAT':'D', 'GAA':'E', 'GAG':'E', 
        'GGA':'G', 'GGC':'G', 'GGG':'G', 'GGT':'G', 
        'TCA':'S', 'TCC':'S', 'TCG':'S', 'TCT':'S', 
        'TTC':'F', 'TTT':'F', 'TTA':'L', 'TTG':'L', 
        'TAC':'Y', 'TAT':'Y', 'TGC':'C', 'TGT':'C', 
        'TGG':'W' 
    } 
    protein = " "
    if len(seq) % 3 == 0:
        for i in range(0, len(seq), 3):
            codon = seq[i:i + 3]
            protein += table[codon].lower()
    return protein
protein = translate(DNA)
print(protein)


def mutation(seq):  
    for i in range(0, len(seq), 3): 
        dna_mut = seq[i:i]
        wt = random.choice(lista_codones)
        mutado = random.choice(codones)
        if wt != mutado:
            dna_mut += seq.replace(wt, mutado, 1)


    return dna_mut

dna_mut = mutation(DNA) 
python bioinformatics
2个回答
1
投票

首先,此行:


0
投票

也许这可以用:

© www.soinside.com 2019 - 2024. All rights reserved.