# 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_1620.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1620 [L=721] Description: # 1739548 # 1741710 # -1 # ID=1_1620;partial=00;start_type=ATG;rbs_motif=GGAGG;rbs_spacer=5-10bp;gc_cont=0.271 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.6e-46 160.1 27.9 1.6e-46 160.1 27.9 2.1 2 TonB_dep_Rec TonB dependent receptor 5.3e-21 75.1 1.8 1.2e-20 73.9 1.8 1.6 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 ? -2.9 9.7 0.51 4.2e+03 50 113 .. 155 215 .. 130 227 .. 0.44 2 ! 160.1 27.9 2e-50 1.6e-46 12 470 .] 222 719 .. 215 719 .. 0.78 Alignments for each domain: == domain 1 score: -2.9 bits; conditional E-value: 0.51 EECCCEEEEEEEEEEEEECEEEECCEEEEEEEEEEEEEEEEEEEEEEEECEEEEEEE.EEEEEEE CS TonB_dep_Rec 50 ysetdnkrlrlgydydlddnlswklslraglnYsstdgdsdynnsglgggrsynrdn.dyyglgn 113 y +++n r +gy++ n+++ +++ag++ +++d d +y + g+r+y+ +n dy++ + 23726_FusoPortal_Gene_1620 155 YKGNNNVRGGVGYQLASFRNNKF--DVSAGTSVGDFDFDINYSKNRKYGYRDYDFTNsDYFS--G 215 33333444455555553233334..4444444443333333333333334444444433333..1 PP == domain 2 score: 160.1 bits; conditional E-value: 2e-50 EEEEEEEEEEEEEETTTCEEEEEEECEEEEEEEEECCC.....EECCCEEEEEEEEEEEEECEEEECCEEEEEEEEEEEEEEE CS TonB_dep_Rec 12 dnllanngdylrytdggdlgsalsrgylgtnplyygdn.....ysetdnkrlrlgydydlddnlswklslraglnYsstdgds 89 + +++ ++y++y+d +++++ l+++ l+ n++++g + +++ ++++++l+y+++ +d+ ++l+ + Y++td s 23726_FusoPortal_Gene_1620 222 NKTSNIAFKYSGYRDKYTYPNFLDQKELDENRRQSGIDkeakeNNRIKKDEFTLTYNTKIGDK----NDLNILGFYQKTDIPS 300 567788899****************77777773333333444566666999999999998777....3888888899888888 PP EEEEEEEEE........................................CEEEEEEEEEEEEEEEEEEEEEEEEECE...... CS TonB_dep_Rec 90 dynnsglgg........................................grsynrdndyyglgnvlaaeldlngdlw...... 126 + + +++ g + +++ + 23726_FusoPortal_Gene_1620 301 ESIEDYTSEykgmlagqtaklrkalrdprlsararaamqnrlnallaelG-----------------STNNVD---Lkkfsqf 363 88888777777788777777777776665554443333333333322221.................111111...0122333 PP .................EECEEEEEEEEEEEEE...EECCCEESCCTCECEEEEEEEEEEEEEEEEEEEEETTEEEEEEEEEE CS TonB_dep_Rec 127 .................fGsnladgldyrklsd...dgydaydplnpsgdasdsvdrdsyglyaqdnikltdrltltlGlRyd 189 +Gsn gl y++++ +d + + d+ + + +++ ++a +++k+ +r+ + GlR++ 23726_FusoPortal_Gene_1620 364 kdtkkaikikdkftydsVGSNVIVGLGYTNNDMlrvSKTDLVGKRTM-ADTKLDLSKKTFEVFALNTFKV-NRFEFIQGLRFE 444 45555555555556566668888889999866644333333333333.5999999999999*********.************ PP EETTESSSS......EEEEEEEEEEEEEEEEEECCTTEEEEEEEEEEEE--THHCHCCCEE.EGGGTEEECCSTT---EEEEE CS TonB_dep_Rec 190 hysvfgkdlrag...afsrsgsafsPsvglryrptddlslrasysrgfraPslgelyaggs.ngaggalsngnpdLkPEtskn 268 +++++g+ +++ + ++s +++ s++++y +d + +++y r+f P g+l + + +++ +n Lk E ++ 23726_FusoPortal_Gene_1620 445 NSKYDGTRKNNDvalDIKKSKDDWAGSLAVNYLYSDTGNAYVKYERAFTSPAPGQLVDKIEiAPRVYTYKVNN--LKSESTNL 525 **********************************************************666766677777777..******** PP EEEEEE.EEETTEEEEEEEEEEEEEEEEEEESSSCTSTTCCTEEEEEEEEEEEEEEEEEEEEEECTTEE.EEEEEEEEEEEEE CS TonB_dep_Rec 269 yelGlk.ldddrlnltlayFridiknfittdpdgpslgglgdfqyqNvgkarikGvElelsyrltpgWg.lglsgggytYyda 349 +e+G++ + +++l l++++F++++k+ i+t +d+++ g f ++N+gk++ +G++l++ + +++ + y++ ++ 23726_FusoPortal_Gene_1620 526 FEVGWNdYLFGSL-LSADVFYSETKDEIATIFDVGGHGFG--FKNTNIGKTKRYGFDLSAEQK----FEkFTFKEA-YSFIET 600 *******9999*9.********************999988..****************99888....776666777.999555 PP ..EECCCCBCTTTSCCCEESEEEEEEEEEESSSCCEEEEECEEEEEEEEECEEEEEECECCEEEE--EEEEEEEEEEECSSSE CS TonB_dep_Rec 350 ..gddsdaddggnplpnvPkhtanlgltYkfpaggWgaglllagarYvavggrrygdaantfrrvpgYtlldlslgYklpknl 430 d +++++g+++++vPkh+++++++Y +++++++ +ga+Y +++ ++ d+an ++++++ +++l+++Yk+++ l 23726_FusoPortal_Gene_1620 601 kiLKDNSNSFEGKHIADVPKHKLVFSVDY-DITSKFT-----VGADYEYRAA-AFIDNANKYGKDKAKSVFNLRANYKITNSL 676 333335779********************.6788999.....8999999999.7777777777******************** PP EEEEEEESTT-TSSCEEESTT...TSCCEEB--SEEEEEEEEE CS TonB_dep_Rec 431 tlslgvnNLfdkryyassgyy...glggypgpgRtfrlglryk 470 ++++g+nN+f +yy+s + + ++ p p ++++g++yk 23726_FusoPortal_Gene_1620 677 NVYAGINNIFGAKYYNSVRGNsrgEKFYDPTPKINYYAGFKYK 719 ***********65555555556977777777789******997 PP >> Plug TonB-dependent Receptor Plug Domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 73.9 1.8 1.5e-24 1.2e-20 2 108 .] 39 147 .. 38 147 .. 0.95 Alignments for each domain: == domain 1 score: 73.9 bits; conditional E-value: 1.5e-24 HGCCSSEEEEEEHHHHHHTTTSSHHHHHTTSTTEECBCCECSTTSEEEEETTB-G....GGSEEEETTEEESC..TCCSCTCC CS Plug 2 kkdtpgsvsvvtaediekqgatsladalqgvpgvnvsssgggggassisiRGlgs....nrvlvlvDGvpian..ggsgagdl 78 ++++++s +++t++die++g ts++++l +vpgvnv + g ++ + +RG g v++lvDGvp++ +++++++ 23726_FusoPortal_Gene_1620 39 LRNSTTSPYIITSKDIETKGYTSVSEILDSVPGVNV----QEGLRPAVDVRGQGFqkakATVQLLVDGVPANMldTSHMNVPI 117 5789999999*************************8....8999999999999777899***********999********** PP TSSBGGGEEEEEEEESCHHHTCSSTCSSEE CS Plug 79 ssidpedierievlkGpasalyGsdaigGv 108 ++++++ ieriev+ G+ ++lyGs+++gGv 23726_FusoPortal_Gene_1620 118 DVVNINEIERIEVIPGGGAVLYGSGTSGGV 147 *****************************8 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (721 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 517 (0.0309359); expected 334.2 (0.02) Passed bias filter: 264 (0.015797); expected 334.2 (0.02) Passed Vit filter: 22 (0.00131642); 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.31u 0.23s 00:00:00.54 Elapsed: 00:00:00.33 # Mc/sec: 6351.42 // [ok]