alignment.write() — write sequences and/or their alignment

write(file, alignment_format='PIR', alignment_features='INDICES CONSERVATION', align_block=0, align_alignment=False)
This command writes the whole alignment to a text file.

file can be either a file name or a writeable file handle (see modfile.File()).

alignment_format selects the format to write the alignment in; see alignment.append().

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 alignment.align2d(), where the corresponding quantities are defined. They refer to the 3D profile defined for the first align_block structures (run alignment.align2d() with fit = False to prepare these structural data without changing the input alignment). Similarly, the 'ACCURACY' option is valid only after the alignment.consensus() command. Options 'PRED_SS' and 'CONF_SS' are best exercised after reading in a ".PSS" file of secondary structure predictions. In the case of multiple sequences, it may be necessary to use the command Sequence.transfer_res_prop() first.

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 = False and align_block is ignored), a '*' indicating a conserved position is printed where all sequences have the same residue type. For blocks (align_alignment = True 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 alignment.align() command.

Example: See alignment.append() command.

Automatic builds 2017-02-17