import os
from collections import OrderedDict
os.chdir('G:/R/Genome/hsa')
sequence = OrderedDict()
with open ("gene.fa", 'r') as f:
for line in f:
line = line.strip()
line = line.upper()
if line.startswith('>'):
idd = line.split(' ')[0][1:]
print(idd)
sequence[idd]=''
else:
sequence[idd]+=line
position = '2:6'
find_chr = position.split(':')[0]
find_pos = position.split(':')[1]
for idd, value in sequence.items():
if find_chr == idd:
print(idd,value[int(find_pos)-1])
用python根据指定染色体及坐标得到位置碱基（第九题）
最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
