# 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_1402.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1402 [L=497] Description: # 1509375 # 1510865 # 1 # ID=1_1402;partial=00;start_type=ATG;rbs_motif=GGAG/GAGG;rbs_spacer=5-10bp;gc_cont=0.250 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.9e-52 177.7 0.8 3.1e-52 176.9 0.8 1.4 1 GATase_3 CobB/CobQ-like glutamine amidotransferase doma 2.2e-22 79.8 0.1 4e-22 78.9 0.1 1.4 1 AAA_26 AAA domain 6.3e-20 71.4 0.1 1.5e-19 70.2 0.1 1.6 1 CbiA CobQ/CobB/MinD/ParA nucleotide binding domain 0.00031 20.5 0.2 0.0014 18.4 0.0 2.0 2 GATase Glutamine amidotransferase class-I 0.0044 16.6 0.7 0.013 15.1 0.1 2.0 3 Peptidase_S51 Peptidase family S51 ------ inclusion threshold ------ 0.026 14.3 2.5 0.19 11.5 0.0 3.1 4 MobB Molybdopterin guanine dinucleotide synthesis p 0.048 13.4 0.1 0.097 12.4 0.1 1.6 1 SNO SNO glutamine amidotransferase family Domain annotation for each model (and alignments): >> GATase_3 CobB/CobQ-like glutamine amidotransferase domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 176.9 0.8 1.3e-55 3.1e-52 1 196 [] 255 442 .. 255 442 .. 0.93 Alignments for each domain: == domain 1 score: 176.9 bits; conditional E-value: 1.3e-55 GATase_3 1 kiavirlpaisnytdldpLeaepevrvrfvplkdeslgedadlvilpGgkdtiqdlallresgldeaikeaaeaggpvlgiCG 83 ki vi+l++isn td+d+L ++v+++fv+ +++++g+ d++i+pG+k+ti+dl++l+esg++e+i + a++++ ++giCG 23726_FusoPortal_Gene_1402 255 KISVIKLKHISNATDIDALSIYNDVEIQFVT-ERSQIGN-EDFIIIPGSKNTIDDLKWLKESGIAEEIIKRARTETIIFGICG 335 699****************************.*****87.9****************************************** PP GATase_3 84 GlqmLGeeiedpegvegeageveglglLdietefekekltgqvegevel......aaegeevkGyEiHsGrtilkekakpllk 160 G+q+LG++++d +++eg+++e++glglLd+et +e+ek++ q +g++ + + ++ e+kGyEiH+G t+ ++k l + 23726_FusoPortal_Gene_1402 336 GFQILGNKVKDLYHIEGDIEELNGLGLLDLETIMENEKTLVQYKGKLVVdngilkTLNDFEIKGYEIHQGITKG--NEKNLTT 416 **********************************************8779**9988888***********7764..5577777 PP GATase_3 161 vkvggnngekdGavkksedgnvlGtYlHglfdnrna 196 e+ +v +++n+++tYlHg+fdn+++ 23726_FusoPortal_Gene_1402 417 D-------ERTIFV---NRDNIIATYLHGIFDNKDF 442 6.......444455...458*************975 PP >> AAA_26 AAA domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 78.9 0.1 1.7e-25 4e-22 2 193 .. 5 228 .. 4 234 .. 0.83 Alignments for each domain: == domain 1 score: 78.9 bits; conditional E-value: 1.7e-25 AAA_26 2 timvtgtdtgvGKTlvslgLvsalkrkgnkvgylKPvqtgle........sdsdaelvqevlrldetheeaeplalseplsph 76 ++m++gt+ g+GK+l ++L++++ + ++kv+ +K + l+ + + a++vq ++ +e++ ++p++l+ ++ 23726_FusoPortal_Gene_1402 5 NLMIVGTSSGAGKSLFVTALCRIFYKDKYKVSPFKSQNMALNsyitkdgkEMGRAQVVQAEASGTEPNVNMNPILLKPSTINK 87 59******************************9998766666666655545669999987777777777899**999888888 PP AAA_26 77 laaaregstvd....lekiv..............tskelekdrdivvvEGagglavpitse..slnadiakslgapvlLVarg 139 ++ g+++ + t++++e+++divv+EGag+ a+ ++e +n+ +a+ +apv+LVa+ 23726_FusoPortal_Gene_1402 88 IQIIVCGKSIGnmsgV---EynqykknlipilkeTYSKIENKNDIVVIEGAGSPAEINIKEedISNFAMARIADAPVILVADI 167 8887778777764443...2444679999**************************7755443689****************** PP AAA_26 140 glgt.inqtllaiealrsr.gvrilgviinkvnep.....envetlrkkipvpvlgavPed 193 ++g+ ++ ++ +i +l+++ ri+g++ink++ + e +++ ++v+ lg++P+ 23726_FusoPortal_Gene_1402 168 DRGGvFASIYGTIMLLHEEdRKRIKGIVINKFRGNkevlkPGFEIIENLTGVKTLGVIPYT 228 99762699999999998751568********87774455559*****************97 PP >> CbiA CobQ/CobB/MinD/ParA nucleotide binding domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 70.2 0.1 6.4e-23 1.5e-19 1 120 [. 6 233 .. 6 239 .. 0.97 Alignments for each domain: == domain 1 score: 70.2 bits; conditional E-value: 6.4e-23 EEEEESSTTSSHHHHHHHHHHHHHHTT--EEEEEEEESSCCCCH.HTTEB-HHHH............................ CS xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx.xxxxxxxxxx............................ RF CbiA 1 iaiastkGGvGKTtlaanLaralarrglrVllvDaDpqannlai.leglendiep............................ 54 ++i++t++G+GK+ ++++L+r++ +++++V +++ +++a n+ i +g+e+ +++ 23726_FusoPortal_Gene_1402 6 LMIVGTSSGAGKSLFVTALCRIFYKDKYKVSPFKSQNMALNSYItKDGKEMGRAQvvqaeasgtepnvnmnpillkpstinki 88 69*****************************************999999999999**************************** PP .............................................................................CCHCCC CS .............................................................................xxxxxx RF CbiA 55 .............................................................................dlkrll 60 d++r + 23726_FusoPortal_Gene_1402 89 qiivcgksignmsgveynqykknlipilketyskienkndivviegagspaeinikeedisnfamariadapvilvaDIDR-G 170 *********************************************************************************.* PP HHHHHHCHCHHHH.......EEEEEEEEESS...-TTCHHHHHHHHHHCTS-EEEEEE--CCHCC CS xxxxxxxxxxxxx.....xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF CbiA 61 ellealgkkkall.....lkiagvvlnrvdgeneetllkealeelkqrygvpvlgviprdeavae 120 +++++++++ ll +i+g+v+n+++g+ ++ lk ++e ++ +gv+ lgvip+++ +e 23726_FusoPortal_Gene_1402 171 GVFASIYGTIMLLheedrKRIKGIVINKFRGN--KEVLKPGFEIIENLTGVKTLGVIPYTDIDIE 233 ***********99*******************..**********99988**********997665 PP >> GATase Glutamine amidotransferase class-I # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -0.7 0.0 0.4 9.5e+02 35 64 .. 118 151 .. 99 162 .. 0.63 2 ! 18.4 0.0 5.8e-07 0.0014 35 91 .. 285 373 .. 246 457 .. 0.65 Alignments for each domain: == domain 1 score: -0.7 bits; conditional E-value: 0.4 HHHHTHHSTSEEEEEEESSCTT....TTTHHHHH CS GATase 35 eilellknpdgilvsgGpgspg....vegkieai 64 e++++ +n + i+v G gsp+ e+ i++ 23726_FusoPortal_Gene_1402 118 ETYSKIENKNDIVVIEGAGSPAeiniKEEDISNF 151 3333567899999999999999844434444443 PP == domain 2 score: 18.4 bits; conditional E-value: 5.8e-07 HHHHTHHSTSEEEEEEESSCTT.......TTTHHHHHHHHCTSTS-EEEETHHHHHHHHHTTT.................... CS GATase 35 eilellknpdgilvsgGpgspg.......vegkieaikearenkiPiLGICLGmQllaiefgr.................... 90 +++++ n d i+++g ++ + + e+ik+ar+ ++ i+GIC G+Q+l++++ 23726_FusoPortal_Gene_1402 285 TERSQIGNEDFIIIPGSKNTIDdlkwlkeSGIAEEIIKRART-ETIIFGICGGFQILGNKVKDlyhiegdieelnglglldle 366 222245678888888888888799999999999*********.777***********99987766666666655555555554 PP ......E CS GATase 91 ......k 91 k 23726_FusoPortal_Gene_1402 367 timeneK 373 4444440 PP >> Peptidase_S51 Peptidase family S51 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -2.4 0.0 1.2 2.8e+03 116 143 .. 5 34 .. 3 46 .. 0.60 2 ? -3.6 0.0 2.8 6.8e+03 67 90 .. 114 134 .. 108 134 .. 0.77 3 ! 15.1 0.1 5.2e-06 0.013 63 158 .. 270 366 .. 251 372 .. 0.83 Alignments for each domain: == domain 1 score: -2.4 bits; conditional E-value: 1.2 T-EEEEETHHH..HTTSSB-GGGTTTS-GG CS Peptidase_S51 116 glvyiGvSAGa..ivagesiettsysdspi 143 l+++G+S Ga + + + + + y+d+ + 23726_FusoPortal_Gene_1402 5 NLMIVGTSSGAgkSLFVTALCRIFYKDKYK 34 57999**97772245555555555555555 PP == domain 2 score: -3.6 bits; conditional E-value: 2.8 -TTTS-HHHHHHHHHHSSEEEE-- CS Peptidase_S51 67 elskdpeeeieeklleadviyvgG 90 ++ k++ ++ie+ ++d+++++G 23726_FusoPortal_Gene_1402 114 PILKETYSKIEN---KNDIVVIEG 134 566777788888...899999988 PP == domain 3 score: 15.1 bits; conditional E-value: 5.2e-06 EEE--TTTS-HHH.HHH..HHHHSSEEEE--S.-HHHHHHHHHHTTHHHHHHHHHHTT-EEEEETHHHHTTSSB-GGGTTTS- CS Peptidase_S51 63 vdeLelskdpeee.iee..klleadviyvgGG.ntfrLlkllketgldkllreavqqglvyiGvSAGaivagesiettsysds 141 +d L++ +d e + ++e ++ ++d+i++ G nt++ lk lke g+++ ++++ + ++++G+ G + g++++ + 23726_FusoPortal_Gene_1402 270 IDALSIYNDVEIQfVTErsQIGNEDFIIIPGSkNTIDDLKWLKESGIAEEIIKRARTETIIFGICGGFQILGNKVKDLYHIEG 352 66666666665554554558999******9976***************************************99988773333 PP GGG-S---B----SSSE CS Peptidase_S51 142 pikapesfeaLgLvkfq 158 e++++LgL++ 23726_FusoPortal_Gene_1402 353 D---IEELNGLGLLDLE 366 3...4889999999865 PP >> MobB Molybdopterin guanine dinucleotide synthesis protein B # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.5 0.0 8e-05 0.19 2 48 .. 6 55 .. 5 70 .. 0.83 2 ? 1.0 0.0 0.15 3.5e+02 79 98 .. 114 134 .. 91 165 .. 0.79 3 ? -1.8 0.1 1 2.5e+03 26 42 .. 252 270 .. 247 273 .. 0.58 4 ? -3.2 0.1 2.8 6.6e+03 91 114 .. 460 484 .. 444 491 .. 0.70 Alignments for each domain: == domain 1 score: 11.5 bits; conditional E-value: 8e-05 MobB 2 lgivgk.sgsGKTtLleklipelkarglrvgtiKhahhgfd..vdkpgkD 48 l ivg+ sg+GK +++ l +++ + +++v+ K ++ +++ + k gk 23726_FusoPortal_Gene_1402 6 LMIVGTsSGAGKSLFVTALCRIFYKDKYKVSPFKSQNMALNsyITKDGKE 55 578877589*************************9998877556666665 PP == domain 2 score: 1.0 bits; conditional E-value: 0.15 MobB 79 lsleellarl.sevDlvlvEG 98 + l+e +++ ++ D+v++EG 23726_FusoPortal_Gene_1402 114 PILKETYSKIeNKNDIVVIEG 134 55788888885569******9 PP == domain 3 score: -1.8 bits; conditional E-value: 1 MobB 26 rglrvgtiK..hahhgfdv 42 +++++++iK h +++ d+ 23726_FusoPortal_Gene_1402 252 NKIKISVIKlkHISNATDI 270 6677777761144455555 PP == domain 4 score: -3.2 bits; conditional E-value: 2.8 MobB 91 vDlvlvEGfkkealpkie.vlrael 114 D + E +k ++++k+e ++r+++ 23726_FusoPortal_Gene_1402 460 NDNISYEEYKLKEFDKLEkLVRENV 484 5777778899999999874555554 PP >> SNO SNO glutamine amidotransferase family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.4 0.1 4.1e-05 0.097 22 87 .. 281 346 .. 259 442 .. 0.80 Alignments for each domain: == domain 1 score: 12.4 bits; conditional E-value: 4.1e-05 EEEE-SGGGGCC-SEEEE-SS-HHHHH.HHHHHTT-HHHHHHHHHTT.--EEEETHHHHHHCCEETT CS SNO 22 nktvkrkedlaevdaliiPGGestams.kliertglyedllefvknsekvvlGtCaGlillskelen 87 + v ++ ++ + d +iiPG t k +++ g+ e++ + ++ e +++G C G+ +l ++++ 23726_FusoPortal_Gene_1402 281 IQFVTERSQIGNEDFIIIPGSKNTIDDlKWLKESGIAEEIIKRART-ETIIFGICGGFQILGNKVKD 346 367889999999*********9996542678999*99999988776.59**********99988754 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (497 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 726 (0.0434418); expected 334.2 (0.02) Passed bias filter: 397 (0.0237554); expected 334.2 (0.02) Passed Vit filter: 52 (0.00311154); expected 16.7 (0.001) Passed Fwd filter: 7 (0.000418861); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 7 [number of targets reported over threshold] # CPU time: 0.33u 0.26s 00:00:00.59 Elapsed: 00:00:00.36 # Mc/sec: 4013.32 // [ok]