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/T22G5.5
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: T22G5.5 CE31996 serine palmitoyltransferase status:Partially_confirmed TR:Q9XVI6 protein_id:CAB03390.1
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
aminotran_1_2 Aminotransferase class I and II 98.8 2.2e-26 1
Dioxygenase_N Catechol dioxygenase N terminus 5.2 6.5 1
DUF507 Protein of unknown function (DUF507) 3.6 5.8 1
DUF98 Protein of unknown function DUF98 2.6 7.2 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
DUF98 1/1 274 291 .. 156 173 .] 2.6 7.2
DUF507 1/1 418 425 .. 176 183 .] 3.6 5.8
aminotran_1_2 1/1 167 495 .. 1 350 [] 98.8 2.2e-26
Dioxygenase_N 1/1 485 503 .. 60 81 .] 5.2 6.5
Alignments of top-scoring domains:
DUF98: domain 1 of 1, from 274 to 291: score 2.6, E = 7.2
*->gkVLieItEiFpverfrv<-*
+kVLi I i+ +e+ v
T22G5.5 274 NKVLIIIEGIYSMEGTIV 291
DUF507: domain 1 of 1, from 418 to 425: score 3.6, E = 5.8
*->ELrkkGll<-*
ELrk+G+l
T22G5.5 418 ELRKRGFL 425
aminotran_1_2: domain 1 of 1, from 167 to 495: score 98.8, E = 2.2e-26
*->klaggtrgpvGlndvrvdleqvlagtGalevaalfisaflapggtkl
gg+r ++G+ + + e ++a+ e a +f ++f+ ++
T22G5.5 167 INCGGSRQEIGNHVAHKSVESTIAQYLNVEDAIVFPMGFATNSM--- 210
devlvpdPtypiysdvlnhaggidGvrlypyplrssnhnDfkaleeaLee
+ + d + i sd lnha++++G+rl+ +++ +++hnD + e++L++
T22G5.5 211 NIPSLVDKGSLILSDRLNHASLVTGCRLSGAHTVVFRHNDASDCERKLRD 260
apeg..........kvvlvanPhNPtGtdltlaelekiadlakehPnill
a+ g +++++++ +kv++++ ++ ++ +t+ +l+++ +++k++ n++l
T22G5.5 261 ALCGvspktgekynKVLIII--EGIYSMEGTIVNLPAFIAVKKKY-NCYL 307
lvDEAYagfvfGgldedgasiaelldeypyanaenkdnllvvqsfSKnfG
+ DEA++ ++ G+ g+++ae+ + + + +++ +G
T22G5.5 308 FLDEAHSIGAVGP---SGRGVAEYWG-------------CNPRDIDIMMG 341
laGkRvGiiidiyangaaggivagsaasadrvssqlrallralysappav
++ k + +agg++ gs+++ d+++ +++ ++ ++pp
T22G5.5 342 TLTKSFA---------SAGGYMGGSKKVIDHIRRYSAGTCYGVTMSPPLI 382
ga.aivalilqdker...nrerwldelkkmlglril.sragfvlwldpsg
++ ++++li++ k+ ++ +r++ ++l++++++r ++gf ++++ +
T22G5.5 383 AQvERAVLIMSGKDGtdiGRQKAIQLLENSRYFRKElRKRGFLVYGNNDS 432
llgkrvpPdlshfekqaglfsyllleeegvvvpgsefptyppgwgRisva
++ + + ++++++ +fs ++l ++ + + + ++p++p ++R++++
T22G5.5 433 PVVPLMT----FYITKVVEFSRRML-KHNIGIVAVGYPATPLLEARVRFC 477
glteetleelleairefl<-*
+ +++t e+l+ ++++
T22G5.5 478 LSADHTKEHLDYILEAVE 495
Dioxygenase_N: domain 1 of 1, from 485 to 503: score 5.2, E = 6.5
*->EhfLDlrvdaadaaaPgreGgT<-*
Eh LD +++a+ ++ g e gT
T22G5.5 485 EH-LDYILEAV-EQV-GMETGT 503