def get_non_backbone_atoms(wtpdb, resids):
    """Return sorted list of atom-numbers that are part of resids and not backbone
    
    WARNING: BACKBONE IS NOT INCLUDING 'CA'

    :param wtpdb: wild-type pdbfile
    :param resids: residue numbers in pdbfile to search
    :type wtpdb: str, path to file
    :type resids: list of integers
    """

    indexes = []
    for line in open(wtpdb):
        if line[:6] == 'ATOM  ' or line[:6] == 'HETATM':
            resid = int(line[22:26].strip())
            if resid in resids:
                if not line[13:17].strip() in ['H', 'N', 'O', 'C']:
                    serial = int(line[6:11])
                    indexes.append(serial)
    return indexes