next up previous contents index
Next: SUPERPOSE superpose Up: Comparison and searching of Previous: MALIGN align   Contents   Index

ALIGN_CONSENSUS -- consensus sequence alignment

GAP_PENALTIES_1D = <real:2> 900 50 gap creation and extension penalties for sequence/sequence alignment
ALIGN_BLOCK = <integer:1> 0 the last sequence in the first block of sequences
STOP_ON_ERROR = <integer:1> 1 whether to stop on error
READ_WEIGHTS = <logical:1> off whether to read the whole NxM weight matrix for ALIGN*
WRITE_WEIGHTS = <logical:1> off whether to write the whole NxM weight matrix for ALIGN*
INPUT_WEIGHTS_FILE = <string:1> ''
OUTPUT_WEIGHTS_FILE = <string:1> ''

Output:
MODELLER_STATUS = <integer:1>

Description:
This command is similar to ALIGN except that a consensus alignment of two blocks of sequences is produced. A consensus alignment is obtained from a consensus similarity matrix using the specified gap penalties and the global dynamic programming method. The consensus similarity matrix is obtained by aligning the two blocks of sequences many times with different parameters and methods and counting how many times each pair was aligned. This command is still experimental and no detailed description is given at this time.

This command also produces the alignment accuracy that can be printed out by the WRITE_ALIGNMENT command in the 'PAP' format (0 inaccurate, 9 accurate). If the gap initiation penalty is 0, the gap extension penalty of say 0.4 means that only those positions will be equivalenced that were aligned in at least 80% of the individual alignments (i.e., 2 times 0.40).

Example:


# Example for: ALIGN_CONSENSUS

# This will read 2 sequences and prepare a consensus alignment
# from many different pairwise alignments.

READ_ALIGNMENT FILE = 'toxin.ali', ALIGN_CODES = '2ctx' '2abx'
ALIGN_CONSENSUS GAP_PENALTIES_1D= 0 0.4, ALIGN_BLOCK = 1
WRITE_ALIGNMENT FILE = 'toxin-seq.pap', ALIGNMENT_FORMAT = 'PAP'


next up previous contents index
Next: SUPERPOSE superpose Up: Comparison and searching of Previous: MALIGN align   Contents   Index
Ben Webb 2004-10-04