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'