Входной текстовый файл в формате fasta:
#!/usr/bin/python
from Bio import AlignIO
seq_file = open('/path/to/fa_alignment_PF00205.txt')
alignment = AlignIO.read(seq_file, "fasta")
Ошибка:
ValueError: Sequences must all be the same length
Входные последовательности не должны быть одинаковой длины, так как на ClustalOmega вы можете выравнивать последовательности разной длины.
Это также не работает... получает ту же ошибку:
alignment = AlignIO.parse(seq_file,"fasta")
for record in alignment:
print(record.id)
Кто-нибудь, кто знаком с BioPython, знает, как обойти это, чтобы выровнять последовательности из файлов fasta?