今天遇到了这个问题,但是搜索了好久也没找到真正的解决办法。后来找到了这个网站,按照里面提供的方法解决了。看到其他小伙伴也有同样困惑。我把网址贴出来,供大家参考。
参考网址:https://warwick.ac.uk/fac/sci/moac/people/students/peter_cock/python/genbank2fasta/
有一点小小的问题给大家提醒一下,我在从gbk文件中提取蛋白质序列时,直接复制了网址中的对应内容,更改文件名之后运行,提示错误。是因为print "Dealing with GenBank record %s" % seq_record.id 这行没加括号。加完括号就成功啦。