576 entries in table, 256 buckets number of buckets with 0 entries: 5 number of buckets with 1 entries: 45 number of buckets with 2 entries: 109 number of buckets with 3 entries: 77 number of buckets with 4 entries: 18 number of buckets with 5 entries: 2 number of buckets with 6 entries: 0 number of buckets with 7 entries: 0 number of buckets with 8 entries: 0 number of buckets with 9 entries: 0 number of buckets with 10 or more entries: 0 average search distance for entry: 1.8 Laa 20 blosum62: 20 20 A 4 4 R 5 9 N 6 6 D 6 6 C 9 9 Q 5 5 E 5 6 G 6 6 H 8 7 I 4 8 L 4 6 K 5 6 M 5 8 F 6 9 P 7 4 S 4 5 T 5 5 W 11 6 Y 7 7 V 4 11 0 : 4 > A 0 : 4 > I 0 : 4 > L 0 : 4 > S 0 : 4 > V 1 : 9 > C 2 : 6 > N 2 : 6 > D 2 : 6 > G 2 : 6 > F 3 : 6 > N 3 : 6 > D 3 : 6 > G 3 : 6 > F 4 : 9 > C 5 : 5 > R 5 : 5 > Q 5 : 5 > E 5 : 5 > K 5 : 5 > M 5 : 5 > T 6 : 6 > N 6 : 6 > D 6 : 6 > G 6 : 6 > F 7 : 6 > N 7 : 6 > D 7 : 6 > G 7 : 6 > F 8 : 7 > P 8 : 7 > Y 9 : 8 > H 10 : 6 > N 10 : 6 > D 10 : 6 > G 10 : 6 > F 11 : 6 > N 11 : 6 > D 11 : 6 > G 11 : 6 > F 12 : 8 > H 13 : 9 > C 14 : 4 > A 14 : 4 > I 14 : 4 > L 14 : 4 > S 14 : 4 > V 15 : 5 > R 15 : 5 > Q 15 : 5 > E 15 : 5 > K 15 : 5 > M 15 : 5 > T 16 : 5 > R 16 : 5 > Q 16 : 5 > E 16 : 5 > K 16 : 5 > M 16 : 5 > T 17 : 6 > N 17 : 6 > D 17 : 6 > G 17 : 6 > F 18 : 7 > P 18 : 7 > Y 19 : 11 > W 
wrong # args: should be "exit code"
    while executing
"exit"
    (in namespace eval "::request" script line 83)
    invoked from within
"namespace eval ::request $script"
    ("::try" body line 12)
OUTPUT BUFFER:
source /home/moumou/ordali2/src/LesProcsDeGscope.tcl
proc loadBLOSUM62Matrix {} {
    set blosum62 [dict create]

    # Provide the path to your BLOSUM62 matrix file
    set matrix_file_path "/home/moumou/ordali2/etc/blosum62.dat"

    set file [open $matrix_file_path "r"]

    set header [gets $file]
    while {[string index $header 0] eq "#"} {
	set header [gets $file]
    }
    regsub -all { +} [string trim $header] " " Laa
    set Laa [split $Laa]

    while {[gets $file line] != -1} {
        set elements [regexp -all -inline {\S+} $line]
        set aa1 [lindex $elements 0]
        set scores [lrange $elements 1 end]

        foreach aa2 $Laa sco $scores {
            dict set blosum62 "$aa1$aa2" $sco
        }
    }
    close $file
    puts "[dict info $blosum62]"
    return $blosum62
}

set blosum62 {
    { 4 -2  0 -1 -2 -1 -1  0 -2 -1 -2 -1 -1 -3 -1 -1 -1 -1 -2 -1 }
    {-2  9 -3 -2 -4  1 -3 -1 -3 -3 -1  0 -2 -3 -3 -2 -2 -3 -2 -3 }
    { 0 -3  6  0 -3  0 -1 -3 -1 -1 -3 -3 -2 -3 -4 -1 -1 -1 -1 -1 }
    {-1 -2  0  6 -3  0  0 -2  0  0 -2 -3 -2 -3 -3  0 -2 -1 -1 -3 }
    {-2 -4 -3 -3  9 -3 -4 -3 -3 -3 -1 -1 -3 -1 -2 -3 -1 -3 -2 -2 }
    {-1  1  0  0 -3  5 -2 -3  1  2 -2 -1  0 -3 -1  0 -1 -2 -1 -2 }
    {-1 -3 -1  0 -4 -2  6 -2 -4 -4 -2 -2 -2 -4 -3 -1 -3 -2 -2 -3 }
    { 0 -1 -3 -2 -3 -3 -2  6 -3  0  0  0 -1 -3 -1  0 -1 -2 -1  3 }
    {-2 -3 -1  0 -3  1 -4 -3  7  2 -3 -2 -1 -2 -4 -2 -3 -1 -1  1 }
    {-1 -3 -1  0 -3  2 -4  0  2  8 -3 -2  0 -3 -1 -1 -2 -3 -1 -2 }
    {-2 -1 -3 -2 -1 -2 -2  0 -3 -3  6  3 -1 -3 -4 -2 -2 -2 -1  0 }
    {-1  0 -3 -3 -1 -1 -2  0 -2 -2  3  6 -1 -2 -3 -1 -1 -2 -1  0 }
    {-1 -2 -2 -2 -3  0 -2 -1 -1  0 -1 -1  8 -3 -3 -1 -1 -2 -2 -2 }
    {-3 -3 -3 -3 -1 -3 -4 -3 -2 -3 -3 -2 -3  9 -3 -1 -1  0 -1 -3 }
    {-1 -3 -4 -3 -2 -1 -3 -1 -4 -1 -4 -3 -3 -3  4 -1 -1 -3 -3 -1 }
    {-1 -2 -1  0 -3  0 -1  0 -2 -1 -2 -1 -1 -1 -1  5  2  0 -1 -2 }
    {-1 -2 -1 -2 -1 -1 -3 -1 -3 -2 -2 -1 -1 -1 -1  2  5 -1 -1 -3 }
    {-1 -3 -1 -1 -3 -2 -2 -2 -1 -3 -2 -2 -2  0 -3  0 -1  6  3 -3 }
    {-2 -2 -1 -1 -2 -1 -2 -1 -1 -1 -1 -1 -2 -1 -3 -1 -1  3  7 -1 }
    {-1 -3 -1 -3 -2 -2 -3  3  1 -2  0  0 -2 -3 -1 -2 -3 -3 -1 11 }
}



set blo1 [loadBLOSUM62Matrix]
#set Laa {A C D E F G H I K L M N P Q R S T V W Y}
set Laa {A R N D C Q E G H I L K M F P S T W Y V}
puts "Laa [llength $Laa]"
puts "blosum62: [llength $blosum62] [llength [lindex $blosum62 0]]"


set i 0
foreach a $Laa {
    set d [dict get $blo1 "$a$a"]
    set m [lindex $blosum62 $i $i]
    puts "$a $d $m"
    incr i
}
for {set i 0} {$i < 20} {incr i} {
    set m [lindex $blosum62 $i $i]
    foreach a $Laa {
	if {[dict get $blo1 "$a$a"] == $m} {
	    puts "$i : $m > $a"
	}
    }
    puts ""
}




exit