# 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_1574.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1574 [L=353] Description: # 1685900 # 1686958 # -1 # ID=1_1574;partial=00;start_type=ATG;rbs_motif=GGAG/GAGG;rbs_spacer=5-10bp;gc_cont=0.252 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.8e-84 283.3 36.6 2.3e-84 283.0 36.6 1.0 1 FecCD FecCD transport family 2.4e-05 23.9 24.0 2.4e-05 23.9 24.0 2.2 2 ABC-3 ABC 3 transport family ------ inclusion threshold ------ 0.014 15.4 0.9 0.014 15.4 0.9 3.9 3 Phage_holin_3_6 Putative Actinobacterial Holin-X, holin supe 0.014 14.6 1.0 0.022 14.0 1.0 1.4 1 Peptidase_M50 Peptidase family M50 Domain annotation for each model (and alignments): >> FecCD FecCD transport family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 283.0 36.6 5.5e-88 2.3e-84 1 310 [. 33 345 .. 33 347 .. 0.97 Alignments for each domain: == domain 1 score: 283.0 bits; conditional E-value: 5.5e-88 FecCD 1 llllllllslllgavkispsevfeallgka.dsleaeivlelrlprvlaallvGaaLavaGallQtvtrNpLadPsilGissg 82 lllll+ll++++g+++is++e+++a++ ++ + ++ i+ ++rlp l+a++vG++L++ G+ +Qt+++Np+a+P++lGi+s+ 23726_FusoPortal_Gene_1574 33 LLLLLFLLNISIGSTSISIKEILKAIFINEgEGNNILIIRKIRLPMSLMAIVVGFSLGIGGCEIQTILKNPIASPYTLGITSA 115 6899********************9886665778888********************************************** PP FecCD 83 aalgavlallll...lgslslfalplaaflgallalllvlllarkrkgasplkliLaGiavsallsaltsllllllneeales 162 a++ga+lal+l l + + a++ af++ +l++ +++l++ +++g+ ++ +iL+Gia+++l++++t++l+++++e++l++ 23726_FusoPortal_Gene_1574 116 ASFGAALALVLNnsiLNLPDTLAVTGNAFFFTFLISTMIYLFS-SQRGIGKTAIILFGIALNFLFTSFTMILQYVADEDKLQN 197 **********9999888999999********************.*************************************** PP FecCD 163 ilfwllGslsgaswskvllllvllllllllllllarkldllalgdetAkslGvnvervrllllllvalltavavalvGpigFv 245 ++fw +Gsl + +w+k+l++l++l++++++l+ + kl++++lgd+ A+s+Gvn+ ++r ++l+v+ll+avav++vG+igF+ 23726_FusoPortal_Gene_1574 198 LIFWNFGSLLKTTWTKFLIVLIVLIICIIFLYKNFWKLTAMTLGDAKAESIGVNPHKLRKQIILIVSLLSAVAVSFVGTIGFI 280 *********************************************************************************** PP FecCD 246 GLivPhiarllvgadhrkllplsallGallllladllarlllsptelpvgivtaliGapffiyll 310 GLi+Phiarl+vg+d r +plsallG+++l l+ +++++++ + lp+g+vt++iG+pffi ++ 23726_FusoPortal_Gene_1574 281 GLIAPHIARLIVGEDQRFFMPLSALLGGFILSLSFFISKVVIPGVILPIGLVTSVIGIPFFISMI 345 **************************************************************887 PP >> ABC-3 ABC 3 transport family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -1.8 1.6 0.38 1.6e+03 127 145 .. 24 39 .. 14 51 .. 0.57 2 ! 23.9 24.0 5.7e-09 2.4e-05 38 235 .. 112 321 .. 79 347 .. 0.79 Alignments for each domain: == domain 1 score: -1.8 bits; conditional E-value: 0.38 ABC-3 127 liliaiieavvllllllfw 145 +i +i+++++lllllf 23726_FusoPortal_Gene_1574 24 VI---VINIIIILLLLLFL 39 33...34444444444433 PP == domain 2 score: 23.9 bits; conditional E-value: 5.7e-09 ABC-3 38 lsHaalpGvalayllgvn.........lllgailvgliaavalellkrktklkedtai..givlsvflalglvlvslikk..s 107 + aa G ala++l+ + ++ a++++++++ ++ l +++ + + ++i gi l + + ++++++ + 23726_FusoPortal_Gene_1574 112 ITSAASFGAALALVLNNSilnlpdtlaVTGNAFFFTFLISTMIYLFSSQRGIGKTAIIlfGIALNFLFTSFTMILQYVADedK 194 55677789999999998888898988666778999999988888888777776665434499998877544455555552244 PP ABC-3 108 kvdldsyLfGniltisqedliliaiieavvllllllfwkelllitfdeelAkveglsvkllkllllalialavvvalqavGvi 190 ++l + fG++l + + +++++i+ +++++ l + +l+++t+ + A+++g++++ l+ ++++++l +va+ +vG i 23726_FusoPortal_Gene_1574 195 LQNLIFWNFGSLLKTTWTKFLIVLIVLIICIIFLYKNFWKLTAMTLGDAKAESIGVNPHKLRKQIILIVSLLSAVAVSFVGTI 277 55566666899999999999999999999999999888999****************************************** PP ABC-3 191 LvialliiPaatArrltrslesllivAillglvsvvaGlllSyyl 235 i li+P ++ + ++++ + + +++llg + + + ++S ++ 23726_FusoPortal_Gene_1574 278 GFIG-LIAPHIARLIVGEDQRFFMPLSALLGGFILSLSFFISKVV 321 9998.589*998888999999999999999998888877777654 PP >> Phage_holin_3_6 Putative Actinobacterial Holin-X, holin superfamily III # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -8.4 9.9 4 1.7e+04 47 47 .. 111 111 .. 20 160 .. 0.62 2 ? 15.4 0.9 3.3e-06 0.014 35 84 .. 166 229 .. 157 246 .. 0.80 3 ? -6.9 8.9 4 1.7e+04 38 60 .. 263 304 .. 253 347 .. 0.48 Alignments for each domain: == domain 1 score: -8.4 bits; conditional E-value: 4 Phage_holin_3_6 47 a 47 + 23726_FusoPortal_Gene_1574 111 G 111 1 PP == domain 2 score: 15.4 bits; conditional E-value: 3.3e-06 Phage_holin_3_6 35 kgagllgaalvlallalllllvala..............lalaa...lwlaalivaavllvlagvla 84 ++++l+g+al +++++++++l +a + w +liv +vl++ + l 23726_FusoPortal_Gene_1574 166 TAIILFGIALNFLFTSFTMILQYVAdedklqnlifwnfgS---LlktTWTKFLIVLIVLIICIIFLY 229 6789********************99*****987654321...2567*******9777766665554 PP == domain 3 score: -6.9 bits; conditional E-value: 4 Phage_holin_3_6 38 gllgaalvlallalllllvala...................l 60 + l+ a+++ +++ + ++ ++a + 23726_FusoPortal_Gene_1574 263 VSLLSAVAVSFVGTIGFIGLIAphiarlivgedqrffmplsA 304 333333444444444444333333334344444443333331 PP >> Peptidase_M50 Peptidase family M50 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 14.0 1.0 5.2e-06 0.022 78 164 .. 140 294 .. 78 344 .. 0.83 Alignments for each domain: == domain 1 score: 14.0 bits; conditional E-value: 5.2e-06 .................XXXXXXXXXXXXXXXXXXXXXXXXXXXXXX.................................... CS x................xxxxxxxxxxxxxxxxxxxxxxxxxxxxxx.................................... RF Peptidase_M50 78 a................flkkspwqrllivlAGplanfilailvfvv.................................... 108 f++++ + + +i+l G++ nf+++ + +++ 23726_FusoPortal_Gene_1574 140 TgnaffftflistmiylFSSQRGIGKTAIILFGIALNFLFTSFTMILqyvadedklqnlifwnfgsllkttwtkflivlivli 222 255666666666667779999************************************************************** PP ..........-EEECE.CTTTSCEE-SEEEC-.....XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX CS ..........xxxxxx.xxxxxxxxxxxxxxx.....xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF Peptidase_M50 109 ..........pkkske.seeekrgviGvspkp.....vkkegllealkegvektlkvvvltvkalgklitgk 164 +k ++ k +iGv+p+ + +ll+a++ ++ t+ + l++ + +li+g+ 23726_FusoPortal_Gene_1574 223 iciiflyknfWKLTAMtLGDAKAESIGVNPHKlrkqiILIVSLLSAVAVSFVGTIGFIGLIAPHIARLIVGE 294 ******9999444444567788899*****9977777777899999**********************9994 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (353 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 1147 (0.0686333); expected 334.2 (0.02) Passed bias filter: 357 (0.0213619); expected 334.2 (0.02) Passed Vit filter: 56 (0.00335089); expected 16.7 (0.001) Passed Fwd filter: 16 (0.000957396); 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.23s 00:00:00.53 Elapsed: 00:00:00.31 # Mc/sec: 3310.27 // [ok]