1 research outputs found
Algorithms for path-constrained sequence alignment
International audienceWe define a novel variation on the constrained sequence alignment problem in which the constraint is given in the form of a regular expression. Given two sequences, an alphabet describing pairwise sequence alignment operations, and a regular expression over , the problem is to compute the highest scoring sequence alignment of the given sequences, such that . Two algorithms are given for solving this problem. The first basic algorithm is general and solves the problem in time and space, where and are the lengths of the two sequences and is the size of the NFA for . The second algorithm is restricted to patterns that do not contain the Kleene-closure star, and exploits this constraint to reduce the NFA size factor in the time complexity to a smaller factor . is compacted by supporting alignment patterns extended by \emph{meta-characters} including general insertion, deletion and match operations, as well as some cases of substitutions. For a regular expression , these time bounds range from to , depending on the meta-characters used in . An additional result obtained along the way is an extension of the algorithm of Fischer and Paterson for String Matching with Wildcards. Our extension allows the input strings to include "negation symbols" (that match all letters but a specific one) while retaining the original time complexity. We implemented both algorithms and applied them to data-mine new miRNA seeding patterns in \textit{C. elegans} Clip-seq experimental data