INPUT:
******
>s1
ACATGCTACACGTATCCGATACCCCGTAACCGATAACGATACACAGACCTCGTACGCTTG
CTACAACGTACTCTATAACCGAGAACGATTGACATGCCTCGTACACATGCTACACGTACT
CCGAT
>s2
ACATGCGACACTACTCCGATACCCCGTAACCGATAACGATACAGAGACCTCGTACGCTTG
CTAATAACCGAGAACGATTGACATTCCTCGTACAGCTACACGTACT
CCGAT
OUTPUT:
********
Scores: match = 1, mismatch = -2, h =-5, g = -2
Sequence 1 = "s1", length = 125 characters
Sequence 2 = "s2", length = 111 characters
s1 1 ACATGCTACACGTATCCGATACCCCGTAACCGATAACGATACACAGACCTCGTACGCTTG 60
|||||| |||| ||||||||||||||||||||||||||||| ||||||||||||||||
s2 1 ACATGCGACACTACTCCGATACCCCGTAACCGATAACGATACAGAGACCTCGTACGCTTG 60
s1 61 CTACAACGTACTCTATAACCGAGAACGATTGACATGCCTCGTACACATGCTACACGTACT 120
||| ||||||||||||||||||||| ||||||||| ||||||||||||
s2 61 CTA-----------ATAACCGAGAACGATTGACATTCCTCGTACA---GCTACACGTACT 106
s1 121 CCGAT 125
|||||
s2 107 CCGAT 111
Report:
Global optimal score = 55
Number of: matches = 105, mismatches = 6, gaps = 14, opening gaps = 2
Identities = 105/125 (84%), Gaps = 14/125 (11%)