Skip to content

Commit 2f291b5

Browse files
author
Hamid Gasmi
committed
#199 refactored
1 parent 30cc9d8 commit 2f291b5

File tree

1 file changed

+4
-3
lines changed

1 file changed

+4
-3
lines changed

09-problems/graph-algorithms-in-genome-sequencing/assemble_phi_X174_error_free_reads.py

+4-3
Original file line numberDiff line numberDiff line change
@@ -74,15 +74,16 @@ def _build_adjacency_list_trie_kmp(self):
7474
for u in range(self._nodes_count):
7575

7676
sharing_kmer_reads = []
77+
suffix_u = self.nodes[u][self._read_size - sharing_kmer_size:self._read_size - 1]
7778
for v in range(self._nodes_count):
7879
if u == v:
7980
continue
8081

81-
if self.nodes[u][self._read_size - sharing_kmer_size:self._read_size - 1] in self.nodes[v]:
82+
if suffix_u in self.nodes[v]:
8283
sharing_kmer_reads.append(v)
83-
84+
8485
max_overlap_nodes = []
85-
max_overlap = 1 # this will prevent adding an edge with 0 weight (if max_overlap(u->v) == 0, then adjacents[u] will be empty)
86+
max_overlap = 1 # this will prevent adding edges of a nul weight (if max_overlap(u->v) == 0, then adjacents[u] will be empty)
8687
for v in sharing_kmer_reads:
8788

8889
overlap = kmp.max_overlap(self.nodes[u], self.nodes[v])

0 commit comments

Comments
 (0)