1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
|
# hmmscan :: search sequence(s) against a profile database
# HMMER 3.0 (March 2010); http://hmmer.org/
# Copyright (C) 2010 Howard Hughes Medical Institute.
# Freely distributed under the GNU General Public License (GPLv3).
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# query sequence file: s04.fasta
# target HMM database: /home/bow/db/hmmer/Pfam-A.hmm
# output directed to file: hmmer_cases/text_hmmscan_s04.out
# per-seq hits tabular output: hmmer_cases/tab_hmmscan_s04.out
# per-dom hits tabular output: hmmer_cases/domtab_hmmscan_s04.out
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: gi|125490392|ref|NP_038661.2| [L=352]
Description: POU domain, class 5, transcription factor 1 isoform 1 [Mus musculus]
Scores for complete sequence (score includes all domains):
--- full sequence --- --- best 1 domain --- -#dom-
E-value score bias E-value score bias exp N Model Description
------- ------ ----- ------- ------ ----- ---- -- -------- -----------
7e-37 124.8 0.5 1.4e-36 123.9 0.3 1.5 1 Pou Pou domain - N-terminal to homeobox domain
2.1e-18 65.5 1.1 4.1e-18 64.6 0.7 1.5 1 Homeobox Homeobox domain
------ inclusion threshold ------
0.012 15.6 0.0 0.16 12.0 0.0 2.2 2 HTH_31 Helix-turn-helix domain
0.039 13.5 0.0 0.095 12.3 0.0 1.6 1 Homeobox_KN Homeobox KN domain
0.14 10.5 0.1 0.26 9.6 0.1 1.4 1 DUF521 Protein of unknown function (DUF521)
Domain annotation for each model (and alignments):
>> Pou Pou domain - N-terminal to homeobox domain
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 123.9 0.3 5e-40 1.4e-36 3 75 .] 133 205 .. 131 205 .. 0.97
Alignments for each domain:
== domain 1 score: 123.9 bits; conditional E-value: 5e-40
Pou 3 eldleeleefakefkqrrikLgltqadvgsalgalyGkefsqttIcrFEalqLslknmckLkpllekWLeeae 75
++ ++ele+fak +kq+ri+Lg+tqadvg +lg+l+Gk+fsqttIcrFEalqLslknmckL+pllekW+eea+
gi|125490392|ref|NP_038661.2| 133 KALQKELEQFAKLLKQKRITLGYTQADVGLTLGVLFGKVFSQTTICRFEALQLSLKNMCKLRPLLEKWVEEAD 205
67899******************************************************************96 PP
>> Homeobox Homeobox domain
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 64.6 0.7 1.5e-21 4.1e-18 1 57 [] 224 280 .. 224 280 .. 0.98
Alignments for each domain:
== domain 1 score: 64.6 bits; conditional E-value: 1.5e-21
SS--SS--HHHHHHHHHHCCTSSS--HHHHHHHHHH----HHHHHHHHHHHHHHHHH CS
Homeobox 1 rrkRttftkeqleeLeelFeknrypsaeereeLAkklgLterqVkvWFqNrRakekk 57
+rkRt++++ Le +F k+++ps ++++++A++lgL++++V+vWF+NrR+k k+
gi|125490392|ref|NP_038661.2| 224 KRKRTSIENRVRWSLETMFLKCPKPSLQQITHIANQLGLEKDVVRVWFCNRRQKGKR 280
79****************************************************997 PP
>> HTH_31 Helix-turn-helix domain
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ? 12.0 0.0 5.7e-05 0.16 1 35 [. 141 181 .. 141 184 .. 0.96
2 ? 0.8 0.0 0.19 5.2e+02 39 62 .. 245 268 .. 243 270 .. 0.86
Alignments for each domain:
== domain 1 score: 12.0 bits; conditional E-value: 5.7e-05
HTH_31 1 aLGarLralReraGLtqeevAerlg......vSastlsrlE 35
+++ +L++ R + G tq++v+ lg +S++t++r E
gi|125490392|ref|NP_038661.2| 141 QFAKLLKQKRITLGYTQADVGLTLGvlfgkvFSQTTICRFE 181
6999***********************************99 PP
== domain 2 score: 0.8 bits; conditional E-value: 0.19
HTH_31 39 rgrpsaavlaalaralgldpaera 62
++ ps+++++ +a+ lgl+ + ++
gi|125490392|ref|NP_038661.2| 245 CPKPSLQQITHIANQLGLEKDVVR 268
678**************9988765 PP
>> Homeobox_KN Homeobox KN domain
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ? 12.3 0.0 3.5e-05 0.095 7 39 .. 244 276 .. 241 277 .. 0.91
Alignments for each domain:
== domain 1 score: 12.3 bits; conditional E-value: 3.5e-05
Homeobox_KN 7 hnPYPskevkeelakqTglsrkqidnWFiNaRr 39
+ P Ps +++ +a+q gl + + WF N R
gi|125490392|ref|NP_038661.2| 244 KCPKPSLQQITHIANQLGLEKDVVRVWFCNRRQ 276
56779*************************996 PP
>> DUF521 Protein of unknown function (DUF521)
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ? 9.6 0.1 9.4e-05 0.26 273 334 .. 221 280 .. 197 294 .. 0.77
Alignments for each domain:
== domain 1 score: 9.6 bits; conditional E-value: 9.4e-05
DUF521 273 adlaavleelnkakkeevdlvvlGcPhlsleeleelaellkgrkkkvsvelvvttsravlsk 334
+ +++ + +++++ +++ ++l cP sl++++++a++l +k v+++ + r+ ++
gi|125490392|ref|NP_038661.2| 221 QARKRKRTSIENRVRWSLETMFLKCPKPSLQQITHIANQLGLEK--DVVRVWFCNRRQKGKR 280
345666667778888899************************99..9999999988876554 PP
Internal pipeline statistics summary:
-------------------------------------
Query sequence(s): 1 (352 residues)
Target model(s): 13672 (2396357 nodes)
Passed MSV filter: 603 (0.0441047); expected 273.4 (0.02)
Passed bias filter: 465 (0.0340111); expected 273.4 (0.02)
Passed Vit filter: 44 (0.00321826); expected 13.7 (0.001)
Passed Fwd filter: 5 (0.000365711); expected 0.1 (1e-05)
Initial search space (Z): 13672 [actual number of targets]
Domain search space (domZ): 5 [number of targets reported over threshold]
# CPU time: 0.45u 0.22s 00:00:00.67 Elapsed: 00:00:00.24
# Mc/sec: 3514.66
//
|
