def rc(dna):
    complements = maketrans('acgtrymkbdhvACGTRYMKBDHV', 'tgcayrkmvhdbTGCAYRKMVHDB')
    return dna.translate(complements)[::-1]
reverse complement
dataset/ETHPy150Open adamewing/bamsurgeon/bamsurgeon/replacereads.py/rc