biopython脚本给出错误'无效语法'

Question

我正在尝试应用biopython代码来分割一个大的fasta文件。代码如下所示：

def batch_iterator():
    entry=True # to make sure the loop run once
    while entry:
        batch=[]
        while len(batch) < batch_size:
            try:
                entry=iterator.next()
            except StopIteration:
                entry=None
            if entry is None:
                #end of file
                break
            batch.append(entry)
        if batch:
            yield batch
record_iter=SeqIO.parse(open('/home/to/file/sorted_sequence.fa', 'fasta')
for i, batch in enumerate (batch_iterator(record_iter, 93)):
    filename='gene_%i.fasta' % (i + 1)
    with open('/home/path/files/', filename, 'w') as ouput_handle:
        count=SeqIO.write(batch, ouput_handle, 'fasta')
    print ('Wrote %i records to %s' % (count, filename))

在这一行中：对于i，枚举中的批处理（批处理迭代器（record_iter，93））：正在给我语法错误：无效语法。但是我看不到错误，有人能帮我找到吗？我从这里取了密码http://biopython.org/wiki/Split_large_file谢谢

Answer 1

这一行缺少括号

record_iter = SeqIO.parse(open('/home/to/file/sorted_sequence.fa', 'fasta')

试试加一个

record_iter = SeqIO.parse(open('/home/to/file/sorted_sequence.fa'), 'fasta')