String Matching with k Mismatches Moshe Lewenstein Bar

String Matching with k Mismatches Moshe Lewenstein Bar Ilan University Modified by Ariel Rosenfeld

String Matching with k Mismatches Landau – Vishkin 1986 Galil – Giancarlo 1986 Abrahamson 1987 Amir - Lewenstein - Porat 2000

Exact String Matching Input: T = t 1. . . P = p 1 … p m Output: All locations tn i of T where P appears Example: P= T= ABCAAB A B C A A B A A…

Exact String Matching Input: T = t 1. . . P = p 1 … p m Output: All locations tn i of T where P appears Example: P= T= ABCAAB A B C A A B A A… 3

Exact String Matching Input: T = t 1. . . P = p 1 … p m Output: All locations tn i of T where P appears Example: P= T= ABCAAB A B C A A B A A… 3 7

Exact String Matching Input: T = t 1. . . P = p 1 … p m Output: All locations tn i of T where P appears Example: P= T= ABCAAB A B C A A B A A… 3 7 11

Exact String Matching Input: T = t 1. . . P = p 1 … p m Output: All locations tn i of T where P appears Example: P= T= ABCAAB A B C A A B A A… Answer: {3, 7, 11, . . }

Exact String Matching Problem: Matching not exact in applications of: • Computational Biology • Musicology • Text Editing • Meteorology • etc. Need other definitions of string matching!

Approximate String Matching Idea: Find all text locations where distance from pattern is sufficiently small. distance metric: Let S = s 1 s 2…sm R = r 1 r 2…rm HAMMING DISTANCE Ham(S, R) = The number of locations j where sj Example: S = ABCABC R = ABBAAC Ham(S, R) = 2 rj

String Matching with Mismatches Input: T = t 1. . . P = p 1 … p m tn Output: For each i in T Ham(P, titi+1…ti+m-1) Example: P= T= ABBAAC A B C A C…

String Matching with Mismatches Input: T = t 1. . . P = p 1 … p m tn Output: For each i in T Ham(P, titi+1…ti+m-1) Example: P= T= ABBAAC A B C A C… 2 Ham(P, T 1) = 2

String Matching with Mismatches Input: T = t 1. . . P = p 1 … p m tn Output: For each i in T Ham(P, titi+1…ti+m-1) Example: P= T= ABBAAC A B C A C… 2, 4 Ham(P, T 2) = 4

String Matching with Mismatches Input: T = t 1. . . P = p 1 … p m tn Output: For each i in T Ham(P, titi+1…ti+m-1) Example: P= T= ABBAAC A B C A C… 2, 4, 6 Ham(P, T 3) = 6

String Matching with Mismatches Input: T = t 1. . . P = p 1 … p m tn Output: For each i in T Ham(P, titi+1…ti+m-1) Example: P= T= ABBAAC A B C A C… 2, 4, 6, 2 Ham(P, T 4) = 2

String Matching with Mismatches Input: T = t 1. . . P = p 1 … p m tn Output: For each i in T Ham(P, titi+1…ti+m-1) Example: P= T= ABBAAC A B C A C… 2, 4, 6, 2, …

String Matching with k Mismatches Input: T = t 1 . . . t n, Output: Every i in T s. t. Example: P= T= P = p 1 … p m Ham(P, titi+1…ti+m-1) k=2 ABBAAC A B C A C… 2, 4, 6, 2, … k

String Matching with k Mismatches Input: T = t 1 . . . t n, Output: Every i in T s. t. Example: P= T= P = p 1 … p m Ham(P, titi+1…ti+m-1) k=2 ABBAAC A B C A C… 2, 4, 6, 2, … k

String Matching with k Mismatches Input: T = t 1 . . . t n, Output: Every i in T s. t. Example: P= T= P = p 1 … p m Ham(P, titi+1…ti+m-1) k=2 ABBAAC A B C A C… 2, 4, 6, 2, … Y, N, N, Y, … k

Naïve Algorithm (for counting mismatches or k-mismatches problem) - Goto each location of text and compute hamming distance of P and Ti Running Time: O(nm) n = |T|, m = |P|

The Kangaroo Method (for k-mismatches) Landau – Vishkin 1986 Galil – Giancarlo 1986

Trie • A tree representing a set of strings. { aeef ad bbfe bbfg c } a e e f d c b b f e g

Trie (Cont) • Assume no string is a prefix of another Each string corresponds to a leaf. a e e f d c b b f e g

Compressed Trie • Compress unary nodes, label edges by strings a e e f d c b b eef f e g c a bbf d e g

