最近看到同学用Perl写了一个提取fasta序列前500bp的脚本,于是自己用Python也写了一个实现相同功能的脚本。 脚本内容如下: """ @author: will @Date: July 2019 @E-mail: willgyw@126.com @Description: """ import sys seq = { } res = open('500bp.fa', 'w') with open(sys