![](/img/trans.png)
[英]How to convert a set of DNA sequences into protein sequences using python programming?
[英]What changes do I need to make in this python code to convert DNA sequences to protein?
我必須首先找到密碼子“ATG”,然后一旦找到它應該使用字典翻譯的密碼子。 我目前沒有 output。 我在 python 方面仍然缺乏經驗,因此無法使用正確的語法編寫代碼
def translate(line, table):
protein = ""
for i in range (0, len(line), 3):
codon = line[i:i+3]
protein += table[codon]
print(protein)
return protein
def find_start(sequence):
for line in sequence:
if line == "ATG":
translate(line)
return
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', 'TAA':'_', 'TAG':'_',
'TGC':'C', 'TGT':'C', 'TGA':'_', 'TGG':'W',
}
sequence = "ggatcatagtcttgttgcattgtaggtgtttctttacagtatccttcttaatgaatcatgattatgttctaagtagaccagatcgattcttacgactacaatattttcttttgagccatcatagattttttttagtttcaaaccagccttgcattgtgttctcctcgatgatgtgttacggattctaggattactagctaacctttctgtatttctaccctccattgaacaatttaac"
sequence = sequence.upper()
find_start(sequence)
將您的查找開始 function 更新為
def find_start(sequence):
# index = sequence.find('ATG') # if yow want start index of codon
if ('ATG' in sequence):
# do some thing
我仍然不是 100% 你想要達到的目標,但希望它有所幫助
for i in range(sequence.find("ATG"), len(sequence) - 2):
possible_codon = sequence[i:i + 3]
if possible_codon not in table:
continue
print(f"key: {possible_codon}, val: {table[possible_codon]}")
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.