#!/usr/bin/perl
use strict;
use warnings;
use Bio::SeqIO;
use CGI;
my $cgi = new CGI;
print $cgi->header();
print <<__EOF;
Sequences with SD sequence or terminator structure:
__EOF
my %hash;
my $name1= "name1";
my $name2= "name2";
$hash{$name1} = "ATATTATCCCCCTATATATGGAGGGAGAGGGGGGGGGGGGGGGGGGGGGGGGGGGAGAGAGGAGATTTTTTTTTTTTTTTT"
$hash{$name2} = "ATATATTATATATATTATATTCGCGCGCGCGGCGCGCGCGGCGCGCGCGTTTTTTTTTTTTTTAGGAGAGAGAGGGAGGAGGAGAGGGGAGT"
foreach my $key (keys %hash) {
if ($hash{$key} =~ s!(aggag)!$1!gi) {
$hash2{$key} = $hash{$key};
}
if ($hash{$key} =~ s!(ttttt)!$1!ig) {
$hash2{$key} = $hash{$key};
}
my %hash2;
foreach my $key (keys %hash2) {
print " $key
";
print "$_
\n" for unpack '(A50)*', $hash2{$key};
}
print <<__EOF;
__EOF