codon_table={"GCU":"A","GCC":"A","GCA":"A","GCG":"A","CGU":"R","CGC":"R","CGA":"R","CGG":"R","AGA":"R","AGG":"R","UCU":"S","UCC":"S","UCA":"S","UCG":"S","AGU":"S","AGC":"I","AUU":"I","AUC":"I","AUA":"I","AUU":"I","AUC":"I","AUA":"I","UUA":"L","UUG":"L","CUU":"L","CUC":"L","CUA":"L","CUG":"L","GGU":"G","GGC":"G","GGA":"G","GGG":"G","GUU":"V","GUC":"V","GUA":"V","GUG":"V","ACU":"T","ACC":"T","ACA":"T","ACG":"T","CCU":"P","CCC":"P","CCA":"P","CCG":"P","AAU":"N","AAC":"N","GAU":"D","GAC":"D","UGU":"C","UGC":"C","CAA":"Q","CAG":"Q","GAA":"E","GAG":"E","CAU":"H","CAC":"H","AAA":"K","AAG":"K","UUU":"F","UUC":"F","UAU":"Y","UAC":"Y","AUG":"M","UGG":"W","UAG":"STOP","UGA":"STOP","UAA":"STOP"}
rna=''
for line in open("A06662-RNA.fasta"):
if not line.startswith(">"):
rna=rna+line.strip()
#translate one frame at a time
for frame in range(3):
prot=''
print 'reading frame '+ str(frame + 1)
for i in range(frame, len(rna), 3):
codon=rna[i:i+3]
if codon in codon_table:
if codon_table[codon] == 'stop':
prot = prot + '*'
else:
prot = prot + codon_table[codon]
else:
#handle too short codons
prot = prot + '+'
#format to blocks of 48 columns
i=0
while i < len(prot):
print prot[i:i + 48]
i= i + 48
將RNA序列翻譯為蛋白質(zhì)序列
最后編輯于 :
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
- 文/潘曉璐 我一進店門麸折,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人粘昨,你說我怎么就攤上這事磕谅。” “怎么了雾棺?”我有些...
- 文/不壞的土叔 我叫張陵,是天一觀的道長衬浑。 經(jīng)常有香客問我捌浩,道長,這世上最難降的妖魔是什么工秩? 我笑而不...
- 正文 為了忘掉前任尸饺,我火速辦了婚禮,結(jié)果婚禮上助币,老公的妹妹穿的比我還像新娘浪听。我一直安慰自己,他們只是感情好眉菱,可當(dāng)我...
- 文/花漫 我一把揭開白布迹栓。 她就那樣靜靜地躺著,像睡著了一般俭缓。 火紅的嫁衣襯著肌膚如雪克伊。 梳的紋絲不亂的頭發(fā)上酥郭,一...
- 文/蒼蘭香墨 我猛地睜開眼坷衍,長吁一口氣:“原來是場噩夢啊……” “哼寝优!你這毒婦竟也來了?” 一聲冷哼從身側(cè)響起惫叛,我...
- 正文 年R本政府宣布,位于F島的核電站吟温,受9級特大地震影響序仙,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜鲁豪,卻給世界環(huán)境...
- 文/蒙蒙 一潘悼、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧爬橡,春花似錦治唤、人聲如沸。這莊子的主人今日做“春日...
- 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至,卻和暖如春辞槐,著一層夾襖步出監(jiān)牢的瞬間掷漱,已是汗流浹背。 一陣腳步聲響...
推薦閱讀更多精彩內(nèi)容
- 瀾媽想問下魅寶們,是否有過這樣的疑惑毅该,為什么同樣的飯博秫,閏蜜吃了沒反應(yīng),自己卻在瘋狂的長肉眶掌?為什么挡育!為什么!朴爬!為什么...
- 陪你見證孩子的每一步成長 觀察是有目的、有計劃痊银、比較持久的知覺[1]抵蚊。兒童期觀察力開始形成并迅速發(fā)展。由無到有,觀...
- 大家好,我叫馬小丹冈闭,她叫馬小花俱尼! 是的,我還是那個馬小丹萎攒,她遇八,也還是那個馬小花。 昨天也榄,2016年的最后一天顽分,我?guī)?..