def eq2(j):
    mval = 0
    for i in range(len(mat)): #előző oszlop maximuma
        if mat[i][j-1]>mval:
            mval=mat[i][j-1]
    return max(mat[0][j-1],mval-T)

def eq3(i,j,d,s):
    return max(mat[i-1][j-1]+int(blosum[s][d]),mat[i-1][j]+int(gapval),mat[i][j-1]+int(gapval),mat[0][j-1])

def fillmat(_seq,_dom):
    for j in range(len(mat[0])):
        for i in range(len(mat)):
            if i == 0:
                mat[i][j] = eq2(j)
            else:
                mat[i][j] = eq3(i,j,alph.index(_dom[i-1]),alph.index(_seq[j-1]))