# 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_1866.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1866 [L=751] Description: # 1994587 # 1996839 # 1 # ID=1_1866;partial=00;start_type=ATG;rbs_motif=GGA/GAG/AGG;rbs_spacer=5-10bp;gc_cont=0.296 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 4.9e-65 218.2 0.2 1.1e-64 217.1 0.2 1.6 1 Transgly Transglycosylase 1.3e-22 80.3 1.3 1.6e-21 76.7 1.3 2.3 1 Transpeptidase Penicillin binding protein transpeptidase dom ------ inclusion threshold ------ 0.039 13.8 0.2 0.13 12.2 0.2 1.8 1 BOFC_N Bypass of Forespore C, N terminal Domain annotation for each model (and alignments): >> Transgly Transglycosylase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 217.1 0.2 2e-68 1.1e-64 3 177 .. 58 231 .. 56 232 .. 0.98 Alignments for each domain: == domain 1 score: 217.1 bits; conditional E-value: 2e-68 -EEEEE-----XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX CS Transgly 3 ellaeleeeenrelvkleeippalkkavvaiEDrrFyehkGvdlkgiarAvlknlksgkssqGgSTitqQlvknlflseertl 85 +++++l + e+re++kle+i+p++k+a+++iED++Fy+h+G++++gi+rA+++++ +g+ +qGgS+itqQl+kn+fl+ ert+ 23726_FusoPortal_Gene_1866 58 NVVDTL-SVESREIAKLENISPYVKEAFLSIEDKQFYSHHGLNFRGITRAIVTTFLKGRPTQGGSSITQQLAKNAFLTPERTF 139 677778.699************************************************************************* PP XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX CS Transgly 86 erKlkeavlAlklekklsKdeILeaYlNsvylgrevyGieaAarayFgkdakdLtlaeaalLagllkaPsrynplenpeeakk 168 +rK+kea+l++++e++++KdeILe YlN++y+g++ yGi++Aa++yF+kd+kdL++aeaalLag++++P++y+p++n+e+a 23726_FusoPortal_Gene_1866 140 SRKVKEAILTYQIERTYTKDEILERYLNEIYFGSGSYGIRNAAEQYFKKDVKDLNIAEAALLAGIPNRPTKYDPNRNLENALY 222 *********************************************************************************** PP XXXXXXXXX CS Transgly 169 rrkvvlerm 177 r+k++l++m 23726_FusoPortal_Gene_1866 223 RQKIILKEM 231 ********9 PP >> Transpeptidase Penicillin binding protein transpeptidase domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 76.7 1.3 3e-25 1.6e-21 1 302 [. 335 602 .. 335 606 .. 0.78 Alignments for each domain: == domain 1 score: 76.7 bits; conditional E-value: 3e-25 EEEEEEETTT-BEEEEEEESE-TTTHHCCTCGCCSTHHHH--BB-GGGGHHHHHHHHHHTTSSTTSEEEEETCCEETTSS.EC CS Transpeptidase 1 gaavvldaktgevlamaskpsydpnakvgkneplanravqsqyePGStfKvvtaaaaleagvikpdetlddsggkiqggqsik 83 ga++ ld+ tg ++++++ ++++ ++ ra+ + + GS fK+++ + al++g ++p++++ + ++ 23726_FusoPortal_Gene_1866 335 GAMITLDPFTGGIVSIVGGKNFKAGNF--------DRATMARRQLGSSFKPFVYLEALQNG-FDPYSVVVNDFVA------FG 402 677889***************988887........469999********************.**99977666554......55 PP TCCCC-...SSSSEEEHHHHHHTT-HHHHHHHHHHH-HHHHHHHHHHTTTTS-SSSSSTT-B--B--SSHHHHHHHHHTTT.S CS Transpeptidase 84 dweqdn...kgsegitlrealekSsntyfvklaqklgadklrkylkkfglGnktgidlpgesaGvptrkwlegatisfgiG.l 162 +w +n + ++ tl ++l S n+ vkl +++ + +++++++ +++l +e++ +t + G + 23726_FusoPortal_Gene_1866 403 KWAPKNfdgRYTYNSTLVNSLNLSLNVPAVKLLDAVTVEAFKEAIGD-------NVKLTSEVKD-LTAA----------LGsV 467 566665543444555***************************99864.......3566667777.3333..........4448 PP BB-HHHHHHHHHHHHTTTEEE--BSEEEEEECTTEEEESS-EEEEE-S-HHHHHHHHHHHHHHHHCTTTTGCG....TCTTS- CS Transpeptidase 163 tvtplqlaqayaaiangGklvkPhlvkaivdpdgvrkvkpkelnklpisketaekvkeglkgvvsrgsgtata....avegyk 241 tp+ +a ++ ++ngG++vkP ++ +i+d+++ + ++ +++++ +++ ++ +lk+vvs+g+++ ++ + + ++ 23726_FusoPortal_Gene_1866 468 DSTPVNVAANFSIFVNGGYIVKPNIIREIRDNQDILIYVAEIEKTKAFDSVDVSVITAMLKTVVSNGTASKARvvdkTGKPIQ 550 99******************************995555555567789********************5555554444444455 PP -EEEEEEEEECSSSTTSEEEEEEEEEEEEEESTTS-SEEEEEEECTCS....SSSSCCHHHHHHH CS Transpeptidase 242 vagKTGTAevagkkgngkkegakiswfvgyapadnpkyavavlieqpk....eaaggkaaapiak 302 gKTGT++ +wfvg++p + a ++ +++k +a+gg+a ap+ + 23726_FusoPortal_Gene_1866 551 QGGKTGTTNEH-----------RTAWFVGITPE--YVTACYIGRDDNKpmygKATGGSAVAPMWA 602 55555554444...........457******97..666666777777788887888888777665 PP >> BOFC_N Bypass of Forespore C, N terminal # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.2 0.2 2.3e-05 0.13 2 23 .. 225 246 .. 225 253 .. 0.91 Alignments for each domain: == domain 1 score: 12.2 bits; conditional E-value: 2.3e-05 EEEEEEE-SSS-EEEEEEEEEE CS BOFC_N 2 tviLervYlDGevseekveEti 23 ++iL+++Y DG++ +e++ E++ 23726_FusoPortal_Gene_1866 225 KIILKEMYTDGRITKEQYDEAL 246 69****************9975 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (751 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 722 (0.0432025); expected 334.2 (0.02) Passed bias filter: 422 (0.0252513); expected 334.2 (0.02) Passed Vit filter: 46 (0.00275251); expected 16.7 (0.001) Passed Fwd filter: 3 (0.000179512); 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.34u 0.24s 00:00:00.58 Elapsed: 00:00:00.33 # Mc/sec: 6615.70 // [ok]