next up previous contents index
Next: EDIT_ALIGNMENT edit Up: Comparison and searching of Previous: SEQUENCE_TO_ALI copy   Contents   Index

WRITE_ALIGNMENT -- write sequences and/or their alignment

FILE = <string:1> 'default' partial or complete filename
OUTPUT_DIRECTORY = <string:1> '' output directory
ATOM_FILES = <string:0> '' complete or partial atom filenames
ALIGN_CODES = <string:0> 'all' codes of proteins in the alignment
ALIGNMENT_FORMAT = <string:1> 'PIR' format of the alignment file: 'PIR' | 'PAP' | 'QUANTA' | 'INSIGHT' | 'FASTA'
ALIGN_BLOCK = <integer:1> 0 the last sequence in the first block of sequences
ALIGN_ALIGNMENT = <logical:1> off writing out an alignment of alignments (for *)
ALIGNMENT_FEATURES = <string:1> 'INDICES CONSERVATION' what alignment features to write out: 'ACCURACY' | 'HELIX' | 'BETA' | 'ACCESSIBILITY' | 'STRAIGHTNESS' | 'CONSERVATION' | 'INDICES' | 'ALL' | 'GAPS'
HETATM_IO = <logical:1> off whether to read HETATM coordinates
WATER_IO = <logical:1> off whether to read water coordinates
HYDROGEN_IO = <logical:1> off whether to read hydrogen coordinates
ATOM_FILES_DIRECTORY = <string:1> './' input atom files directory list (e.g., 'dir1:dir2:dir3:./:/')

Description:
This command writes the whole alignment to a text file.

The 'PAP' format, which corresponds to a relatively nice looking alignment, has several additional formatting options that can be selected by the ALIGNMENT_FEATURES variable. This scalar variable can contain any combination of the following keywords:

Options 'HELIX', 'BETA', 'ACCESSIBILITY', and 'STRAIGHTNESS' are valid only after executing command ALIGN2D, where the corresponding quantities are defined. They refer to the 3D profile defined for the first ALIGN_BLOCK structures (run ALIGN2D with FIT = off to prepare these structural data without changing the input alignment). Similarly, the 'ACCURACY' option is valid only after the CONSENSUS_ALIGNMENT command.

ALIGN_ALIGNMENT and ALIGN_BLOCK are used to ensure correct indication of identical alignment positions, depending on whether sequences or two blocks of sequences were aligned: For sequences (ALIGN_ALIGNMENT = off and ALIGN_BLOCK is ignored), a '*' indicating a conserved position is printed where all sequences have the same residue type. For blocks (ALIGN_ALIGNMENT = on and ALIGN_BLOCK indicates the last sequence of the first block), a '*' is printed only where the two blocks have the same order of residue types (there has to be the same number of sequences in both blocks). The blocks option is useful when comparing two alignments, possibly aligned by the ALIGN command.

Example: See READ_ALIGNMENT command.


next up previous contents index
Next: EDIT_ALIGNMENT edit Up: Comparison and searching of Previous: SEQUENCE_TO_ALI copy   Contents   Index
Ben Webb 2004-10-04