Дополнительное задание

В данном задании необходимо было разобраться с альтернативным способом построения дерева по выравниванию. С помощью программы fprotdist из пакета EMBOSS, получающей на вход файл с выровненными последовательностями, был получен файл с расширением .fprotdist, в котором в виде таблицы были записана расстояния между последовательностями белков, которые затем будут использованы другой программой - fneighbor.
Команда:

fprotdist al.fasta al.fprotdist

Программа fneighbor получает на вход файл al.fprotdist и с помощью алгоритма Neighbor-joining строит дерево, которое записывается в файл al.treefile в виде скобочной формулы.
Команда:

fneighbor al.fprotdist

С помощью программы MEGA эта формула была показана графически. На Рис.1 и 2 можно сравнить дерево, построенное алгоритмом Neighbor-joining в программе MEGA и программой fneighbor. Можно заметить, что они являются топологически идентичными, хотя и по-разному укорененными.

Рис. 1. Дерево, построенное алгоритмом
Neighbor-joining в программе MEGA.

Рис. 2. Дерево, построенное алгоритмом
Neighbor-joining в программе fneighbor.

Далее с помощью программы fdrawtree было построено неукорененное дерево из пакета PHYLIP, выдача программы происходит в файл в формате .ps, поэтому для его визуализации необходио было конвертировать полученный файл, например, в pdf с помощью программы ps2pdf.
Команда:

fdrawtree -intreefile al.treefile -plotfile al.ps

Результат работы программы представлен на Рис. 3, можно убедиться, что это дерево так же топологически идентично двум предыдущим.

Рис. 3. Неукорененное дерево.

Работа с программой fdrawgram не принесла положительных результатов, так как выдавала ошибку при открытии начального файла.