use strict;
use warnings;
use Term::ANSIColor;
use autodie;

#Program to find motif site in a given protein sequence using files

my $motif = "AGGGGG";
open( my $read, "<dna.txt" );
my @e = <$read>;
$_ = join( " ", @e );
s/\s+//g;
my @c;
push @c, pos( ) - length( $motif ) + 1 while /$motif/g;
s/$motif/color( 'bold green' ) . $motif . color( 'black' )/eg;
print $_, "\n";
print "Number of sites the motif (AGGGGG) is present: ", scalar @c, "\n";
print "And the positions in the string are: ", join( ',', @c ), "\n\n";