# 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_1912.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1912 [L=365] Description: # 2052572 # 2053666 # -1 # ID=1_1912;partial=00;start_type=ATG;rbs_motif=AGGA;rbs_spacer=5-10bp;gc_cont=0.320 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 3.6e-29 100.7 1.6 8.6e-29 99.5 1.6 1.6 1 Peptidase_M23 Peptidase family M23 1e-16 60.6 0.7 1e-07 31.8 0.0 2.4 2 LysM LysM domain ------ inclusion threshold ------ 0.32 11.1 2.7 3.3 7.9 0.5 2.9 3 DUF5055 Domain of unknown function (DUF5055) 1.5 7.9 7.1 7.1 5.7 0.9 2.9 3 Ribosomal_S2 Ribosomal protein S2 Domain annotation for each model (and alignments): >> Peptidase_M23 Peptidase family M23 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 99.5 1.6 2.1e-32 8.6e-29 3 96 .] 267 359 .. 266 359 .. 0.98 Alignments for each domain: == domain 1 score: 99.5 bits; conditional E-value: 2.1e-32 BSSEEEEETTT-EEEESSSEEEEEEEEETTTEEEEEEEECTCEEEEEEEESEESS-TTEEE-TTCEEEEEBSTSSCSSSEEEE CS Peptidase_M23 3 hkGidiaapkgtpVlAaadGtVvfagksegyGnvviidhgdgystlYahlskisvkegqkVkaGqviGtvgstgrstgphLhf 85 h+G+d+ a+ p +A ++G+V+fag+++gyG+++ii+h++gy+t+Yahls is ++g++V++G++iG++g++gr+tg hLhf 23726_FusoPortal_Gene_1912 267 HTGVDLVAK-YVPLRASKAGVVTFAGNMSGYGKIIIIKHDNGYETRYAHLSVISTNVGEHVNKGDLIGKTGNSGRTTGAHLHF 348 99****999.699**************99****************************************************** PP EEEETTEEE-C CS Peptidase_M23 86 eirkngkavdp 96 eir+ng + +p 23726_FusoPortal_Gene_1912 349 EIRHNGVPKNP 359 ******99988 PP >> LysM LysM domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 31.8 0.0 2.4e-11 1e-07 1 44 [] 85 128 .. 85 128 .. 0.97 2 ! 27.1 0.1 7.1e-10 3e-06 1 42 [. 135 175 .. 135 177 .. 0.96 Alignments for each domain: == domain 1 score: 31.8 bits; conditional E-value: 2.4e-11 EEE-TT--HHHHHHHCTCHHHHHHHHCTCSSCGG--TTEEEEEE CS LysM 1 YtVkkGDTLwsIAkrygitveklleaNnlsspdkifpGqkLkiP 44 Y+V+k DT+ sIAk++g+ e+++ N+ ++k+++G +L+ P 23726_FusoPortal_Gene_1912 85 YKVQKKDTVQSIAKKFGVKPETIMINNQTAMDNKLKVGEVLTFP 128 9****************************9999********977 PP == domain 2 score: 27.1 bits; conditional E-value: 7.1e-10 EEE-TT--HHHHHHHCTCHHHHHHHHCTCSSCGG--TTEEEE CS LysM 1 YtVkkGDTLwsIAkrygitveklleaNnlsspdkifpGqkLk 42 Y+ +k++ L ++Ak+yg+ v ++ ++Nn+++ + +++G++L 23726_FusoPortal_Gene_1912 135 YKLQKNEMLAKVAKKYGVKVVDIVDYNNINPKK-LKAGTTLF 175 899******************************.******96 PP >> DUF5055 Domain of unknown function (DUF5055) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 7.9 0.5 0.00078 3.3 14 78 .. 42 107 .. 32 124 .. 0.85 2 ? 1.0 0.1 0.11 4.5e+02 64 94 .. 130 159 .. 127 169 .. 0.81 3 ? -0.5 0.0 0.32 1.3e+03 57 81 .. 311 335 .. 307 339 .. 0.83 Alignments for each domain: == domain 1 score: 7.9 bits; conditional E-value: 0.00078 DUF5055 14 ytleftrrtvkemekk.gfvaadieekpmttlpalfagaflahhrfekkdvineiyskltkkeeli 78 le+t + +ek+ fv +++ ek p + a + ++ +kkd ++ i k+ k e i 23726_FusoPortal_Gene_1912 42 GGLELTTSNFTTFEKEyNFVKEEVVEKKEEKPPVQQKRAEKITYKVQKKDTVQSIAKKFGVKPETI 107 469999999999999757*************************************99999988765 PP == domain 2 score: 1.0 bits; conditional E-value: 0.11 DUF5055 64 ineiyskltkkeeligklaemynepilalve 94 i+ +y kl +k+e++ k+a+ y ++ +v+ 23726_FusoPortal_Gene_1912 130 IDGLYYKL-QKNEMLAKVAKKYGVKVVDIVD 159 67778887.578999*******999988886 PP == domain 3 score: -0.5 bits; conditional E-value: 0.32 DUF5055 57 rfekkdvineiyskltkkeeligkl 81 r+ + vi + +k +ligk 23726_FusoPortal_Gene_1912 311 RYAHLSVISTNVGEHVNKGDLIGKT 335 6777788888889999*******96 PP >> Ribosomal_S2 Ribosomal protein S2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 5.7 0.9 0.0017 7.1 93 159 .. 48 114 .. 24 130 .. 0.72 2 ? 4.6 0.0 0.0036 15 157 185 .. 137 166 .. 133 170 .. 0.84 3 ? -1.2 0.3 0.22 9.2e+02 93 126 .. 181 206 .. 174 216 .. 0.44 Alignments for each domain: == domain 1 score: 5.7 bits; conditional E-value: 0.0017 Ribosomal_S2 93 LTNfktikksikklkeleeleedgeldltkkealklkkekekleknlggikkmkklPdlvvvldpkk 159 Nf+t +k+ + +ke ++++++ +++k+a k++ + +k + + kk +P+ +++ +++ 23726_FusoPortal_Gene_1912 48 TSNFTTFEKEYNFVKEEVVEKKEEKPPVQQKRAEKITYKVQKKDTVQSIAKKFGVKPETIMINNQTA 114 56999999995555533333233333366777777888888866666667888999****9998765 PP == domain 2 score: 4.6 bits; conditional E-value: 0.0036 Ribosomal_S2 157 pkkekiavkEakklgipviaivDtn.vdpd 185 +k++++ k akk+g++v+ ivD n ++p+ 23726_FusoPortal_Gene_1912 137 LQKNEMLAKVAKKYGVKVVDIVDYNnINPK 166 68999999***************8625665 PP == domain 3 score: -1.2 bits; conditional E-value: 0.22 Ribosomal_S2 93 LTNfktikksikklkeleeleedgeldltkkeal 126 L +k++++++ ++ +e++ kke++ 23726_FusoPortal_Gene_1912 181 LKKYKEVEQRLIAAQQAKEEQ--------KKEKA 206 444455555533333333332........22222 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (365 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 1439 (0.0861058); expected 334.2 (0.02) Passed bias filter: 559 (0.033449); expected 334.2 (0.02) Passed Vit filter: 50 (0.00299186); 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): 4 [number of targets reported over threshold] # CPU time: 0.30u 0.25s 00:00:00.55 Elapsed: 00:00:00.33 # Mc/sec: 3215.35 // [ok]