# 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_913.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_913 [L=277] Description: # 1017699 # 1018529 # 1 # ID=1_913;partial=00;start_type=ATG;rbs_motif=AGGAG;rbs_spacer=5-10bp;gc_cont=0.285 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 6.1e-23 81.7 0.0 8.5e-23 81.3 0.0 1.2 1 ABC_tran ABC transporter 0.00076 19.8 0.0 0.0012 19.2 0.0 1.3 1 AAA_16 AAA ATPase domain 0.0025 17.5 0.0 0.004 16.8 0.0 1.4 1 PRK Phosphoribulokinase / Uridine kinase family 0.0027 17.0 0.0 0.0046 16.3 0.0 1.4 1 Mg_chelatase Magnesium chelatase, subunit ChlI 0.0093 15.6 0.1 0.018 14.7 0.1 1.4 1 AAA_29 P-loop containing region of AAA domain ------ inclusion threshold ------ 0.019 13.9 0.0 0.026 13.5 0.0 1.3 1 T2SSE Type II/IV secretion system protein 0.021 14.2 0.0 0.047 13.0 0.0 1.6 2 PhoH PhoH-like protein 0.04 14.1 0.1 0.68 10.1 0.0 2.1 2 AAA_22 AAA domain 0.046 13.0 0.1 0.1 11.9 0.0 1.6 2 ATPase KaiC 0.077 12.8 0.1 0.18 11.6 0.0 1.7 2 RsgA_GTPase RsgA GTPase 0.095 11.4 0.0 0.14 10.9 0.0 1.2 1 DUF2791 P-loop Domain of unknown function (DUF2791) Domain annotation for each model (and alignments): >> ABC_tran ABC transporter # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 81.3 0.0 5.6e-26 8.5e-23 3 136 .. 30 170 .. 28 171 .. 0.90 Alignments for each domain: == domain 1 score: 81.3 bits; conditional E-value: 5.6e-26 EEEEEEETTSEEEEEESTTSSHHHHHHHHTTSS--...SEEEEEETTEEGGG.B-HHHCCHCEEEE-SS.--S-TTSBHHHH.. CS xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx...xxxxxxxxxxxxxx.xxxxxxxxxxxxxxxx.xxxxxxxxxxxx.. RF ABC_tran 3 nvslklkegekvaivGenGaGKStLlkllagllkp...teGeilldgkdlke.qeleslrkeigvlpqe.pqlfpeltvren.. 79 n+++++kege++ai G++G+GKS ++ +++g l + + ei+++++ ++e + + kei ++pq + l+p +t++++ 23726_FusoPortal_Gene_913 30 NLNVTIKEGEILAIAGSSGSGKSLMAHAILGILPKnavVSAEIKFKNEIVDEnRLYQLRGKEITFVPQSiAYLDPLMTIEDQlm 113 89******************************988777889999999999994455666799*******99***********99 PP HHHTCHHHHHHCSTTGGTTCBGESTTSGGS-HHHHHHHHHHHHHHHTSSEEEEESTT CS xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF ABC_tran 80 esdeeiekalsklglkelkdtvvksspssLSgGqkqrvalarallkkpkllllDEpt 136 +d + +++ + ++ +++++ +++ p +LSgG+ +rv +a+++l+k++l+++DEpt 23726_FusoPortal_Gene_913 114 RKDIDRQDFFKVMDTLGFTKSDLGKYPFQLSGGMARRVLIANTILSKADLIIADEPT 170 666667777777777778889999********************************9 PP >> AAA_16 AAA ATPase domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 19.2 0.0 7.8e-07 0.0012 19 144 .. 33 202 .. 21 216 .. 0.54 Alignments for each domain: == domain 1 score: 19.2 bits; conditional E-value: 7.8e-07 AAA_16 19 vrsgrpelvllsGesGvGKsslvrellralekdggif.lrgkcdql..esa................ipysa.laqalareell 82 v ++e++ + G sG+GKs ++ ++l l k+ +++ ++ k++ +++ i+y+ l ++ ++l+ 23726_FusoPortal_Gene_913 33 VTIKEGEILAIAGSSGSGKSLMAHAILGILPKN-AVVsAEIKFKNEivDENrlyqlrgkeitfvpqsIAYLDpLMTIE--DQLM 113 4455669********************888888.44444444444433333233333333111111144433122222..2333 PP AAA_16 83 rqllde.........aesaeleawrekleeala..................................avpelagaeasrllell 123 r+ d+ + + +++e ++ ++ +l+ ++ 23726_FusoPortal_Gene_913 114 RKDIDRqdffkvmdtL----------------GftksdlgkypfqlsggmarrvliantilskadliIADEPTPGLSLDLAVEV 181 3333332333333330................2333334444444444444777777776666666655555555556666666 PP AAA_16 124 krllallsrderplvlvlDdl 144 ++++ +++d++ ++l+ d+ 23726_FusoPortal_Gene_913 182 LNHFRNMANDGKGILLISHDI 202 666666666666666665555 PP >> PRK Phosphoribulokinase / Uridine kinase family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 16.8 0.0 2.6e-06 0.004 1 39 [. 40 78 .. 40 104 .. 0.86 Alignments for each domain: == domain 1 score: 16.8 bits; conditional E-value: 2.6e-06 PRK 1 vIgvagsSGaGkttvarklekifgrekvkaaeiegdsfl 39 ++++agsSG+Gk+ +a+++ i+ +++v +aei+ + ++ 23726_FusoPortal_Gene_913 40 ILAIAGSSGSGKSLMAHAILGILPKNAVVSAEIKFKNEI 78 689***************************999766544 PP >> Mg_chelatase Magnesium chelatase, subunit ChlI # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 16.3 0.0 3e-06 0.0046 15 50 .. 31 66 .. 21 90 .. 0.83 Alignments for each domain: == domain 1 score: 16.3 bits; conditional E-value: 3e-06 Mg_chelatase 15 leiaaaGgHnlllvGppGsGKtmlakrlpsllppls 50 l++++ g l ++G++GsGK+++a+++ ++lp+ 23726_FusoPortal_Gene_913 31 LNVTIKEGEILAIAGSSGSGKSLMAHAILGILPKNA 66 6677777778899********************864 PP >> AAA_29 P-loop containing region of AAA domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 14.7 0.1 1.2e-05 0.018 17 46 .. 33 62 .. 26 70 .. 0.82 Alignments for each domain: == domain 1 score: 14.7 bits; conditional E-value: 1.2e-05 AAA_29 17 ieldlkgvtlltGpsGsGKSTllDAiqalL 46 +++++ ++ +++G sGsGKS + Ai +L 23726_FusoPortal_Gene_913 33 VTIKEGEILAIAGSSGSGKSLMAHAILGIL 62 677766699***********9988887666 PP >> T2SSE Type II/IV secretion system protein # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.5 0.0 1.7e-05 0.026 121 164 .. 30 71 .. 5 116 .. 0.80 Alignments for each domain: == domain 1 score: 13.5 bits; conditional E-value: 1.7e-05 T2SSE 121 vlrkavrerglilvsGgTgsGKttllyallkeinksderivtvE 164 l+ ++e ++ ++G gsGK+ ++ a+l +++k+ +v +E 23726_FusoPortal_Gene_913 30 NLNVTIKEGEILAIAGSSGSGKSLMAHAILGILPKN--AVVSAE 71 3677788888999*********************98..555555 PP >> PhoH PhoH-like protein # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.0 0.0 3.1e-05 0.047 15 80 .. 34 106 .. 18 131 .. 0.64 2 ? -2.8 0.0 2.2 3.4e+03 186 203 .. 230 247 .. 223 248 .. 0.81 Alignments for each domain: == domain 1 score: 13.0 bits; conditional E-value: 3.1e-05 HHHH-SEEEEE--TTSSTTHHHHHHHHHHHHTTS-SEEEEEE-S--T........T----SS---------TTT CS PhoH 15 alkkndivfaiGpaGtGKtylavalavealkkkkvkriiltrpavea........geelGflpgdleekvdpyl 80 ++k+ +i+ ++G++G+GK+++a a+ + k+ v+ i +++ + g+e++f+p ++dp + 23726_FusoPortal_Gene_913 34 TIKEGEILAIAGSSGSGKSLMAHAILGILPKNAVVSAEIKFKNEIVDenrlyqlrGKEITFVPQS-IAYLDPLM 106 466889****************9875555555556655555555543123333466777777664.34555544 PP == domain 2 score: -2.8 bits; conditional E-value: 2.2 --GGG----HHHHHHHHH CS PhoH 186 lkkkdvvRhplvkkivka 203 lk ++ +Rhpl+k++ ka 23726_FusoPortal_Gene_913 230 LKGEKYIRHPLTKAFWKA 247 6778999***99987665 PP >> AAA_22 AAA domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 10.1 0.0 0.00045 0.68 4 29 .. 37 62 .. 33 91 .. 0.84 2 ? 1.9 0.0 0.16 2.4e+02 82 118 .. 153 191 .. 112 205 .. 0.77 Alignments for each domain: == domain 1 score: 10.1 bits; conditional E-value: 0.00045 AAA_22 4 tagilvvtGesGtGKTtllrrlleql 29 + il++ G sG+GK++ + +l l 23726_FusoPortal_Gene_913 37 EGEILAIAGSSGSGKSLMAHAILGIL 62 6789**************99887666 PP == domain 2 score: 1.9 bits; conditional E-value: 0.16 AAA_22 82 leqalraravielliiDEaqh.lrs....eealeeLrsllnl 118 +++ + ++a+ l+i DE l s e+l+++r+++n+ 23726_FusoPortal_Gene_913 153 IANTILSKAD--LIIADEPTPgL-SldlaVEVLNHFRNMAND 191 4455555555..99999999989.555669999999999994 PP >> ATPase KaiC # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.9 0.0 6.6e-05 0.1 15 37 .. 34 56 .. 14 68 .. 0.89 2 ? -3.0 0.0 2.4 3.6e+03 118 149 .. 95 127 .. 75 132 .. 0.64 Alignments for each domain: == domain 1 score: 11.9 bits; conditional E-value: 6.6e-05 ATPase 15 GipegrvvlitGepGtGKtifal 37 i eg++ i+G++G+GK+++a 23726_FusoPortal_Gene_913 34 TIKEGEILAIAGSSGSGKSLMAH 56 5889*****************96 PP == domain 2 score: -3.0 bits; conditional E-value: 2.4 ATPase 118 ikakrvviDsitgl.ayeekeavaReilrqlvr 149 + ++ +D++ ++ ++ ++++ R+ +++++ 23726_FusoPortal_Gene_913 95 VPQSIAYLDPLMTIeDQLMRKDIDRQDFFKVMD 127 455555667777775556666666666555544 PP >> RsgA_GTPase RsgA GTPase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.6 0.0 0.00012 0.18 82 121 .. 20 60 .. 3 89 .. 0.80 2 ? -3.9 0.0 6.8 1e+04 52 63 .. 157 168 .. 154 178 .. 0.67 Alignments for each domain: == domain 1 score: 11.6 bits; conditional E-value: 0.00012 --TTTHTTHHHHHHHCTT.SEEEEECSTTSSHHHHHHHHCT CS RsgA_GTPase 82 vsaktkegleelkellkd.ktsvlvGqSGvGKStLinallp 121 +++k ++ + +l+ +k+ ++ +++G+SG GKS + a+l+ 23726_FusoPortal_Gene_913 20 LNQKLHQMIFNLNVTIKEgEILAIAGSSGSGKSLMAHAILG 60 5666678888899999988**************99888876 PP == domain 2 score: -3.9 bits; conditional E-value: 6.8 EEE-GGG--HHH CS RsgA_GTPase 52 vlnKiDlleeee 63 +l+K+Dl+ +e 23726_FusoPortal_Gene_913 157 ILSKADLIIADE 168 567777765444 PP >> DUF2791 P-loop Domain of unknown function (DUF2791) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 10.9 0.0 8.9e-05 0.14 44 88 .. 33 77 .. 6 95 .. 0.76 Alignments for each domain: == domain 1 score: 10.9 bits; conditional E-value: 8.9e-05 DUF2791 44 vaeggaafrfirGeYGsGKtFllqllreralergfvvaeadlspe 88 v+ ++++ +i G GsGK+ +++++ + ++ +v ae+++ +e 23726_FusoPortal_Gene_913 33 VTIKEGEILAIAGSSGSGKSLMAHAILGILPKNAVVSAEIKFKNE 77 6667788999****************************9987654 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (277 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 427 (0.0255505); expected 334.2 (0.02) Passed bias filter: 385 (0.0230373); expected 334.2 (0.02) Passed Vit filter: 64 (0.00382958); expected 16.7 (0.001) Passed Fwd filter: 11 (0.00065821); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 11 [number of targets reported over threshold] # CPU time: 0.25u 0.23s 00:00:00.48 Elapsed: 00:00:00.31 # Mc/sec: 2597.57 // [ok]