package Sequence;

sub new
	{
	my $class = shift;
	my $self = {
		_id		=> shift,
		_name	=> shift,
		_seq	=> shift,
	};
	bless $self, $class;
	return $self;
}

sub randomseq
{
my ($sel, $lenght) = @_;
my @nucleotides = ('a','c','g','t');
	$sel->{_seq} = '';
	for (my $i = 0; $i < $lenght; $i++)
	{
		my $base = $nucleotides[int(rand(4))];
		$sel->{_seq} .= $base;
	}
   return $sel->{_seq};
}

sub reversecomplement
	{
	my $self = shift;
	my $seqv = $self->{_seq};
	$self->{_seq} =~ tr/actg/tgac/;
	return $self->{_seq};
	}
1;