今天遇到了這個問題朴上,但是搜索了好久也沒找到真正的解決辦法垒棋。后來找到了這個網(wǎng)站,按照里面提供的方法解決了痪宰〉鸺埽看到其他小伙伴也有同樣困惑。我把網(wǎng)址貼出來酵镜,供大家參考碉碉。
參考網(wǎng)址:https://warwick.ac.uk/fac/sci/moac/people/students/peter_cock/python/genbank2fasta/
有一點小小的問題給大家提醒一下垢粮,我在從gbk文件中提取蛋白質(zhì)序列時,直接復(fù)制了網(wǎng)址中的對應(yīng)內(nèi)容靠粪,更改文件名之后運行蜡吧,提示錯誤毫蚓。是因為print "Dealing with GenBank record %s" % seq_record.id 這行沒加括號。加完括號就成功啦昔善。