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条)

报告相同问题?

悬赏问题

  • ¥30 Matlab打开默认名称带有/的光谱数据
  • ¥50 easyExcel模板 动态单元格合并列
  • ¥15 res.rows如何取值使用
  • ¥15 在odoo17开发环境中,怎么实现库存管理系统,或独立模块设计与AGV小车对接?开发方面应如何设计和开发?请详细解释MES或WMS在与AGV小车对接时需完成的设计和开发
  • ¥15 CSP算法实现EEG特征提取,哪一步错了?
  • ¥15 游戏盾如何溯源服务器真实ip?需要30个字。后面的字是凑数的
  • ¥15 vue3前端取消收藏的不会引用collectId
  • ¥15 delphi7 HMAC_SHA256方式加密
  • ¥15 关于#qt#的问题:我想实现qcustomplot完成坐标轴
  • ¥15 下列c语言代码为何输出了多余的空格