fix 2

        for m, mnext in zip(moltypes[:-1], moltypes[1:]):

            sect = self.sections[m]

            molname = [x.split()[0] for x in sect if len(x.split()) > 1

                       and not x.startswith('[') and not x.startswith(';')][0]

                       and not x.lstrip().startswith('[') and not x.lstrip().startswith(';')][0]

            mols[molname] = [x for x in bondeds if m < x < mnext]

            mols[molname].sort()

        return mols