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/F31D5.4
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  F31D5.4  CE31799    status:Partially_confirmed TR:Q19930 protein_id:AAM98026.1

Scores for sequence family classification (score includes all domains):
Model          Description                              Score    E-value  N
--------       -----------                              -----    ------- ---
GPS            Latrophilin/CL-1-like GPS domain          67.6    2.6e-17   1
Glyco_hydro_48 Glycosyl hydrolase family 48               4.8        1.1   1
FixO           Cytochrome C oxidase, mono-heme subuni     3.0        6.2   1

Parsed for domains:
Model          Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
--------       ------- ----- -----    ----- -----      -----  -------
Glyco_hydro_48   1/1      45    67 ..   614   636 .]     4.8      1.1
FixO             1/1     430   448 ..   218   236 .]     3.0      6.2
GPS              1/1     625   674 ..     1    54 []    67.6  2.6e-17

Alignments of top-scoring domains:
Glyco_hydro_48: domain 1 of 1, from 45 to 67: score 4.8, E = 1.1
                   *->frYHRFWaqsdiAiAngvYaiLF<-*
                      f YH FWa++d  + +  +   F
     F31D5.4    45    FAYHLFWAECDSSLRLSLFNVGF    67

FixO: domain 1 of 1, from 430 to 448: score 3.0, E = 6.2
                   *->LqvLGTmVDfstynpkanr<-*
                      L++LGTmVD++  +  +
     F31D5.4   430    LNSLGTMVDNTVDQRIEYL    448

GPS: domain 1 of 1, from 625 to 674: score 67.6, E = 2.6e-17
                   *->snpiCvfWdeselslgvWstdrGCelvetskpshttCsCnHLTsFAv
                      +n+++++Wd ++   + W++dr+Ce+v +s+ ++ +++C+HLT+FA+
     F31D5.4   625    GNLKVTWWDIGR---NAWAEDRQCEIVSESD-GILEARCTHLTDFAI 667

                   Lmdvspn<-*
                   ++d+  n
     F31D5.4   668 IVDAALN    674



[DB home][top]