use strict;
my $filename1 = "sequences_with_upstream_stuff.fasta";

open(sequences,'sequences_with_upstream_stuff.fasta');
my @sequences = <sequences>; 
my @mainseq = ();
chomp @sequences;
foreach my $seqline(@sequences) 
	{ 
      my @temp = split("\t",$seqline);
	  push(@mainseq, \@temp);
    }
 #print $mainseq[16][0];
 
my $rvscomp = (); 
my $i = (); 
 
foreach (@mainseq)
{$rvscomp = reverse $mainseq[$i][1];
$rvscomp =~ tr/ACGT/TGCA/; # to get the reverse complement strand
print "$mainseq[$i][0]\n\nForward:\n\n$mainseq[$i][1]\n\nReverse:\n\n$rvscomp";
++$i;
}