# 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_1001.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1001 [L=172] Description: # 1108267 # 1108782 # 1 # ID=1_1001;partial=00;start_type=TTG;rbs_motif=None;rbs_spacer=None;gc_cont=0.264 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2.4e-06 27.7 4.3 2.4e-06 27.7 4.3 1.5 2 OMP_b-brl Outer membrane protein beta-barrel domain 0.0027 17.4 2.4 0.011 15.4 1.8 2.0 2 OmpA_membrane OmpA-like transmembrane domain 0.0084 15.9 0.1 0.01 15.6 0.1 1.1 1 LamB LamB porin Domain annotation for each model (and alignments): >> OMP_b-brl Outer membrane protein beta-barrel domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 1.5 1.0 0.049 2.7e+02 110 165 .. 18 69 .. 9 72 .. 0.45 2 ! 27.7 4.3 4.3e-10 2.4e-06 20 141 .. 45 170 .. 36 171 .. 0.78 Alignments for each domain: == domain 1 score: 1.5 bits; conditional E-value: 0.049 OMP_b-brl 110 gytkvevdvsgngasdsdtdfglglga..Gveyalad..nwslkaeYrytdlksdadndl 165 +y ++e+d++++ + +++ g++l a +++ ++s +eY+y ++ n++ 23726_FusoPortal_Gene_1001 18 TYIEFETDINKV---FGRKKHGYSLLAslKGTSNIGYgvQFSNVLEYEYLNY-----NQY 69 455566666555...3333333333221112222222124444666666666.....444 PP == domain 2 score: 27.7 bits; conditional E-value: 4.3e-10 OMP_b-brl 20 ggfyvggglglskvkv......dgv..gsdsgfgfgvgagyefndnvgvegglsysgdkdysgsgktdikinlls...vpvna 91 g ++g g+++s+v +++ + ++++ ++++ye+n+n+++++++++ k++++ + +n+l + ++ 23726_FusoPortal_Gene_1001 45 GTSNIGYGVQFSNVLEyeylnyNQYknAHKTKWETVLRWTYEINENIAFSPEVTFKVEKYNNS----K--ENYLIessAGPYV 121 55678888888877765555443332244458***********************99999888....3..3444445588888 PP OMP_b-brl 92 gyafnlsdsvlpyggaGlg.ytkvevdvsgngasdsdtdfglglgaGveya 141 + n++d +++yg++G+ + k e + g +s+++++ ++ +G+ey 23726_FusoPortal_Gene_1001 122 LFTKNINDDLRIYGKVGVPvFRKDESKAEGY--RYSKSQTSAYGKIGFEYI 170 888999************9766666777777..****************96 PP >> OmpA_membrane OmpA-like transmembrane domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 1.3 0.5 0.041 2.3e+02 146 164 .. 59 77 .. 14 116 .. 0.57 2 ! 15.4 1.8 2e-06 0.011 87 136 .. 125 169 .. 49 171 .. 0.73 Alignments for each domain: == domain 1 score: 1.3 bits; conditional E-value: 0.041 EEEEEES..SSSSBTSTTS.. CS OmpA_membrane 146 leYqyln..nvgdlrkaakky 164 leY yln + +++k+ k+ 23726_FusoPortal_Gene_1001 59 LEYEYLNynQYKNAHKT--KW 77 33333330022222222..22 PP == domain 2 score: 15.4 bits; conditional E-value: 2e-06 EESSSSEEEEEEEEEEEEEEEEESSSSB.....EEEEEEEEEEEEEEEEE CS OmpA_membrane 87 yellddldvYakvGvalvradykvkedkkskekknhdtkvslllavGvey 136 +++ddl +Y+kvGv + r+d +k++ + ++ ++s +G+ey 23726_FusoPortal_Gene_1001 125 KNINDDLRIYGKVGVPVFRKDE-----SKAEGYRYSKSQTSAYGKIGFEY 169 479******************9.....55555677888888888888888 PP >> LamB LamB porin # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 15.6 0.1 1.9e-06 0.01 244 314 .. 28 103 .. 3 125 .. 0.80 Alignments for each domain: == domain 1 score: 15.6 bits; conditional E-value: 1.9e-06 ........EEEEEEEEEEEEEETTTEEEEEEEEEEEEEE..TTS--EEEEEEEEEEEEECCTTEEEEEEEEEEEEE CS LamB 244 e.....kdakslrlleegvvqlsdkfslayavvyqksddrydngddrkwlsagvRpvykltdnfklllElGydrvk 314 + k ++sl +g+ +++ ++++ ++ y++ + ++ ++ ++ +++++R +y++++n++ E++++ k 23726_FusoPortal_Gene_1001 28 KvfgrkKHGYSLLASLKGTSNIGYGVQFSNVLEYEYLNYNQYKNAHKTKWETVLRWTYEINENIAFSPEVTFKVEK 103 1222336699999999**********************999999999999********************998887 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (172 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 438 (0.0262087); expected 334.2 (0.02) Passed bias filter: 204 (0.0122068); expected 334.2 (0.02) Passed Vit filter: 23 (0.00137626); 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.25u 0.28s 00:00:00.53 Elapsed: 00:00:00.33 # Mc/sec: 1515.18 // [ok]