# 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_1495.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1495 [L=661] Description: # 1592617 # 1594599 # -1 # ID=1_1495;partial=00;start_type=ATG;rbs_motif=GGAG/GAGG;rbs_spacer=5-10bp;gc_cont=0.287 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1.3e-39 137.3 47.6 3e-39 136.1 47.4 1.6 1 TonB_dep_Rec TonB dependent receptor 1e-16 61.3 2.2 1.9e-16 60.5 2.2 1.5 1 Plug TonB-dependent Receptor Plug Domain 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 ! 136.1 47.4 3.6e-43 3e-39 23 470 .] 182 659 .. 160 659 .. 0.69 Alignments for each domain: == domain 1 score: 136.1 bits; conditional E-value: 3.6e-43 EEETTTCEEEEEEECEEEEEEEEECCCEECCCEEEEEEEEEEEEECEEEECCEEEEEEEEEEEEEEEEEEEEEEEE....... CS TonB_dep_Rec 23 rytdggdlgsalsrgylgtnplyygdnysetdnkrlrlgydydlddnlswklslraglnYsstdgdsdynnsglgg....... 98 + td ++ ++sr+++++ +++++++ ++++++++ y++d+ +l++ +Y++++ + + + + 23726_FusoPortal_Gene_1495 182 KATDKLTFKVDYSRDRYRSERDSDENGKVVSRSQEVSIDIKYKFDNA-----DLTV--KYTRNEKHRADGGDLEEAdyyknrk 257 55556666666666666666655555566666777777777776444.....3332..2332222111111111112222244 PP ...................................CEEEEEEE.....EEEE...EEEEEEEEEEEEEECEEE.CEEEEEEEE CS TonB_dep_Rec 99 ...................................grsynrdn.....dyyg...lgnvlaaeldlngdlwfG.snladgldy 137 +++ + + + ++++++a+l+l +++ + + + g+dy 23726_FusoPortal_Gene_1495 258 mvtwvardftrsndwyvnyrqnigentelltyidlYDTR---EnddisK--TldrDYSRKTAKLQLKHKY-LNsHYFIIGTDY 334 555552222222222222222222222222222111111...1222321..0224445555555555553.555666666999 PP EEEEEEECCCEESCCTCECEEEEEEEEEEEEEEEEEEEEETTEEEEEEEEEEEETTE.....SSSS............EEEEE CS TonB_dep_Rec 138 rklsddgydaydplnpsgdasdsvdrdsyglyaqdnikltdrltltlGlRydhysvf.....gkdlrag.........afsrs 206 k + ++n+ ++ +++ + yg++a +++k+ +++t+ GlRy +++ + 23726_FusoPortal_Gene_1495 335 MKEKL----KGLGYNGE-YTGKNTEKTDYGIFAMNELKF-GKFTFAQGLRYNKAEYDfywrnK-----YpvprnirgeHGEQE 406 95444....44666666.888899***************.**********9877777333321.....133334555545555 PP EEEEEEEEEEEEECCTTEEEEEEEEEEEE--THHCHCCCEEEGGGTEEECCSTT---EEEEEEEEEEE.EEETTEEEEEEEEE CS TonB_dep_Rec 207 gsafsPsvglryrptddlslrasysrgfraPslgelyaggsngaggalsngnpdLkPEtsknyelGlk.ldddrlnltlayFr 288 ++++ ++ lry+ +d ++++++sr fr P e+y+ +l + Lk++t +++e+G k + + + ++l+ F+ 23726_FusoPortal_Gene_1495 407 YKNYAANLELRYDYSDTGMVYGKWSRDFRTPIAREMYY--------TLEGSK--LKAQTQNTFEIGAKdYIAGTY-ISLSTFY 478 699999999*****************************........455544..**************8776666.******* PP EEEEEEEEEESSSCTSTTCCT...EEEEEEEEEEEEEEEEEEEEEECTTEEEEEEEEEEEEEEEEECCCCBCTTTSCCCEESE CS TonB_dep_Rec 289 idiknfittdpdgpslgglgd...fqyqNvgkarikGvElelsyrltpgWglglsgggytYydagddsdaddggnplpnvPkh 368 +++++ i+++++++ ++ + f+y+N g++r G+El l+ + ++++++ s + y + d+++ +++++ +p vP + 23726_FusoPortal_Gene_1495 479 KKTNGEIYYQGTPNKESTRPGavvFPYYNMGDTRRLGIEL-LTEQYVKNFTFTESIS-YLNHKIVDSDFESRKNKEIPMVPNW 559 *************99999988*******************.6556667788877777.4424444448889************ PP EEEEEEEEESSSCCEEEEECEEEEEEEEECEEEEEECECCEE..EE--EEEEEEEEEEECSSSEEEEEEEESTT-TSSCEEES CS TonB_dep_Rec 369 tanlgltYkfpaggWgaglllagarYvavggrrygdaantfr..rvpgYtlldlslgYklpknltlslgvnNLfdkryyassg 449 ++ +g+ Ykf +++++ ++a+ v++g+++++d+ ++r + Y+++ ls++Yk+ +++ +++ vnNLfdk+y ++ g 23726_FusoPortal_Gene_1495 560 KLAFGVGYKF-NNKLN-----VNADIVYYGKFYDSDDPENVRpkDRGNYATVSLSANYKFENGFAINARVNNLFDKKYEDYVG 636 ********55.56666.....566666669977777777777899************************************** PP TT..TSCCEEB--SEEEEEEEEE CS TonB_dep_Rec 450 yy..glggypgpgRtfrlglryk 470 y+ + ++ p++g ++++g++y+ 23726_FusoPortal_Gene_1495 637 YWdgTRQYSPAAGKYYSIGASYT 659 **9888888888*********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 ! 60.5 2.2 2.3e-20 1.9e-16 2 108 .] 40 142 .. 39 142 .. 0.95 Alignments for each domain: == domain 1 score: 60.5 bits; conditional E-value: 2.3e-20 HGCCSSEEEEEEHHHHHHTTTSSHHHHHTTSTTEECBCCECSTTSEEEEETTB-G....GGSEEEETTEEESCTCCSCTCCTS CS Plug 2 kkdtpgsvsvvtaediekqgatsladalqgvpgvnvsssgggggassisiRGlgs....nrvlvlvDGvpianggsgagdlss 80 +++t+++v++vta+die++ga++l dal+ vpg+ v ++ g+ + i G +s ++ ++ DGv+i +a + + 23726_FusoPortal_Gene_1495 40 VRNTASNVTIVTAKDIEEKGAQNLVDALRMVPGIMV---KNYYGNIAFDIGGYSSvhaeRNSIITYDGVRI-----SAKEATN 114 689********************************9...7778888999****99***999999******9.....6****** PP SBGGGEEEEEEEESCHHHTCSSTCSSEE CS Plug 81 idpedierievlkGpasalyGsdaigGv 108 i++++ier+ev+ + lyG++a gGv 23726_FusoPortal_Gene_1495 115 IPISAIERVEVIPNGGGILYGDGANGGV 142 ***************************8 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (661 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 731 (0.043741); expected 334.2 (0.02) Passed bias filter: 359 (0.0214816); expected 334.2 (0.02) Passed Vit filter: 29 (0.00173528); expected 16.7 (0.001) Passed Fwd filter: 2 (0.000119674); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.39u 0.26s 00:00:00.65 Elapsed: 00:00:00.40 # Mc/sec: 4803.87 // [ok]