my $mismatch = get_mismatches (\@dna, \@complement); sub get_mismatches { my ($dna, $complement) = @_; print "
TESTING TESTING TESTING SEGMENT::::: @$dna
COMP::::::: @$complement
"; my $i; my %good_pairs = ( 'A' => 'T', 'C' => 'G', 'G' => 'C', 'T' => 'A', ); foreach (my $i =0; $i < @$dna; $i++) { if ($dna[$i] ne $good_pairs {$complement[$i]}) { print "SUB - MISMATCH: $dna[$i-1]$dna[$i]/$complement[$i-1]$complement[$i]\n"; push @seg_mismatches,"$dna[$i-1]$dna[$i]/$complement[$i-1]$complement[$i]\n"; print "SUB - MISMATCH: $dna[$i]$dna[$i+1]/$complement[$i]$complement[$i+1]\n"; push @seg_mismatches, "$dna[$i]$dna[$i+1]/$complement[$i]$complement[$i+1]\n"; push @seg_mis, "$dna[$i-1]$dna[$i] $dna[$i]$dna[$i+1]\n"; } } return (@seg_mismatches); }