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

Для данного задания я выбрал три белка, принадлежащие организму Saccharomyces cerevisiae: Succinate dehydrogenase assembly factor 2, mitochondrial (SDHF2_YEAST); Helicase SWR1 (SWR1_YEAST), tRNA (guanine(9)-N1)-methyltransferase (TRM10_YEAST). Поисковой запрос был следующий:
taxonomy:ascomycota existence:"Evidence at protein level [1]" AND reviewed:yes

Последовательности белков получены с помощью команды (вместо name - название белка):

seqret sw:"NAME_YEAST" name.fasta

Затем была создана локальная база данных, и были выполнены поиски гомологов:

makeblastdb -in X5.fasta -dbtype nucl
tblastn -query name.fasta -db X5.fasta -out name.out

Результаты:

• Succinate dehydrogenase assembly factor 2, mitochondrial - получено 2 скэффолда с Е-value порядка 1e-20.
• Helicase SWR1 - найдено множество скэффолдов с низким E-value (до 1e-68), разбросанных по всему белку, что также может свидетельствовать о гомологии.
• tRNA (guanine(9)-N1)-methyltransferase - 2 скэффолда с E-value порядка 1е-32. Вероятна гомология участков белка.