Entrez Direct, BLAST+, EMBOSS

Поиск гомологов белков в неаннотированном геноме

Для получения последовательностей белков организма относительно близкого к Amoeboaphelidium protococcarum (Saccharomyces cerevisiae) использовалась команда:

seqret sw:*** out.fasta (где вместо звездочек АС белка)

В результате я получил три белка:

Актин

Миозин

Тубулин-В

Сравнение белков с геномом было осуществлено с помощью команды:

tblastn -query NAME.fasta -db X5.fasta > NAME.txt

Результат:

Актин

Миозин

Тубулин-В

Все три белка являются гомологами, о чем свидетельствует e-value 0. И достаточно хорошие проценты идентичности: 89%, 44% и 75% соответственно.