def get_rfam_non_microrna_families_matching_mirbase(rfam_accs, rfam_microrna_accs, rfam_cmscan_accs):
    """
    Get a list of non-microRNA families that match miRBase sequences.
    """
    s1 = set(rfam_accs)
    s2 = set([x['rfam_acc'] for x in rfam_microrna_accs])
    rfam_non_microrna_accs = s1 - s2
    diff = rfam_non_microrna_accs.intersection(rfam_cmscan_accs)
    print('Found {} Rfam non-microRNA families matching miRBase sequences'.format(len(diff)))
    return diff