Ticket #1969: bug1969.py

File bug1969.py, 367 bytes (added by Greg Couch, 6 years ago)
Line 
1from collections import Counter
2m = session.models[0]
3ms = session.models[1]
4counts = Counter(pb.atoms[0].residue.chain_id for pb in ms.pseudobonds)
5print("# of missing segments:", len(ms.pseudobonds), counts)
6for c in m.chains:
7 print('Chain:', c.chain_id, len(c.residues), 'residues')
8 print(' missing:', [i for i, r in enumerate(c.residues) if r is None])
9