Example #1
0
#!/usr/bin/env python3
# Usage: ./msa_approx.py score_matrix 5 test.fa
# Requires Alignment.py in the directory

from Alignment import Alignment
from Alignment import GetArguments
import sys

arguments = GetArguments(sys.argv)
test = Alignment(arguments.seqs, arguments.score_matrix, arguments.gapcost)
alignm = test.multiple_align()

for i in range(len(alignm)):
    print(">", arguments.heads[test.seqOrder[i]])
    print(test.num_to_sequence(alignm[i]), "\n")