Было - правильно.
Это Вы меня извините, Игорь.
Было неправильно - стало еще неправильнее. =)
У меня ccode /99 67.87;, а надо ccode /99 66.86;. Правильно?
Правда первый вариант был безобиднее.
Правильно будет:
ccode /99 67.86;Вашей ошибки нет. Это мне, надо лучше помнить, сколько у меня маркеров в расчете.
По скрипту: было бы очень замечательно ввести, в теле скрипта, переменную/флаг, при установке которой, tree-tags прописывались бы только на узлах дерева, а на листьях, там где сейчас пишется только кол-во мутаций, все оставалось бы чистым.
Поясню. При перегонке результатов, скриптом, в ньюик, теги на листьях, создают проблемы, которые еще не нашел как обойти. Хотя в будущем, возможно, и напишу решение.
И вопрос: возможно восстановить гаплотип конкретного узла дерева? Пусть и в формате матрицы состояний. Было бы очень хорошо, если можно "вытащить" из дерева предковые гаплотипы, на которые укореняются наши бранчи. И в дальнейшем вести расчет не всего дерева r1a, а только одного кластера, подсунув в качестве предкового гаплотипа, тот который получили ранее.
Спасибо.