def reverse_complement(seq):
    """Get reverse complement of a given sequence"""
    comp_dict = {'A':'T','T':'A','G':'C','C':'G','N':'N'}
    return "".join([comp_dict[base] for base in reversed(seq)])