# 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_1496.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1496 [L=287] Description: # 1594802 # 1595662 # 1 # ID=1_1496;partial=00;start_type=ATG;rbs_motif=AGGAG;rbs_spacer=5-10bp;gc_cont=0.272 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.7e-28 99.7 0.4 2.2e-28 99.3 0.4 1.1 1 Peripla_BP_2 Periplasmic binding protein ------ inclusion threshold ------ 0.034 14.1 0.6 0.12 12.4 0.4 2.0 2 Creatinase_N_2 Creatinase/Prolidase N-terminal domain Domain annotation for each model (and alignments): >> Peripla_BP_2 Periplasmic binding protein # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 99.3 0.4 2.7e-32 2.2e-28 3 234 .. 25 258 .. 24 260 .. 0.92 Alignments for each domain: == domain 1 score: 99.3 bits; conditional E-value: 2.7e-32 E-SHHHHHHHHHTT-GGGEEEB-GGG--GGGGGGT..TS-B-BESSSB-HHHHHHH--SEEEE-GGG......GGGHHHHHHH CS Peripla_BP_2 3 slspsvteillalGatdkvvgvdtttsdpldadkleeaivevgadgeidvEklaalkPDlvilskgl...gltdedeellksg 82 s + +teill +Ga++++vg++ ++ l++ k++ +++ v + ++++ E+ a+kPD + + + + ++el+k+g 23726_FusoPortal_Gene_1496 25 SAAHFSTEILLSIGAEKQMVGTAYPDNEILPSLKEKYDKIPVLSMKNPTKEQFYAVKPDFLTGWDSTvldKNLGPIKELEKNG 107 667789**************************98877************************987777789999********** PP S-EEE----TTS..HHHHHHHHHHHHHTT-HHHHHHHHHHHHHHHHHHHHT--TT.-EEEEEEEETTE..EEEE-TTSTTTHH CS Peripla_BP_2 83 iptvivessntg..eslkeqiralgellgledeaeklvaeidsklaaaksaidsekpkpvlvfggadeggyvvagsnsyigdl 163 ++++i+ s + + e i ++g++++le++a+k+v +++s laa+++ +++ + +v+++ ++d+ +++v g +s ++ + 23726_FusoPortal_Gene_1496 108 VQVYIMKSLHSSdiNLVFEDILTYGKIFNLENNAKKVVGKMKSDLAAVQKQLPKN-KVKVFTYDSGDK-APFVVGGDSIGNTI 188 *******9877677*****************************************.667777666666.77777778899999 PP HHTS--EES-TT--STT-EEE-HHHHHHH--SEEEEEEHHHHT-S--HHHHT-SHHHHTSHHHHTT-EEEE CS Peripla_BP_2 164 ldalgieniaaeskkseyeeisaEqilasdpdvifvsgraeadnktgleflkenpllaelpAvkngrVyal 234 + +g++ni+++ k+++++ ++E +l ++pd+i++ + + ++ + +flke ++l+Avkn++++ + 23726_FusoPortal_Gene_1496 189 ITLAGGDNIFKDI-KKAWADGNWEKVLVENPDIIVIIDYGDQSAESKIKFLKEKSPIKDLKAVKNNKFVVI 258 999**********.***********************999999************************9875 PP >> Creatinase_N_2 Creatinase/Prolidase N-terminal domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -0.4 0.0 0.12 1e+03 127 127 .. 147 147 .. 92 175 .. 0.51 2 ? 12.4 0.4 1.4e-05 0.12 101 159 .. 190 253 .. 139 254 .. 0.68 Alignments for each domain: == domain 1 score: -0.4 bits; conditional E-value: 0.12 G CS Creatinase_N_2 127 p 127 + 23726_FusoPortal_Gene_1496 147 K 147 1 PP == domain 2 score: 12.4 bits; conditional E-value: 1.4e-05 EEEEGGGHHHHHHHCCH.CCT-BHCCCGCCCEHHH.HHHC-.....GGGEEE-S-CHHHHHHC--H CS Creatinase_N_2 101 eikpyeevvsdlkalas.askkkvlldpsktsyal.asalk.....ekkvvekespiallKavKNe 159 ++ +++++d+k+ + ++ +kvl+++ + + + + + + k ++++spi++lKavKN+ 23726_FusoPortal_Gene_1496 190 TLAGGDNIFKDIKKAWAdGNWEKVLVENPDIIVIIdYG--DqsaesKIKFLKEKSPIKDLKAVKNN 253 45555666666665544444677787766655555222..12343455568889***********7 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (287 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 578 (0.0345859); expected 334.2 (0.02) Passed bias filter: 386 (0.0230972); expected 334.2 (0.02) Passed Vit filter: 28 (0.00167544); 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.27u 0.24s 00:00:00.51 Elapsed: 00:00:00.33 # Mc/sec: 2528.24 // [ok]