function [matrix, matrixInfo] = pam250cdi
%PAM250CDI substitution matrix in 1/10 Bit Units, 
%   Expected score = -2.81, Entropy = 0.354 bits
%   Lowest score = -26, Highest score = 57
%
%   Order:
%   A  R  N  D  C  Q  E  G  H  I  L  K  M  F  P  S  T  W  Y  V  B  Z  X  *
%
%   [MATRIX,MATRIXINFO] = PAM250CDI returns a structure of information about
%   the matrix with fields Name, Scale, Entropy, ExpectedScore, LowestScore,
%   HighestScore and Order.

% Source:  ftp://ftp.ncbi.nih.gov/blast/matrices/
%   $Revision: 1.5 $ 

matrix = [...
    6  -5   1   1  -7  -1   1   4  -4  -2  -6  -4  -4 -12   4   4   4 -19 -12   1   1   0  -1 -26;...
   -5  20   0  -4 -12   4  -4  -8   5  -7 -10  11  -1 -15  -1  -1  -3   7 -14  -8  -2   0  -3 -26;...
    1   0   7   7 -12   3   5   1   5  -6 -10   3  -6 -12  -2   2   1 -13  -7  -6   7   4  -1 -26;...
    1  -4   7  13 -17   5  11   2   2  -8 -13   0  -9 -19  -3   1   0 -22 -14  -7  10   9  -2 -26;...
   -7 -12 -12 -17  40 -18 -18 -11 -11  -8 -20 -18 -17 -14  -9   0  -7 -26   1  -6 -15 -18 -10 -26;...
   -1   4   3   5 -18  13   8  -4  10  -7  -6   2  -3 -15   1  -2  -3 -16 -13  -6   4  11  -2 -26;...
    1  -4   5  11 -18   8  13   1   2  -7 -11   0  -7 -18  -2   0  -1 -23 -14  -6   9  11  -2 -26;...
    4  -8   1   2 -11  -4   1  16  -7  -8 -13  -6  -9 -16  -2   4   0 -23 -17  -4   2  -1  -3 -26;...
   -4   5   5   2 -11  10   2  -7  22  -8  -7   0  -7  -6  -1  -3  -4  -9   0  -7   4   6  -2 -26;...
   -2  -7  -6  -8  -8  -7  -7  -8  -8  15   8  -6   7   3  -7  -5   0 -17  -3  12  -7  -7  -2 -26;...
   -6 -10 -10 -13 -20  -6 -11 -13  -7   8  20  -9  12   6  -8  -9  -6  -6  -3   6 -11  -9  -4 -26;...
   -4  11   3   0 -18   2   0  -6   0  -6  -9  15   1 -17  -4  -1   0 -12 -15  -8   2   1  -3 -26;...
   -4  -1  -6  -9 -17  -3  -7  -9  -7   7  12   1  21   1  -7  -5  -2 -14  -8   6  -7  -5  -2 -26;...
  -12 -15 -12 -19 -14 -15 -18 -16  -6   3   6 -17   1  30 -15 -11 -10   1  23  -4 -15 -17  -8 -26;...
    4  -1  -2  -3  -9   1  -2  -2  -1  -7  -8  -4  -7 -15  19   3   1 -19 -16  -4  -2  -1  -2 -26;...
    4  -1   2   1   0  -2   0   4  -3  -5  -9  -1  -5 -11   3   5   4  -8  -9  -3   2  -1  -1 -26;...
    4  -3   1   0  -7  -3  -1   0  -4   0  -6   0  -2 -10   1   4   9 -17  -9   1   0  -2  -1 -26;...
  -19   7 -13 -22 -26 -16 -23 -23  -9 -17  -6 -12 -14   1 -19  -8 -17  57   0 -21 -18 -19 -13 -26;...
  -12 -14  -7 -14   1 -13 -14 -17   0  -3  -3 -15  -8  23 -16  -9  -9   0  34  -8 -10 -14  -8 -26;...
    1  -8  -6  -7  -6  -6  -6  -4  -7  12   6  -8   6  -4  -4  -3   1 -21  -8  14  -6  -6  -2 -26;...
    1  -2   7  10 -15   4   9   2   4  -7 -11   2  -7 -15  -2   2   0 -18 -10  -6   9   7  -2 -26;...
    0   0   4   9 -18  11  11  -1   6  -7  -9   1  -5 -17  -1  -1  -2 -19 -14  -6   7  11  -2 -26;...
   -1  -3  -1  -2 -10  -2  -2  -3  -2  -2  -4  -3  -2  -8  -2  -1  -1 -13  -8  -2  -2  -2  -3 -26;...
  -26 -26 -26 -26 -26 -26 -26 -26 -26 -26 -26 -26 -26 -26 -26 -26 -26 -26 -26 -26 -26 -26 -26   1;...
    ];

if nargout >1
    matrixInfo.Name = 'PAM250CDI';
    matrixInfo.Scale = 0.1;
    matrixInfo.Entropy = 0.354 ;
    matrixInfo.ExpectedScore = -2.81;
    matrixInfo.LowestScore = -26;
    matrixInfo.HighestScore = 57;
    matrixInfo.Order = 'ARNDCQEGHILKMFPSTWYVBZX*';
end

