OESequenceAlignment OEGetAlignment(OEChem::OEMolBase &mol1, OEChem::OEMolBase &mol2,
unsigned int assume = OEAssumption::Default,
unsigned int method = OESeqAlignmentMethod::PAM250,
int gap=-10, int extend=-2)
Returns the highest scored alignment of all pairwise chain alignments calculated by
- mol1, mol2
- The two proteins of which alignment is being calculated.
- A bitmask that indicates whether the algorithm can assume that the given proteins
have perceived residues, have PDB ordered atoms, or have bonded residues.
See the OEAssumption namespace for more information.
- This value has to be from the OESeqAlignmentMethod
namespace and it defines the algorithm that is used to gerenerate the
- gap, extent
- The gap and extension penalty to be used in the alignment.
See also Gap penalty in Wikipedia.