# MySeqTools.pm - my collection of useful sequence utilities
package MySeqTools;

# subroutine for working out the
# reverse complement of a sequence
sub reverse_complement {
        my $seq = shift @_;
        my $rev = reverse $seq;
        $rev =~ tr/GATC/CTAG/;
        return $rev;
}

# we could include more subroutine definitions here...

return 1;  # end of package