sub reverse_complement {
    local $_ = shift;
    y/ACGT/TGCA/;
    reverse;
}