Suffix tree of string s: a compressed trie of all suffixes of s Prefix-free: add a special character, say $, at the end of s

Suffix tree (Example) Let s = abab, a suffix tree of s is a compressed trie of all suffixes of s=abab$ { $ b$ ab$ bab$ abab$ } $ b a b $ $ $ a b $

$ Suffix Tree properties a b 5 b $ a a b $ $ b $ 3 1 - Succint in space - O(n). - Can be built in O(n) time. Mc. Creight, Weiner, Ukkonen, Farach-Colton 2 4

Exact string matching $ s=abab$ a b $ 1 b $ 3 5 $ a b $ 4 2 Given a pattern P = ab we traverse the tree according to the pattern.

Exact string matching $ s=abab$ 1 3 a b $ 1 b $ 3 5 $ a b $ 4 2 Leaves correspond to locations of appearance!

Exact string matching $ s=abab$ 1 3 a b $ 1 b $ 3 5 $ a b $ 4 2 Prepare Tree: O(n) time Find matches: O(m + occ) time occ = # of matches

Lowest common ancestors A lot more can be gained from the suffix tree if we preprocess it so that we can answer LCA queries on it

Why? The LCA of two leaves represents the longest common prefix (LCP) of these 2 suffixes $ s = abbaab$ a a $b b a a b 4 $ 1 b b $ 5 7 $ a a b $ b a 3 6 a b $ 2

Why? The LCA of two leaves represents the longest common prefix (LCP) of these 2 suffixes $ s = abbaab$ a a $b b a a b 4 $ 1 b b $ 5 7 $ a a b $ b a 3 6 a b $ 2

Why? The LCA of two leaves represents the longest common prefix (LCP) of these 2 suffixes $ s = abbaab$ a a $b b a a b 4 $ 1 b b $ 5 7 $ a a b $ b a 3 6 a b $ 2

Why? The LCA of two leaves represents the longest common prefix (LCP) of these 2 suffixes $ s = abbaab$ a a $b b a a b 4 $ 1 b b $ 5 7 $ a a b $ b a 3 6 a b $ 2

LCA/LCP properties $ b $ 4 a b b b a a b $ 1 7 a$ $ a a b $ 5 Preprocesssing time : O(n) Query Time: O(1) Harel & Tarjan 1984, Schieber & Vishkin 1988, Berkman & Vishkin 1993 3 b 6 a a b $ 2

The Kangaroo Method (for k-mismatches) - Create suffix tree for: s = P#T -Check P at each location i of T by kangrooing Example: P = ABABAABACAB T = ABBACABABABCABCA… i

The Kangaroo Method (for k-mismatches) - Create suffix tree for: s = P#T -Check P at each location i of T by kangrooing Example: P = ABABAABACAB T = ABBACABABABCABCA… i

The Kangaroo Method (for k-mismatches) - Create suffix tree for: s = P#T -Check P at each location i of T by kangrooing Example: P = ABABAABACAB T = ABBACABABABCABCA… i

The Kangaroo Method (for k-mismatches) - Create suffix tree for: s = P#T -Check P at each location i of T by kangrooing Example: P = ABABAABACAB T = ABBACABABABCABCA… i

The Kangaroo Method (for k-mismatches) - Create suffix tree for: s = P#T -Check P at each location i of T by kangrooing Example: P = ABABAABACAB T = ABBACABABABCABCA… i

The Kangaroo Method (for k-mismatches) - Create suffix tree for: s = P#T -Check P at each location i of T by kangrooing Example: P = ABABAABACAB T = ABBACABABABCABCA… i

The Kangaroo Method (for k-mismatches) - Create suffix tree for: s = P#T -Check P at each location i of T by kangrooing Example: P = ABABAABACAB T = ABBACABABABCABCA… i

The Kangaroo Method (for k-mismatches) - Create suffix tree for: s = P#T -Check P at each location i of T by kangrooing Example: P = ABABAABACAB T = ABBACABABABCABCA… i

The Kangaroo Method (for k-mismatches) - Create suffix tree for: s = P#T -Check P at each location i of T by kangrooing Example: P = ABABAABACAB T = ABBACABABABCABCA… i

The Kangaroo Method (for k-mismatches) - Create suffix tree for: s = P#T - Do up to k LCP queries for every text location Example: P = ABABAABACAB T = ABBACABABABCABCA… i

The Kangaroo Method (for k-mismatches) Preprocess: Build suffix tree of both P and T - O(n+m) time LCA preprocessing - O(n+m) time Check P at given text location Kangroo jump till next mismatch - O(k) time Overall time: O(nk)