Align methionine synthase; EC 2.1.1.13 (characterized)
to candidate Ac3H11_3274 5-methyltetrahydrofolate--homocysteine methyltransferase (EC 2.1.1.13)
Query= CharProtDB::CH_090726 (1227 letters) >FitnessBrowser__acidovorax_3H11:Ac3H11_3274 Length = 903 Score = 1052 bits (2721), Expect = 0.0 Identities = 566/908 (62%), Positives = 666/908 (73%), Gaps = 32/908 (3%) Query: 341 RLSGLEPLNIGEDSLFVNVGERTNVTGSAKFKRLIKEEKYSEALDVARQQVENGAQIIDI 400 +LSGLEP++IGE +LFVN+GERTNVTGS F R+I +Y EAL VARQQVENGAQ+IDI Sbjct: 2 KLSGLEPVSIGEGTLFVNIGERTNVTGSKAFARMILNGQYEEALAVARQQVENGAQVIDI 61 Query: 401 NMDEGMLDAEAAMVRFLNLIAGEPDIARVPIMIDSSKWDVIEKGLKCIQGKGIVNSISMK 460 NMDE MLD++AAMVRFL LIA EPDIARVPIM+DSSKWDVIE GL+CIQGKGIVNSISMK Sbjct: 62 NMDEAMLDSKAAMVRFLQLIASEPDIARVPIMVDSSKWDVIEAGLRCIQGKGIVNSISMK 121 Query: 461 EGVDAFIHHAKLLRRYGAAVVVMAFDEQGQADTRARKIEICRRAYKILTEEVGFPPEDII 520 EGV+ F H A+L+RRYGAA VVMAFDE GQADT ARKIEIC RAY++L +EVGFPPEDII Sbjct: 122 EGVEKFKHEARLVRRYGAAAVVMAFDEVGQADTYARKIEICERAYRVLVDEVGFPPEDII 181 Query: 521 FDPNIFAVATGIEEHNNYAQDFIGACEDIKRELPHALISGGVSNVSFSFRGNDPVREAIH 580 FDPNIFAVATGIEEHNNYA DFI A IK+ LP A +SGGVSNVSFSFRGNDPVREAIH Sbjct: 182 FDPNIFAVATGIEEHNNYAVDFIEATRWIKQNLPGAKVSGGVSNVSFSFRGNDPVREAIH 241 Query: 581 AVFLYYAIRNGMDMGIVNAGQLAIYDDLPAELRDAVEDVILNRRDDGTERLLELAEKYRG 640 VFLY+AI+ GMDMGIVNAG + +YDDL LR+ VEDV+LNRR D ERL+E+AE + Sbjct: 242 TVFLYHAIQAGMDMGIVNAGMVGVYDDLEPVLRERVEDVVLNRRADAGERLVEIAETAKS 301 Query: 641 SKTDDTANAQQAEWRSW-----EVNKRLEYSLVKGITEFIEQDTEEARQ----QATRPIE 691 D+ +++ EWR V +RL ++LV GIT+FI +DTEEA Q + RP+ Sbjct: 302 GAKDE---SKKLEWRGTPEHPKTVGERLSHALVHGITDFIVEDTEEAYQAIVAKGGRPLH 358 Query: 692 VIEGPLMDGMNVVGDLFGEGKMFLPQVVKSARVMKQAVAYLEPFIEASKEQG-------K 744 VIEGPLMDGMNVVGDLFG GKMFLPQVVKSARVMK AVA+L P+IE K Q + Sbjct: 359 VIEGPLMDGMNVVGDLFGAGKMFLPQVVKSARVMKSAVAHLIPYIEEEKRQDELAGRDVR 418 Query: 745 TNGKMVIATVKGDVHDIGKNIVGVVLQCNNYEIVDLGVMVPAEKILRTAKEVNADLIGLS 804 + GK++IATVKGDVHDIGKNIV VVLQCNN+E+V++GVMVP +IL AK AD++GLS Sbjct: 419 SKGKIIIATVKGDVHDIGKNIVTVVLQCNNFEVVNMGVMVPCHEILARAKVEGADIVGLS 478 Query: 805 GLITPSLDEMVNVAKEMERQGF----TIPLLIGGATTSKAHTAVKIEQNYSGPTVYVQNA 860 GLITPSL+EM VA EM++ IPLLIGGATTS+ HTAVKI +Y GP VYV +A Sbjct: 479 GLITPSLEEMQYVAGEMQKDDHFRIKKIPLLIGGATTSRVHTAVKIAPHYEGPVVYVPDA 538 Query: 861 SRTVGVVAALLSDTQRDDFVARTRKEYETVRIQHGRKKPRTPPVTLEAARDNDFAFDWQA 920 SR+V V +LL D + ++V +YE VR QH KK +TP +LE AR N ++ Sbjct: 539 SRSVSVAQSLLGDGVQ-NYVQELNADYEKVRTQHANKK-QTPMWSLEKARANRTPVNFST 596 Query: 921 YTPPVAHRLGVQEVE-ASIETLRNYIDWTPFFMTWSLAGKYPRILEDEVVGVEAQRLFKD 979 Y PP LG + + + L YIDW PFF TW LAG YP IL DEVVGVEA R+F D Sbjct: 597 YVPPTPRLLGRRVFKNFDLAELAKYIDWGPFFQTWDLAGPYPAILTDEVVGVEATRVFAD 656 Query: 980 ANDMLDKLSAEKTLNPRGVVGLFPANRVGDDIEIYRDETRTHVINVSHHLRQQTEKTGF- 1038 ML KL + L+ GV+ L+PAN VGDDIE Y DETRT V + +RQQTEK Sbjct: 657 GQAMLKKLIEGRWLSASGVMALYPANSVGDDIEFYTDETRTQVAMTWYGMRQQTEKQMID 716 Query: 1039 ----ANYCLADFVAPKLSGKADYIGAFAVTGGLEEDALADAFEAQHDDYNKIMVKALADR 1094 + CLADFVAPK SG ADY G FAVT GL + AF DDY+ IM+K+LADR Sbjct: 717 GVMRPSRCLADFVAPKDSGIADYAGLFAVTAGLGVEKKEQAFVDALDDYSAIMLKSLADR 776 Query: 1095 LAEAFAEYLHERVRKVYWGYAPNENLSNEELIRENYQGIRPAPGYPACPEHTEKATIWEL 1154 LAEAFAE LH+RVR WGYA E LSN+E+I E Y GIRPAPGYPACP+H+ K ++ + Sbjct: 777 LAEAFAECLHQRVRTDLWGYAAGEQLSNDEMIAEKYAGIRPAPGYPACPDHSAKQELFRV 836 Query: 1155 LEVEKHTGMKLTESFAMWPGASVSGWYFSHPDSKYYAVAQIQRDQVEDYARRKGMSVTEV 1214 L+ E+ GM LTES AM P ASVSG+Y HPDS Y+ V +I DQ+ D A R+GM + Sbjct: 837 LQCEE-IGMGLTESLAMTPAASVSGFYIGHPDSVYFNVGKIGEDQLHDMAERRGMDEKVL 895 Query: 1215 ERWLAPNL 1222 ER LAPNL Sbjct: 896 ERLLAPNL 903 Lambda K H 0.318 0.134 0.391 Gapped Lambda K H 0.267 0.0410 0.140 Matrix: BLOSUM62 Gap Penalties: Existence: 11, Extension: 1 Number of Sequences: 1 Number of Hits to DB: 2674 Number of extensions: 118 Number of successful extensions: 11 Number of sequences better than 1.0e-02: 1 Number of HSP's gapped: 2 Number of HSP's successfully gapped: 1 Length of query: 1227 Length of database: 903 Length adjustment: 45 Effective length of query: 1182 Effective length of database: 858 Effective search space: 1014156 Effective search space used: 1014156 Neighboring words threshold: 11 Window for multiple hits: 40 X1: 16 ( 7.3 bits) X2: 38 (14.6 bits) X3: 64 (24.7 bits) S1: 41 (21.7 bits) S2: 58 (26.9 bits)
Align methionine synthase; EC 2.1.1.13 (characterized)
to candidate Ac3H11_3297 5-methyltetrahydrofolate--homocysteine methyltransferase (EC 2.1.1.13)
Query= CharProtDB::CH_090726 (1227 letters) >FitnessBrowser__acidovorax_3H11:Ac3H11_3297 Length = 354 Score = 371 bits (953), Expect = e-107 Identities = 186/344 (54%), Positives = 253/344 (73%), Gaps = 15/344 (4%) Query: 3 SKVEQLRAQLNERILVLDGGMGTMIQSYRLNEADFRGE----------RFADWPCDLKGN 52 ++ QL L +RI++LDG MGTMIQ ++L EA +RGE RF D+ D+KGN Sbjct: 7 TRAAQLPDILAQRIVILDGAMGTMIQRFKLGEAQYRGEGYTGPGKAGDRFKDFAYDVKGN 66 Query: 53 NDLLVLSKPEVIAAIHNAYFEAGADIIETNTFNSTTIAMADYQMESLSAEINFAAAKLAR 112 N+LL L++P+VI+ IH Y AGAD+IETNTF +TTIA DY M +L+ E+N +A+LAR Sbjct: 67 NELLSLTRPDVISDIHERYLAAGADLIETNTFGATTIAQEDYHMANLAREMNLRSAELAR 126 Query: 113 ACADEWTARTPEKPRYVAGVLGPTNRTASISPDVNDPAFRNITFDGLVAAYRESTKALVE 172 A D+++ TP+KPR+VAG LGPT +TASISPDVNDP RN+ F+ L AAY E T+ALVE Sbjct: 127 AACDKYS--TPDKPRFVAGALGPTPKTASISPDVNDPGARNVDFEALRAAYYEQTEALVE 184 Query: 173 GGADLILIETVFDTLNAKAAVFAVKTEFEALGVELPIMISGTITDASGRTLSGQTTEAFY 232 GG+D++L+ET+FDTLNAKAA+FA+ FE G LP++ISGT+TDASGR LSGQT AF+ Sbjct: 185 GGSDVLLVETIFDTLNAKAALFAIDEYFENSGQRLPLIISGTVTDASGRILSGQTVTAFW 244 Query: 233 NSLRHAEALTFGLNCALGPDELRQYVQELSRIAE-CYVTAHPNAGLPNAFGE--YDLDAD 289 +S+RH+ L GLNCALG +R Y+QEL+R+AE +++ +PNAGLPN + +D + Sbjct: 245 HSVRHSRPLAIGLNCALGATLMRPYIQELNRVAEDTFISCYPNAGLPNPMSDTGFDETPE 304 Query: 290 TMAKQIREWAQAGFLNIVGGCCGTTPQHIAAMSRAVEGLAPRKL 333 ++ + E+A G +NI+GGCCGTTP HI A+++AVE + RK+ Sbjct: 305 ITSRLVHEFAAEGLVNILGGCCGTTPDHIGAIAKAVESVPTRKM 348 Lambda K H 0.318 0.134 0.391 Gapped Lambda K H 0.267 0.0410 0.140 Matrix: BLOSUM62 Gap Penalties: Existence: 11, Extension: 1 Number of Sequences: 1 Number of Hits to DB: 1007 Number of extensions: 32 Number of successful extensions: 5 Number of sequences better than 1.0e-02: 1 Number of HSP's gapped: 1 Number of HSP's successfully gapped: 1 Length of query: 1227 Length of database: 354 Length adjustment: 38 Effective length of query: 1189 Effective length of database: 316 Effective search space: 375724 Effective search space used: 375724 Neighboring words threshold: 11 Window for multiple hits: 40 X1: 16 ( 7.3 bits) X2: 38 (14.6 bits) X3: 64 (24.7 bits) S1: 41 (21.7 bits) S2: 54 (25.4 bits)
Align candidate Ac3H11_3274 (5-methyltetrahydrofolate--homocysteine methyltransferase (EC 2.1.1.13))
to HMM TIGR02082 (metH: methionine synthase (EC 2.1.1.13))
# hmmsearch :: search profile(s) against a sequence database # HMMER 3.3.1 (Jul 2020); http://hmmer.org/ # Copyright (C) 2020 Howard Hughes Medical Institute. # Freely distributed under the BSD open source license. # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - # query HMM file: ../tmp/path.aa/TIGR02082.hmm # target sequence database: /tmp/gapView.26825.genome.faa # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: TIGR02082 [M=1182] Accession: TIGR02082 Description: metH: methionine synthase Scores for complete sequences (score includes all domains): --- full sequence --- --- best 1 domain --- -#dom- E-value score bias E-value score bias exp N Sequence Description ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 0 1191.4 0.0 0 1191.2 0.0 1.0 1 lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274 5-methyltetrahydrofolate--homocy Domain annotation for each sequence (and alignments): >> lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274 5-methyltetrahydrofolate--homocysteine methyltransferase (EC 2.1.1.1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 1191.2 0.0 0 0 328 1182 .] 2 873 .. 1 873 [. 0.96 Alignments for each domain: == domain 1 score: 1191.2 bits; conditional E-value: 0 TIGR02082 328 vlsglealkiaqessfvniGeRtnvaGskkfrklikaedyeealkiakqqveeGaqilDi 387 lsgle+++i++ + fvniGeRtnv+Gsk f+++i +++yeeal +a+qqve+Gaq++Di lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274 2 KLSGLEPVSIGEGTLFVNIGERTNVTGSKAFARMILNGQYEEALAVARQQVENGAQVIDI 61 68********************************************************** PP TIGR02082 388 nvDevllDgeadmkkllsllasepdiakvPlmlDssefevleaGLkviqGkaivnsislk 447 n+De++lD++a+m+++l+l+asepdia+vP+m+Dss+++v+eaGL++iqGk+ivnsis+k lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274 62 NMDEAMLDSKAAMVRFLQLIASEPDIARVPIMVDSSKWDVIEAGLRCIQGKGIVNSISMK 121 ************************************************************ PP TIGR02082 448 dGeerFlekaklikeyGaavvvmafDeeGqartadkkieiakRayklltekvgfppedii 507 +G+e+F ++a+l+++yGaa+vvmafDe Gqa+t+++kiei++Ray++l+++vgfppedii lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274 122 EGVEKFKHEARLVRRYGAAAVVMAFDEVGQADTYARKIEICERAYRVLVDEVGFPPEDII 181 ************************************************************ PP TIGR02082 508 fDpniltiatGieehdryaidfieaireikeelPdakisgGvsnvsFslrgndavRealh 567 fDpni+++atGieeh++ya+dfiea+r+ik++lP ak+sgGvsnvsFs+rgnd+vRea+h lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274 182 FDPNIFAVATGIEEHNNYAVDFIEATRWIKQNLPGAKVSGGVSNVSFSFRGNDPVREAIH 241 ************************************************************ PP TIGR02082 568 svFLyeaikaGlDmgivnagklavyddidkelrevvedlildrrreatekLlelaelykg 627 +vFLy+ai+aG+Dmgivnag++ vydd+++ lre ved++l+rr +a e+L+e+ae+ k lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274 242 TVFLYHAIQAGMDMGIVNAGMVGVYDDLEPVLRERVEDVVLNRRADAGERLVEIAETAKS 301 *********************************************************776 PP TIGR02082 628 tkeksskeaqeaewrnlp.....veeRLeralvkGeregieedleear....kklkaple 678 +++ +++ ewr++p v eRL++alv+G++++i ed+eea+ k +pl+ lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274 302 GAKDE---SKKLEWRGTPehpktVGERLSHALVHGITDFIVEDTEEAYqaivAKGGRPLH 358 66664...44666666543333389***********************444446689*** PP TIGR02082 679 iiegpLldGmkvvGdLFGsGkmfLPqvvksarvmkkavayLePylekekeed.......k 731 +iegpL+dGm+vvGdLFG+GkmfLPqvvksarvmk ava+L+Py+e+ek +d + lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274 359 VIEGPLMDGMNVVGDLFGAGKMFLPQVVKSARVMKSAVAHLIPYIEEEKRQDelagrdvR 418 *************************************************9999******* PP TIGR02082 732 skGkivlatvkGDvhDiGknivdvvLscngyevvdlGvkvPvekileaakkkkaDvigls 791 skGki++atvkGDvhDiGkniv vvL+cn++evv++Gv+vP+++il ak + aD++gls lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274 419 SKGKIIIATVKGDVHDIGKNIVTVVLQCNNFEVVNMGVMVPCHEILARAKVEGADIVGLS 478 ************************************************************ PP TIGR02082 792 GLivksldemvevaeemerrg.vk...iPlllGGaalskahvavkiaekYkgevvyvkda 847 GLi++sl+em++va em++ ++ iPll+GGa++s+ h+avkia++Y+g+vvyv da lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274 479 GLITPSLEEMQYVAGEMQKDDhFRikkIPLLIGGATTSRVHTAVKIAPHYEGPVVYVPDA 538 *****************9975333446********************************* PP TIGR02082 848 seavkvvdkllsekkkaeelekikeeyeeirekfgekkeklialsekaarkevfaldrse 907 s++v+v+++ll + ++++++++++ye++r+++ + k+++ s+++ar ++ ++ s lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274 539 SRSVSVAQSLLGDGV-QNYVQELNADYEKVRTQHAN-KKQTPMWSLEKARANRTPVNFS- 595 *************96.89****************98.7789999**********99999. PP TIGR02082 908 dlevpapkflGtkvleas.ieellkyiDwkalFvqWelrgkypkilkdeleglearklfk 966 + +p+p+ lG++v+++ ++el kyiDw ++F +W+l+g yp il+de++g+ea+++f+ lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274 596 TYVPPTPRLLGRRVFKNFdLAELAKYIDWGPFFQTWDLAGPYPAILTDEVVGVEATRVFA 655 9*********************************************************** PP TIGR02082 967 dakelldklsaekllrargvvGlfPaqsvgddieiytdetvsqetkpiatvrekleqlrq 1026 d +++l+kl++ + l a+gv+ l+Pa+svgddie ytdet++ ++ t + ++q+++ lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274 656 DGQAMLKKLIEGRWLSASGVMALYPANSVGDDIEFYTDETRT---QVAMTWYGMRQQTEK 712 ****************************************94...445555555555555 PP TIGR02082 1027 qsdr.....ylclaDfiaskesGikDylgallvtaglgaeelakkleakeddydsilvka 1081 q claDf+a+k+sGi+Dy g+++vtaglg+e+ +++ + ddy++i++k+ lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274 713 QMIDgvmrpSRCLADFVAPKDSGIADYAGLFAVTAGLGVEKKEQAFVDALDDYSAIMLKS 772 54444667899************************************************* PP TIGR02082 1082 ladrlaealaellhervRkelwgyaeeenldkedllkerYrGirpafGYpacPdhtekat 1141 ladrlaea+ae lh+rvR +lwgya+ e+l++++++ e+Y Girpa+GYpacPdh+ k++ lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274 773 LADRLAEAFAECLHQRVRTDLWGYAAGEQLSNDEMIAEKYAGIRPAPGYPACPDHSAKQE 832 ************************************************************ PP TIGR02082 1142 lleLleaeriGlklteslalaPeasvsglyfahpeakYfav 1182 l+++l+ e+iG+ ltesla++P+asvsg+y+ hp++ Yf+v lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274 833 LFRVLQCEEIGMGLTESLAMTPAASVSGFYIGHPDSVYFNV 873 ***************************************98 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (1182 nodes) Target sequences: 1 (903 residues searched) Passed MSV filter: 1 (1); expected 0.0 (0.02) Passed bias filter: 1 (1); expected 0.0 (0.02) Passed Vit filter: 1 (1); expected 0.0 (0.001) Passed Fwd filter: 1 (1); expected 0.0 (1e-05) Initial search space (Z): 1 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.08u 0.03s 00:00:00.11 Elapsed: 00:00:00.09 # Mc/sec: 11.27 // [ok]
This GapMind analysis is from Apr 09 2024. The underlying query database was built on Apr 09 2024.
Each pathway is defined by a set of rules based on individual steps or genes. Candidates for each step are identified by using ublast (a fast alternative to protein BLAST) against a database of manually-curated proteins (most of which are experimentally characterized) or by using HMMer with enzyme models (usually from TIGRFam). Ublast hits may be split across two different proteins.
A candidate for a step is "high confidence" if either:
Otherwise, a candidate is "medium confidence" if either:
Other blast hits with at least 50% coverage are "low confidence."
Steps with no high- or medium-confidence candidates may be considered "gaps." For the typical bacterium that can make all 20 amino acids, there are 1-2 gaps in amino acid biosynthesis pathways. For diverse bacteria and archaea that can utilize a carbon source, there is a complete high-confidence catabolic pathway (including a transporter) just 38% of the time, and there is a complete medium-confidence pathway 63% of the time. Gaps may be due to:
GapMind relies on the predicted proteins in the genome and does not search the six-frame translation. In most cases, you can search the six-frame translation by clicking on links to Curated BLAST for each step definition (in the per-step page).
For more information, see:
If you notice any errors or omissions in the step descriptions, or any questionable results, please let us know
by Morgan Price, Arkin group, Lawrence Berkeley National Laboratory