def subseq_la_search(target, data, matrix, gap_cost, min_score, first_only):
# Substitution matrix
sub_matrix = SubMatrix.SubMatrix(matrix)
# The maximum score for given target
max_score = alignment.calculate_max_score(target, sub_matrix)
match_list = list()
for model in data.keys():
for chain in data[model].keys():
sequence = data[model][chain]['sequence']
sw = SmithWaterman.SmithWaterman(target, sequence, gap_cost,
sub_matrix)
# Skip if alignment best score is less than minimum passing score
if float(sw.get_best_score()) / max_score * 100 < min_score:
continue
for i, j in sw.get_coordinates():
aligned_target, aligned_sequence, start_i, start_j =\
sw.get_traceback(i, j)
start_pos = start_j - 1
for _ in range(0, len(aligned_sequence.replace('-', ''))):
resi = data[model][chain]['ids'][start_pos]
match_list.append((model, chain, resi))
start_pos += 1
alignment_string, identities, gaps, mismatches = \
alignment.create_alignment_string(aligned_target, aligned_sequence)
alignment.print_alignment(model, chain, target, sequence,
sub_matrix.get_name(), gap_cost,
sw.get_best_score(), max_score,
identities, mismatches, gaps,
aligned_target, aligned_sequence,
alignment_string, start_i, start_j,
data[model][chain]['ids'])
if first_only:
break
else:
continue
break
return match_list if len(match_list) != 0 else None
def setUp(self):
self.SW = SmithWaterman('', '')
def setUp(self):
self.seq1 = "AGC"
self.seq2 = "ACA"
self.SW = SmithWaterman(self.seq1, self.seq2)
self.function_toBeTested = self.SW._score_diag
def setUp(self):
self.seq1 = 'AG'
self.seq2 = 'AC'
self.SW = SmithWaterman(self.seq1, self.seq2)
self.test_function = self.SW._score_diag
def setUp(self):
self.seq1 = "TTACCGGCCAACTAA"
self.seq2 = "ACCGTGTCACTAC"
self.SW = SmithWaterman(self.seq1, self.seq2)
def setUp(self):
self.seq1 = "AGCACACA"
self.seq2 = "ACACACTA"
self.SW = SmithWaterman(self.seq1, self.seq2)
def setUp(self):
self.seq1 = "AGCA"
self.seq2 = "ACAC"
self.SW = SmithWaterman(self.seq1, self.seq2)
self.function_toBeTested = self.SW._build_scoreMatrix
self.SW._build_scoreMatrix()
def setUp(self):
self.seq1 = "AG"
self.seq2 = "AC"
self.SW = SmithWaterman(self.seq1, self.seq2)
self.function_toBeTested = self.SW._score_horizontal
humanAPOE = "CTACTCAGCCCCAGCGGAGGTGAAGGACGTCCTTCCCCAGGAGCCGGTGAGAAGCGCAGTCGGGGGCACG GGGATGAGCTCAGGGGCCTCTAGAAAGAGCTGGGACCCTGGGAACCCCTGGCCTCCAGGTAGTCTCAGGA GAGCTACTCGGGGTCGGGCTTGGGGAGAGGAGGAGCGGGGGTGAGGCAAGCAGCAGGGGACTGGACCTGG GAAGGGCTGGGCAGCAGAGACGACCCGACCCGCTAGAAGGTGGGGTGGGGAGAGCAGCTGGACTGGGATG TAAGCCATAGCAGGACTCCACGAGTTGTCACTATCATTTATCGAGCACCTACTGGGTGTCCCCAGTGTCC TCAGATCTCCATAACTGGGGAGCCAGGGGCAGCGACACGGTAGCTAGCCGTCGATTGGAGAACTTTAAAA TGAGGACTGAATTAGCTCATAAATGGAACACGGCGCTTAACTGTGAGGTTGGAGCTTAGAATGTGAAGGG AGAATGAGGAATGCGAGACTGGGACTGAGATGGAACCGGCGGTGGGGAGGGGGTGGGGGGATGGAATTTG AACCCCGGGAGAGGAAGATGGAATTTTCTATGGAGGCCGACCTGGGGATGGGGAGATAAGAGAAGACCAG GAGGGAGTTAAATAGGGAATGGGTTGGGGGCGGCTTGGTAAATGTGCTGGGATTAGGCTGTTGCAGATAA TGCAACAAGGCTTGGAAGGCTAACCTGGGGTGAGGCCGGGTTGGGGCCGGGCTGGGGGTGGGAGGAGTCC TCACTGGCGGTTGATTGACAGTTTCTCCTTCCCCAGACTGGCCAATCACAGGCAGGAAGATGAAGGTTCT GTGGGCTGCGTTGCTGGTCACATTCCTGGCAGGTATGGGGGCGGGGCTTGCTCGGTTCCCCCCGCTCCTC CCCCTCTCATCCTCACCTCAACCTCCTGGCCCCATTCAGGCAGACCCTGGGCCCCCTCTTCTGAGGCTTC TGTGCTGCTTCCTGGCTCTGAACAGCGATTTGACGCTCTCTGGGCCTCGGTTTCCCCCATCCTTGAGATA GGAGTTAGAAGTTGTTTTGTTGTTGTTGTTTGTTGTTGTTGTTTTGTTTTTTTGAGATGAAGTCTCGCTC TGTCGCCCAGGCTGGAGTGCAGTGGCGGGATCTCGGCTCACTGCAAGCTCCGCCTCCCAGGTCCACGCCA TTCTCCTGCCTCAGCCTCCCAAGTAGCTGGGACTACAGGCACATGCCACCACACCCGACTAACTTTTTTG TATTTTCAGTAGAGACGGGGTTTCACCATGTTGGCCAGGCTGGTCTGGAACTCCTGACCTCAGGTGATCT GCCCGTTTCGATCTCCCAAAGTGCTGGGATTACAGGCGTGAGCCACCGCACCTGGCTGGGAGTTAGAGGT TTCTAATGCATTGCAGGCAGATAGTGAATACCAGACACGGGGCAGCTGTGATCTTTATTCTCCATCACCC CCACACAGCCCTGCCTGGGGCACACAAGGACACTCAATACATGCTTTTCCGCTGGGCGCGGTGGCTCACC CCTGTAATCCCAGCACTTTGGGAGGCCAAGGTGGGAGGATCACTTGAGCCCAGGAGTTCAACACCAGCCT GGGCAACATAGTGAGACCCTGTCTCTACTAAAAATACAAAAATTAGCCAGGCATGGTGCCACACACCTGT GCTCTCAGCTACTCAGGAGGCTGAGGCAGGAGGATCGCTTGAGCCCAGAAGGTCAAGGTTGCAGTGAACC ATGTTCAGGCCGCTGCACTCCAGCCTGGGTGACAGAGCAAGACCCTGTTTATAAATACATAATGCTTTCC AAGTGATTAAACCGACTCCCCCCTCACCCTGCCCACCATGGCTCCAAAGAAGCATTTGTGGAGCACCTTC TGTGTGCCCCTAGGTACTAGATGCCTGGACGGGGTCAGAAGGACCCTGACCCACCTTGAACTTGTTCCAC ACAGGATGCCAGGCCAAGGTGGAGCAAGCGGTGGAGACAGAGCCGGAGCCCGAGCTGCGCCAGCAGACCG AGTGGCAGAGCGGCCAGCGCTGGGAACTGGCACTGGGTCGCTTTTGGGATTACCTGCGCTGGGTGCAGAC ACTGTCTGAGCAGGTGCAGGAGGAGCTGCTCAGCTCCCAGGTCACCCAGGAACTGAGGTGAGTGTCCCCA TCCTGGCCCTTGACCCTCCTGGTGGGCGGCTATACCTCCCCAGGTCCAGGTTTCATTCTGCCCCTGTCGC TAAGTCTTGGGGGGCCTGGGTCTCTGCTGGTTCTAGCTTCCTCTTCCCATTTCTGACTCCTGGCTTTAGC TCTCTGGAATTCTCTCTCTCAGCTTTGTCTCTCTCTCTTCCCTTCTGACTCAGTCTCTCACACTCGTCCT GGCTCTGTCTCTGTCCTTCCCTAGCTCTTTTATATAGAGACAGAGAGATGGGGTCTCACTGTGTTGCCCA GGCTGGTCTTGAACTTCTGGGCTCAAGCGATCCTCCCGCCTCGGCCTCCCAAAGTGCTGGGATTAGAGGC ATGAGCCACCTTGCCCGGCCTCCTAGCTCCTTCTTCGTCTCTGCCTCTGCCCTCTGCATCTGCTCTCTGC ATCTGTCTCTGTCTCCTTCTCTCGGCCTCTGCCCCGTTCCTTCTCTCCCTCTTGGGTCTCTCTGGCTCAT CCCCATCTCGCCCGCCCCATCCCAGCCCTTCTCCCCGCCTCCCACTGTGCGACACCCTCCCGCCCTCTCG GCCGCAGGGCGCTGATGGACGAGACCATGAAGGAGTTGAAGGCCTACAAATCGGAACTGGAGGAACAACT GACCCCGGTGGCGGAGGAGACGCGGGCACGGCTGTCCAAGGAGCTGCAGGCGGCGCAGGCCCGGCTGGGC GCGGACATGGAGGACGTGTGCGGCCGCCTGGTGCAGTACCGCGGCGAGGTGCAGGCCATGCTCGGCCAGA GCACCGAGGAGCTGCGGGTGCGCCTCGCCTCCCACCTGCGCAAGCTGCGTAAGCGGCTCCTCCGCGATGC CGATGACCTGCAGAAGCGCCTGGCAGTGTACCAGGCCGGGGCCCGCGAGGGCGCCGAGCGCGGCCTCAGC GCCATCCGCGAGCGCCTGGGGCCCCTGGTGGAACAGGGCCGCGTGCGGGCCGCCACTGTGGGCTCCCTGG CCGGCCAGCCGCTACAGGAGCGGGCCCAGGCCTGGGGCGAGCGGCTGCGCGCGCGGATGGAGGAGATGGG CAGCCGGACCCGCGACCGCCTGGACGAGGTGAAGGAGCAGGTGGCGGAGGTGCGCGCCAAGCTGGAGGAG CAGGCCCAGCAGATACGCCTGCAGGCCGAGGCCTTCCAGGCCCGCCTCAAGAGCTGGTTCGAGCCCCTGG TGGAAGACATGCAGCGCCAGTGGGCCGGGCTGGTGGAGAAGGTGCAGGCTGCCGTGGGCACCAGCGCCGC CCCTGTGCCCAGCGACAATCACTGAACGCCGAAGCCTGCAGCCATGCGACCCCACGCCACCCCGTGCCTC CTGCCTCCGCGCAGCCTGCAGCGGGAGACCCTGTCCCCGCCCCAGCCGTCCTCCTGGGGTGGACCCTAGT TTAATAAAGATTCACCAAGTTTCACGCA"
monkeyAPOE = "ATGAGCTCAGGGGCCTCTAGAAAGATGTAGCTGGGACCTCGGGAAGCCCTGGCCTCCAGGTAGTCTCAGG AGAGCTACTCGAGGTCGGGCTTGGGGATAGGAGGAGCGGGGGTGAGGCCAGCAGCAGGGGACTGGACCTG GTAAGGGCTGGGCAGCAGAGACGACCCGACCCGCTAGAAGGTGGGGTGGGCAGAGCATGTGGACTAGGAG CTAAGCCACAGCAGGACCCCCACGAGTTGTCACTGTCATTTATCGAGCACCTACTGGGTGTCCCCAGTGT CCTCAGATCTCCATAACAGGGAAGCCAGGGGCAGCGACACGGTAGCTAGCCGTCGATTGGAGAACTTTAA AATGAGGACTGAATTAGCTCATAAATGGAAAACGGCGCTTAAATGTGAGGTTAGAGCTTAGAATGTGAAG GGAGAATGAGGAATGCGAGACTGGGACTGAGATGGAACCGGCGGTGGGGAGGGGGAGGGGGTGTGGAATT TGAACCCCGGGAGAGAAAGATGGAATTTTGACTATGGAGGCCGACCTGGGGATGGGGAAATAAGAGAAGA CCAGGAGGGAGTTAAATAGGGAATGGGTTGGGGGCGGCTTGGTAACTGTTTGTGCTAGGATTAGGCTGTT GCAGATAATGGAACTAGGCTTGGAAGGCTAACCTGGGGTGGGGCCGGGTTGGGGTCGGGCTGGGAGTGGG AGGAGTCCTTACTGGCGGTTGATTGACTGTTTCTCCCTCCCCAGACTGGCCAATCACAGGCAGGAAGATG AAGGTTCTGTGGGCTGCGTTGCTGGTCACATTCCTGGCAGGTATGGGGGCGGGGCTTGCTCGGTTTCCCT GCTCCTCCCCCTCTCATCCTCACCTCAACCTCCTGGCCCCATTCAGGCAGACTTCGGGCCCCCTTTTCTT CTGCTGGTCTGTCTTCCCCTTGAGGGGAAAGCCCAGGTCTGAGGCTTCTATGCTGCTTTCTGGCTCAGAA CAGCGATTTGACGCTCTGTGAGCCTCGTTTTCCTGCCCCCGCTTTTTTTTTTTTTTTTTTTTTTTGAGCC AGAGTCTCACTCTGTCGCCCGGCTGGAGTGCAGTGGCACAATCTCAGCTCACTGCAAGCTCCGCCTCCCG GGTTCACGCTATTCTCCCGCCTCAGCCTCCCGAGTAGCTGGGACTACAGGCGCCCGCCACCATGCCCGGC TAATTTTTTGTACTTTGAGTAGAGAAGGGGTTTCACTGTATTATCCAGGATGGTCTCTATCTCCTGACCT CGTGATCTGCCCGCCTTGGCCTCCCAAAGTGCTGGAATTACAGGCGTGAGCCACCGCGCCCGGCCTCCCC ATCCTTAATATAGGAGTTAGAAGTTTTTGTTTGTTTTTGTTTTGTTTTGAGATGAAGTCCCTCTGTCGCC CAGGCTGGAGTGCGGTGGCTCCCAGGCTGGAGTTCAGTGGCAGGATCTCAGCTCACTGCAAGCTCCCCCT CCCAGGTTCATGCCATTCTCCTGCCTCAGCCTCCGGAGTAGCTGGGACTACAGGAACATGCCACCACACC TGACTAACTTTTTTTGTATTTTTAGTAGAGACGGGGTTTCACCATGTTGGCCAGGCTGGTCTGGAACTCC TGACCTCAGGTGATCTGCCTGCTTCAACCTCCCAAAGTGCTGGGATTACAGGCGTGGGCCACCGCGCCCA GCTGGGAGTTAAGGGGCTTCTAATGCATTGCATTAGAATACCAGACACGGGACAGCTGTGATCTTTATTC TCTATCACCCCACACAGCCCTGCCTGGGACACACAAGGACACTCAATACATGCTTTTCCGCTAGGCACGG TGGCTCACCCCTGTAATCCCAGCATTTTGGGAGGCCAAGGTGGGAGGATCACTTGAGCCCAGGAGTTCAA CACCAGACTCGGCAACATAGTGAGACTCTGTCTCTACTAAAAATACAAAAATTAGCCAGGCCTGGTGCCA CACACCTGTGGTCCCAGCTACTCAGGAGGCTGAGGCAGGAGGATTGCTTGAGCCCAGAAGGTCAAGGTTG CAGTGAACCATGTTCAGGCCACTGCAATCCAGCCTGGGTGACAGAGCAAGACCCTGTCTGTAAATAAATA ACGCTTTTCAAGTGATTAAACAGACTCCCCCCTCACCCTGCCCACCATGGCTCCAAAGCAGCATTTGTGG AGCACCTTCTGTGTGCCCCTAGGTACTAGGTGCCTGGACGGGGTCAGAAGGAACCTGAACCACCTTCAAC TTGTTCCACACAGGATGCCAGGCCAAGGTGGAGCAACCGGTGGAGCCAGAGACGGAACCCGAGCTTCGCC AGCAGGCTGAGGGGCAGAGCGGCCAGCCCTGGGAGCTGGCACTGGGTCGCTTTTGGGATTACCTGCGCTG GGTGCAGACACTGTCTGAGCAGGTGCAGGAGGAGCTGCTCAGCCCCCAGGTCACCCAGGAACTGACGTGA GTGTCCCCATCCCGGCCCTTGACCCTTCTGGTGGGCGGCTATACCTCCCCAGGTCCAGGTTTCATTCTGC CCCTGCCACTAAGTCTTGGGAGTCCTGGGTCTCTGCTGGTTCTAGCTTCCTCTTCCCATTTCTGACTCCT GGCTTTAGCTCTCTGGAATTCTCTCTCTCAGTTCTGTTTCTCCCTCTTCCCTTCTGACTCAGCCTCTCAC ACTCGTCCTGGCTCTGTCTCTGTCCTTCACTAGCTCTTTTATATAGAGACAGAGAGATGGGGTCTCACTG TGTTGCCCAGGCTGGTCTTGAACTTCTGGGCTCAAGCGATCCTCCCACCTCGCCTCCCAAAGTGCTGGGA ATAAAGACATGAGCCACCTTGCCCGGCCTCCTAGCTCTTTCTTCGTCTCTGCCTCTGCTCTCTGCGTCTG TCTTTGTCTCCTCTCTGCCTCTGTCCCGTTCCTTCTCTCTTGGTTCACTGCCCTTCTGTCTCTCCCTGTT CTCCTTAGGAGACTCTCCTCTCTTCCTTCTCGGGTCTCTCTGGCTGATCCCCATCTCACCCACACCTATC CCAGCCCTTCTCGCCTCCCCCTGTGCGCACACCCTCCTGCTCTTTCGGCTGCAGGACGCTGATGGATGAG ACCATGAAGGAGTTGAAGGCCTACAAATCGGAACTGGAGGAACAGCTGAGCCCGGTGGCGGAGGAGACGC GGGCACGGCTGTCCAAGGAGCTACAGGCGGCGCAGGCCCGGCTGGGTGCCGACATGGAGGACGTGCGCAG CCGCCTGGTGCAGTACCGCAGCGAGGTGCAGGCCATGCTGGGCCAGAGTACCGAGGAGCTGCGGGCGCGC CTCGCCTCCCACCTGCGCAAGCTGCGCAAGCGGCTCCTCCGCGATGCTGATGACCTGCAGAAGCGCCTGG CAGTGTATCAGGCCGGGGCCCGCGAGGGCGCCGAGCGCGGGGTCAGCGCCATCCGCGAGCGCCTGGGACC CCTGGTGGAGCAGGGCCGCGTGCGGGCCGCCACTGTGGGCTCCCTGGCCAGCCAGCCGCTTCAGGAGCGG GCCCAGGCCTTGGGTGAGCGGCTTCGCGCACGGATGGAGGAGATGGGCAGCCGGACCCGCGACCGCCTGG ACGAGGTGAAGGAGCAGGTGGCGGAGGTGCGCGCCAAGCTGGAGGAACAGGCCCAGCAGATAAGCCTGCA GGCCGAGGCCTTCCAGGCCCGCCTCAAGAGCTGGTTCGAGCCCCTGGTGGAAGATATGCAGCGCCAGTGG GCTGGGCTGGTGGAGAAGGTGCAGGCTGCCGTGGGCGCCAGCACCGCCCCTGTGCCCAGCGACAATCACT GAACGCCCAGGCCTACAGCCATGCGACCCGACTCCACCCCATGCCTCCTCTCTCCGCTCAGCCTGCAGCG GGAGACCCTGTCCCCACCCCAGCCGTCCTCCAGGGGTGGGCCCTAGTTTAATAAAGATTCGCCAAGTTTC ACCGCA"
humanInsulin = "AGCCCTCCAGGACAGGCTGCATCAGAAGAGGCCATCAAGCAGGTCTGTTCCAAGGGCCTTTGCGTCAGGT GGGCTCAGGATTCCAGGGTGGCTGGACCCCAGGCCCCAGCTCTGCAGCAGGGAGGACGTGGCTGGGCTCG TGAAGCATGTGGGGGTGAGCCCAGGGGCCCCAAGGCAGGGCACCTGGCCTTCAGCCTGCCTCAGCCCTGC CTGTCTCCCAGATCACTGTCCTTCTGCCATGGCCCTGTGGATGCGCCTCCTGCCCCTGCTGGCGCTGCTG GCCCTCTGGGGACCTGACCCAGCCGCAGCCTTTGTGAACCAACACCTGTGCGGCTCACACCTGGTGGAAG CTCTCTACCTAGTGTGCGGGGAACGAGGCTTCTTCTACACACCCAAGACCCGCCGGGAGGCAGAGGACCT GCAGGGTGAGCCAACTGCCCATTGCTGCCCCTGGCCGCCCCCAGCCACCCCCTGCTCCTGGCGCTCCCAC CCAGCATGGGCAGAAGGGGGCAGGAGGCTGCCACCCAGCAGGGGGTCAGGTGCACTTTTTTAAAAAGAAG TTCTCTTGGTCACGTCCTAAAAGTGACCAGCTCCCTGTGGCCCAGTCAGAATCTCAGCCTGAGGACGGTG TTGGCTTCGGCAGCCCCGAGATACATCAGAGGGTGGGCACGCTCCTCCCTCCACTCGCCCCTCAAACAAA TGCCCCGCAGCCCATTTCTCCACCCTCATTTGATGACCGCAGATTCAAGTGTTTTGTTAAGTAAAGTCCT GGGTGACCTGGGGTCACAGGGTGCCCCACGCTGCCTGCCTCTGGGCGAACACCCCATCACGCCCGGAGGA GGGCGTGGCTGCCTGCCTGAGTGGGCCAGACCCCTGTCGCCAGGCCTCACGGCAGCTCCATAGTCAGGAG ATGGGGAAGATGCTGGGGACAGGCCCTGGGGAGAAGTACTGGGATCACCTGTTCAGGCTCCCACTGTGAC GCTGCCCCGGGGCGGGGGAAGGAGGTGGGACATGTGGGCGTTGGGGCCTGTAGGTCCACACCCAGTGTGG GTGACCCTCCCTCTAACCTGGGTCCAGCCCGGCTGGAGATGGGTGGGAGTGCGACCTAGGGCTGGCGGGC AGGCGGGCACTGTGTCTCCCTGACTGTGTCCTCCTGTGTCCCTCTGCCTCGCCGCTGTTCCGGAACCTGC TCTGCGCGGCACGTCCTGGCAGTGGGGCAGGTGGAGCTGGGCGGGGGCCCTGGTGCAGGCAGCCTGCAGC CCTTGGCCCTGGAGGGGTCCCTGCAGAAGCGTGGCATTGTGGAACAATGCTGTACCAGCATCTGCTCCCT CTACCAGCTGGAGAACTACTGCAACTAGACGCAGCCCGCAGGCAGCCCCACACCCGCCGCCTCCTGCACC GAGAGAGATGGAATAAAGCCCTTGAACCAGC"
hamsterInsulin = "ATGGCCCTGTGGATGCGCCTTCTGCCCCTGTTGGCCCTGCTGGCCCTCTGGGAGCCGAACCCTGCCCAGG CTTTTGTCAACCAGCACCTTTGTGGCTCCCACCTTGTGGAGGCACTCTACCTGGTGTGTGGGGAGCGTGG CTTCTTCTACACACCCAAATCCCGTCGTGGAGTGGAGGACCCACAAGGTGAGTTCTGCCCCTGAATTCTG TCCCCAGTGCTTGCCACCCTGGTTTTCCTTGCCCACAGGACCTCACAGATTATGTCCTGGGTGTGGAGGG TCTCAGAGGAACTGGGCAGGGGCACATTTCCGTGGGAAGCTAGACATAGCTAAACACAGCGGCTGCTGGG AATGAATGAGAATCCTGCCTTGAGGCTCCTAGGTGCAGACATGTGGGCAGGCCCCAGGATAGGCACCTAT TTGGGGCCGCCATAAAACACTAGGGGTTGGTGGCAGGATGCGTAGGCTTTAGAGCTCTTTGTGTCCATGC CCGGTGACTTGTCCCACATACTGACTTAGCAGGAGAGACAAGGTGAGAGGAAGCCTGGGGTAGGCAGGAG GCTGCTCAGCTGCCCCTGACTGGATTGTCCCATGTGTCTTTGCTTCTATGTTGCTGACACTCTGGCTTGC TCTGACACTGCCTCCCTGGCAGTGACACAGCTGGAGCTGGGTGGCGGCCCTGGAGCAGGTGACCTTCAGA CCTTGGCACTGGAGGTGGCCCAGCAGAAGCGCGGCATTGTGGATCAGTGCTGCACCAGCATCTGCTCGCT CTACCAGCTAGAGAACTACTGCAACTAG"
scores = [
[10, -1, -1],
[2, -6, -2],
[2, -2, -5]
]
for score in scores:
wunsch = nw.NeedlemanWunsch(score[0], score[1], score[2])
smith = sw.SmithWaterman(score[0], score[1], score[2])
print(score)
print("Needleman-Wunsch")
print("\n APOE \n")
wunsch.solve(humanAPOE, monkeyAPOE)
print("\n INS \n")
wunsch.solve(humanInsulin, hamsterInsulin)
print("\nSmith-Waterman")
print("\n APOE \n")
smith.solve(humanAPOE, monkeyAPOE)
print("\n INS \n")
smith.solve(humanInsulin, hamsterInsulin)
