#!/usr/bin/perl use strict; use warnings; my $S=" MSATPLTQEQKKA"; my $T=" MNATLTQQTKA"; my $n=length($S); my $m=length($T); my @seq1=split("",$S); my @seq2=split("",$T); print"\n\t "; print"@seq1\n"; print"\t "; for (my$j=0;$j<$n;$j=$j+1){ print" $j"; } print"\n\t 0"; print"\n\t"; $initialvalue=0; if($seq1[$i]==$seq2[$j]){ $match=0 if (($seq1[$i-1],$seq2[$j])>($seq1[$i],$seq2[$j-1])||($seq1[$i-1],$seq2[$j-1]) $add=1 print"$initialvalue+$match+$add\n" } else($seq1[$i]!=$seq2[$j]{ $initialvalue= $initialvalue+1 if (($seq1[$i-1],$seq2[$j])>($seq1[$i],$seq2[$j-1])||($seq1[$i-1],$seq2[$j-1]) $add=1 print"$match+$add+$initialvalue\n" } print"\n"; for (my$k=0;$k<$m;$k=$k+1){ print"$seq2[$k] $k\n"; }