Я пытаюсь оценить уже выровненные последовательности. Скажем,
seq1 = 'PAVKDLGAEG-ASDKGT--SHVVY----------TI-QLASTFE'
seq2 = 'PAVEDLGATG-ANDKGT--LYNIYARNTEGHPRSTV-QLGSTFE'
с заданными параметрами
substitution matrix : blosum62
gap open penalty : -5
gap extension penalty : -1
Я просмотрел кулинарную книгу биопайтона, но все, что я могу получить, это матрица подстановки blogsum62, но я чувствую, что в ней должен быть кто-то уже реализовал такую библиотеку.
Так что кто-нибудь может предложить любую библиотеки или кратчайший код, который может решить мою проблему?
Спасибо заранее