import re, sys, os
# 构建序列名-序列字典
with open("protein_rm_enter.faa") as f:
Dict = {}
keys = []
for line in f:
if line[0] == ">":
key = line.strip()
keys.append(key)
Dict[key] = []
else:
Dict[key].append(line.strip())
# 遍历目标序列名,提取目标序列
with open("protein2.target", 'w') as o:
with open("apd.target") as target:
for tar in target:
tar = tar.strip()
tar = ">{}".format(tar)
if tar in keys:
o.write("{}\t{}\n".format(tar, ''.join(Dict[tar])))
#o.write("{}\n{}\n".format(tar, ''.join(Dict[tar])))
python: pick fna sequence
最后编辑于 :
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
- 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
推荐阅读更多精彩内容
- 主要内容源自解读《Fluent Python》,理解如有错误敬请指正:-) Python标准库提供了大量使用C来实...
- 这是由于我们在抓取网页的时候采用的是UTF-8编码,而存储时没有指定编码,在存储到文件的过程中就会报错。 处理方法...
- 错误 UnicodeEncodeError: 'gbk' codec can't encode character...
- 【背景】http://www.tuicool.com/articles/nEjiEv 问题参见: python2....
- Document Link: http://www.openbookproject.net/books/bpp4a...