Python dna补体中的关键错误
这里的错误是文件Python dna补体中的关键错误,python,python-2.7,python-3.x,bioinformatics,dna-sequence,Python,Python 2.7,Python 3.x,Bioinformatics,Dna Sequence,这里的错误是文件 import string import os,sys file=open("C:\Python27\\New Text Document.txt",'r')\ seq =file.readlines() basecomplement = {'A': 'T', 'C': 'G', 'G': 'C', 'T': 'A'} def translate(seq): aaseq = [] for str in seq: aaseq.append(bas
import string
import os,sys
file=open("C:\Python27\\New Text Document.txt",'r')\
seq =file.readlines() basecomplement = {'A': 'T', 'C': 'G', 'G': 'C', 'T': 'A'}
def translate(seq):
aaseq = []
for str in seq:
aaseq.append(basecomplement[str])
return aaseq
for frame in range(1):
rseqn= (''.join(item.split('|')[0] for item in translate(seq[frame:])))
rseqn = rseqn[::-1]
print(rseqn)
print(seq)
如果我使用简单字符串,我也会获得序列,但我试图将文件作为输入读取,如果我理解您试图获取的内容,则会显示以下错误
"C:\Users\ram\Desktop\pythonhw\dnacomp.py", line 13, in <module>
rseqn= (''.join(item.split('|')[0] for item in
translate(seq[frame:]))) File
"C:\Users\ram\Desktop\pythonhw\dnacomp.py", line 10, in translate
aaseq.append(basecomplement[str]) KeyError:
'AGTCTGGCATACCAGTACAGACTATCA'
测试
def translate(seq):
basecomplement = {'A': 'T', 'C': 'G', 'G': 'C', 'T': 'A'}
return ''.join(basecomplement[letter] for letter in seq)
要翻译文件,您可以执行以下操作
s = 'ATCGGCATAATAC'
>>> translate(s)
'TAGCCGTATTATG'
请清理代码!!!我知道我可以用这个函数作为文件输入吗?你能举个例子说明你的文件是什么样子吗?dna文件快速格式如果你想让我在这里粘贴数据,我也可以,我必须在shell或程序中使用这个代码
with open('file.txt', 'r') as f:
for line in f:
print(translate(line))