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: hmm2/seq/F37A4.4
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: F37A4.4 CE24953 status:Partially_confirmed SW:P41882 protein_id:AAA50632.2
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
WSN Domain of unknown function 106.1 2.9e-28 1
BRCT BRCA1 C Terminus (BRCT) domain 59.2 8.9e-17 2
ank Ankyrin repeat 33.1 1.3e-07 1
RelA_SpoT Region found in RelA / SpoT proteins 9.1 0.31 1
fusion_gly Fusion glycoprotein F0 5.3 1.5 1
Mnd1 Mnd1 family 5.3 2 1
mce mce related protein 4.6 6.1 1
vwa von Willebrand factor type A domain 3.7 3.1 1
UPF0128 Uncharacterised protein family (UPF0128) 2.9 9.8 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
WSN 1/1 27 95 .. 1 72 [] 106.1 2.9e-28
UPF0128 1/1 226 244 .. 206 224 .] 2.9 9.8
vwa 1/1 440 474 .. 151 189 .] 3.7 3.1
Mnd1 1/1 687 696 .. 196 205 .] 5.3 2
fusion_gly 1/1 753 780 .. 489 516 .] 5.3 1.5
mce 1/1 762 795 .. 1 37 [. 4.6 6.1
ank 1/1 856 888 .. 1 33 [] 33.1 1.3e-07
RelA_SpoT 1/1 910 925 .. 1 16 [. 9.1 0.31
BRCT 1/2 932 1026 .. 1 99 [] 52.2 1.2e-14
BRCT 2/2 1051 1154 .. 1 99 [] 6.9 0.83
Alignments of top-scoring domains:
WSN: domain 1 of 1, from 27 to 95: score 106.1, E = 2.9e-28
*->skLqtvienlsilaRVtNAIsLqaglingsiksddvIaELLniGtin
s+Lq++++ lsil+RVtNAI+Lqa+++++++k+++vI+ELL+ ++
F37A4.4 27 SQLQQIYDELSILSRVTNAIALQAAALSKTVKIREVITELLK---VD 70
psqlselinvdlkkiveglngLkei<-*
++++s+l+++d++++v++l++L+++
F37A4.4 71 NGNFSNLLSLDPAHLVKNLDELHKK 95
UPF0128: domain 1 of 1, from 226 to 244: score 2.9, E = 9.8
*->aLFDvDIAekvdkiLgiLk<-*
aL +DIA+k+ L +k
F37A4.4 226 ALNELDIANKIPNMLQKFK 244
vwa: domain 1 of 1, from 440 to 474: score 3.7, E = 3.1
*->dadeeeLr.iasepeseghvfyvedfddlwlqlkdlqeal<-*
d+ +eeL++++ + + + v ++++fdd + ++++++
F37A4.4 440 DVYFEELEkLTPSN-IKEQVEKITNFDDFL----KIINKF 474
Mnd1: domain 1 of 1, from 687 to 696: score 5.3, E = 2
*->GIPEdFdelt<-*
GIPE+F +++
F37A4.4 687 GIPEVFSYID 696
fusion_gly: domain 1 of 1, from 753 to 780: score 5.3, E = 1.5
*->svgvsktttsvliiiivlvvllliiill<-*
+v+v++ +ts +++++++v+++l++++l
F37A4.4 753 KVSVQSESTSPIVVVLICVAIVLVLVIL 780
mce: domain 1 of 1, from 762 to 795: score 4.6, E = 6.1
*->lglfllvalvvllvvllwfglrslPlprfsgtgpgry<-*
++++++ v++++vl+ ++++ +++f+++g+++y
F37A4.4 762 SPIVVVLICVAIVLVLVILAIV---GYGFTSNGRNQY 795
ank: domain 1 of 1, from 856 to 888: score 33.1, E = 1.3e-07
*->dGnTPLHlAakrgnlevvklLlqrGAdlnaqdd<-*
+GnT+LH+A +rg+ ++v+ L+++GAd +
F37A4.4 856 YGNTALHVATRRGYQNLVEILIKHGADRSFLNP 888
RelA_SpoT: domain 1 of 1, from 910 to 925: score 9.1, E = 0.31
*->gRvKhiySIyrKmqrk<-*
R K+i+SIy+K ++k
F37A4.4 910 ERFKSIESIYNKYRKK 925
BRCT: domain 1 of 2, from 932 to 1026: score 52.2, E = 1.2e-14
*->eklfsgltfvltgtknaaviieamldseer.elkelieklGGkvt.s
+ +f +f+++ d ++ + e+ k+ +++
F37A4.4 932 PEKFPVSSFHIYI-----------EDRTDDnVTNEFTTKFQSITSdE 967
lskktthvivgenagsklekalk....aialgipivteeWlldcikegkl
+ +tthv+v+++++ le+ +++ +i +g +iv+ W++dc++ kl
F37A4.4 968 AMITTTHVVVKTTEDGILETDDLnlliWIFHGSIIVRDTWMVDCLRDEKL 1017
lpe.dpyll<-*
++ +++yl+
F37A4.4 1018 IEKdCDYLV 1026
BRCT: domain 2 of 2, from 1051 to 1154: score 6.9, E = 0.83
*->eklfsgltfvltgtknaaviieamldseer.elkelieklGGkvt.s
++++g+++ l+ + + +l ++i+ GG+ +
F37A4.4 1051 TPFLYGVHVALCM-----------KNCPYLaSLTAIIQGQGGTMLdK 1086
l.........skktthvivgenagsklekalk.....aialgipivteeW
+++++ +++s+ h +g+ + +++ ++++++ + + te+
F37A4.4 1087 FpdkdafnkgSHPYLHKNLGPIFLLHDGTGDLdlyrsDPDKMFTLFTEQQ 1136
lldcikegkllpedpyll<-*
+ d++ ++++ ++ ++
F37A4.4 1137 FMDLLFKREINKDTNPKI 1154
[output cut off at A = 10 top alignments]