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/T23G11.2
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: T23G11.2 CE16458 locus:gna-2 Acetyltransferase (GNAT) family (2 domains) status:Partially_confirmed TR:O45811 protein_id:CAB03416.1
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
Acetyltransf Acetyltransferase (GNAT) family 83.5 4.1e-28 2
Transglutamin_N Transglutaminase family 3.7 7 2
Tropomyosin Tropomyosin 2.5 9.8 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
Tropomyosin 1/1 99 128 .. 208 237 .] 2.5 9.8
Transglutamin_N 1/2 129 141 .. 119 131 .] 3.3 8.6
Acetyltransf 1/2 70 152 .. 1 80 [] 46.3 5e-15
Transglutamin_N 2/2 306 317 .. 120 131 .] 0.3 69
Acetyltransf 2/2 246 328 .. 1 80 [] 38.0 4.2e-12
Alignments of top-scoring domains:
Tropomyosin: domain 1 of 1, from 99 to 128: score 2.5, E = 9.8
*->dLEDkLaneKekykaiseeLDqTlleLnnm<-*
+ ED ++ eK + k++ L + l+e+ +
T23G11.2 99 RVEDVVVDEKMRGKKLGALLNEVLVEMAKT 128
Transglutamin_N: domain 1 of 2, from 129 to 141: score 3.3, E = 8.6
*->iGkYrLtvevstp<-*
iG+Y+L++e +t+
T23G11.2 129 IGVYKLSLECKTE 141
Acetyltransf: domain 1 of 2, from 70 to 152: score 46.3, E = 5e-15
*->lvaeedgelvGfaglspide.......vaeieglaVdpeyRgkGiGt
l e +++G a+l ++ ++ ++++ +e+++Vd+ Rgk +G
T23G11.2 70 LEHIESSKIIGAATLLIEFKyiheagtRGRVEDVVVDEKMRGKKLGA 116
aLlealeeyareelGlkrieleveedNeaAialYeklGFk<-*
L e l+e a+ G+++++le+ + ++Y k+G+
T23G11.2 117 LLNEVLVEMAKT-IGVYKLSLECKTELI---PFYNKFGYS 152
Transglutamin_N: domain 2 of 2, from 306 to 317: score 0.3, E = 69
*->GkYrLtvevstp<-*
G+Y+L++e +t+
T23G11.2 306 GVYKLSLECKTE 317
Acetyltransf: domain 2 of 2, from 246 to 328: score 38.0, E = 4.2e-12
*->lvaeedgelvGfaglspide.......vaeieglaVdpeyRgkGiGt
l + ++vG a+l+++ ++ + ++ +e+++Vd Rgk +G
T23G11.2 246 LEDVNSSKIVGAATLVVELKyihecglRGRVEDVVVDLTMRGKRLGI 292
aLlealeeyareelGlkrieleveedNeaAialYeklGFk<-*
eal++ ar+ lG+++++le+ + ++Y k+G+
T23G11.2 293 LINEALVKMARS-LGVYKLSLECKTELI---PFYNKFGYN 328