A B C D E F . . * * * * . . . * * * . . . . * *
msbar my.fasta ABCDE.fasta -point 4 -count 293 -auto msbar my.fasta F.fasta -point 4 -count 1171 -auto msbar ABCDE.fasta ABCD.fasta -point 4 -count 68 -auto msbar ABCDE.fasta E.fasta -point 4 -count 878 -auto msbar ABCD.fasta ABC.fasta -point 4 -count 195 -auto msbar ABCD.fasta D.fasta -point 4 -count 810 -auto msbar ABC.fasta AB.fasta -point 4 -count 195 -auto msbar ABC.fasta C.fasta -point 4 -count 615 -auto msbar AB.fasta B.fasta -point 4 -count 420 -auto msbar AB.fasta A.fasta -point 4 -count 420 -auto cat A.fasta >> mut1.fasta cat B.fasta >> mut1.fasta cat C.fasta >> mut1.fasta cat D.fasta >> mut1.fasta cat E.fasta >> mut1.fasta cat F.fasta >> mut1.fasta
fdnaml mut1.fasta -ttratio 0.5 -auto Параметр -ttratio равен 0.5, поскольку программа msbar делает все замены с одинаковой вероятностью и отношение числа транзиций к числу трансверсий равно 1/2. +----------B | | +---------------------------------------------F | +----4 | | | +---------------------------------E 1-------2 +--3 | | +--------------------D | | | +------------------C | +-------------A
Дерево по алгоритму UPGMA (построено по матрице эволюционных расстояний) fneighbor mut1.fdnadist -treetype u -auto -outfile UPGMA.fneighbor +-----------A +------1 +---2 +-----------B ! ! +------3 +------------------C ! ! +------4 +---------------------D ! ! --5 +----------------------------E ! +-----------------------------------F
Дерево по алгоритму ближайших соседей (построено по матрице эволюционных расстояний). fneighbor mut1.fdnadist -auto -outfile NJ.fneighbor +--------B ! ! +----------------C 1--------2 ! ! +----------------------------------E ! +-----3 ! ! +-------------------D ! +-4 ! +-------------------------------------------F ! +--------------A
Для сравнение предложено было сделать таблицу, в левой части которой приведены (в виде точек и звёздочек) все ветви,
встреченные во всех деревьях (исходном и трёх реконструкциях), а в правой добавлены четыре столбца, соостветствующие
четырём деревьям. Знаком + отмечено, в каких деревьях встретилась каждая из ветвей.
A | B | C | D | E | F | правильное дерево | 1ое | 2ое | 3е |
. | . | * | * | * | * | + | + | + | + |
. | . | . | * | * | * | + | + | + | + |
. | . | . | . | * | * | + | - | + | - |
. | . | . | * | . | * | - | - | - | + |
. | . | . | * | * | . | - | + | - | - |
Species in order: 1. C 2. F 3. D 4. E 5. B 6. A Sets included in the consensus tree Set (species in order) How many times out of 100.00 ....** 83.00 .***.. 51.00 .**... 45.00 Sets NOT included in consensus tree: Set (species in order) How many times out of 100.00 .*.*.. 22.00 И полученное при этом дерево: +-------------E +-51.0-| | | +------F | +-45.0-| +------| +------D | | | | +------A | +--------83.0-| | +------B | +---------------------------CПолучили довольно интересный результат - ветка E все так же "не спокойна" и стремится к A,B,C. Это вызвано, по-видимому, тем что у нас есть короткие (а потому несовсем достоверные) ветки, мы видим после надписи "Sets NOT included in consensus tree" возможные альтернативные ветви, которые могли возникнуть с несколько меньшими вероятностями, обратим внимание на то, что разделение на AB и CDEF происходит в большом количестве случаев (83 из 100), разделение на ABC и DEF, происходит уже с меньшей вероятностью (51 процент), ветвь FD, которая не встречалась в реальном дереве существует в 45 из 100 случаев, смотрим далее - ниже приведены альтернативные варианты размещения ветвей, и наш вариант расположения, что существует ветвь FE, такая ветвь существует в 22 случаях (что конечно тоже не мало, но можно было ожидать и лучшего).
©Метелев Михаил