def get_reciprocal_hits():
    """Checks if candidate is reciprocal best blast hit with gene

    returns: list with reciprocal hits"""
    reciprocal = {}
    proccesed = set()
    for line in open('tmp/dataset_diamond.res'):
        sline = line.split('\t')
        full_name = sline[0]
        sequence, gene = full_name.split('@')
        sequence = sequence[:-4]
        if sequence not in proccesed:
            proccesed.add(sequence)
            hit_gene = sline[1]
            reciprocal[sequence] = hit_gene
    return reciprocal
