# hmmscan :: search sequence(s) against a profile database # HMMER 3.1b2 (February 2015); http://hmmer.org/ # Copyright (C) 2015 Howard Hughes Medical Institute. # Freely distributed under the GNU General Public License (GPLv3). # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - # query sequence file: AA/23726_FusoPortal_Gene_597.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_597 [L=658] Description: # 661428 # 663401 # 1 # ID=1_597;partial=00;start_type=ATG;rbs_motif=AGGAGG;rbs_spacer=3-4bp;gc_cont=0.302 Scores for complete sequence (score includes all domains): --- full sequence --- --- best 1 domain --- -#dom- E-value score bias E-value score bias exp N Model Description ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1e-32 114.5 52.2 1.3e-32 114.2 52.2 1.0 1 TonB_dep_Rec TonB dependent receptor 5.4e-24 84.7 2.1 7.2e-24 84.3 0.8 2.0 2 Plug TonB-dependent Receptor Plug Domain ------ inclusion threshold ------ 0.13 12.3 0.0 0.3 11.1 0.0 1.6 1 AAA_2 AAA domain (Cdc48 subfamily) Domain annotation for each model (and alignments): >> TonB_dep_Rec TonB dependent receptor # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 114.2 52.2 2.4e-36 1.3e-32 9 470 .] 184 656 .. 177 656 .. 0.75 Alignments for each domain: == domain 1 score: 114.2 bits; conditional E-value: 2.4e-36 TSXEEEEEEEEEEEEEETTTCEEEEEEECEEEEEEEEECCCEECCCEEEEEEEEEEEEECEEEECCEEEEEEEEEEEEEEEEEE CS TonB_dep_Rec 9 lsddnllanngdylrytdggdlgsalsrgylgtnplyygdnysetdnkrlrlgydydlddnlswklslraglnYsstdgdsdyn 92 ++d+ ll+ + +++ d +d ++++ +++ +++ + ++ + dn + ++y+++ +d++ +++sl+ ++++ ++ ++ 23726_FusoPortal_Gene_597 184 IGDKFLLNASYSGNSSKDYRDRSPQYENKKDKRDS-LWLRGKYLLDNGSIAINYNHS-EDKDYYTGSLS-KEQFDKNPRQVGSW 264 56677788889999999999999999999999999.*****9***************.77777777777.44455444444444 PP EEEEEE......................CEEEEEEEEEEEEEEEEEEEEEEEEECE........EE..CEEEEEEEEEE.EEEE CS TonB_dep_Rec 93 nsglgg......................grsynrdndyyglgnvlaaeldlngdlw........fG..snladgldyrk.lsdd 143 ++++ g + y+ + + ++ ++++++ + + s + gldyr+ +d+ 23726_FusoPortal_Gene_597 265 SGYTYGindiinakynqkindkidifltS-GYYHNKNKFQKNSTSEY--------FirpevkltYAkdSYVTLGLDYRDgKRDF 339 44444444444334333444444444422.22222222222222222........1333333333334444444****644442 PP ECCCEESCCTCECEEEEEEEEEEEEEEEEEEEEETTEEEEEEEEEEEETTESSSS.............EEEEEEEEE.EEEEEE CS TonB_dep_Rec 144 gydaydplnpsgdasdsvdrdsyglyaqdnikltdrltltlGlRydhysvfgkdlrag..........afsrsgsaf.sPsvgl 216 d + + ++ ++++r+s++ y+ + ++ +++++t G+R++ +++ ++ + + s ++s+ s +g+ 23726_FusoPortal_Gene_597 340 KDDVFVNGV--NQKAPDDKRESFAGYVMNKTTF-GNWQFTQGYRREKVKYEYSSKVYDpmtwqlkeikPKSADYSNNdSFEFGV 420 333333333..56777779****9*********.**********9988888777777777777888874444444444****** PP EEECCTTEEEEEEEEEEEE--THHCHCCCEEEGGGTEEECCSTT---EEEEEEEEEEE..EEETTEEEEEEEEEEEEEEEEEEE CS TonB_dep_Rec 217 ryrptddlslrasysrgfraPslgelyaggsngaggalsngnpdLkPEtsknyelGlk..ldddrlnltlayFridiknfittd 298 +y +d +++++y+r+ r P++++ a+ +g k ++ + +e+Gl+ + + + ++ + F+ + kn i++d 23726_FusoPortal_Gene_597 421 NYLYSDTGNVFFNYARALRTPTIQDAGAW----------YGP--VKTQKNDIFEIGLRdaYKNTS--ISTSIFYVNSKNEIYYD 490 ****************************4..........576..8*************9665444..***************** PP SSSCTSTTCCTEEEEEE.EEEEEEEEEEEEEEEECTTEE.EEEEEEEEEEEEEEECCCCBCTTTSCCCEESEEEEEEEEEESSS CS TonB_dep_Rec 299 pdgpslgglgdfqyqNv.gkarikGvElelsyrltpgWg.lglsgggytYydagddsdaddggnplpnvPkhtanlgltYkfpa 380 ++p + + +qN gk+r G+ l+l+++ ++ l+l+ ++Y +++ ++ ++g+ +++v+++tan+g+tY +++ 23726_FusoPortal_Gene_597 491 KTNPFS-----SNNQNFdGKVRRIGAQLSLAHY----FDkLTLRER-VSYIVPKVT-SGIYDGKEFAGVSRWTANVGATY-NIT 562 ***998.....788899****************....997666666.777666444.88999******************.889 PP CCEEEEECEEEEEEEEECEEEEEECE..CCEEEE--EEEEEEEEEEECSSSEEEEEEEESTT-TSSCEEESTT........TSC CS TonB_dep_Rec 381 ggWgaglllagarYvavggrrygdaa..ntfrrvpgYtlldlslgYklpknltlslgvnNLfdkryyassgyy........glg 454 g++a++ +++Y ++ +y++++ n f++ Y+++d+sl+Y + +++ l+ gv NLfdk+y++ ++ 23726_FusoPortal_Gene_597 563 KGLTANI---DGYY---QSNAYAEDDfdNYFSKGNNYVTVDASLSYAFENGIELYTGVSNLFDKKYANAVTSTrstfgagpRKV 640 99**888...9999...443555555546666888***********999**************877777777789999**99** PP CEEB--SEEEEEEEEE CS TonB_dep_Rec 455 gypgpgRtfrlglryk 470 +yp+ gR++++g++y+ 23726_FusoPortal_Gene_597 641 YYPANGRSVYAGIKYT 656 **************97 PP >> Plug TonB-dependent Receptor Plug Domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 84.3 0.8 1.3e-27 7.2e-24 2 108 .] 42 147 .. 41 147 .. 0.97 2 ? -3.4 0.0 2.5 1.4e+04 60 73 .. 342 355 .. 337 366 .. 0.71 Alignments for each domain: == domain 1 score: 84.3 bits; conditional E-value: 1.3e-27 HGCCSSEEEEEEHHHHHHTTTSSHHHHHTTSTTEECBCCECSTTSEEEEETTB-G...GGSEEEETTEEESCTCCSCTCCTSSB CS Plug 2 kkdtpgsvsvvtaediekqgatsladalqgvpgvnvsssgggggassisiRGlgs...nrvlvlvDGvpianggsgagdlssid 82 +++t+++++v++a++i+++ga+++adal+gvpgv +++ g ++ + +RG g+ ++ +l+DG+p++ g ++++l++++ 23726_FusoPortal_Gene_597 42 VRETAKNMTVINAKEIKEKGAKTIADALRGVPGVV--VRQMDGASPMLDLRGSGAtsqFNTVILLDGIPVS--GLAGFNLNTVP 121 6899*******************************..7***************9999999**********6..9********** PP GGGEEEEEEEESCHHHTCSSTCSSEE CS Plug 83 pedierievlkGpasalyGsdaigGv 108 ++ ie+ievl G+ +++yG++aigGv 23726_FusoPortal_Gene_597 122 VDEIEKIEVLQGAGAVMYGDGAIGGV 147 *************************8 PP == domain 2 score: -3.4 bits; conditional E-value: 2.5 EEEETTEEESCTCC CS Plug 60 lvlvDGvpianggs 73 v+v+Gv++ +++ 23726_FusoPortal_Gene_597 342 DVFVNGVNQKAPDD 355 57899998876444 PP >> AAA_2 AAA domain (Cdc48 subfamily) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.1 0.0 5.4e-05 0.3 30 82 .. 44 130 .. 19 152 .. 0.84 Alignments for each domain: == domain 1 score: 11.1 bits; conditional E-value: 5.4e-05 ....-EEEEEGGG-------THHHHHHHHHTTT-HHHHTTS...................................EEEEETGG CS AAA_2 30 kkekkliradaselleeeevehlvsrligappgyvgaeegG...................................ivllDEid 78 ++ k++++++a+e+ e + ++++ + l g p v++ +g v +DEi+ 23726_FusoPortal_Gene_597 44 ETAKNMTVINAKEIKE-KGAKTIADALRGVPGVVVRQMDGAspmldlrgsgatsqfntvilldgipvsglagfnlnTVPVDEIE 126 6789************.778888888888888888777777*****************99999999999999999999999999 PP G--- CS AAA_2 79 Kiap 82 Ki 23726_FusoPortal_Gene_597 127 KIEV 130 9976 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (658 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 637 (0.0381163); expected 334.2 (0.02) Passed bias filter: 337 (0.0201652); expected 334.2 (0.02) Passed Vit filter: 28 (0.00167544); expected 16.7 (0.001) Passed Fwd filter: 4 (0.000239349); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 3 [number of targets reported over threshold] # CPU time: 0.30u 0.22s 00:00:00.52 Elapsed: 00:00:00.32 # Mc/sec: 5977.58 // [ok]