Skip to content
Permalink
Branch: master
Find file Copy path
Find file Copy path
Fetching contributors…
Cannot retrieve contributors at this time
21 lines (15 sloc) 397 Bytes
from LocalAlignment import *
from scores import *
from fasta import *
s=readFASTA("rosalind_laff.txt")
def pam250_penalty(x,y):
return -pam250_score(x,y)
def blosum62_penalty(x,y):
return -blosum62_score(x,y)
def agap(j):
return 10+j
X=LocalAlignment(s.values()[0],s.values()[1],gap=agap,match_score=blosum62_penalty)
(R,S)=X.aligned()
print X.score()
print R
print S
You can’t perform that action at this time.