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/C31E10.8
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: C31E10.8 CE31423 S.pombe protein C4G8.04 like status:Partially_confirmed TR:Q93316 protein_id:CAB01733.2
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
TBC TBC domain 165.4 5e-47 1
PH PH domain 6.4 2.1 1
DUF593 Protein of unknown function, DUF593 6.3 2 1
Transp_Tc5_C Tc5 transposase C-terminal domain 4.1 9.4 1
tyrosinase Common central domain of tyrosinase 3.2 4.8 1
DBR1 Lariat debranching enzyme, C-terminal do 2.8 9 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
PH 1/1 16 41 .. 60 85 .] 6.4 2.1
DUF593 1/1 182 205 .. 1 24 [. 6.3 2
DBR1 1/1 185 207 .. 121 144 .. 2.8 9
tyrosinase 1/1 548 588 .. 336 379 .] 3.2 4.8
TBC 1/1 372 596 .. 1 341 [] 165.4 5e-47
Transp_Tc5_C 1/1 1137 1147 .. 1 11 [. 4.1 9.4
Alignments of top-scoring domains:
PH: domain 1 of 1, from 16 to 41: score 6.4, E = 2.1
*->drtlllqaeseeerkeWvkaiqsair<-*
+r++++ a++ee+++ W++ +++ +
C31E10.8 16 GRSIEFEADNEESCRAWLHTFCNRYS 41
DUF593: domain 1 of 1, from 182 to 205: score 6.3, E = 2
*->velLretvseerkalqdLyeELEe<-*
++lLre +++k ++d+++E ++
C31E10.8 182 LDLLREITNNQKKKIEDMKHEMNA 205
DBR1: domain 1 of 1, from 185 to 207: score 2.8, E = 9
*->dkeitesfkklfvDvsklrlNPqs<-*
++eit ++kk++ + +k+ N ++
C31E10.8 185 LREITNNQKKKI-EDMKHEMNAMT 207
tyrosinase: domain 1 of 1, from 548 to 588: score 3.2, E = 4.8
*->dwlishfrfydenpfppkrvkvrDcfdtsrpklGYdYedqptpw<-*
+w is f+ d +pf + + +Dcf + pk + +++++ +
C31E10.8 548 NWFISLFF--DAVPFNT-LLRIWDCFLLEGPKVLFRFAIVLIGK 588
TBC: domain 1 of 1, from 372 to 596: score 165.4, E = 5e-47
*->vrqgvpsslRgkVWklllgaqelnnclltdnfkgldlfglvpvllla
r+g+p+slR+ VWk l++ q ++ l+
C31E10.8 372 CRKGIPNSLRATVWKILINQQVED-------------------LKNV 399
dkdeYeellnk....nkektvqdqneKssvgirrldyVEAVEKHPLSDDN
+++Y ++l +++++ekt d+ K
C31E10.8 400 YGKYYYRNLCNiqggEDEKTYSDVHQK----------------------- 426
DKTKGSLekgsdekalklredldkIekDlsRTfpdeiffqtrlaeqqlkk
+I +Dl RT+p++++f +
C31E10.8 427 -----------------------QINLDLLRTMPNNVHFMSA-------- 445
dqdldaydkDEfddeddkneppsikqLrrlLvaYswknpqehlgYvQGMn
++ qL ++L a+ ++n q +gY+QGMn
C31E10.8 446 -------------------NSKGVTQLLQVLHAFCLHNSQ--IGYCQGMN 474
vilspLLlflkhgvdldeideeqAFwclvkLmdnylpqkyflndlsglne
++++ Llf e+AFw+l + +++y+ yf ++l g+++
C31E10.8 475 FLAATALLFV---------GPEDAFWFLIAVTERYFDKTYFDSNLTGAQA 515
dlrvLdslvkeslPeLyshlkkkenktgsgkKknllaldltllifafpwf
d+ vL+ l+++ P++ +hlk+ d++ +f+++wf
C31E10.8 516 DQEVLKGLLEVQHPKIMKHLKSL---------------DIDVASFTLNWF 550
LtlFarelPleivlrIwDilftyYlgshflifvalAiLkllkskllkh<-
+lF++ +P +++lrIwD+++++ g+ l+++a+ ++ +++++++
C31E10.8 551 ISLFFDAVPFNTLLRIWDCFLLE--GPKVLFRFAIVLIGKHEEEIISR 596
*
C31E10.8 - -
Transp_Tc5_C: domain 1 of 1, from 1137 to 1147: score 4.1, E = 9.4
*->dshdeFlTPsq<-*
+sh F TPs+
C31E10.8 1137 TSHKRFATPSK 1147