#!/opt/local/bin/python import sys seqsPerFile = int(sys.argv[1]) filename = sys.argv[2] if len(sys.argv) > 3 : stubfilename = sys.argv[3] else : stubfilename = filename print seqsPerFile print filename file = open(filename, 'r') data = file.readlines() index = 0 fileindex = 0 whole = "" for line in data: if line[0] == '>': index = index + 1 #print index #print seqsPerFile #print index > seqsPerFile if index > seqsPerFile: newFileName = stubfilename +".part"+ str(fileindex) print newFileName newfile = open(newFileName, 'w') newfile.write(whole) newfile.close() whole = "" index = 1 fileindex = fileindex + 1 whole = whole + line newFileName = stubfilename +".part"+ str(fileindex) print newFileName newfile = open(newFileName, 'w') newfile.write(whole) newfile.close() file.close()