Biopython翻译错误：TypeError：'generator'对象没有属性'__getitem__'

Question

我的代码有问题，这让我感到沮丧。

#!usr/bin/python
import os,sys
import Bio
from Bio import SeqIO
from Bio.Seq import Seq
from Bio.SeqRecord import SeqRecord
from Bio.Alphabet import IUPAC
fastaname = sys.argv[1]
in_put = open(fastaname,"rU")
out_put = open(fastaname[:-6] + "pro.fa","w")
m_record = SeqIO.parse(in_put, "fasta")
x = str(m_record)
start = x.find("ATG")
n = m_record[start:]
pro_records = (rec.translate(id = rec.id, to_stop = True) for rec in n)
SeqIO.write(pro_records, out_put, "fasta")
in_put.close()
out_put.close()

运行它

$ python script test.fasta
Traceback (most recent call last):
File "translation.py", line 14, in <module>
n = m_record[start:]
TypeError: 'generator' object has no attribute '__getitem__'

和测试文件test.fasta

>GRMZM5G800101_T01  
ATGATCTGGCATGTACAGAATTAG

Answer 1

像错误消息状态一样， m_record是生成器。

如果要像执行列表一样使用它，则应首先从中创建一个列表。

m_record = list(SeqIO.parse(in_put, "fasta"))