A Python implementation of the global sequence-alignment algorithm that optimises DNA / protein alignments with dynamic programming. It supports custom scoring matrices (e.g. BLOSUM 62) or simple match/mismatch schemes, calculates detailed statistics, and exports publication-ready alignments.
O(mn) time/space Needleman-Wunsch matrix with traceback.
Imports any text-matrix (BLOSUM, PAM, DNA identity) or simple match/mismatch values.
Alignment blocks wrapped at 80 chars, plus score, matches, mismatches, gaps & similarity %.
GATTACA-CTG ||| | | | GA-T-CAACTG
Score: 23 Matches: 7 Mismatches: 2 Gaps: 2 Similarity: 77.8%