Эволюционная модель.
Моделирование эволюциигена trpC из Escherichia coli по заданной скобочной структуре дерева.
Пограмы с которыми я работал, выполня это задание:distmat, msbar, emma, ednadist, ednaml
Сначала была получена последовательность гена trpC из документа EMBL V00366.
Далее взята правильная скобочная структура дерева:
((A:50,B:50):50, ((C:50,D:50):5,(E:50,F:50):5):45)
Расчет числа мутаций на полную последовательность:
Длина последовательности trpC – 1359 нуклеотидных пар.
В скобочной структуре даны расстояния в виде числа мутаций на 100 нуклеотидных пар.
Стало быть, чтобы вычислить, сколко мутаций необходимо сделатьь в полной последовательности, необходимо число мутаций на 100 нуклеотидов умножить на 1359\100. Скобочная структура принимает вид:
((A:680,B:680):680, ((C:680,D:680):68,(E:680,F:680):68):612)
Графическая структура дерева:

Скрипт UNIX для получения мутантных последовательностей при помощи программы msbar:
msbar trpC.gene trpc2.gene -point 4 -count 680 -auto
msbar trpC.gene trpc3.gene -point 4 -count 612 -auto
msbar trpc2.gene trpca.gene -point 4 -count 680 -auto
msbar trpc2.gene trpcb.gene -point 4 -count 680 -auto
msbar trpc3.gene trpc4.gene -point 4 -count 68 -auto
msbar trpc3.gene trpc5.gene -point 4 -count 68 -auto
msbar trpc4.gene trpcc.gene -point 4 -count 680 -auto
msbar trpc4.gene trpcd.gene -point 4 -count 680 -auto
msbar trpc5.gene trpce.gene -point 4 -count 680 -auto
msbar trpc5.gene trpcf.gene -point 4 -count 680 -auto
Сравнение методов реконструкции деревьев.
Определение эволюционного расстояния между нуклеотидными последовательностями.