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/D2013.1
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: D2013.1 CE30340 locus:rab-39 GTP binding protein (YPT1 subfamily) status:Partially_confirmed TR:Q18969 protein_id:CAA87774.2
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
ras Ras family 225.5 6.7e-66 1
arf ADP-ribosylation factor family 3.7 6.1 2
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
arf 1/2 67 90 .. 62 85 .. 1.7 23
arf 2/2 164 186 .. 169 191 .] 2.0 19
ras 1/1 19 229 .] 1 198 [] 225.5 6.7e-66
Alignments of top-scoring domains:
arf: domain 1 of 2, from 67 to 90: score 1.7, E = 23
*->ikFtvWDvGGQeslRPlWrnYfpn<-*
+k++ WD GQe+ R+ ++Y++n
D2013.1 67 VKLQLWDTAGQEKFRSITKSYYRN 90
arf: domain 2 of 2, from 164 to 186: score 2.0, E = 19
*->csAvkGeGLyEGldWLsnnikkr<-*
+s ++G+++ E++ +++ i++r
D2013.1 164 TSSRTGDNVNEAFHMIAQEIQNR 186
ras: domain 1 of 1, from 19 to 229: score 225.5, E = 6.7e-66
*->KlVliGDsGVGKSsLliRFtdnkFveeYipTIGvDFytktvevD.Gk
l++iGDs VGKSsLl Ft++k e ++pT+GvDFy + +e+ +G
D2013.1 19 RLIVIGDSTVGKSSLLRYFTEGKMAEISDPTVGVDFYARMIELRpGY 65
tvkLqIWDTAGQERFralrpaYYRgAqGfLLVYDITsrdSFenvkkWlee
+vkLq+WDTAGQE+Fr+++ +YYR+++G+L YD T+r+SFe+v +W++e
D2013.1 66 RVKLQLWDTAGQEKFRSITKSYYRNSVGVLAIYDTTNRESFEHVENWVKE 115
ilrhad..kdenvpivLVGNKcDLeddedleltegqkRvVsteeGealAk
+++++ ++ LVG+K+D+ +R+V +eeGe +Ak
D2013.1 116 AALNLGgpSPSKCVFQLVGTKSDMDS----------QRQVNYEEGEYFAK 155
elgAlpFmETSAKtntNVeeaFeelareilkkvse...............
++ ++F+ETS t++NV+eaF ++a+ei ++v++++ ++ ++ ++ +++
D2013.1 156 YHK-MKFIETSSRTGDNVNEAFHMIAQEIQNRVDDgelrpvdgweglktg 204
......vnvnldqpakkkkskCcil<-*
++++v +++++ +++++++ c +
D2013.1 205 imrsqsVCLSERSFPQNSSAGACGC 229