sub reverse_complement { local $_ = shift; y/ACGT/TGCA/; reverse; } #### while () { chomp; if (/^>\s*(.*)/) { # it's a sequence label $name = $1; } else { # it's sequence data } }