Elsa镇魂女孩 2019-06-27 21:38 采纳率: 62.5%
浏览 1837
已采纳

如何用python写一个循环,或者其他方式,完成批量选取内容并保存。

下面是将文件 split__1.fasta中选取一段序列,如下421至480. 存于selected_split_1.fasta文件中__。
如何写一个循环,可以从split__1.fasta 至 split__68.fasta,中分别选取421至480. 并存于selected_split_1.fasta至selected_split_68.fasta文件中文件中

from Bio import SeqIO

fin = open('split_1.fasta', 'r')
fout = open('selected_split_1.fasta', 'w')

with open("selected_split_1.fasta","w") as f:
        for seq_record in SeqIO.parse("split_1.fasta", "fasta"):
                f.write(">")
                f.write(str(seq_record.id) + "\n")
                f.write(str(seq_record.seq[421:480]) + "\n")  #start 421 to end 480 base positions

fin.close()
fout.close()
  • 写回答

3条回答 默认 最新

  • threenewbee 2019-06-28 05:18
    关注

    from Bio import SeqIO

    for xx in range(1, 68):
    xn = "split__" + str(xx) + ".fasta"
    yn = "selected_split_" + str(xx) + ".fasta"

    fin = open(xn, 'r')
    fout = open(yn, 'w')

    with open(yn,"w") as f:
    for seq_record in SeqIO.parse(xn, "fasta"):
    f.write(">")
    f.write(str(seq_record.id) + "\n")
    f.write(str(seq_record.seq[421:480]) + "\n") #start 421 to end 480 base positions

    fin.close()
    fout.close()

    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论
查看更多回答(2条)

报告相同问题?

悬赏问题

  • ¥15 多尺度特征图提取和融合
  • ¥15 微信小程序:渲染收货地址时页面不显示
  • ¥20 win7 64位DirectShow提示初始化失败如何解决?
  • ¥20 小规模孤立词识别系统设计
  • ¥15 关于Java对接海康威视车牌识别一体机SDK是否需要固定外网的IP?
  • ¥15 Linux扩容时,格式化卡住了:vgdispaly查看卷组信息,没有输出
  • ¥18 关于#ubuntu#的问题:使用背景-工作职责内有七八台ubuntu系统的电脑,平时需要互相调取资料,想实现把这几台电脑用交换机组成一个局域网,来实现指定文件夹的互相调取和写入
  • ¥20 求一个简易射频信号综测仪
  • ¥15 esp8266 tally灯 接收端改为发射端
  • ¥30 Labview代码调用access 数据库,相同代码其中一个调用不出来是为什么