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/T18D3.4
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  T18D3.4  CE31619  locus:myo-2 myosin heavy chain C status:Confirmed TR:Q20439 protein_id:CAA92183.2

Scores for sequence family classification (score includes all domains):
Model           Description                             Score    E-value  N
--------        -----------                             -----    ------- ---
Myosin_tail     Myosin tail                            1508.6          0   2
myosin_head     Myosin head (motor domain)             1451.3          0   1
Myosin_N        Myosin N-terminal SH3-like domain        66.3      8e-19   1
UPF0042         Uncharacterised P-loop ATPase protein     7.5        0.3   1
ABC_tran        ABC transporter                           4.8        3.1   1
DUF641          Plant protein of unknown function (DU     4.5        5.2   1
RICH            RICH domain                               4.0        7.2   1
Phage_min_cap   Phage minor capsid protein                3.7        3.3   1
Hanta_nucleocap Hantavirus nucleocapsid protein           3.1        1.9   1
DUF511          Protein of unknown function (DUF511)      2.6        7.5   1
aconitase       Aconitase family (aconitate hydratase     1.6        8.7   1

Parsed for domains:
Model           Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
--------        ------- ----- -----    ----- -----      -----  -------
Myosin_N          1/1      30    73 ..     1    48 []    66.3    8e-19
ABC_tran          1/1     167   181 ..     1    15 [.     4.8      3.1
UPF0042           1/1     168   181 ..     1    15 [.     7.5      0.3
aconitase         1/1     552   562 ..     1    11 [.     1.6      8.7
myosin_head       1/1      83   782 ..     1   734 []  1451.3        0
Hanta_nucleocap   1/1     912   945 ..     1    35 [.     3.1      1.9
Myosin_tail       1/2     979  1083 ..   175   281 ..   -15.9  5.9e+04
Phage_min_cap     1/1    1078  1095 ..   477   494 .]     3.7      3.3
RICH              1/1    1076  1122 ..    43    85 .]     4.0      7.2
DUF641            1/1    1096  1124 ..   113   141 .]     4.5      5.2
DUF511            1/1    1399  1410 ..   309   320 .]     2.6      7.5
Myosin_tail       2/2    1084  1943 ..     1   864 []  1534.1        0

Alignments of top-scoring domains:
Myosin_N: domain 1 of 1, from 30 to 73: score 66.3, E = 8e-19
                   *->KknVWVpDekegfikgeIkesteGdkVtVktedtrsGktvtvkkDdv
                      KknVW+pD+++g+i+g I+  t+Gd+VtV       G ++tvkkD v
     T18D3.4    30    KKNVWIPDSEDGYIEGVIT-KTAGDNVTVSIGQ---GAEKTVKKDVV 72

                   q<-*
                   q
     T18D3.4    73 Q    73

ABC_tran: domain 1 of 1, from 167 to 181: score 4.8, E = 3.1
                   *->GevlalvGpNGaGKS<-*
                      ++   ++G++GaGK+
     T18D3.4   167    NQSMLITGESGAGKT    181

UPF0042: domain 1 of 1, from 168 to 181: score 7.5, E = 0.3
                   *->ielvivtGrSGAGKs<-*
                       + +++tG SGAGK+
     T18D3.4   168    -QSMLITGESGAGKT    181

aconitase: domain 1 of 1, from 552 to 562: score 1.6, E = 8.7
                   *->TlyeKildaHl<-*
                      Tl+ K++d+Hl
     T18D3.4   552    TLAQKLIDQHL    562

myosin_head: domain 1 of 1, from 83 to 782: score 1451.3, E = 0
                   *->vEDmveLtyLnEpsvlhNLKkRYksdlIYTYsGlvLvsvNPYkrLpq
                      +EDm++Lt+Ln++svl+NLK RY+++lIYTYsGl++v++NPYkrLp
     T18D3.4    83    TEDMSNLTFLNDASVLYNLKARYAAMLIYTYSGLFCVVINPYKRLP- 128

                   iYteeiiakYrGKrryElPPHiFAiADeAYRsMlsdkeNQsilISGESGA
                   iYt+++  +++GKrr+E+PPH+FA++DeAYR+Ml+ +eNQs+lI+GESGA
     T18D3.4   129 IYTDSVARMFMGKRRTEMPPHLFAVSDEAYRNMLQNHENQSMLITGESGA 178

                   GKTEntKkvmqYlAaVsggnsgngeev.......psvkvgrvEdqILqsN
                   GKTEntKkv++Y+AaV++++++  ++ +  ++++ ++k+ ++EdqI+q+N
     T18D3.4   179 GKTENTKKVISYFAAVGAAQQETFGAKkaateedKNKKKVTLEDQIVQTN 228

                   PiLEAFGNAKTtRNNNSSRFGKyieIqFdktGkivGakIenYLLEKSRVv
                   P+LEAFGNAKT+RNNNSSRFGK+i+I+F k G+ ++++Ie+YLLEKSRV+
     T18D3.4   229 PVLEAFGNAKTVRNNNSSRFGKFIRIHFSKQGRVASCDIEHYLLEKSRVI 278

                   yQtegERNFHIFYQLLaGasqqnlkkeLkLtndpedYhYLnqggevkpcy
                   +Q++gER++HIFYQ  +   +  lkk+L+L + ++dY+++ q++     +
     T18D3.4   279 RQAPGERCYHIFYQVFSDYLPN-LKKDLLLNKPVKDYWFIAQAE-----L 322

                   tvdGiDDsegnveeFketrkAmdilGftdeeqrsIFrivAaILhlGNikF
                   ++dGi+D+    ee ++t++A+dil ft +e+++++r+vAa +h+GN+kF
     T18D3.4   323 IIDGINDK----EEHQLTDEAFDILKFTPTEKMECYRLVAAMMHMGNMKF 368

                   kqrrkeeaaipddnnadtkalekaaeLlGvdatelekALlsrriktGteg
                   kqr++ee+a+pd+    t ++e+aa ++G+d++e++kAL+++r+k+G+e+
     T18D3.4   369 KQRPREEQAEPDG----TDDAERAAKCFGIDSEEFLKALTRPRVKVGNEW 414

                   rkStvtkpqnveQAsyARDALAKalYsRLFdWIVnrINktLdfkakegqd
                       v+k+qn+eQ+++A++A+AK lYsR F+W+V+++N+tLd+k +++
     T18D3.4   415 ----VNKGQNIEQVNWAVGAMAKGLYSRIFNWLVKKCNQTLDQKGISR-- 458

                   asfIGVLDIyGFEIFekNSFEQLCINYvNEKLQQfFNhhmFklEQEEYkr
                   ++fIGVLDI+GFEIF++NSFEQL+IN+vNEKLQQfFNhhmF+lEQEEY r
     T18D3.4   459 DHFIGVLDIAGFEIFDFNSFEQLWINFVNEKLQQFFNHHMFVLEQEEYAR 508

                   EGIeWtfIdFgdNLQpcIDLIEkKsPpGILsLLDEeClfPkaqSGtDqtF
                   EGI+WtfIdFg++LQ+cI+LIEk+  +GI+++LDEeC++Pka   tD t+
     T18D3.4   509 EGIQWTFIDFGLDLQACIELIEKP--LGIIAMLDEECIVPKA---TDLTL 553

                   ldKLystfskhpahfekfsPrfrqkksg.ahFiikHYAGdVeYnvegFle
                   ++KL++++++++++fek  P+++++k+ +ahF+++HYAG+V+Yn+ ++le
     T18D3.4   554 AQKLIDQHLGKHPNFEK--PKPPKGKQAeAHFAMRHYAGTVRYNCLNWLE 601

                   KNKDpLfddlisllksSs.npllaeLFpdeetlagpfeadpsslskkrks
                   KNKDpL+d+++ ++k+S+++ l++e ++d++t++++++a+ + +  ++
     T18D3.4   602 KNKDPLNDTVVTVMKASKeHALIVEVWQDYTTQEEAAAAAAKGTAGAK-- 649

                   gskNkstgkktkksnfiTvGaqfKeslneLMktLsstnLPHFvRCIkPNe
                           kk+k+++f+Tv+++++esln+LM++L st+ PHF+RCI+PNe
     T18D3.4   650 --------KKGKSGSFMTVSMLYRESLNKLMTMLHSTH-PHFIRCIIPNE 690

                   kKkagvfDaslVlhQLrclGVLEgiRIrRaGFPnRitfdeFlqRYriLap
                   kK++gv+Da lVl+QL+c+GVLEgiRI+R+GFPnR ++++F+qRY++La+
     T18D3.4   691 KKASGVIDAGLVLNQLTCNGVLEGIRICRKGFPNRTLHPDFVQRYALLAA 740

                   ktwPkwsgdakkgeknEIvaceklLqsLnlDkg...eeyrfGkTKIFFR<
                   ++   ++ dakkg       + ++L +L  +k+ ++e +r+G TK+FF+
     T18D3.4   741 DESIIGKTDAKKG-------SALMLARLVKEKKleeENFRVGLTKVFFK  782

                   -*

     T18D3.4     -    -

Hanta_nucleocap: domain 1 of 1, from 912 to 945: score 3.1, E = 1.9
                   *->msnlkElqeeitaHEqqLviARqKlkdAEkavEvD<-*
                          +E+ e +ta  +q v+   Kl dA k++Ev+
     T18D3.4   912    -GSTREVEERMTAMNEQKVALEGKLADASKKLEVE    945

Myosin_tail: domain 1 of 2, from 979 to 1083: score -15.9, E = 5.9e+04
                   *->KlnaEKkakqlEsQlsElqvKldElqRqlnDltsqKsRLqsENsdLt
                      K +aEK ak  E+Q+  lq  ++ +  +++ l+  ++     N  Lt
     T18D3.4   979    KVEAEKNAK--EHQIRALQDEMRQQDENISKLNKERKNQEEQNKKLT 1023

                   rqleEaEaqvsqlsklKsqlesQLEeAkRslEEEsReRanLqaqlrnleh
                    +l++aE q+ +++klK+ l + LE+  +++E E R+Ra      r  e
     T18D3.4  1024 EDLQAAEEQNLAANKLKAKLMQSLEDSEQTMEREKRNRADMDKNKRKAEG 1073

                   DlDslrEqlE<-*
                   +l  ++E+lE
     T18D3.4  1074 ELKIAQETLE    1083

Phage_min_cap: domain 1 of 1, from 1078 to 1095: score 3.7, E = 3.3
                   *->AdEwlAkIakEkqteipe<-*
                      A+E l++++k k+++++
     T18D3.4  1078    AQETLEELNKSKSDAENA    1095

RICH: domain 1 of 1, from 1076 to 1122: score 4.0, E = 7.2
                   *->aIKtkYLyeLdvlkskseKs.....aelqelvsktkaeLDaAFekFK
                      +I +++L eL+++ks  e + +++++el +l +k   e  aA+ k +
     T18D3.4  1076    KIAQETLEELNKSKSDAENAlrrkeTELHTLGMKLEDE-QAAVAKLQ 1121

                   K<-*
                   K
     T18D3.4  1122 K    1122

DUF641: domain 1 of 1, from 1096 to 1124: score 4.5, E = 5.2
                   *->vkaKdseIdsLkekleealrankklekrl<-*
                      ++ K++e + L +kle +++a  kl k +
     T18D3.4  1096    LRRKETELHTLGMKLEDEQAAVAKLQKGI    1124

	[output cut off at A = 10 top alignments]


[DB home][top]