Hello All

Some genome sequence share a counter. I mapped each genome sequence to the reference of the update counter, so there is some dereferencing.

Chris

#!/usr/bin/perl
use strict;
use warnings;

my @genome = ('a', 'a');
my $i = 0;

my @seq = qw/ aa ag ac at ta tg tc tt ca cg cc ct ga gg gc gt /;
my ($tt, $ag, $ac, $at, $ta, $tg, $ga, $cg, $cc, $gc) = (0) x 10;
my %counter;
@counter{@seq} = \($tt, $ag, $ac, $at, $ta, $tg, $ga, $tt, 
                  $tg, $cg, $cc, $ag, $ga, $cc, $gc, $ac);            
+      

if (exists $counter{ $genome[$i] . $genome[$i+1] }) {                 
+                   
    ++${ $counter{ $genome[$i] . $genome[$i+1] }} 
}
[download]