0.95 0.02 0.07 0.07 #A
0.03 0.01 0.06 0.83 #C
0.01 0.02 0.80 0.09 #G
0.01 0.95 0.07 0.01 #T

##</code><code>##

sub readdt2 {
  my $ifn = shift;
  open(my $IFH, "<$ifn") or die "cannot open file $ifn\n";
  my $line;
  my @nt = ("A","C","G","T");
  my %ret;
  my @tmp;
  for my $j(@nt) {
    $ret{$j} = [];
    $line = <$IFH>;
    chomp($line);
    @tmp = split(/\s+/,$line);
    for (my $i=0;$i<=$#tmp;$i++) {
      $ret{$j}[$i]= $tmp[$i] +0.0001;
    }
  }
  close($IFH);
  return(\%ret);
}