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/F26D2.10
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: F26D2.10 CE18640 Ank repeat, BRCA1 C Terminus (BRCT) domain status:Partially_confirmed TR:Q9XV62 protein_id:CAB04180.1
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
WSN Domain of unknown function 97.6 7.2e-26 2
BRCT BRCA1 C Terminus (BRCT) domain 51.9 1.4e-14 1
ank Ankyrin repeat 31.8 2.9e-07 1
Endonuclease DNA/RNA non-specific endonuclease 5.6 2.2 1
Neuregulin Neuregulin family 3.4 1.5 1
CO_dh CO dehydrogenase beta subunit/acetyl-CoA 3.1 7.5 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
WSN 1/2 31 99 .. 1 72 [] 102.2 3.5e-27
WSN 2/2 136 160 .. 48 72 .] -0.2 2.5e+02
Neuregulin 1/1 824 849 .. 15 40 .. 3.4 1.5
ank 1/1 911 943 .. 1 33 [] 31.8 2.9e-07
CO_dh 1/1 957 974 .. 163 180 .] 3.1 7.5
Endonuclease 1/1 945 982 .. 218 256 .] 5.6 2.2
BRCT 1/1 987 1081 .. 1 99 [] 51.9 1.4e-14
Alignments of top-scoring domains:
WSN: domain 1 of 2, from 31 to 99: score 102.2, E = 3.5e-27
*->skLqtvienlsilaRVtNAIsLqaglingsiksddvIaELLniGtin
+Lq+++e+++il+RVtNAI+Lqag+i i+++++I+ELL+ ++
F26D2.10 31 ANLQLISEKFRILSRVTNAIYLQAGAIKNEINARELISELLH---MP 74
psqlselinvdlkkiveglngLkei<-*
+ ql++++nv++ + ++l +++++
F26D2.10 75 SNQLNTVVNVNVEVASSELRRMHQE 99
WSN: domain 2 of 2, from 136 to 160: score -0.2, E = 2.5e+02
*->psqlselinvdlkkiveglngLkei<-*
+s lse+ ++ +k++++++++++++
F26D2.10 136 SSLLSEVSSISTKNLLNETKSIAKD 160
Neuregulin: domain 1 of 1, from 824 to 849: score 3.4, E = 1.5
*->iCi.dLLvvGimCvvdAYCKtKKqRkK<-*
iCi+ +L +G Cvv AY tK R K
F26D2.10 824 ICIgSILLIGLFCVV-AYGFTKNGRAK 849
ank: domain 1 of 1, from 911 to 943: score 31.8, E = 2.9e-07
*->dGnTPLHlAakrgnlevvklLlqrGAdlnaqdd<-*
GnT+LH A k+ +e+vk L+++GAdl + +
F26D2.10 911 FGNTALHAATKSALPEIVKILIEHGADLTLLNT 943
CO_dh: domain 1 of 1, from 957 to 974: score 3.1, E = 7.5
*->pNlskeeeeylemLeEml<-*
p+ sk++ e+le++++++
F26D2.10 957 PDNSKDKIEKLEKIRKLY 974
Endonuclease: domain 1 of 1, from 945 to 982: score 5.6, E = 2.2
*->nepidvaelpltnfdvsvdeiekltglvffdnldpdvqk<-*
n +++ + +pl+++d s+d+iekl ++ ++n ++ +++
F26D2.10 945 NKTPE-QMIPLDSPDNSKDKIEKLEKIRKLYNKANGKKY 982
BRCT: domain 1 of 1, from 987 to 1081: score 51.9, E = 1.4e-14
*->eklfsgltfvltgtknaaviieamldseer.elkelieklGGkvt.s
+ +f +f++ + ++ ++ +++ek+ + +++
F26D2.10 987 PQKFPATSFHIWM-----------DEITDNkITNDFMEKFQAITSdE 1022
lskktthvivgenagsklekalk....aialgipivteeWlldcikegkl
++tthviv+++++ +l + ++ + +g +iv+e W++dc++ kl
F26D2.10 1023 AAPNTTHVIVKTDENGVLDTDNLdllfWMFNGAIIVKEKWMIDCLNDEKL 1072
lpe.dpyll<-*
++++++yl+
F26D2.10 1073 IEQdSKYLV 1081