mdl = <model> | model to read sequence from | |
align_codes = <str:0> | 'all' | codes of proteins in the alignment |
atom_files = <str:0> | '' | complete or partial atom filenames |
You should additionally set align_codes and atom_files to the PDB ID and file name, respectively. This information is added to the alignment with the new sequence. (Alternatively, you can set this information later by setting alignment[x].code and alignment[x].atom_file, where 'x' is the index of the sequence in the alignment, starting at zero.)
Example: See alignment.compare_with() command.