sf

Python code posted by vdfb
created at 06 Aug 00:42

Edit | Back 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
 
kmer_size = 12
def split_dna(dna, kmer_size):
    kmers = []
    for start in range(0,len(dna)-(kmer_size-1),1):
        kmer = dna[start:start+kmer_size]
        kmers.append(kmer)
    return kmers

kmer_counts = {}
dna_file = open('rosalind_1bba.txt','r')
for line in dna_file:
    dna = line.rstrip("\n")

    for kmer in split_dna(dna, kmer_size):
            current_count = kmer_counts.get(kmer, 0)
            new_count = current_count + 1
            kmer_counts[kmer] = new_count

max_kmers = ([key for key,val in kmer_counts.items() if val == max(kmer_counts.values())])
for kmer in max_kmers:
    print(kmer)
639 Bytes in 2 ms with coderay