Next: APPENDICES Up: No Title Previous: Generating Alignments of

Generating multiple pattern alignments with each sequence

In the previous examples, only the best scoring alignment of the pattern with the sequence entries is output. However, a pattern may match equally well at more than one position in a sequence (for example, if a repeat is present). The following example, illustrates the principal of generating multiple pattern alignments, though the results are of limited interest.

(file bash_ge_4_scan1_top20_patt2.com)


output_file=bash_ge_4_scan1_top20.junk2
mode = scan
block_file/pattern=bash_ge_4.bloc,1
matrix_file=md.mat
pattern_level=5
database=bash_ge_4_scan1.top20seq
print_horizontal/pattern=bash_ge_4_scan1_top20.patt2

The only new command necessary is the PATTERN_LEVEL=N command. This specifies that the top N ways of aligning the pattern to EACH sequence will be output.

As before, the program PATT may be used to compress the results and store them in the file bash_ge_4_scan1_top20.patt2out. Inspection of this file shows that in this example, the alternative pattern alignments have identical scores and differ only in the positioning of the first group of pattern elements.


gjb@bioch.ox.ac.uk