# we only want to show HBD/HBA
if contacted_lig_at.type == "N" and contacted_res_at.type == "O" or \
contacted_lig_at.type == "O" and contacted_res_at.type == "N":
contacts.append([contacted_lig_at, contacted_res_at])
In util is a bit fragile as it depends on atom types, instead should query atom elements, as types can be anything.
In util is a bit fragile as it depends on atom
types
, instead should query atom elements, as types can be anything.