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 delta降尺度计算的一些细节,有偿
  • ¥15 Arduino红外遥控代码有问题
  • ¥15 数值计算离散正交多项式
  • ¥30 数值计算均差系数编程
  • ¥15 redis-full-check比较 两个集群的数据出错
  • ¥15 Matlab编程问题
  • ¥15 训练的多模态特征融合模型准确度很低怎么办
  • ¥15 kylin启动报错log4j类冲突
  • ¥15 超声波模块测距控制点灯,灯的闪烁很不稳定,经过调试发现测的距离偏大
  • ¥15 import arcpy出现importing _arcgisscripting 找不到相关程序