Dear PyMol users!

I'm looking for the possible python script which using the pymol source would introduce selected mutations in the defined PDB file and produce PDB output containing such protein with the selected substitution residues. It would be also good if rotamers for mutation residues would be backbone dependent or taken from the existing rotamer libraries (although its not essence because I can run minimization on the mutants)

Previously I've done such tasks with selectivity mutations via pymol-GUI but now I'd like to obtain big series of the mutant of studied protein for further examinations of such mutants by means of molecular dynamics simulation.
Could someone provide me with such script if it could not be very complicated to make it?

Thanks for help,