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/F32A5.2a
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: F32A5.2a CE28288 Peroxidase status:Partially_confirmed TR:Q95QH6 protein_id:AAK68356.1
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
An_peroxidase Animal haem peroxidase 732.3 2.2e-216 1
ShTK ShTK domain 41.1 8.2e-11 1
Androgen_recep Androgen receptor 2.4 2 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
ShTK 1/1 40 76 .. 1 44 [] 41.1 8.2e-11
Androgen_recep 1/1 682 698 .. 444 460 .] 2.4 2
An_peroxidase 1/1 463 971 .. 1 596 [] 732.3 2.2e-216
Alignments of top-scoring domains:
ShTK: domain 1 of 1, from 40 to 76: score 41.1, E = 8.2e-11
*->tCvDyvDpasdCaawaslGefCnnppYyqdfmkenCpkTCGfCs<-*
+C D++ +C++was Ge C++++ d+m +nC+ +C++C+
F32A5.2a 40 GCC---DHHEWCRFWASVGE-CKTNA---DWMVDNCQLACNTCT 76
Androgen_recep: domain 1 of 1, from 682 to 698: score 2.4, E = 2
*->ssWHtLFtAEEGqLYGP<-*
s+WHt+Ft E L
F32A5.2a 682 SAWHTIFTEEHNRLVTA 698
An_peroxidase: domain 1 of 1, from 463 to 971: score 732.3, E = 2.2e-216
*->rYRtlDGtCNNle.nPtwGaanspfrRllpPsYdDGidaPRsrsvdg
+ RtlDGtCNN ++ P Ga+++p++Rllp Yd +++P++ s
F32A5.2a 463 MFRTLDGTCNNMKgEPLRGASYRPYTRLLPTIYDNEVSEPVG-SLFT 508
spLPsaRvvSrkllsdrdfpipdpkfgtTlmvmqfGQFmaHDitktssst
+Ps+R++ r+l s++ ++ ++p + +++ mqfGQF++HD+ kt+ +
F32A5.2a 509 DARPSPREITRRLTSSQASV-ESPDY--NALIMQFGQFISHDMAKTTLVP 555
sPiEengealeCckcdspvspqlshpeCfPIpvppdDPffgqlgvrClpF
s ++ C +c q+++++C+++p+++dD +++++ +C+++
F32A5.2a 556 S---SK-----CNVC------QNITSRCMSVPITFDDSNANFRQAQCIRV 591
vRSlpaprtgpqcrlGgpReQlNqaTayiDaSpiYGssaeearklRaFqg
RS p+++ ++l +pR+QlN++T+yiDaSpiYGss +++k+R++
F32A5.2a 592 SRSSPICGS---GNL-KPRQQLNENTGYIDASPIYGSSVHDSKKFRDGNS 637
GlLktspdnggpvrlgrdllPldentqdcgsCvyprsaagsgCfiaGDdR
G+Lk ++n g + +lP+d+n +c+ ++ +s f+aGD R
F32A5.2a 638 GFLKLPMFN-G-----KAFLPFDQN--KCR-----NRGQCSVIFTAGDSR 674
vnlqpgLaaLHtLflREHNRiAdaLkalNPhWdDErlFQEARkIviAqmQ
vnl++gL+a Ht+f EHNR+ +a k+lNPhWd Erl+QEARk ++Aq+Q
F32A5.2a 675 VNLFVGLSAWHTIFTEEHNRLVTAFKRLNPHWDGERLYQEARKMIGAQVQ 724
hItynEyLPilLGrelmkkMrANWyGLLGfgLvygplksGyVgyyngYdp
I+y E+LP++LG ++ + Vg y gYd+
F32A5.2a 725 AIVYREWLPKVLGASFATV----------------------VGDYRGYDS 752
nvQYnpsianEFatAAFRfGHslipgkfplvdendkiyfgkpisnlsgsv
+v + ++anEF++AAFRfGH++i+++++++d + +++is g+
F32A5.2a 753 DV--DSTVANEFTSAAFRFGHGMIQEFYQRLDNS-----FRNISF--GAL 793
rLsdlffnpselylilaeggvdqllrGlatqpaGrmLnrvDrfitdevtn
+++ + +++ + eggvd+l+rG+ +q++ +++++ t vt+
F32A5.2a 794 PFQKGTLHSDV---LVNEGGVDPLIRGMFSQNV-----KRPQRVTTTVTE 835
kLFgtrgeppfglDLasinIqRGRDhGlpsYNeyRefCGLpratsFeDLt
++Fg + DL inIqRGRDhG p Y +yRe+CG+++a +Fe L+
F32A5.2a 836 NMFG-------STDLSTINIQRGRDHGHPAYVKYRELCGMGTAFNFEHLS 878
dlispdeenlarLkslYendvdDvDLyVGgllEkpvpGglvGpTfaCIia
+i + + ++++L+++Y+ +vd +DL+VG+llE+p+ glvGpT aCIi+
F32A5.2a 879 REILN-TGTRNKLQEIYG-SVDKIDLWVGALLEDPIIRGLVGPTVACIIG 926
eQFqRlrrGDRFwYEngrgkpsgFTpeVQLaeiRKtsLarliCdNtdvlt
+QF+R+r+GDRF+YEn p++F++ QL eiRK+sL+r+iCdNt+ +
F32A5.2a 927 PQFKRTRDGDRFYYEN----PGVFSRR-QLVEIRKSSLSRIICDNTNTIS 971
<-*
F32A5.2a - -