--- /dev/null
+
+char *amino_acid_order = "ABCDEFGHIKLMNPQRSTVWXYZ";
+
+//Jones Taylor Thornton, tm matrix,
+//H JOND940101
+//D The 250 PAM transmembrane protein exchange matrix (Jones et al., 1994)
+//R LIT:2006072 PMID:8112466
+//A Jones, D.T., Taylor, W.R. and Thornton, J.M.
+//T A mutation data matrix for transmembrane proteins
+//J FEBS Lett. 339, 269-275 (1994)
+int jtttm250mt[]={
+ 2,
+ 0, 0,
+ 0, 0, 6,
+ 0, 0, -3, 12,
+ 0, 0, -3, 8, 13,
+ -2, 0, 1, -6, -6, 5,
+ 1, 0, -1, 3, 3, -4, 6,
+ -3, 0, -1, 3, 2, -3, -3, 11,
+ 0, 0, -1, -3, -4, -1, -2, -4, 2,
+ -2, 0, -3, 3, 1, -5, -1, 4, -4, 12,
+ -2, 0, -1, -5, -5, 1, -4, -4, 1, -4, 3,
+ -1, 0, -1, -3, -3, 0, -3, -3, 1, -1, 1, 3,
+ -1, 0, -1, 6, 1, -4, -2, 3, -3, 5, -4, -2, 11,
+ 0, 0, -4, -2, -3, -4, -2, -4, -3, -4, -1, -3, -2, 11,
+ -2, 0, -3, 2, 7, -4, -1, 7, -4, 6, -2, -2, 3, 0, 11,
+ -1, 0, -1, 1, 2, -4, 0, 5, -3, 9, -3, 0, 2, -3, 6, 7,
+ 2, 0, 1, 0, 0, -1, 1, -2, -1, -1, -2, -2, 2, -1, -1, -1, 3,
+ 1, 0, 0, 0, -1, -2, 0, -2, 0, -2, -1, 0, 1, -1, -2, -1, 2, 3,
+ 0, 0, 0, -3, -2, -1, -1, -4, 2, -4, 0, 1, -3, -3, -4, -2, -1, 0, 2,
+ -4, 0, 1, -4, -3, -3, -2, -1, -3, 3, -2, -2, -3, -6, 0, 5, -3, -4, -2, 12,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+ -3, 0, 3, -2, -5, 2, -5, 6, -4, 1, -3, -3, -1, -5, 0, -1, 0, -3, -4, -2, 0, 10,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
+
+
+int blosum30mt[]={
+ 4,
+ 0, 5,
+ -3, -2, 17,
+ 0, 5, -3, 9,
+ 0, 0, 1, 1, 6,
+ -2, -3, -3, -5, -4, 10,
+ 0, 0, -4, -1, -2, -3, 8,
+ -2, -2, -5, -2, 0, -3, -3, 14,
+ 0, -2, -2, -4, -3, 0, -1, -2, 6,
+ 0, 0, -3, 0, 2, -1, -1, -2, -2, 4,
+ -1, -1, 0, -1, -1, 2, -2, -1, 2, -2, 4,
+ 1, -2, -2, -3, -1, -2, -2, 2, 1, 2, 2, 6,
+ 0, 4, -1, 1, -1, -1, 0, -1, 0, 0, -2, 0, 8,
+ -1, -2, -3, -1, 1, -4, -1, 1, -3, 1, -3, -4, -3, 11,
+ 1, -1, -2, -1, 2, -3, -2, 0, -2, 0, -2, -1, -1, 0, 8,
+ -1, -2, -2, -1, -1, -1, -2, -1, -3, 1, -2, 0, -2, -1, 3, 8,
+ 1, 0, -2, 0, 0, -1, 0, -1, -1, 0, -2, -2, 0, -1, -1, -1, 4,
+ 1, 0, -2, -1, -2, -2, -2, -2, 0, -1, 0, 0, 1, 0, 0, -3, 2, 5,
+ 1, -2, -2, -2, -3, 1, -3, -3, 4, -2, 1, 0, -2, -4, -3, -1, -1, 1, 5,
+ -5, -5, -2, -4, -1, 1, 1, -5, -3, -2, -2, -3, -7, -3, -1, 0, -3, -5, -3, 20,
+ 0, -1, -2, -1, -1, -1, -1, -1, 0, 0, 0, 0, 0, -1, 0, -1, 0, 0, 0, -2, -1,
+ -4, -3, -6, -1, -2, 3, -3, 0, -1, -1, 3, -1, -4, -2, -1, 0, -2, -1, 1, 5, -1, 9,
+ 0, 0, 0, 0, 5, -4, -2, 0, -3, 1, -1, -1, -1, 0, 4, 0, -1, -1, -3, -1, 0, -2, 4};
+
+int blosum40mt[]={
+ 5,
+ -1, 5,
+ -2, -2, 16,
+ -1, 6, -2, 9,
+ -1, 1, -2, 2, 7,
+ -3, -3, -2, -4, -3, 9,
+ 1, -1, -3, -2, -3, -3, 8,
+ -2, 0, -4, 0, 0, -2, -2, 13,
+ -1, -3, -4, -4, -4, 1, -4, -3, 6,
+ -1, 0, -3, 0, 1, -3, -2, -1, -3, 6,
+ -2, -3, -2, -3, -2, 2, -4, -2, 2, -2, 6,
+ -1, -3, -3, -3, -2, 0, -2, 1, 1, -1, 3, 7,
+ -1, 4, -2, 2, -1, -3, 0, 1, -2, 0, -3, -2, 8,
+ -2, -2, -5, -2, 0, -4, -1, -2, -2, -1, -4, -2, -2, 11,
+ 0, 0, -4, -1, 2, -4, -2, 0, -3, 1, -2, -1, 1, -2, 8,
+ -2, -1, -3, -1, -1, -2, -3, 0, -3, 3, -2, -1, 0, -3, 2, 9,
+ 1, 0, -1, 0, 0, -2, 0, -1, -2, 0, -3, -2, 1, -1, 1, -1, 5,
+ 0, 0, -1, -1, -1, -1, -2, -2, -1, 0, -1, -1, 0, 0, -1, -2, 2, 6,
+ 0, -3, -2, -3, -3, 0, -4, -4, 4, -2, 2, 1, -3, -3, -3, -2, -1, 1, 5,
+ -3, -4, -6, -5, -2, 1, -2, -5, -3, -2, -1, -2, -4, -4, -1, -2, -5, -4, -3, 19,
+ 0, -1, -2, -1, -1, -1, -1, -1, -1, -1, -1, 0, -1, -2, -1, -1, 0, 0, -1, -2, -1,
+ -2, -3, -4, -3, -2, 4, -3, 2, 0, -1, 0, 1, -2, -3, -1, -1, -2, -1, -1, 3, -1, 9,
+ -1, 2, -3, 1, 5, -4, -2, 0, -4, 1, -2, -2, 0, -1, 4, 0, 0, -1, -3, -2, -1, -2, 5};
+
+int blosum45mt[]={
+ 5,
+ -1, 4,
+ -1, -2, 12,
+ -2, 5, -3, 7,
+ -1, 1, -3, 2, 6,
+ -2, -3, -2, -4, -3, 8,
+ 0, -1, -3, -1, -2, -3, 7,
+ -2, 0, -3, 0, 0, -2, -2, 10,
+ -1, -3, -3, -4, -3, 0, -4, -3, 5,
+ -1, 0, -3, 0, 1, -3, -2, -1, -3, 5,
+ -1, -3, -2, -3, -2, 1, -3, -2, 2, -3, 5,
+ -1, -2, -2, -3, -2, 0, -2, 0, 2, -1, 2, 6,
+ -1, 4, -2, 2, 0, -2, 0, 1, -2, 0, -3, -2, 6,
+ -1, -2, -4, -1, 0, -3, -2, -2, -2, -1, -3, -2, -2, 9,
+ -1, 0, -3, 0, 2, -4, -2, 1, -2, 1, -2, 0, 0, -1, 6,
+ -2, -1, -3, -1, 0, -2, -2, 0, -3, 3, -2, -1, 0, -2, 1, 7,
+ 1, 0, -1, 0, 0, -2, 0, -1, -2, -1, -3, -2, 1, -1, 0, -1, 4,
+ 0, 0, -1, -1, -1, -1, -2, -2, -1, -1, -1, -1, 0, -1, -1, -1, 2, 5,
+ 0, -3, -1, -3, -3, 0, -3, -3, 3, -2, 1, 1, -3, -3, -3, -2, -1, 0, 5,
+ -2, -4, -5, -4, -3, 1, -2, -3, -2, -2, -2, -2, -4, -3, -2, -2, -4, -3, -3, 15,
+ 0, -1, -2, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 0, 0, -1, -2, -1,
+ -2, -2, -3, -2, -2, 3, -3, 2, 0, -1, 0, 0, -2, -3, -1, -1, -2, -1, -1, 3, -1, 8,
+ -1, 2, -3, 1, 4, -3, -2, 0, -3, 1, -2, -1, 0, -1, 4, 0, 0, -1, -3, -2, -1, -2, 4};
+
+int blosum50mt[] = {
+ 5,
+ -2, 5,
+ -1, -3, 13,
+ -2, 5, -4, 8,
+ -1, 1, -3, 2, 6,
+ -3, -4, -2, -5, -3, 8,
+ 0, -1, -3, -1, -3, -4, 8,
+ -2, 0, -3, -1, 0, -1, -2, 10,
+ -1, -4, -2, -4, -4, 0, -4, -4, 5,
+ -1, 0, -3, -1, 1, -4, -2, 0, -3, 6,
+ -2, -4, -2, -4, -3, 1, -4, -3, 2, -3, 5,
+ -1, -3, -2, -4, -2, 0, -3, -1, 2, -2, 3, 7,
+ -1, 4, -2, 2, 0, -4, 0, 1, -3, 0, -4, -2, 7,
+ -1, -2, -4, -1, -1, -4, -2, -2, -3, -1, -4, -3, -2, 10,
+ -1, 0, -3, 0, 2, -4, -2, 1, -3, 2, -2, 0, 0, -1, 7,
+ -2, -1, -4, -2, 0, -3, -3, 0, -4, 3, -3, -2, -1, -3, 1, 7,
+ 1, 0, -1, 0, -1, -3, 0, -1, -3, 0, -3, -2, 1, -1, 0, -1, 5,
+ 0, 0, -1, -1, -1, -2, -2, -2, -1, -1, -1, -1, 0, -1, -1, -1, 2, 5,
+ 0, -4, -1, -4, -3, -1, -4, -4, 4, -3, 1, 1, -3, -3, -3, -3, -2, 0, 5,
+ -3, -5, -5, -5, -3, 1, -3, -3, -3, -3, -2, -1, -4, -4, -1, -3, -4, -3, -3, 15,
+ -1, -1, -2, -1, -1, -2, -2, -1, -1, -1, -1, -1, -1, -2, -1, -1, -1, 0, -1, -3, -1,
+ -2, -3, -3, -3, -2, 4, -3, 2, -1, -2, -1, 0, -2, -3, -1, -1, -2, -2, -1, 2, -1, 8,
+ -1, 2, -3, 1, 5, -4, -2, 0, -3, 1, -3, -1, 0, -1, 4, 0, 0, -1, -3, -2, -1, -2, 5};
+
+
+int blosum55mt[]={
+ 5,
+ -2, 5,
+ 0, -4, 13,
+ -2, 5, -4, 8,
+ -1, 1, -4, 2, 7,
+ -3, -5, -3, -5, -4, 9,
+ 0, -1, -3, -2, -3, -4, 8,
+ -2, 0, -4, -1, -1, -1, -2, 11,
+ -2, -4, -2, -4, -4, 0, -5, -4, 6,
+ -1, 0, -4, -1, 1, -4, -2, 0, -4, 6,
+ -2, -4, -2, -5, -4, 1, -5, -3, 2, -3, 6,
+ -1, -3, -2, -4, -3, 0, -3, -2, 2, -2, 3, 8,
+ -2, 4, -3, 2, 0, -4, 0, 1, -4, 0, -4, -3, 8,
+ -1, -2, -3, -2, -1, -5, -3, -3, -3, -1, -4, -3, -2, 10,
+ -1, 0, -4, 0, 2, -4, -2, 1, -4, 2, -3, 0, 0, -1, 7,
+ -2, -1, -4, -2, 0, -3, -3, 0, -4, 3, -3, -2, -1, -3, 1, 8,
+ 2, 0, -1, 0, 0, -3, 0, -1, -3, 0, -3, -2, 1, -1, 0, -1, 5,
+ 0, -1, -1, -1, -1, -3, -2, -2, -1, -1, -2, -1, 0, -1, -1, -1, 2, 6,
+ 0, -4, -1, -4, -3, -1, -4, -4, 4, -3, 1, 1, -4, -3, -3, -3, -2, 0, 5,
+ -4, -5, -4, -5, -3, 2, -3, -3, -3, -4, -3, -2, -5, -5, -2, -3, -4, -3, -4, 15,
+ -1, -1, -2, -2, -1, -2, -2, -1, -1, -1, -1, -1, -1, -2, -1, -1, -1, -1, -1, -3, -1,
+ -2, -3, -3, -3, -2, 4, -4, 2, -1, -2, -1, -1, -2, -4, -1, -2, -2, -2, -2, 3, -1, 9,
+ -1, 2, -4, 1, 5, -4, -3, 0, -4, 1, -3, -2, 0, -1, 4, 0, 0, -1, -3, -3, -1, -2, 5};
+
+int blosum62mt[]={
+ 4,
+ -2, 4,
+ 0, -3, 9,
+ -2, 4, -3, 6,
+ -1, 1, -4, 2, 5,
+ -2, -3, -2, -3, -3, 6,
+ 0, -1, -3, -1, -2, -3, 6,
+ -2, 0, -3, -1, 0, -1, -2, 8,
+ -1, -3, -1, -3, -3, 0, -4, -3, 4,
+ -1, 0, -3, -1, 1, -3, -2, -1, -3, 5,
+ -1, -4, -1, -4, -3, 0, -4, -3, 2, -2, 4,
+ -1, -3, -1, -3, -2, 0, -3, -2, 1, -1, 2, 5,
+ -2, 3, -3, 1, 0, -3, 0, 1, -3, 0, -3, -2, 6,
+ -1, -2, -3, -1, -1, -4, -2, -2, -3, -1, -3, -2, -2, 7,
+ -1, 0, -3, 0, 2, -3, -2, 0, -3, 1, -2, 0, 0, -1, 5,
+ -1, -1, -3, -2, 0, -3, -2, 0, -3, 2, -2, -1, 0, -2, 1, 5,
+ 1, 0, -1, 0, 0, -2, 0, -1, -2, 0, -2, -1, 1, -1, 0, -1, 4,
+ 0, -1, -1, -1, -1, -2, -2, -2, -1, -1, -1, -1, 0, -1, -1, -1, 1, 5,
+ 0, -3, -1, -3, -2, -1, -3, -3, 3, -2, 1, 1, -3, -2, -2, -3, -2, 0, 4,
+ -3, -4, -2, -4, -3, 1, -2, -2, -3, -3, -2, -1, -4, -4, -2, -3, -3, -2, -3, 11,
+ 0, -1, -2, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -2, -1, -1, 0, 0, -1, -2, -1,
+ -2, -3, -2, -3, -2, 3, -3, 2, -1, -2, -1, -1, -2, -3, -1, -2, -2, -2, -1, 2, -1, 7,
+ -1, 1, -3, 1, 4, -3, -2, 0, -3, 1, -3, -1, 0, -1, 3, 0, 0, -1, -2, -3, -1, -2, 4};
+
+int blosum62mt3[]={
+ 0,
+ -2, 0,
+ 0, -3, 0,
+ -2, 4, -3, 0,
+ -1, 1, -4, 2, 0,
+ -2, -3, -2, -3, -3, 0,
+ 0, -1, -3, -1, -2, -3, 0,
+ -2, 0, -3, -1, 0, -1, -2, 0,
+ -1, -3, -1, -3, -3, 0, -4, -3, 0,
+ -1, 0, -3, -1, 1, -3, -2, -1, -3, 0,
+ -1, -4, -1, -4, -3, 0, -4, -3, 2, -2, 0,
+ -1, -3, -1, -3, -2, 0, -3, -2, 1, -1, 2, 0,
+ -2, 3, -3, 1, 0, -3, 0, 1, -3, 0, -3, -2, 0,
+ -1, -2, -3, -1, -1, -4, -2, -2, -3, -1, -3, -2, -2, 0,
+ -1, 0, -3, 0, 2, -3, -2, 0, -3, 1, -2, 0, 0, -1, 0,
+ -1, -1, -3, -2, 0, -3, -2, 0, -3, 2, -2, -1, 0, -2, 1, 0,
+ 1, 0, -1, 0, 0, -2, 0, -1, -2, 0, -2, -1, 1, -1, 0, -1, 0,
+ 0, -1, -1, -1, -1, -2, -2, -2, -1, -1, -1, -1, 0, -1, -1, -1, 1, 0,
+ 0, -3, -1, -3, -2, -1, -3, -3, 3, -2, 1, 1, -3, -2, -2, -3, -2, 0, 0,
+ -3, -4, -2, -4, -3, 1, -2, -2, -3, -3, -2, -1, -4, -4, -2, -3, -3, -2, -3, 0,
+ 0, -1, -2, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -2, -1, -1, 0, 0, -1, -2, 0,
+ -2, -3, -2, -3, -2, 3, -3, 2, -1, -2, -1, -1, -2, -3, -1, -2, -2, -2, -1, 2, -1, 0,
+ -1, 1, -3, 1, 4, -3, -2, 0, -3, 1, -3, -1, 0, -1, 3, 0, 0, -1, -2, -3, -1, -2, 0};
+
+int blosum62mt2[]={
+ 0,
+ -2, 4,
+ 0, -3, 0,
+ -2, 4, -3, 2,
+ -1, 1, -4, 2, 2,
+ -2, -3, -2, -3, -3, 3,
+ 0, -1, -3, -1, -2, -3, 0,
+ -2, 0, -3, -1, 0, -1, -2, 2,
+ -1, -3, -1, -3, -3, 0, -4, -3, 2,
+ -1, 0, -3, -1, 1, -3, -2, -1, -3, 2,
+ -1, -4, -1, -4, -3, 0, -4, -3, 2, -2, 1,
+ -1, -3, -1, -3, -2, 0, -3, -2, 1, -1, 2, 0,
+ -2, 3, -3, 1, 0, -3, 0, 1, -3, 0, -3, -2, 0,
+ -1, -2, -3, -1, -1, -4, -2, -2, -3, -1, -3, -2, -2, 0,
+ -1, 0, -3, 0, 2, -3, -2, 0, -3, 1, -2, 0, 0, -1, 2,
+ -1, -1, -3, -2, 0, -3, -2, 0, -3, 2, -2, -1, 0, -2, 1, 0,
+ 1, 0, -1, 0, 0, -2, 0, -1, -2, 0, -2, -1, 1, -1, 0, -1, 1,
+ 0, -1, -1, -1, -1, -2, -2, -2, -1, -1, -1, -1, 0, -1, -1, -1, 1, 1,
+ 0, -3, -1, -3, -2, -1, -3, -3, 3, -2, 1, 1, -3, -2, -2, -3, -2, 0, 1,
+ -3, -4, -2, -4, -3, 1, -2, -2, -3, -3, -2, -1, -4, -4, -2, -3, -3, -2, -3, 1,
+ 0, -1, -2, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -2, -1, -1, 0, 0, -1, -2, -1,
+ -2, -3, -2, -3, -2, 3, -3, 2, -1, -2, -1, -1, -2, -3, -1, -2, -2, -2, -1, 2, -1, 3,
+ -1, 1, -3, 1, 4, -3, -2, 0, -3, 1, -3, -1, 0, -1, 3, 0, 0, -1, -2, -3, -1, -2, 4};
+
+
+int blosum80mt[]={
+ 7,
+ -3, 6,
+ -1, -6, 13,
+ -3, 6, -7, 10,
+ -2, 1, -7, 2, 8,
+ -4, -6, -4, -6, -6, 10,
+ 0, -2, -6, -3, -4, -6, 9,
+ -3, -1, -7, -2, 0, -2, -4, 12,
+ -3, -6, -2, -7, -6, -1, -7, -6, 7,
+ -1, -1, -6, -2, 1, -5, -3, -1, -5, 8,
+ -3, -7, -3, -7, -6, 0, -7, -5, 2, -4, 6,
+ -2, -5, -3, -6, -4, 0, -5, -4, 2, -3, 3, 9,
+ -3, 5, -5, 2, -1, -6, -1, 1, -6, 0, -6, -4, 9,
+ -1, -4, -6, -3, -2, -6, -5, -4, -5, -2, -5, -4, -4, 12,
+ -2, -1, -5, -1, 3, -5, -4, 1, -5, 2, -4, -1, 0, -3, 9,
+ -3, -2, -6, -3, -1, -5, -4, 0, -5, 3, -4, -3, -1, -3, 1, 9,
+ 2, 0, -2, -1, -1, -4, -1, -2, -4, -1, -4, -3, 1, -2, -1, -2, 7,
+ 0, -1, -2, -2, -2, -4, -3, -3, -2, -1, -3, -1, 0, -3, -1, -2, 2, 8,
+ -1, -6, -2, -6, -4, -2, -6, -5, 4, -4, 1, 1, -5, -4, -4, -4, -3, 0, 7,
+ -5, -8, -5, -8, -6, 0, -6, -4, -5, -6, -4, -3, -7, -7, -4, -5, -6, -5, -5, 16,
+ -1, -3, -4, -3, -2, -3, -3, -2, -2, -2, -2, -2, -2, -3, -2, -2, -1, -1, -2, -5, -2,
+ -4, -5, -5, -6, -5, 4, -6, 3, -3, -4, -2, -3, -4, -6, -3, -4, -3, -3, -3, 3, -3, 11,
+ -2, 0, -7, 1, 6, -6, -4, 0, -6, 1, -5, -3, -1, -2, 5, 0, -1, -2, -4, -5, -1, -4, 6};
+
+int pam120mt[]={
+ 3,
+ 0, 4,
+ -3, -6, 9,
+ 0, 4, -7, 5,
+ 0, 3, -7, 3, 5,
+ -4, -5, -6, -7, -7, 8,
+ 1, 0, -4, 0, -1, -5, 5,
+ -3, 1, -4, 0, -1, -3, -4, 7,
+ -1, -3, -3, -3, -3, 0, -4, -4, 6,
+ -2, 0, -7, -1, -1, -7, -3, -2, -3, 5,
+ -3, -4, -7, -5, -4, 0, -5, -3, 1, -4, 5,
+ -2, -4, -6, -4, -3, -1, -4, -4, 1, 0, 3, 8,
+ -1, 3, -5, 2, 1, -4, 0, 2, -2, 1, -4, -3, 4,
+ 1, -2, -4, -3, -2, -5, -2, -1, -3, -2, -3, -3, -2, 6,
+ -1, 0, -7, 1, 2, -6, -3, 3, -3, 0, -2, -1, 0, 0, 6,
+ -3, -2, -4, -3, -3, -5, -4, 1, -2, 2, -4, -1, -1, -1, 1, 6,
+ 1, 0, 0, 0, -1, -3, 1, -2, -2, -1, -4, -2, 1, 1, -2, -1, 3,
+ 1, 0, -3, -1, -2, -4, -1, -3, 0, -1, -3, -1, 0, -1, -2, -2, 2, 4,
+ 0, -3, -3, -3, -3, -3, -2, -3, 3, -4, 1, 1, -3, -2, -3, -3, -2, 0, 5,
+ -7, -6, -8, -8, -8, -1, -8, -3, -6, -5, -3, -6, -4, -7, -6, 1, -2, -6, -8, 12,
+ -1, -1, -4, -2, -1, -3, -2, -2, -1, -2, -2, -2, -1, -2, -1, -2, -1, -1, -1, -5, -2,
+ -4, -3, -1, -5, -5, 4, -6, -1, -2, -5, -2, -4, -2, -6, -5, -5, -3, -3, -3, -2, -3, 8,
+ -1, 2, -7, 3, 4, -6, -2, 1, -3, -1, -3, -2, 0, -1, 4, -1, -1, -2, -3, -7, -1, -5, 4};
+
+int pam160mt[]={
+ 2,
+ 0, 3,
+ -2, -4, 9,
+ 0, 3, -5, 4,
+ 0, 2, -5, 3, 4,
+ -3, -4, -5, -6, -5, 7,
+ 1, 0, -3, 0, 0, -4, 4,
+ -2, 1, -3, 0, 0, -2, -3, 6,
+ -1, -2, -2, -3, -2, 0, -3, -3, 5,
+ -2, 0, -5, 0, -1, -5, -2, -1, -2, 4,
+ -2, -4, -6, -4, -3, 1, -4, -2, 2, -3, 5,
+ -1, -3, -5, -3, -2, 0, -3, -3, 2, 0, 3, 7,
+ 0, 2, -4, 2, 1, -3, 0, 2, -2, 1, -3, -2, 3,
+ 1, -1, -3, -2, -1, -4, -1, -1, -2, -2, -3, -2, -1, 5,
+ -1, 1, -5, 1, 2, -5, -2, 2, -2, 0, -2, -1, 0, 0, 5,
+ -2, -1, -3, -2, -2, -4, -3, 1, -2, 3, -3, -1, -1, -1, 1, 6,
+ 1, 0, 0, 0, 0, -3, 1, -1, -2, -1, -3, -2, 1, 1, -1, -1, 2,
+ 1, 0, -2, -1, -1, -3, -1, -2, 0, 0, -2, -1, 0, 0, -1, -1, 1, 3,
+ 0, -2, -2, -3, -2, -2, -2, -2, 3, -3, 1, 1, -2, -2, -2, -3, -1, 0, 4,
+ -5, -5, -7, -6, -7, -1, -7, -3, -5, -4, -2, -4, -4, -5, -5, 1, -2, -5, -6, 12,
+ 0, -1, -3, -1, -1, -3, -1, -1, -1, -1, -2, -1, 0, -1, -1, -1, 0, 0, -1, -4, -1,
+ -3, -3, 0, -4, -4, 5, -5, 0, -2, -4, -2, -3, -2, -5, -4, -4, -3, -3, -3, -1, -3, 8,
+ 0, 2, -5, 2, 3, -5, -1, 1, -2, 0, -3, -2, 1, -1, 3, 0, -1, -1, -2, -6, -1, -4, 3};
+
+int pam250mt[]={
+ 2,
+ 0, 3,
+ -2, -4, 12,
+ 0, 3, -5, 4,
+ 0, 3, -5, 3, 4,
+ -3, -4, -4, -6, -5, 9,
+ 1, 0, -3, 1, 0, -5, 5,
+ -1, 1, -3, 1, 1, -2, -2, 6,
+ -1, -2, -2, -2, -2, 1, -3, -2, 5,
+ -1, 1, -5, 0, 0, -5, -2, 0, -2, 5,
+ -2, -3, -6, -4, -3, 2, -4, -2, 2, -3, 6,
+ -1, -2, -5, -3, -2, 0, -3, -2, 2, 0, 4, 6,
+ 0, 2, -4, 2, 1, -3, 0, 2, -2, 1, -3, -2, 2,
+ 1, -1, -3, -1, -1, -5, 0, 0, -2, -1, -3, -2, 0, 6,
+ 0, 1, -5, 2, 2, -5, -1, 3, -2, 1, -2, -1, 1, 0, 4,
+ -2, -1, -4, -1, -1, -4, -3, 2, -2, 3, -3, 0, 0, 0, 1, 6,
+ 1, 0, 0, 0, 0, -3, 1, -1, -1, 0, -3, -2, 1, 1, -1, 0, 2,
+ 1, 0, -2, 0, 0, -3, 0, -1, 0, 0, -2, -1, 0, 0, -1, -1, 1, 3,
+ 0, -2, -2, -2, -2, -1, -1, -2, 4, -2, 2, 2, -2, -1, -2, -2, -1, 0, 4,
+ -6, -5, -8, -7, -7, 0, -7, -3, -5, -3, -2, -4, -4, -6, -5, 2, -2, -5, -6, 17,
+ 0, -1, -3, -1, -1, -2, -1, -1, -1, -1, -1, -1, 0, -1, -1, -1, 0, 0, -1, -4, -1,
+ -3, -3, 0, -4, -4, 7, -5, 0, -1, -4, -1, -2, -2, -5, -4, -4, -3, -3, -2, 0, -2, 10,
+ 0, 2, -5, 3, 3, -5, 0, 2, -2, 0, -3, -2, 1, 0, 3, 0, 0, -1, -2, -6, -1, -4, 3};
+
+int pam350mt[]={
+ 2,
+ 1, 3,
+ -2, -5, 18,
+ 1, 3, -6, 4,
+ 1, 3, -6, 4, 4,
+ -4, -5, -5, -6, -6, 13,
+ 2, 1, -4, 1, 1, -6, 5,
+ -1, 1, -4, 1, 1, -2, -2, 7,
+ 0, -2, -3, -2, -2, 2, -2, -2, 5,
+ -1, 1, -6, 1, 0, -6, -1, 1, -2, 5,
+ -2, -4, -7, -4, -4, 3, -4, -2, 4, -3, 8,
+ -1, -2, -6, -3, -2, 1, -3, -2, 3, 0, 5, 6,
+ 0, 2, -4, 2, 2, -4, 1, 2, -2, 1, -3, -2, 2,
+ 1, 0, -3, 0, 0, -5, 0, 0, -2, -1, -3, -2, 0, 6,
+ 0, 2, -6, 2, 3, -5, -1, 3, -2, 1, -2, -1, 1, 1, 4,
+ -1, 0, -4, -1, 0, -5, -2, 2, -2, 4, -3, 0, 1, 0, 2, 7,
+ 1, 1, 0, 1, 0, -4, 1, -1, -1, 0, -3, -2, 1, 1, 0, 0, 1,
+ 1, 0, -2, 0, 0, -3, 1, -1, 0, 0, -2, -1, 1, 1, 0, -1, 1, 2,
+ 0, -2, -2, -2, -2, -1, -1, -2, 4, -2, 3, 2, -2, -1, -2, -3, -1, 0, 5,
+ -7, -6,-10, -8, -8, 1, -8, -3, -6, -4, -2, -5, -5, -7, -5, 4, -3, -6, -7, 27,
+ 0, 0, -3, -1, 0, -2, -1, 0, 0, -1, -1, 0, 0, 0, 0, -1, 0, 0, 0, -5, -1,
+ -4, -4, 1, -5, -5, 11, -6, 0, 0, -5, 0, -2, -3, -6, -5, -5, -3, -3, -2, 1, -2, 14,
+ 0, 2, -6, 3, 3, -6, 0, 2, -2, 1, -3, -2, 2, 0, 3, 1, 0, 0, -2, -7, 0, -5, 3};
+
+int md_40mt[]={
+ 9,
+ 0, 0,
+ -7, 0, 16,
+ -6, 0,-13, 11,
+ -5, 0,-15, 3, 11,
+-11, 0, -5,-15,-16, 13,
+ -3, 0, -7, -4, -4,-15, 10,
+ -9, 0, -6, -4, -8, -7,-10, 14,
+ -6, 0,-11,-12,-12, -5,-13,-11, 11,
+ -8, 0,-12, -8, -3,-16, -9, -6,-11, 11,
+ -9, 0,-10,-14,-13, -1,-14, -7, -1,-12, 9,
+ -6, 0, -9,-12,-11, -7,-12, -9, 1, -7, 1, 14,
+ -6, 0, -8, 1, -5,-12, -5, 0, -8, -1,-12, -9, 12,
+ -2, 0,-11,-11,-11,-11, -9, -4,-11,-10, -5,-10, -9, 12,
+ -7, 0,-12, -6, 0,-14, -9, 2,-12, -1, -6, -8, -5, -3, 12,
+ -7, 0, -5,-10, -8,-15, -4, 0,-10, 3, -9, -8, -6, -6, 0, 11,
+ 0, 0, -2, -6, -8, -6, -2, -6, -8, -7, -7, -8, 1, -1, -7, -5, 9,
+ 1, 0, -7, -8, -8,-11, -7, -7, -2, -5, -9, -2, -2, -4, -7, -6, 1, 10,
+ -1, 0, -7, -9, -8, -6, -8,-12, 4,-12, -2, 0,-10, -9,-11,-11, -7, -4, 10,
+-14, 0, -4,-15,-15, -7, -7,-13,-13,-13, -8,-11,-14,-14,-11, -4, -9,-12,-10, 18,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+-13, 0, -2, -8,-14, 2,-13, 2, -9,-13, -9,-11, -6,-13, -9,-10, -7,-10,-11, -6, 0, 14,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
+
+int md_120mt[]={
+ 6,
+ 0, 0,
+ -3, 0, 14,
+ -2, 0, -7, 8,
+ -2, 0, -8, 5, 8,
+ -6, 0, -2, -9,-10, 11,
+ 0, 0, -3, 0, -1, -9, 8,
+ -4, 0, -2, -1, -3, -2, -4, 11,
+ -1, 0, -5, -7, -7, -1, -6, -6, 7,
+ -4, 0, -6, -2, 0, -9, -4, -1, -6, 8,
+ -4, 0, -5, -8, -8, 2, -8, -4, 2, -6, 7,
+ -2, 0, -5, -7, -6, -2, -6, -5, 3, -4, 3, 10,
+ -1, 0, -3, 3, -1, -6, -1, 2, -4, 1, -6, -5, 8,
+ 0, 0, -5, -5, -5, -5, -4, -1, -5, -4, -2, -5, -3, 9,
+ -3, 0, -6, -1, 2, -7, -4, 4, -6, 2, -3, -4, -1, 0, 9,
+ -3, 0, -2, -4, -3, -8, -1, 2, -6, 4, -5, -4, -2, -2, 2, 8,
+ 2, 0, 0, -2, -3, -3, 0, -2, -3, -3, -3, -3, 2, 1, -3, -2, 5,
+ 2, 0, -3, -3, -4, -6, -2, -3, 0, -2, -4, 0, 1, 0, -3, -3, 2, 6,
+ 1, 0, -3, -5, -5, -2, -4, -6, 5, -6, 1, 2, -5, -4, -6, -6, -3, 0, 7,
+ -8, 0, 0, -9, -9, -3, -3, -6, -7, -6, -4, -6, -8, -8, -6, -1, -5, -7, -6, 17,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+ -7, 0, 2, -4, -7, 5, -8, 4, -5, -7, -4, -6, -2, -7, -4, -5, -3, -6, -6, -2, 0, 12,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
+
+int md_250mt[]={
+ 2,
+ 0, 0,
+ -1, 0, 11,
+ -1, 0, -3, 5,
+ -1, 0, -4, 4, 5,
+ -3, 0, 0, -5, -5, 8,
+ 1, 0, -1, 1, 1, -5, 5,
+ -2, 0, 0, 0, 0, 0, -2, 6,
+ 0, 0, -2, -3, -3, 0, -3, -3, 4,
+ -1, 0, -3, 0, 1, -5, -1, 1, -3, 5,
+ -1, 0, -2, -4, -4, 2, -4, -2, 2, -3, 5,
+ 0, 0, -2, -3, -3, 0, -3, -2, 3, -2, 3, 6,
+ 0, 0, -1, 2, 1, -3, 0, 1, -2, 1, -3, -2, 3,
+ 1, 0, -2, -2, -2, -2, -1, 0, -2, -1, 0, -2, -1, 6,
+ -1, 0, -3, 0, 2, -4, -1, 3, -3, 2, -2, -2, 0, 0, 5,
+ -1, 0, -1, -1, 0, -4, 0, 2, -3, 4, -3, -2, 0, -1, 2, 5,
+ 1, 0, 1, 0, -1, -2, 1, -1, -1, -1, -2, -1, 1, 1, -1, -1, 2,
+ 2, 0, -1, -1, -1, -2, 0, -1, 1, -1, -1, 0, 1, 1, -1, -1, 1, 2,
+ 1, 0, -2, -3, -2, 0, -2, -3, 4, -3, 2, 2, -2, -1, -3, -3, -1, 0, 4,
+ -4, 0, 1, -5, -5, -1, -1, -3, -4, -3, -2, -3, -4, -4, -3, 0, -3, -4, -3, 15,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+ -3, 0, 2, -2, -4, 5, -4, 4, -2, -3, -1, -3, -1, -3, -2, -2, -1, -3, -3, 0, 0, 9,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
+
+int md_350mt[]={
+ 1,
+ 0, 0,
+ 0, 0, 9,
+ 0, 0, -2, 3,
+ 0, 0, -2, 3, 3,
+ -2, 0, 1, -3, -4, 6,
+ 1, 0, 0, 1, 1, -3, 4,
+ -1, 0, 0, 0, 0, 0, -1, 3,
+ 0, 0, -1, -2, -2, 1, -2, -2, 3,
+ -1, 0, -1, 0, 1, -3, 0, 1, -2, 3,
+ -1, 0, -1, -3, -3, 2, -2, -1, 2, -2, 3,
+ 0, 0, -1, -2, -2, 1, -2, -1, 2, -2, 2, 3,
+ 0, 0, -1, 1, 1, -2, 0, 1, -1, 1, -2, -1, 2,
+ 1, 0, -1, -1, -1, -2, -1, 0, -1, -1, 0, -1, 0, 4,
+ -1, 0, -2, 1, 1, -2, 0, 2, -2, 2, -1, -1, 0, 0, 3,
+ -1, 0, 0, 0, 0, -3, 0, 1, -2, 3, -2, -1, 0, 0, 2, 3,
+ 1, 0, 0, 0, 0, -1, 1, 0, -1, 0, -1, -1, 1, 1, 0, 0, 1,
+ 1, 0, 0, 0, -1, -1, 0, -1, 0, 0, -1, 0, 0, 1, -1, 0, 1, 1,
+ 0, 0, -1, -2, -2, 0, -1, -2, 2, -2, 1, 2, -1, -1, -2, -2, 0, 0, 2,
+ -3, 0, 1, -4, -3, 0, -1, -2, -3, -2, -1, -2, -3, -3, -2, 0, -2, -3, -2, 14,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+ -2, 0, 2, -2, -2, 5, -3, 3, -1, -2, 0, -1, -1, -2, -1, -1, -1, -2, -2, 0, 0, 7,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
+
+int idmat[]={
+10,
+ 0, 10,
+ 0, 0, 10,
+ 0, 0, 0, 10,
+ 0, 0, 0, 0, 10,
+ 0, 0, 0, 0, 0, 10,
+ 0, 0, 0, 0, 0, 0, 10,
+ 0, 0, 0, 0, 0, 0, 0, 10,
+ 0, 0, 0, 0, 0, 0, 0, 0, 10,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 10,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 10,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 10,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 10,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 10,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 10,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 10,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 10,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 10,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,10,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,10,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,10,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,10,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,10};
+
+int dna_idmat[]={
+10,
+ -1, 10,
+ -1, -1, 10,
+ -1, -1, -1, 10,
+ -1, -1, -1, -1, 10,
+ -1, -1, -1, -1, -1, 10,
+ -1, -1, -1, -1, -1, -1, 10,
+ -1, -1, -1, -1, -1, -1, -1, 10,
+ -1, -1, -1, -1, -1, -1, -1, -1, 10,
+ -1, -1, -1, -1, -1, -1, -1, -1, -1, 10,
+ -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 10,
+ -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 10,
+ -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 10,
+ -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 10,
+ -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 10,
+ -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 10,
+ -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 10,
+ -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 10,
+ -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,10,
+ -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,10,
+ -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,10,
+ -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,10,
+ -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,10};
+
+int est_idmat[]={
+10,
+ -10, 10,
+ -10, -10, 10,
+ -10, -10, -10, 10,
+ -10, -10, -10, -10, 10,
+ -10, -10, -10, -10, -10, 10,
+ -10, -10, -10, -10, -10, -10, 10,
+ -10, -10, -10, -10, -10, -10, -10, 10,
+ -10, -10, -10, -10, -10, -10, -10, -10, 10,
+ -10, -10, -10, -10, -10, -10, -10, -10, -10, 10,
+ -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, 10,
+ -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, 10,
+ -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, 10,
+ -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, 10,
+ -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, 10,
+ -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, 10,
+ -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, 10,
+ -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, 10,
+ -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10,10,
+ -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10,10,
+ -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10,10,
+ -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10,10,
+ -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10,10};
+
+/*These are the three structure specific matrices descibed by Luthy, R., McLachlan, A.D. and Eisenberg, D.
+ in Proteins 10, 229-239 (1991), taken from http://www.genome.ad.jp/dbget/dbget2.html
+ ID: coil :LUTR910107
+ ID: alpha:LUTR910108
+ ID: beta :LUTR910109
+*/
+int coil_mat []={
+ 11,
+ 0, 0,
+ 3, 0, 108,
+ 2, 0, -5, 27,
+ 2, 0, 7, -9, 15,
+-15, 0, 3, -11, -4, 83,
+ 1, 0, 4, -15, -4, 2, 41,
+ 5, 0, 4, 108, -6, -13, -3, 55,
+ 27, 0, -11, 3, -2, -14, 0, -1, 49,
+ -9, 0, 3, 4, -21, -2, -3, -4, -7, 38,
+ -5, 0, 14, 7, -5, 2, 2, -6, -3, 1, 48,
+ -9, 0, 3, 4, -6, -16, -23, -9, -7, -1, 7, 52,
+ 28, 0, 1, 2, 27, 7, -5, -9, -9, 108, -15, 3, 10,
+108, 0, 8, 14, 41, 0, -2, -4, -10, -3, 1, -15, 4, 58,
+ 10, 0, 3, -9, 3, -4, 3, 3, 8, 6, -15, 15, 5, -6, 14,
+ 1, 0, 10, 3, 1, 108, 5, 27, -5, -9, -9, -15, 2, 3, 2, 28,
+ 3, 0, -15, 3, -3, -8, -12, -5, -6, -9, -19, -14, 4, 4, -2, 7, 19,
+ 7, 0, 6, 3, 0, -7, 55, -18, -8, -7, 6, -20, -11, 2, -5, 4, 3, 31,
+-11, 0, -6, -15, -23, -6, -4, -3, 19, -15, -6, -22, 3, -11, -4, 14, -1, -7, 37,
+ 4, 0, 15, 3, 2, -8, -1, -12, -22, -2, 5, -19, 14, 4, -21, 4, 0, -9, -1, 129,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+ 4, 0, -4, 8, -3, -30, -6, 49, -10, -7, -12, -18, 3, 3, -6, -11, 3, 2, 2, -5, 0, 59,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 };
+
+int alpha_mat []={
+ 23,
+ 0, 0,
+ 5, 0, 51,
+ 1, 0, 4, 28,
+ 0, 0, 5, 1, 23,
+ -4, 0, 5, 1, 3, 88,
+ 3, 0, 5, -4, 2, 4, 37,
+ 7, 0, 6, 51, 3, 5, -6, 52,
+ 28, 0, 1, 1, 6, -21, 2, -6, 38,
+ 1, 0, 0, 5, 4, -2, 3, -7, -10, 34,
+ 4, 0, 11, 5, 4, 3, 5, -6, -19, -6, 45,
+ 1, 0, 4, 6, 3, -23, -5, -12, -14, -4, 9, 40,
+ 19, 0, 3, 0, 28, 5, 4, 1, 1, 51, -4, 1, 13,
+ 51, 0, 10, 11, 37, 1, 1, -6, -8, -5, 2, -22, 5, 59,
+ 13, 0, 1, 1, 0, 2, 4, 5, 10, 7, -5, 23, 7, 3, 11,
+ -1, 0, 13, 5, 3, 51, 7, 28, 4, 1, 1, -4, 1, 1, 0, 19,
+ 1, 0, -5, 0, -6, 2, -1, -10, -14, -9, -6, -22, 6, 3, 6, 5, 17,
+ 5, 0, 7, 4, 2, 3, 52, -8, -18, 3, 4, -24, 1, 4, 4, 5, 2, 14,
+ 1, 0, 3, -5, -5, 5, -7, -19, 7, -6, -2, -28, 4, 7, 3, 11, 2, -38, 40,
+ 5, 0, 23, 5, 5, 5, -6, -9, -19, -2, 3, -4, 11, 5, 4, 6, 2, -43, -8, 162,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+ 6, 0, 2, 10, 3, -2, -6, 38, -14, -7, -1, -19, 0, 5, 3, 1, 4, -29, -6, -19, 0, 97,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 };
+
+int beta_mat []={
+ 23,
+ 0, 0,
+ 4, 0, 97,
+ 1, 0, -5, 49,
+ -1, 0, 9, -5, 31,
+-16, 0, 8, -11, 6, 57,
+ 2, 0, 7, -16, 12, 8, 54,
+ 13, 0, 9, 97, 2, -30, -4, 55,
+ 49, 0, -11, 0, 0, -22, 2, -7, 28,
+-19, 0, 0, 7, -18, 2, 17, -12, -10, 34,
+ -5, 0, 31, 9, 1, 7, 13, -9, -7, -3, 34,
+ -5, 0, 2, 9, 4, -9, -12, -17, -13, 1, 14, 25,
+ 40, 0, 2, -1, 49, 9, -5, -19, -5, 97, -16, 0, 15,
+ 97, 0, 17, 31, 54, 5, 5, -13, -34, -2, 8, -11, 7, 78,
+ 15, 0, 0, -19, 0, 12, 2, 8, 17, 11, -21, 31, 13, 2, 31,
+ -2, 0, 15, 4, 2, 97, 13, 49, -5, -5, -19, -16, 1, 0, -1, 40,
+ 0, 0, -21, 0, -4, -6, -5, -13, -12, -5, -23, -9, 9, 7, 0, 9, 26,
+ 9, 0, 11, 2, 2, 2, 55, -11, -14, -28, 12, -21, -11, 1, 1, 7, 4, 21,
+-11, 0, 2, -21, -12, -2, -12, -7, 13, 2, -1, -13, 2, -10, 6, 31, 2, -22, 23,
+ 7, 0, 31, 8, 13, -7, -7, -15, -15, -3, 10, -31, 31, 6, -18, 9, 6, -17, -5, 133,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+ 9, 0, 12, 17, 17, -12, -9, 28, -13, -3, 0, -12, 0, 8, 4, -11, 2, -5, -5, -5, 0, 47,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 };
+/*EXPERIMENTAL MATRICES: Blosum Series*/
+/*Blosum Series rescaled on blosum62mt,
+Lambda Blosum30mt=17
+Lambda Blosum62mt=10
+Lambda Blosum80mt=4.5
+*/
+
+int rblosum30mt []={
+ 52,
+ 33, 56,
+ 18, 23, 113,
+ 33, 56, 18, 75,
+ 33, 33, 37, 37, 61,
+ 23, 18, 18, 9, 14, 80,
+ 33, 33, 14, 28, 23, 18, 71,
+ 23, 23, 9, 23, 33, 18, 18, 99,
+ 33, 23, 23, 14, 18, 33, 28, 23, 61,
+ 33, 33, 18, 33, 42, 28, 28, 23, 23, 52,
+ 28, 28, 33, 28, 28, 42, 23, 28, 42, 23, 52,
+ 37, 23, 23, 18, 28, 23, 23, 42, 37, 42, 42, 61,
+ 33, 52, 28, 37, 28, 28, 33, 28, 33, 33, 23, 33, 71,
+ 28, 23, 18, 28, 37, 14, 28, 37, 18, 37, 18, 14, 18, 85,
+ 37, 28, 23, 28, 42, 18, 23, 33, 23, 33, 23, 28, 28, 33, 71,
+ 28, 23, 23, 28, 28, 28, 23, 28, 18, 37, 23, 33, 23, 28, 47, 71,
+ 37, 33, 23, 33, 33, 28, 33, 28, 28, 33, 23, 23, 33, 28, 28, 28, 52,
+ 37, 33, 23, 28, 23, 23, 23, 23, 33, 28, 33, 33, 37, 33, 33, 18, 42, 56,
+ 37, 23, 23, 23, 18, 37, 18, 18, 52, 23, 37, 33, 23, 14, 18, 28, 28, 37, 56,
+ 9, 9, 23, 14, 28, 37, 37, 9, 18, 23, 23, 18, 0, 18, 28, 33, 18, 9, 18, 128,
+ 0, -1, -2, -1, -1, -1, -1, -1, 0, 0, 0, 0, 0, -1, 0, -1, 0, 0, 0, -2, -1,
+ -4, -3, -6, -1, -2, 3, -3, 0, -1, -1, 3, -1, -4, -2, -1, 0, -2, -1, 1, 5, -1, 9,
+ 0, 0, 0, 0, 5, -4, -2, 0, -3, 1, -1, -1, -1, 0, 4, 0, -1, -1, -3, -1, 0, -2, 4,
+
+
+};
+
+
+
+int rblosum62mt[]={
+ 80,
+ 20, 80,
+ 40, 10, 130,
+ 20, 80, 10, 100,
+ 30, 50, 0, 60, 90,
+ 20, 10, 20, 10, 10, 100,
+ 40, 30, 10, 30, 20, 10, 100,
+ 20, 40, 10, 30, 40, 30, 20, 120,
+ 30, 10, 30, 10, 10, 40, 0, 10, 80,
+ 30, 40, 10, 30, 50, 10, 20, 30, 10, 90,
+ 30, 0, 30, 0, 10, 40, 0, 10, 60, 20, 80,
+ 30, 10, 30, 10, 20, 40, 10, 20, 50, 30, 60, 90,
+ 20, 70, 10, 50, 40, 10, 40, 50, 10, 40, 10, 20, 100,
+ 30, 20, 10, 30, 30, 0, 20, 20, 10, 30, 10, 20, 20, 110,
+ 30, 40, 10, 40, 60, 10, 20, 40, 10, 50, 20, 40, 40, 30, 90,
+ 30, 30, 10, 20, 40, 10, 20, 40, 10, 60, 20, 30, 40, 20, 50, 90,
+ 50, 40, 30, 40, 40, 20, 40, 30, 20, 40, 20, 30, 50, 30, 40, 30, 80,
+ 40, 30, 30, 30, 30, 20, 20, 20, 30, 30, 30, 30, 40, 30, 30, 30, 50, 90,
+ 40, 10, 30, 10, 20, 30, 10, 10, 70, 20, 50, 50, 10, 20, 20, 10, 20, 40, 80,
+ 10, 0, 20, 0, 10, 50, 20, 20, 10, 10, 20, 30, 0, 0, 20, 10, 10, 20, 10, 150,
+ 0, -1, -2, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -2, -1, -1, 0, 0, -1, -2, -1,
+ -2, -3, -2, -3, -2, 3, -3, 2, -1, -2, -1, -1, -2, -3, -1, -2, -2, -2, -1, 2, -1, 7,
+ -1, 1, -3, 1, 4, -3, -2, 0, -3, 1, -3, -1, 0, -1, 3, 0, 0, -1, -2, -3, -1, -2, 4,
+
+ };
+int rblosum80mt[]={
+77,
+ 25, 72,
+ 36, 10, 108,
+ 25, 72, 5, 93,
+ 31, 46, 5, 51, 82,
+ 20, 10, 20, 10, 10, 93,
+ 41, 31, 10, 25, 20, 10, 87,
+ 25, 36, 5, 31, 41, 31, 20, 103,
+ 25, 10, 31, 5, 10, 36, 5, 10, 77,
+ 36, 36, 10, 31, 46, 15, 25, 36, 15, 82,
+ 25, 5, 25, 5, 10, 41, 5, 15, 51, 20, 72,
+ 31, 15, 25, 10, 20, 41, 15, 20, 51, 25, 56, 87,
+ 25, 67, 15, 51, 36, 10, 36, 46, 10, 41, 10, 20, 87,
+ 36, 20, 10, 25, 31, 10, 15, 20, 15, 31, 15, 20, 20, 103,
+ 31, 36, 15, 36, 56, 15, 20, 46, 15, 51, 20, 36, 41, 25, 87,
+ 25, 31, 10, 25, 36, 15, 20, 41, 15, 56, 20, 25, 36, 25, 46, 87,
+ 51, 41, 31, 36, 36, 20, 36, 31, 20, 36, 20, 25, 46, 31, 36, 31, 77,
+ 41, 36, 31, 31, 31, 20, 25, 25, 31, 36, 25, 36, 41, 25, 36, 31, 51, 82,
+ 36, 10, 31, 10, 20, 31, 10, 15, 62, 20, 46, 46, 15, 20, 20, 20, 25, 41, 77,
+ 15, 0, 15, 0, 10, 41, 10, 20, 15, 10, 20, 25, 5, 5, 20, 15, 10, 15, 15, 124,
+ -1, -3, -4, -3, -2, -3, -3, -2, -2, -2, -2, -2, -2, -3, -2, -2, -1, -1, -2, -5, -2,
+ -4, -5, -5, -6, -5, 4, -6, 3, -3, -4, -2, -3, -4, -6, -3, -4, -3, -3, -3, 3, -3, 11,
+ -2, 0, -7, 1, 6, -6, -4, 0, -6, 1, -5, -3, -1, -2, 5, 0, -1, -2, -4, -5, -1, -4, 6,
+};
+int rpam120mt[]={
+ 52,
+ 38, 57,
+ 24, 9, 81,
+ 38, 57, 4, 62,
+ 38, 52, 4, 52, 62,
+ 19, 14, 9, 4, 4, 77,
+ 43, 38, 19, 38, 33, 14, 62,
+ 24, 43, 19, 38, 33, 24, 19, 72,
+ 33, 24, 24, 24, 24, 38, 19, 19, 67,
+ 28, 38, 4, 33, 33, 4, 24, 28, 24, 62,
+ 24, 19, 4, 14, 19, 38, 14, 24, 43, 19, 62,
+ 28, 19, 9, 19, 24, 33, 19, 19, 43, 38, 52, 77,
+ 33, 52, 14, 48, 43, 19, 38, 48, 28, 43, 19, 24, 57,
+ 43, 28, 19, 24, 28, 14, 28, 33, 24, 28, 24, 24, 28, 67,
+ 33, 38, 4, 43, 48, 9, 24, 52, 24, 38, 28, 33, 38, 38, 67,
+ 24, 28, 19, 24, 24, 14, 19, 43, 28, 48, 19, 33, 33, 33, 43, 67,
+ 43, 38, 38, 38, 33, 24, 43, 28, 28, 33, 19, 28, 43, 43, 28, 33, 52,
+ 43, 38, 24, 33, 28, 19, 33, 24, 38, 33, 24, 33, 38, 33, 28, 28, 48, 57,
+ 38, 24, 24, 24, 24, 24, 28, 24, 52, 19, 43, 43, 24, 28, 24, 24, 28, 38, 62,
+ 4, 9, 0, 0, 0, 33, 0, 24, 9, 14, 24, 9, 19, 4, 9, 43, 28, 9, 0, 96,
+ -1, -1, -4, -2, -1, -3, -2, -2, -1, -2, -2, -2, -1, -2, -1, -2, -1, -1, -1, -5, -2,
+ -4, -3, -1, -5, -5, 4, -6, -1, -2, -5, -2, -4, -2, -6, -5, -5, -3, -3, -3, -2, -3, 8,
+ -1, 2, -7, 3, 4, -6, -2, 1, -3, -1, -3, -2, 0, -1, 4, -1, -1, -2, -3, -7, -1, -5, 4,
+};
+int rpam160mt[]={
+ 47,
+ 36, 52,
+ 26, 15, 83,
+ 36, 52, 10, 57,
+ 36, 47, 10, 52, 57,
+ 20, 15, 10, 5, 10, 73,
+ 41, 36, 20, 36, 36, 15, 57,
+ 26, 41, 20, 36, 36, 26, 20, 68,
+ 31, 26, 26, 20, 26, 36, 20, 20, 62,
+ 26, 36, 10, 36, 31, 10, 26, 31, 26, 57,
+ 26, 15, 5, 15, 20, 41, 15, 26, 47, 20, 62,
+ 31, 20, 10, 20, 26, 36, 20, 20, 47, 36, 52, 73,
+ 36, 47, 15, 47, 41, 20, 36, 47, 26, 41, 20, 26, 52,
+ 41, 31, 20, 26, 31, 15, 31, 31, 26, 26, 20, 26, 31, 62,
+ 31, 41, 10, 41, 47, 10, 26, 47, 26, 36, 26, 31, 36, 36, 62,
+ 26, 31, 20, 26, 26, 15, 20, 41, 26, 52, 20, 31, 31, 31, 41, 68,
+ 41, 36, 36, 36, 36, 20, 41, 31, 26, 31, 20, 26, 41, 41, 31, 31, 47,
+ 41, 36, 26, 31, 31, 20, 31, 26, 36, 36, 26, 31, 36, 36, 31, 31, 41, 52,
+ 36, 26, 26, 20, 26, 26, 26, 26, 52, 20, 41, 41, 26, 26, 26, 20, 31, 36, 57,
+ 10, 10, 0, 5, 0, 31, 0, 20, 10, 15, 26, 15, 15, 10, 10, 41, 26, 10, 5, 99,
+ 0, -1, -3, -1, -1, -3, -1, -1, -1, -1, -2, -1, 0, -1, -1, -1, 0, 0, -1, -4, -1,
+ -3, -3, 0, -4, -4, 5, -5, 0, -2, -4, -2, -3, -2, -5, -4, -4, -3, -3, -3, -1, -3, 8,
+ 0, 2, -5, 2, 3, -5, -1, 1, -2, 0, -3, -2, 1, -1, 3, 0, -1, -1, -2, -6, -1, -4, 3,
+};
+int rpam250mt[]={
+ 42,
+ 34, 47,
+ 25, 17, 85,
+ 34, 47, 12, 51,
+ 34, 47, 12, 47, 51,
+ 21, 17, 17, 8, 12, 72,
+ 38, 34, 21, 38, 34, 12, 55,
+ 30, 38, 21, 38, 38, 25, 25, 60,
+ 30, 25, 25, 25, 25, 38, 21, 25, 55,
+ 30, 38, 12, 34, 34, 12, 25, 34, 25, 55,
+ 25, 21, 8, 17, 21, 42, 17, 25, 42, 21, 60,
+ 30, 25, 12, 21, 25, 34, 21, 25, 42, 34, 51, 60,
+ 34, 42, 17, 42, 38, 21, 34, 42, 25, 38, 21, 25, 42,
+ 38, 30, 21, 30, 30, 12, 34, 34, 25, 30, 21, 25, 34, 60,
+ 34, 38, 12, 42, 42, 12, 30, 47, 25, 38, 25, 30, 38, 34, 51,
+ 25, 30, 17, 30, 30, 17, 21, 42, 25, 47, 21, 34, 34, 34, 38, 60,
+ 38, 34, 34, 34, 34, 21, 38, 30, 30, 34, 21, 25, 38, 38, 30, 34, 42,
+ 38, 34, 25, 34, 34, 21, 34, 30, 34, 34, 25, 30, 34, 34, 30, 30, 38, 47,
+ 34, 25, 25, 25, 25, 30, 30, 25, 51, 25, 42, 42, 25, 30, 25, 25, 30, 34, 51,
+ 8, 12, 0, 4, 4, 34, 4, 21, 12, 21, 25, 17, 17, 8, 12, 42, 25, 12, 8, 107,
+ 0, -1, -3, -1, -1, -2, -1, -1, -1, -1, -1, -1, 0, -1, -1, -1, 0, 0, -1, -4, -1,
+ -3, -3, 0, -4, -4, 7, -5, 0, -1, -4, -1, -2, -2, -5, -4, -4, -3, -3, -2, 0, -2, 10,
+ 0, 2, -5, 3, 3, -5, 0, 2, -2, 0, -3, -2, 1, 0, 3, 0, 0, -1, -2, -6, -1, -4, 3,
+};
+int rpam350mt[]={
+ 39,
+ 36, 43,
+ 26, 16, 92,
+ 36, 43, 13, 46,
+ 36, 43, 13, 46, 46,
+ 19, 16, 16, 13, 13, 76,
+ 39, 36, 19, 36, 36, 13, 49,
+ 29, 36, 19, 36, 36, 26, 26, 56,
+ 33, 26, 23, 26, 26, 39, 26, 26, 49,
+ 29, 36, 13, 36, 33, 13, 29, 36, 26, 49,
+ 26, 19, 9, 19, 19, 43, 19, 26, 46, 23, 59,
+ 29, 26, 13, 23, 26, 36, 23, 26, 43, 33, 49, 53,
+ 33, 39, 19, 39, 39, 19, 36, 39, 26, 36, 23, 26, 39,
+ 36, 33, 23, 33, 33, 16, 33, 33, 26, 29, 23, 26, 33, 53,
+ 33, 39, 13, 39, 43, 16, 29, 43, 26, 36, 26, 29, 36, 36, 46,
+ 29, 33, 19, 29, 33, 16, 26, 39, 26, 46, 23, 33, 36, 33, 39, 56,
+ 36, 36, 33, 36, 33, 19, 36, 29, 29, 33, 23, 26, 36, 36, 33, 33, 36,
+ 36, 33, 26, 33, 33, 23, 36, 29, 33, 33, 26, 29, 36, 36, 33, 29, 36, 39,
+ 33, 26, 26, 26, 26, 29, 29, 26, 46, 26, 43, 39, 26, 29, 26, 23, 29, 33, 49,
+ 9, 13, 0, 6, 6, 36, 6, 23, 13, 19, 26, 16, 16, 9, 16, 46, 23, 13, 9, 122,
+ 0, 0, -3, -1, 0, -2, -1, 0, 0, -1, -1, 0, 0, 0, 0, -1, 0, 0, 0, -5, -1,
+ -4, -4, 1, -5, -5, 11, -6, 0, 0, -5, 0, -2, -3, -6, -5, -5, -3, -3, -2, 1, -2, 14,
+ 0, 2, -6, 3, 3, -6, 0, 2, -2, 1, -3, -2, 2, 0, 3, 1, 0, 0, -2, -7, 0, -5, 3,
+};
+int ralpha_mat []={
+ 42,
+ 27, 27,
+ 31, 27, 60,
+ 28, 27, 30, 46,
+ 27, 27, 31, 28, 42,
+ 25, 27, 31, 28, 29, 84,
+ 29, 27, 31, 25, 29, 30, 51,
+ 32, 27, 31, 60, 29, 31, 23, 61,
+ 46, 27, 28, 28, 31, 14, 29, 23, 52,
+ 28, 27, 27, 31, 30, 26, 29, 23, 21, 49,
+ 30, 27, 35, 31, 30, 29, 31, 23, 15, 23, 57,
+ 28, 27, 30, 31, 29, 12, 24, 20, 18, 25, 33, 53,
+ 40, 27, 29, 27, 46, 31, 30, 28, 28, 60, 25, 28, 36,
+ 60, 27, 34, 35, 51, 28, 28, 23, 22, 24, 29, 13, 31, 66,
+ 36, 27, 28, 28, 27, 29, 30, 31, 34, 32, 24, 42, 32, 29, 35,
+ 27, 27, 36, 31, 29, 60, 32, 46, 30, 28, 28, 25, 28, 28, 27, 40,
+ 28, 27, 24, 27, 23, 29, 27, 21, 18, 22, 23, 13, 31, 29, 31, 31, 38,
+ 31, 27, 32, 30, 29, 29, 61, 22, 16, 29, 30, 12, 28, 30, 30, 31, 29, 36,
+ 28, 27, 29, 24, 24, 31, 23, 15, 32, 23, 26, 9, 30, 32, 29, 35, 29, 3, 53,
+ 31, 27, 42, 31, 31, 31, 23, 22, 15, 26, 29, 25, 35, 31, 30, 31, 29, 0, 22, 132,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+ 6, 0, 2, 10, 3, -2, -6, 38, -14, -7, -1, -19, 0, 5, 3, 1, 4, -29, -6, -19, 0, 97,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+};
+
+int rbeta_mat[]={
+ 44,
+ 26, 26,
+ 29, 26, 102,
+ 27, 26, 22, 64,
+ 25, 26, 33, 22, 50,
+ 14, 26, 32, 17, 31, 70,
+ 28, 26, 31, 14, 35, 32, 68,
+ 36, 26, 33, 102, 28, 3, 23, 69,
+ 64, 26, 17, 26, 26, 9, 28, 21, 48,
+ 11, 26, 26, 31, 12, 28, 39, 17, 18, 52,
+ 22, 26, 50, 33, 27, 31, 36, 19, 21, 24, 52,
+ 22, 26, 28, 33, 29, 19, 17, 13, 16, 27, 37, 45,
+ 57, 26, 28, 25, 64, 33, 22, 11, 22, 102, 14, 26, 38,
+102, 26, 39, 50, 68, 30, 30, 16, 0, 24, 32, 17, 31, 87,
+ 38, 26, 26, 11, 26, 35, 28, 32, 39, 35, 10, 50, 36, 28, 50,
+ 24, 26, 38, 29, 28, 102, 36, 64, 22, 22, 11, 14, 27, 26, 25, 57,
+ 26, 26, 10, 26, 23, 21, 22, 16, 17, 22, 8, 19, 33, 31, 26, 33, 46,
+ 33, 26, 35, 28, 28, 28, 69, 17, 15, 4, 35, 10, 17, 27, 27, 31, 29, 42,
+ 17, 26, 28, 10, 17, 24, 17, 21, 36, 28, 25, 16, 28, 18, 31, 50, 28, 9, 44,
+ 31, 26, 50, 32, 36, 21, 21, 14, 14, 24, 34, 2, 50, 31, 12, 33, 31, 13, 22, 130,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+ 9, 0, 12, 17, 17, -12, -9, 28, -13, -3, 0, -12, 0, 8, 4, -11, 2, -5, -5, -5, 0, 47,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+};
+int rcoil_mat[]={
+ 36,
+ 27, 27,
+ 29, 27, 124,
+ 28, 27, 22, 51,
+ 28, 27, 33, 18, 40,
+ 13, 27, 29, 17, 23, 101,
+ 27, 27, 30, 13, 23, 28, 64,
+ 31, 27, 30, 124, 21, 15, 24, 76,
+ 51, 27, 17, 29, 25, 14, 27, 26, 71,
+ 18, 27, 29, 30, 8, 25, 24, 23, 20, 61,
+ 22, 27, 39, 33, 22, 28, 28, 21, 24, 27, 70,
+ 18, 27, 29, 30, 21, 12, 6, 18, 20, 26, 33, 73,
+ 52, 27, 27, 28, 51, 33, 22, 18, 18, 124, 13, 29, 36,
+124, 27, 34, 39, 64, 27, 25, 23, 18, 24, 27, 13, 30, 79,
+ 36, 27, 29, 18, 29, 23, 29, 29, 34, 32, 13, 40, 31, 21, 39,
+ 27, 27, 36, 29, 27, 124, 31, 51, 22, 18, 18, 13, 28, 29, 28, 52,
+ 29, 27, 13, 29, 24, 19, 16, 22, 21, 18, 9, 14, 30, 30, 25, 33, 44,
+ 33, 27, 32, 29, 27, 20, 76, 10, 19, 20, 32, 9, 17, 28, 22, 30, 29, 55,
+ 17, 27, 21, 13, 6, 21, 23, 24, 44, 13, 21, 7, 29, 17, 23, 39, 26, 20, 60,
+ 30, 27, 40, 29, 28, 19, 26, 16, 7, 25, 31, 9, 39, 30, 8, 30, 27, 18, 26, 143,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+ 4, 0, -4, 8, -3, -30, -6, 49, -10, -7, -12, -18, 3, 3, -6, -11, 3, 2, 2, -5, 0, 59,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+};
+int tmpam250mt[]={
+ 2,
+ 0, 0,
+ 2, 0, 6,
+ -1, 0, 0, 12,
+ 0, 0, 6, -1, 13,
+ -3, 0, 1, -3, 3, 5,
+ 1, 0, 3, -3, 1, 7, 6,
+ 6, 0, 2, 6, 0, -1, -3, 11,
+ 12, 0, -3, -2, -2, -5, 5, 0, 2,
+ -1, 0, 0, 3, -1, 0, 3, -3, -3, 12,
+ 0, 0, 11, 6, 3, 5, 3, -3, -2, -2, 3,
+ -1, 0, 2, 2, -1, -2, -1, -3, -4, -1, 9, 3,
+ 7, 0, 1, 0, 12, 6, 0, -1, -1, 6, -3, -2, 11,
+ 6, 0, 8, 11, 6, 0, 2, -4, -1, -2, 3, -4, 3, 11,
+ 11, 0, -2, -1, 0, 1, 2, 1, 8, 7, -3, 13, 6, 0, 11,
+ -1, 0, 11, 2, 1, 6, 6, 12, 0, -1, -1, -3, -1, -2, 0, 7,
+ -2, 0, -3, 0, -3, -3, -3, -4, -2, -3, -3, -4, 2, 2, -2, 6, 3,
+ 6, 0, 7, 2, 5, -2, 11, -2, -5, -1, 6, -6, -3, -1, 3, 3, 1, 3,
+ -3, 0, 0, -3, -1, 0, -3, -2, 1, -3, 4, -6, 2, 1, 3, 11, 0, -3, 2,
+ 3, 0, 13, 1, 3, -2, 0, -4, -4, -2, 1, 1, 11, 2, -1, 2, -1, -4, -1, 12,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+ 2, 0, 1, 8, 3, -4, -3, 2, -4, -3, -1, -4, 0, 0, -1, -3, 1, 5, 0, -3, 0, 10,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, };
+int rtmpam250mt[]={
+ 32,
+ 24, 24,
+ 32, 24, 49,
+ 20, 24, 24, 73,
+ 24, 24, 49, 20, 77,
+ 12, 24, 28, 12, 36, 45,
+ 28, 24, 36, 12, 28, 53, 49,
+ 49, 24, 32, 49, 24, 20, 12, 69,
+ 73, 24, 12, 16, 16, 4, 45, 24, 32,
+ 20, 24, 24, 36, 20, 24, 36, 12, 12, 73,
+ 24, 24, 69, 49, 36, 45, 36, 12, 16, 16, 36,
+ 20, 24, 32, 32, 20, 16, 20, 12, 8, 20, 61, 36,
+ 53, 24, 28, 24, 73, 49, 24, 20, 20, 49, 12, 16, 69,
+ 49, 24, 57, 69, 49, 24, 32, 8, 20, 16, 36, 8, 36, 69,
+ 69, 24, 16, 20, 24, 28, 32, 28, 57, 53, 12, 77, 49, 24, 69,
+ 20, 24, 69, 32, 28, 49, 49, 73, 24, 20, 20, 12, 20, 16, 24, 53,
+ 16, 24, 12, 24, 12, 12, 12, 8, 16, 12, 12, 8, 32, 32, 16, 49, 36,
+ 49, 24, 53, 32, 45, 16, 69, 16, 4, 20, 49, 0, 12, 20, 36, 36, 28, 36,
+ 12, 24, 24, 12, 20, 24, 12, 16, 28, 12, 40, 0, 32, 28, 36, 69, 24, 12, 32,
+ 36, 24, 77, 28, 36, 16, 24, 8, 8, 16, 28, 28, 69, 32, 20, 32, 20, 8, 20, 73,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+ 2, 0, 1, 8, 3, -4, -3, 2, -4, -3, -1, -4, 0, 0, -1, -3, 1, 5, 0, -3, 0, 10,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+};
+int promoter_tf1[]={
+ 0,
+ 0, 0,
+ 0, 0, 0,
+ 0, 0, 0, 6,
+ 0, 0, 0, 1, 6,
+ 0, 0, 0, 2, 2, 5,
+ 0, 0, 0, 0, 0, 0, 0,
+ 0, 0, 0, 2, 2, 1, 0, 6,
+ 0, 0, 0, 1, -1, -1, 0, -2, 6,
+ 0, 0, 0, -3, 1, -2, 0, -3, 1, 6,
+ 0, 0, 0, -2, -1, 2, 0, -2, 2, 1, 6,
+ 0, 0, 0, -1, -1, -1, 0, 1, 2, 2, 2, 6,
+ 0, 0, 0, 1, -1, 0, 0, -2, 1, -2, -1, -1, 6,
+ 0, 0, 0, -2, 2, -1, 0, -2, -2, 1, -2, -1, 2, 6,
+ 0, 0, 0, -3, -3, 2, 0, -4, -3, -3, 0, -3, 1, 0, 4,
+ 0, 0, 0, -2, -1, -2, 0, 2, -2, -2, -2, 2, 1, 2, 0, 6,
+ 0, 0, 0, 1, -1, -1, 0, -1, 1, -2, -2, -1, 1, -2, -3, -1, 6,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+ 0, 0, 0, -2, 2, -2, 0, -2, -1, 2, -2, 0, -2, 1, -3, -1, 1, 0, 6,
+ 0, 0, 0, -2, -2, 2, 0, -2, -2, -2, 2, -1, -1, -2, 0, -1, 2, 0, 1, 5,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+ 0, 0, 0, -1, -2, -2, 0, 2, -2, -2, -2, 2, -2, -3, -4, 1, 2, 0, 2, 1, 0, 6,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+};
+int blosumR[]={
+ 0,
+ 0, 0,
+ 0, 0, 0,
+ 0, 0, 0, 6,
+ 0, 0, 0, 1, 6,
+ 0, 0, 0, 2, 2, 5,
+ 0, 0, 0, 0, 0, 0, 0,
+ 0, 0, 0, 2, 2, 1, 0, 6,
+ 0, 0, 0, 1, -1, -1, 0, -2, 6,
+ 0, 0, 0, -3, 1, -2, 0, -3, 1, 6,
+ 0, 0, 0, -2, -1, 2, 0, -2, 2, 1, 6,
+ 0, 0, 0, -1, -1, -1, 0, 1, 2, 2, 2, 6,
+ 0, 0, 0, 1, -1, 0, 0, -2, 1, -2, -1, -1, 6,
+ 0, 0, 0, -2, 2, -1, 0, -2, -2, 1, -2, -1, 2, 6,
+ 0, 0, 0, -3, -3, 2, 0, -4, -3, -3, 0, -3, 1, 0, 4,
+ 0, 0, 0, -2, -1, -2, 0, 2, -2, -2, -2, 2, 1, 2, 0, 6,
+ 0, 0, 0, 1, -1, -1, 0, -1, 1, -2, -2, -1, 1, -2, -3, -1, 6,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+ 0, 0, 0, -2, 2, -2, 0, -2, -1, 2, -2, 0, -2, 1, -3, -1, 1, 0, 6,
+ 0, 0, 0, -2, -2, 2, 0, -2, -2, -2, 2, -1, -1, -2, 0, -1, 2, 0, 1, 5,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+ 0, 0, 0, -1, -2, -2, 0, 2, -2, -2, -2, 2, -2, -3, -4, 1, 2, 0, 2, 1, 0, 6,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,};
+/******************************COPYRIGHT NOTICE*******************************/
+/*© Centro de Regulacio Genomica */
+/*and */
+/*Cedric Notredame */
+/*Fri Feb 18 08:27:45 CET 2011 - Revision 596. */
+/*All rights reserved.*/
+/*This file is part of T-COFFEE.*/
+/**/
+/* T-COFFEE is free software; you can redistribute it and/or modify*/
+/* it under the terms of the GNU General Public License as published by*/
+/* the Free Software Foundation; either version 2 of the License, or*/
+/* (at your option) any later version.*/
+/**/
+/* T-COFFEE is distributed in the hope that it will be useful,*/
+/* but WITHOUT ANY WARRANTY; without even the implied warranty of*/
+/* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the*/
+/* GNU General Public License for more details.*/
+/**/
+/* You should have received a copy of the GNU General Public License*/
+/* along with Foobar; if not, write to the Free Software*/
+/* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA*/
+/*............................................... |*/
+/* If you need some more information*/
+/* cedric.notredame@europe.com*/
+/*............................................... |*/
+/**/
+/**/
+/* */
+/******************************COPYRIGHT NOTICE*******************************/