def getExclusiveReactionsMutants(project, ref_id, muts=set()):
    """
    Given a bunch of organisms, get two dicts
    mreacts: org_id -->  set(re_id, ...) (wild-type)
    mix: org_id --> set(re_id, ...) (mutated but also in wild-type)
    ereacts: org_id -->  set(re_id, ...) (exclusive mutated)
    """
    from ductape.storage.SQLite.database import Kegg
    kegg = Kegg(project)
    ereacts = kegg.getExclusiveReactionsMutants(ref_id, muts)
    mreacts = {}
    mix = {}
    for mut_id in muts:
        mut = set()
        for oR in kegg.getOrgReact(mut_id):
            mut.add(oR.re_id)
        mreacts[mut_id] = set()
        mix[mut_id] = set()
        for oR in kegg.getReferenceReact(mut_id, ref_id):
            if oR.re_id in mut:
                mix[mut_id].add(oR.re_id)
            else:
                mreacts[mut_id].add(oR.re_id)
    return mreacts, mix, ereacts
