Motif confirmation by HMMer
hmmpfam - search a single seq against HMM database
HMMER 2.1.1 (Dec 1998)
Copyright (C) 1992-1998 Washington University School of Medicine
HMMER is freely distributed under the GNU General Public License (GPL).
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
HMM file: /home/niguts/usr02/tshini/ykclst/db/Pfam_fs
Sequence file: hmm1/seq/C18H2.1
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: C18H2.1 CE29670 status:Partially_confirmed TR:Q18102 protein_id:AAA50664.2
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
WSN Domain of unknown function 189.8 4.4e-53 2
ank Ankyrin repeat 39.9 1.4e-09 1
BRCT BRCA1 C Terminus (BRCT) domain 19.2 0.00014 2
DUF278 DUF278 15.7 0.0054 1
synaptobrevin Synaptobrevin 5.3 5.9 2
MtN3_slv MtN3/saliva family 4.7 7.6 1
Baculo_E56 Baculoviral E56 protein, specific to OD 3.9 3.3 1
LacD Tagatose 1,6-diphosphate aldolase, (Lac 3.3 3.7 1
Fragilysin Fragilysin metallopeptidase (M10C) ente 2.6 6 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
WSN 1/2 230 298 .. 1 72 [] 92.4 2.1e-24
LacD 1/1 499 522 .. 309 332 .] 3.3 3.7
synaptobrevin 1/2 963 981 .. 71 89 .] -6.4 1e+04
Baculo_E56 1/1 966 990 .. 297 321 .] 3.9 3.3
ank 1/1 1053 1085 .. 1 33 [] 39.9 1.4e-09
BRCT 1/2 1168 1223 .. 49 99 .] 16.4 0.001
BRCT 2/2 1248 1287 .. 1 49 [. 2.9 14
WSN 2/2 1385 1453 .. 1 72 [] 101.0 8e-27
synaptobrevin 2/2 1447 1472 .. 33 58 .. 1.5 66
DUF278 1/1 1487 1523 .. 148 184 .. 15.7 0.0054
Fragilysin 1/1 1555 1564 .. 373 382 .] 2.6 6
MtN3_slv 1/1 1729 1744 .. 74 89 .] 4.7 7.6
Alignments of top-scoring domains:
WSN: domain 1 of 2, from 230 to 298: score 92.4, E = 2.1e-24
*->skLqtvienlsilaRVtNAIsLqaglingsiksddvIaELLniGtin
+L +++++l+i+ RVtNAI Lqa+ i++ +k++dv+aEL++ ++
C18H2.1 230 PALVLAYNKLAIITRVTNAIVLQAEGIRKLVKVRDVVAELIR---SP 273
psqlselinvdlkkiveglngLkei<-*
++++++l++vd+++i ln+++e+
C18H2.1 274 NASFYDLMDVDPTNIFPVLNEIIEA 298
LacD: domain 1 of 1, from 499 to 522: score 3.3, E = 3.7
*->tvgfknldeLnkVldetaspwtek<-*
t gf+n+ ++ kV ++ s+w +
C18H2.1 499 TLGFPNIGDIAKVSEDLKSSWFKS 522
synaptobrevin: domain 1 of 2, from 963 to 981: score -6.4, E = 1e+04
*->ivglvvvlllliiivylcg<-*
++gl+++++++++iv++++
C18H2.1 963 AIGLCIGIFIILFIVVFVV 981
Baculo_E56: domain 1 of 1, from 966 to 990: score 3.9, E = 3.3
*->LviGgVlllilIgyvvFKvlnkkrg<-*
L iG+ + l ++++vv+++ ++r+
C18H2.1 966 LCIGIFIILFIVVFVVYGRTPSGRQ 990
ank: domain 1 of 1, from 1053 to 1085: score 39.9, E = 1.4e-09
*->dGnTPLHlAakrgnlevvklLlqrGAdlnaqdd<-*
+GnTPLH+A+k+g++e+v+ L+++GAd + +
C18H2.1 1053 YGNTPLHVASKSGYPEIVEILIKNGADRTLLNA 1085
BRCT: domain 1 of 2, from 1168 to 1223: score 16.4, E = 0.001
*->ktthvivgenagsklekalk....aialgipivteeWlldcikegkl
++th+iv+++++ l + + + ++ g+ i++ +W++dc+ + l
C18H2.1 1168 SITHFIVKTDENGTLVTDAFqllsWVLSGVLIMKNSWMKDCLGNESL 1214
lpe.dpyll<-*
+ +++y +
C18H2.1 1215 IKKdYNYFV 1223
BRCT: domain 2 of 2, from 1248 to 1287: score 2.9, E = 14
*->eklfsgltfvltgtknaaviieamldseer.elkelieklGGkvt.s
+ +++g+ v+ + d ++ + + +++ ++GG+ ++
C18H2.1 1248 MLFLHGVYAVVVS-----------IDCSNIgAITAIVTNNGGVMCeK 1283
lskk<-*
+ +k
C18H2.1 1284 FPEK 1287
WSN: domain 2 of 2, from 1385 to 1453: score 101.0, E = 8e-27
*->skLqtvienlsilaRVtNAIsLqaglingsiksddvIaELLniGtin
s L++v + l ++aRVtNAIsL+ag+i++++k+ dv++E L i+
C18H2.1 1385 SSLEKVRDVLTMVARVTNAISLEAGAIRKTLKLTDVVSEVLD---ID 1428
psqlselinvdlkkiveglngLkei<-*
p++l+el+n+d+++++++l++++++
C18H2.1 1429 PGKLDELLNLDSAGLNAHLEKMIHR 1453
synaptobrevin: domain 2 of 2, from 1447 to 1472: score 1.5, E = 66
*->LdeLvDktdnLqassdqFrksakklr<-*
L+++ ++ ++L + +qF+++++kl+
C18H2.1 1447 LEKMIHRYQKLTQDTEQFKDKLTKLK 1472
DUF278: domain 1 of 1, from 1487 to 1523: score 15.7, E = 0.0054
*->ikkMDRaatEaeLwkWEfYphFrYYRetLGqFEAklR<-*
+kk+ t+ eL kWE+Yp++ YYR L+ E +R
C18H2.1 1487 DKKLNEEDTDRELTKWEYYPFIKYYRVNLADWERRKR 1523
[output cut off at A = 10 top alignments]