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: hmm4/seq/T13C5.5a
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: T13C5.5a CE04946 locus:bca-1 carbonic anhydrase status:Confirmed TR:Q22460 protein_id:AAK39291.1
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
Pro_CA Carbonic anhydrase 77.4 2e-21 2
DUF62 Protein of unknown function DUF62 4.5 2.9 1
Isochorismatase Isochorismatase family 3.6 4.6 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
DUF62 1/1 78 101 .. 254 279 .] 4.5 2.9
Isochorismatase 1/1 72 104 .. 122 153 .. 3.6 4.6
Pro_CA 1/2 23 111 .. 13 95 .. 57.9 9.8e-16
Pro_CA 2/2 182 223 .. 140 182 .] 19.5 0.00013
Alignments of top-scoring domains:
DUF62: domain 1 of 1, from 78 to 101: score 4.5, E = 2.9
*->EiavnmDnAaeklnvKeGddeievil<-*
E++vn Aa l vK+G+ i+ i+
T13C5.5a 78 EVSVNTEPAALELAVKRGG--IRHIV 101
Isochorismatase: domain 1 of 1, from 72 to 104: score 3.6, E = 4.6
*->iYSAFfahvcvlnTd..LdllLRelGidqLiitG<-*
+Y+AF v v ++ L+l++++ Gi++++++G
T13C5.5a 72 NYGAFS-EVSVNTEPaaLELAVKRGGIRHIVVCG 104
Pro_CA: domain 1 of 2, from 23 to 111: score 57.9, E = 9.8e-16
*->FkeLAkgQsPkyLviaCSDSRVcPslvtqlePGElFViRNvANlVP.
F+e + sP ++C DSR P tq++ G+ FV+RN++N++P+
T13C5.5a 23 FEEIKNNPSPTAVMFTCMDSRMLPTRFTQSQVGDMFVVRNAGNMIPd 69
.....pydknkyaGvGAaiEYAVaaLKVehIVVcGHSnCGGI<-*
++ +++ + + Aa E AV + hIVVcGHS C +I
T13C5.5a 70 apnygAFSEVSVNTEPAALELAVKRGGIRHIVVCGHSDCKAI 111
Pro_CA: domain 2 of 2, from 182 to 223: score 19.5, E = 0.00013
*->talekEnVnvqLaNLlTyPfVReglenGqkLaLhGWvYDvvsG<-*
+ l +nV+ qL N + + f +e le G +L +hG D+ +G
T13C5.5a 182 DKLSQINVLQQLINICSHEFLKEYLESG-RLHIHGMWFDIYKG 223