my ($a,$c,$g,$t) = (0, 0, 0, 0); #### my $a = 0; my $c = 0; my $g = 0; my $t = 0; #### sub comNucCount { my($string1, $string2) = @_; $string1 = uc $string1; $string2 = uc $string2; my %counts = map { $_ => 0 } qw( A C G T ); for my $position (0 .. length($string1)-1) { my $char = substr($string1,$position,1); $counts{$char}++ if substr($string2, $position, 1) eq $char; } return "A=$counts{A}, C=$counts{C}, G=$counts{G}, T=$counts{T}"; } #### sub comNucCount { my($string1, $string2) = @_; $string1 = uc $string1; $string2 = uc $string2; my @letters = qw( A C G T ) my %counts = map { $_ => 0 } @letters; for my $position (0 .. length($string1)-1) { my $char = substr($string1,$position,1); $counts{$char}++ if substr($string2, $position, 1) eq $char; } } return join ', ', map { "$_=$counts{$_}" } @letters; }