Hello all, Coming back after lunch... I managed to load a matrix using this: from Bio import pairwise2 from Bio.SubsMat import MatrixInfo as m #print dir(m) matrix = m.blosum60 pairwise2.align.localdx(seqA, seqB, matrix) Thanks a lot for the help, it was simple after all, just a bit hard to start with..