python筛选fasta中特定位置的基因

最新推荐文章于 2025-02-10 21:15:20 发布

panda-生物

最新推荐文章于 2025-02-10 21:15:20 发布

阅读量675

点赞数 1

文章标签： python

本文链接：https://blog.youkuaiyun.com/qq_45862222/article/details/110114386

版权

from Bio import SeqIO
fasta=input('请输入你的fasta文件：')
for seq_record in SeqIO.parse(fasta, "fasta"):
    name = input('输入你的基因名字（>xxx格式）：')
    star = eval(input('请输入起始位置（注意软件索引从0开始，且是半开区间）:'))
    end = eval(input('请输入终止位置（注意软件索引从0开始，且是半开区间）:'))
    print(name+'\n'+seq_record.seq[star:end])