Эволюционная модель

Была дана скобочная структура дерева, которая описывала ход эволюции гена LacY_Ecoli, где эволюционные расстояния представляют собой число мутаций на 100 пар нуклеотидов. Их надо было перевести в число мутаций на полную последовательность гена ( длину последовательности надо умножить на соответстующее число мутаций на 100 н.п. и разделить на 100).

Внимание! В дереве ошибка и все результаты сделаны на основе ошибочного дерева!(в последующем ошибки будут исправлены)
Скобочная структура дерева
(число мутаций на 100 п.н.)
Скобочная структура дерева
(число мутаций
на всю последовательность)
(
A:90,
(
((B:25,C:25):10,
D:35):5,
(E:30, F:30):10):50
);
(
A:1128.6,
(
((B:313.5,C:313.5):125.4,
D:438.9):62.7,
(E:376.2, F:376.2):125.4):627
);
Соответствующее дерево Соответствующее дерево

Далее с помощь программы msbar пакета EMBOSS были получены мутантные последовательности, соответствующие всем узлам дерева.

1 - файл исходной последовательности LacY_nucl.txt
2 - файл pointA.txt
3 - файл point50.txt
4 - файл point10_n.txt
5 - файл point10_v.txt
6 - файл point5.txt
7 - файл pointF.txt
8 - файл pointE.txt
9 - файл pointD.txt
10 - файл pointC.txt
11 - файл pointB.txt

При этом был составлен следующий скрипт:
msbar lacY_nucl.txt pointA.txt -point 4 -count 1129 -auto
msbar lacY_nucl.txt point50.txt -point 4 -count 627 -auto
msbar point50.txt point10_v.txt -point 4 -count 125 -auto
msbar point50.txt point10_n.txt -point 4 -count 125 -auto
msbar point10_v.txt point5.txt -point 4 -count 63 -auto
msbar point10_v.txt pointD.txt -point 4 -count 439 -auto
msbar point10_n.txt pointE.txt -point 4 -count 376 -auto
msbar point10_n.txt pointF.txt -point 4 -count 376 -auto
msbar point5.txt pointB.txt -point 4 -count 314 -auto
msbar point5.txt pointC.txt -point 4 -count 314 -auto

Анализ модели