def METHOD_NAME(self):
    # Is there a place to read from mol2 file?
    top = Topology.load(get_fn("ethane.mol2"))
    for site in top.sites:
        site.group = "ethane"
    assert np.all(np.array([site.group for site in top.sites]) == "ethane")