有时我在bash shell中使用一组命令非常密集地工作,然后我花了几天时间在一个非常不同的工作.问题是,当我想继续执行第一组任务时,shell历史记录已经被最后一个完全覆盖了,我花了一些时间搞清楚我在做什么(我缺乏长期记忆).
我通常会保留一个日记,所以我可以回去继续我最后一次停止的地方,但我想知道是否有一种简单的方法可以通过保存最近的shell历史(例如,最后100个命令)到文件并使用此文件覆盖将来的历史记录.
今天是个好日子,
我开始在 Biopython 中编程,我想知道如何从具有所有特征坐标的基因组 GenBank 文件 (*.gb) 中提取基因序列和蛋白质标识符。
我的想法是创建一个包含蛋白质标识符、基因坐标和基因序列的文本文件。
如果您有任何想法,我将不胜感激。
到目前为止我试过这个:
for seq_record in seq_record.features:
if seq_record.type == 'CDS':
x=seq_record.qualifiers['protein_id']
i=seq_record.location._start.position
f=seq_record.location._end.position
sq = seq_record.seq
FEAT_LIST.append('START END STRAND ID')
FEAT_LIST.append(str(((i, f), s, x, sq)))
print(FEAT_LIST)
Run Code Online (Sandbox Code Playgroud)
但是,我收到此错误: sq = seq_record.seq AttributeError: 'SeqFeature' object has no attribute 'seq'
感谢您的帮助。