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: hmm3/seq/K01A6.1
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: K01A6.1 CE06034 protein-tyrosine-phosphatase like status:Partially_confirmed TR:Q21075 protein_id:CAA92965.1
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
PDZ PDZ domain (Also known as DHR or GLGF) 234.3 1.7e-66 5
WW WW domain 72.5 1.2e-18 3
ATP-synt_D ATP synthase subunit D 3.2 4.8 1
FH2 Formin Homology 2 Domain 2.5 9.4 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
WW 1/3 133 162 .. 1 30 [] 52.8 4.4e-13
WW 2/3 165 182 .. 1 18 [. 13.6 0.051
WW 3/3 270 281 .. 19 30 .] 12.7 0.088
ATP-synt_D 1/1 273 292 .. 104 123 .. 3.2 4.8
PDZ 1/5 382 438 .. 1 55 [. 30.6 3.1e-07
PDZ 2/5 495 572 .. 1 84 [] 15.1 0.008
FH2 1/1 571 611 .. 443 486 .] 2.5 9.4
PDZ 3/5 677 760 .. 1 84 [] 64.9 4.9e-17
PDZ 4/5 781 862 .. 1 84 [] 56.8 1e-14
PDZ 5/5 928 1009 .. 1 84 [] 85.9 5.1e-23
Alignments of top-scoring domains:
WW: domain 1 of 3, from 133 to 162: score 52.8, E = 4.4e-13
*->lppgWeerkdpdGrpYYyNhnTktTqWekP<-*
lpp+We +++++G +Y+++hnT+tT+W++P
K01A6.1 133 LPPNWETAYTENGDKYFIDHNTGTTTWDDP 162
WW: domain 2 of 3, from 165 to 182: score 13.6, E = 0.051
*->lppgWeerkdpdGrpYYy<-*
lppgWe++ d++ ++Y+
K01A6.1 165 LPPGWEQVDDQNYGTFYV 182
WW: domain 3 of 3, from 270 to 281: score 12.7, E = 0.088
*->NhnTktTqWekP<-*
Nh++++Tq+e+P
K01A6.1 270 NHINRKTQYERP 281
ATP-synt_D: domain 1 of 1, from 273 to 292: score 3.2, E = 4.8
*->verkpdergYgllgtsawvD<-*
++++ +er+Yg++g+sa +D
K01A6.1 273 NRKTQYERPYGFGGSSATID 292
PDZ: domain 1 of 5, from 382 to 438: score 30.6, E = 3.1e-07
*->evtlekevkrgglGfsikggsdk.....givvsevlpGsgaAeagGr
++++ k + glGf++ g+++ +++++ i v++vl G g+A+a G
K01A6.1 382 STKIVK--GAKGLGFTLIGNDSSsrgdeFIQVKSVLSG-GPAAANGV 425
LkeGDvIlsvNGq<-*
L+ GD+++ vNG+
K01A6.1 426 LRSGDILVRVNGR 438
PDZ: domain 2 of 5, from 495 to 572: score 15.1, E = 0.008
*->evtlekevkrgglGfsikggsdkgivvsevlpGsgaAeagGrLkeGD
e++++k +++g+Gf+i+ + + g ++++l + ++ + L eGD
K01A6.1 495 EIDIFK--GSEGFGFTIADNLN-GQRIKKIL---FPSQCPN-LMEGD 534
vIlsvNGqdvenmsheravlaikgsggsev....tLtvlRd<-*
I++ +G++v+ ++h + v+++++ +++ +++L+v+R+
K01A6.1 535 TIVELDGRNVRPIPHTQLVDMLRER---PIgyrgKLVVKRG 572
FH2: domain 1 of 1, from 571 to 611: score 2.5, E = 9.4
*->qskskqkeknplmdfeaepverdfkgvldsalleklrtgaafrr<-*
++++k+ ++p+ +f + + +++++dsa + ++r++ +r
K01A6.1 571 RGSPKTRSRTPSAAFRY---GEPQTNMMDSAAPLPVRSKTPAER 611
PDZ: domain 3 of 5, from 677 to 760: score 64.9, E = 4.9e-17
*->evtlekevkrgglGfsikggsdk..givvsevlpGsgaAeagGrLke
+v+l + k+ g+Gf++ gg + ++++ v ++ G gaAe +GrL+e
K01A6.1 677 TVSLIR--KPVGFGFRLLGGVESktPLSVGQIVIG-GAAEEDGRLQE 720
GDvIlsvNGqdvenmsheravlaikgsgg.sevtLtvlRd<-*
GD I++++G +ve+ sh +av ++ + ++++v+L v+R+
K01A6.1 721 GDEIVEIDGHNVEGASHSEAVVLLEAAAQnKHVKLIVRRP 760
PDZ: domain 4 of 5, from 781 to 862: score 56.8, E = 1e-14
*->evtlekevkrgglGfsikggsdk.givvsevlpGsgaAeagGrLkeG
+v l+++ ++ g+Gf + + + k+g v ++pG ++A++ GrL +G
K01A6.1 781 DVLLHRN-ENDGFGFVLMSSQHKnGSTVGQIQPG-SPAARCGRLSVG 825
DvIlsvNGqdvenmsheravlaikgsggsevtLtvlRd<-*
D++++vNG d+ ++sh + + +ik sg +v+Lt++++
K01A6.1 826 DRVIAVNGIDILSLSHPDTISLIKDSGL-SVRLTIAPP 862
PDZ: domain 5 of 5, from 928 to 1009: score 85.9, E = 5.1e-23
*->evtlekevkrgglGfsikggsdk...givvsevlpGsgaAeagGrLk
+vtle+ + g+Gfsi gg++ ++ +++v ++ + g+A+a+GrL+
K01A6.1 928 DVTLER--GTKGFGFSIRGGQEFgsmPLFVLRIADD-GPAKADGRLQ 971
eGDvIlsvNGqdvenmsheravlaikgsggsevtLtvlRd<-*
+GD+++ +NGq++++msh++a+ +ik+ v LtvlR+
K01A6.1 972 VGDQLTTINGQSTKGMSHDDAIRIIKQHT--MVNLTVLRN 1009
[output cut off at A = 10 top alignments]