| 1 | import chimera
|
|---|
| 2 | incomplete = []
|
|---|
| 3 | for m in chimera.openModels.list(modelTypes=[chimera.Molecule]):
|
|---|
| 4 | for r in m.residues:
|
|---|
| 5 | tmplRes = chimera.restmplFindResidue(r.type, False, False)
|
|---|
| 6 | if not tmplRes:
|
|---|
| 7 | continue
|
|---|
| 8 | for aname in tmplRes.atomsMap.keys():
|
|---|
| 9 | if aname[0] == "H":
|
|---|
| 10 | continue
|
|---|
| 11 | if aname not in r.atomsMap:
|
|---|
| 12 | incomplete.append(r)
|
|---|
| 13 | break
|
|---|
| 14 | if incomplete:
|
|---|
| 15 | chimera.replyobj.info("Incomplete residues: "
|
|---|
| 16 | + ", ".join([str(r) for r in incomplete]) + "\n")
|
|---|
| 17 | chimera.selection.setCurrent(incomplete)
|
|---|
| 18 | else:
|
|---|
| 19 | chimera.replyobj.info("No incomplete residues\n")
|
|---|