Skip to content
Permalink
Branch: master
Find file Copy path
Find file Copy path
Fetching contributors…
Cannot retrieve contributors at this time
30 lines (25 sloc) 906 Bytes
import os
# Input Hint
filename = input("Input filename: ")
k = int(input("how many files you want to split to: "))
# Load the original file
f = open(filename)
file = f.readlines()
length = len(file) # Total length of the file
length_seq = length//2 # Total number of seq(reads)
length_seq_num = length_seq//int(k) # Number of reads in each splitted file
# Splitting the file
for i in range(k-1):
temp = file[i*length_seq_num*2:i*length_seq_num*2+length_seq_num*2]
os.system("mkdir splitted"+str(i+1))
f = open("splitted"+str(i+1)+"/splitted"+".fa", 'w')
for L in temp:
f.write(L)
f.close()
# For the Last splitted file, save all the remained reads
os.system("mkdir splitted"+str(k))
temp = file[(k-1)*length_seq_num:]
f = open("splitted"+str(k)+"/splitted"+".fa", 'w')
for L in temp:
f.write(L)
f.close()
You can’t perform that action at this time.