Align Methionine synthase; EC 2.1.1.13; 5-methyltetrahydrofolate--homocysteine methyltransferase; Methionine synthase, vitamin-B12 dependent; MS (uncharacterized)
to candidate WP_012470188.1 GLOV_RS10575 methionine synthase
Query= curated2:O33259 (1192 letters) >NCBI__GCF_000020385.1:WP_012470188.1 Length = 1169 Score = 771 bits (1992), Expect = 0.0 Identities = 459/1193 (38%), Positives = 691/1193 (57%), Gaps = 60/1193 (5%) Query: 13 LLDVLSQRVMVGDGAMGTQLQAADLTLDDFRG--LEGCNEILNETRPDVLETIHRNYFEA 70 +L LSQ +++ DGAMGTQLQA +LT DF G EGCNE+L TRPDV+E +H+ Y EA Sbjct: 6 VLQALSQCILILDGAMGTQLQARNLTAADFGGEAYEGCNEMLVLTRPDVIEDVHKAYLEA 65 Query: 71 GADAVETNTFGCNLSNLGDYDIADRIRDLSQKGTAIARRVADELGSPDRKRYVLGSMGPG 130 GAD VET +FG L +Y +AD++ +L++ +AR+ D +PD+ R+V GSMGP Sbjct: 66 GADIVETCSFGSTDIVLAEYGLADKVFELNKAAALVARKACDAYSTPDKPRFVAGSMGPT 125 Query: 131 TK-LPTLGHTEYAVIRDAYTEAALGMLDGGADAILVETCQDLLQLKAAVLGSRRAMTRAG 189 T+ + G + + A+ + +G+L GG D +L+ET QD L LKA G R A + G Sbjct: 126 TRTISVTGGVTFEQLVTAFRDQTIGLLAGGVDLLLLETAQDTLNLKAGAEGIRLAFEQTG 185 Query: 190 RHIPVFAHVTVETTGTMLLGSEIGAALTAV----EPLGVDMIGLNCATGPAEMSEHLRHL 245 + +P+ T+E TGTML G + A ++ E LG+ IGLNCATGP M++HLR L Sbjct: 186 QRVPLMVSGTIEPTGTMLAGQNVEALYASLTHLEENLGLISIGLNCATGPEFMTDHLRAL 245 Query: 246 SRHARIPVSVMPNAGLPVLGAKGAEYPLLPDELAEALAGFIAEFGLSLVGGCCGTTPAHI 305 S A +SV PNAGLP + Y PD LA+ L+ F+ E L+++GGCCGT+PAHI Sbjct: 246 SELATCHISVYPNAGLP---DENGNYAESPDSLAQKLSRFVDEGWLNIIGGCCGTSPAHI 302 Query: 306 REVAAAVANIKRPERQVSYEPSVSSLYTAIPFAQDASVLVIGERTNANGSKGFREAMIAE 365 +A A ++P + V+ + S + +D +++GERTN GS+ F+ ++A Sbjct: 303 AAIANMAAG-RQPRKPVTIRRRLVSGIEPLFIEEDNRPILVGERTNVIGSRKFKNMIVAG 361 Query: 366 DYQKCLDIAKDQTRDGAHLLDLCVDYVGRDGVADMKALASRLATSSTLPIMLDSTETAVL 425 +++ +IA+ Q + GA ++D+CV RD +ADM+ + L PIM+DST+ V Sbjct: 362 QFEEAAEIARAQVKTGAQVIDICVANPDRDELADMEQMLVYLPKKVRAPIMIDSTDAKVT 421 Query: 426 QAGLEHLGGRCAINSVNYEDGDGPESRFAKTMALVAEHGAAVVALTIDEE---GQARTAQ 482 + L+ L G+C INS+N EDG E RFAK L+ +G AVV IDE+ G A T Q Sbjct: 422 ELALQRLQGKCVINSINLEDG---EERFAKVAPLLRHYGGAVVVGCIDEDPANGMAVTRQ 478 Query: 483 KKVEIAERLINDITGNWGVDESSILIDTLTFTIATGQEESRRDGIETIEAIRELKKRHPD 542 +K+E+A+R + + +G+ ++ D L F + +G + +ETIE +R + + P Sbjct: 479 RKLEVAQRSYDLLVNKYGLRPEDLIFDPLVFPVGSGDDNYIGSAVETIEGVRLITETFPQ 538 Query: 543 VQTTLGLSNISFGLNPAARQVLNSVFLHECQEAGLDSAIVHASKILPMNRIPEEQRNVAL 602 T LG+SN+SFGL A R+VLN+VFL+ C +AGL AIV+ K+ IPE +R +A Sbjct: 539 CSTILGISNVSFGLPAAGREVLNAVFLYHCVKAGLGYAIVNTEKLERYASIPEAERRLAE 598 Query: 603 DLVYDRRREDYDPLQELMRLFEGVSAASSKEDRLAELAGLPLFERLAQRIVDGERNGLDA 662 DL++ R DP+ F S A A LPL ERL I++G ++GL A Sbjct: 599 DLIFWR---GADPVAAFAAAFRDKKPVSH-----APAAELPLDERLPLYIIEGSKDGLTA 650 Query: 663 DLDEAMTQ-KPPLQIINEHLLAGMKTVGELFGSGQMQLPFVLQSAEVMKAAVAYLEPHME 721 DLD A+T+ PL+IIN L+AGM VG LF Q+ + VLQSAE MKAAV++LEPH+ Sbjct: 651 DLDAALTRGDKPLEIINGPLMAGMAEVGRLFNDNQLIVAEVLQSAEAMKAAVSHLEPHLS 710 Query: 722 RSDDDSGKGRIVLATVKGDVHDIGKNLVDIILSNNGYEVVNIGIKQPIATILEVAEDKSA 781 + D+ S KG+++LATVKGDVHDIGKNLV+IILSNNG+EV+N+GIK ++ A+ ++ Sbjct: 711 K-DETSSKGKMLLATVKGDVHDIGKNLVEIILSNNGFEVINLGIKVGPEELIAAAKKENP 769 Query: 782 DVVGMSGLLVKSTVVMKENLEEMNTRGVAEKFPVLLGGAALTRSYVENDLAEIYQGEVHY 841 D +G+SGLLVKS + M ++ G+ P+L+GGAAL+R++ + + Y G V Y Sbjct: 770 DFIGLSGLLVKSALQMIVTAADLKAAGI--DAPLLVGGAALSRAFADTRITPEYNGPVLY 827 Query: 842 ARDAFEGLKLMDTIMS-AKRGEAPDENSPEAIKAREKEAERKARHQRSKRIAAQRKAAEE 900 A+DA GL+L + ++ A R + + AR++EA S +IAA + A + Sbjct: 828 AKDAMAGLELANQLVDPALRQQLMLD------LARQQEA--------SAKIAAAKAAGQS 873 Query: 901 PVEV-PERSDVAADIEVPAPPFWGSRIVKGLAVADYTGLLDERALFLGQWGLRGQRGGEG 959 V S ++++ + P I++ + V L+ + L+ GL G + Sbjct: 874 TVATGSTHSAISSNAPILVAPDLEQHILRDIPVGQIIPYLNRQMLYTKHLGLTGS-VDKL 932 Query: 960 PSYEDLVETEGRPRLRYWLDRLSTDGILAHAAVVYGYFPAVSEGNDIVVLTEPKPDAPVR 1019 + +D T+ + L+R+ +G++ A +Y ++ A +GND+++ + +A Sbjct: 933 LAGQDEKATKLHLTVEAMLERVLQEGLIKPQA-IYRFYQANGDGNDLILFNQDGSEAT-- 989 Query: 1020 YRFHFPRQQRGRFLCIADFIRSRELAAERGEVDVLPFQLVTMGQPIADFANELFASNAYR 1079 RF PRQ+ G LC+ADF+R E D + VT GQ + + + + A Y Sbjct: 990 -RFTLPRQKSGEQLCVADFVR----PLSGTEKDTMALFAVTCGQGVRELSEQWKAEGDYL 1044 Query: 1080 DYLEVHGIGVQLTEALAEYWHRRIREELKFSGDRAMAAEDPEAKEDYFKLGYRGARFAFG 1139 + + + +++ EA AEY H+RIR D+ + + F Y+G R +FG Sbjct: 1045 NSHLLQALALEMAEATAEYLHKRIRTSWGIVDDQTL------TMKQLFNAEYQGIRVSFG 1098 Query: 1140 YGACPDLEDRAKMMALLEPERIGVTLSEELQLHPEQSTDAFVLHHPEAKYFNV 1192 Y ACP+L+D+ K+ +LL+PE+IG+ L+E + PE S A V HHPE KYF++ Sbjct: 1099 YPACPNLDDQKKLFSLLKPEQIGINLTEGDMMDPEASVSALVFHHPEGKYFDL 1151 Lambda K H 0.317 0.135 0.386 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: 2978 Number of extensions: 121 Number of successful extensions: 17 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: 1192 Length of database: 1169 Length adjustment: 47 Effective length of query: 1145 Effective length of database: 1122 Effective search space: 1284690 Effective search space used: 1284690 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 candidate WP_012470188.1 GLOV_RS10575 (methionine synthase)
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.5930.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 1227.0 0.0 0 1226.8 0.0 1.0 1 lcl|NCBI__GCF_000020385.1:WP_012470188.1 GLOV_RS10575 methionine synthase Domain annotation for each sequence (and alignments): >> lcl|NCBI__GCF_000020385.1:WP_012470188.1 GLOV_RS10575 methionine synthase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 1226.8 0.0 0 0 2 1181 .. 11 1150 .. 10 1151 .. 0.95 Alignments for each domain: == domain 1 score: 1226.8 bits; conditional E-value: 0 TIGR02082 2 nkrilvlDGamGtqlqsanLteadFrgeeadlarelkGnndlLnltkPeviaaihrayfeaGaDive 68 ++ il+lDGamGtqlq+ nLt+adF ge ++G+n++L+lt+P+vi+++h+ay+eaGaDive lcl|NCBI__GCF_000020385.1:WP_012470188.1 11 SQCILILDGAMGTQLQARNLTAADFGGE------AYEGCNEMLVLTRPDVIEDVHKAYLEAGADIVE 71 788************************7......9******************************** PP TIGR02082 69 tntFnsteialadYdledkayelnkkaaklarevadeftltpekkRfvaGslGPtnklatlspdver 135 t F+st i+la+Y+l+dk+ elnk+aa +ar+++d ++ tp+k+RfvaGs+GPt+++ +++ lcl|NCBI__GCF_000020385.1:WP_012470188.1 72 TCSFGSTDIVLAEYGLADKVFELNKAAALVARKACDAYS-TPDKPRFVAGSMGPTTRTISVT----- 132 ***************************************.*********************9..... PP TIGR02082 136 pefrnvtydelvdaYkeqvkglldGGvDllLietvfDtlnakaalfaveevfeekgrelPilisgvi 202 vt+++lv a+++q+ gll GGvDllL+et +Dtln+ka +++ +fe+ g+++P+++sg+i lcl|NCBI__GCF_000020385.1:WP_012470188.1 133 ---GGVTFEQLVTAFRDQTIGLLAGGVDLLLLETAQDTLNLKAGAEGIRLAFEQTGQRVPLMVSGTI 196 ...789**********************************************************876 PP TIGR02082 203 vdksGrtLsGqtleaflaslehae....ililGLnCalGadelrefvkelsetaealvsviPnaGLP 265 + +G++L Gq++ea +asl h e ++++GLnCa+G++ + +++++lse a++ +sv+PnaGLP lcl|NCBI__GCF_000020385.1:WP_012470188.1 197 -EPTGTMLAGQNVEALYASLTHLEenlgLISIGLNCATGPEFMTDHLRALSELATCHISVYPNAGLP 262 .899***************987542222899************************************ PP TIGR02082 266 nalgeYdltpeelakalkefaeegllnivGGCCGttPehiraiaeavkdikprkrqeleeksvlsgl 332 ++ g+Y + p++la++l f+ eg lni+GGCCGt P+hi+aia+++++ +prk ++ ++sg+ lcl|NCBI__GCF_000020385.1:WP_012470188.1 263 DENGNYAESPDSLAQKLSRFVDEGWLNIIGGCCGTSPAHIAAIANMAAGRQPRKPVTIRR-RLVSGI 328 *****************************************************9887765.5689** PP TIGR02082 333 ealkiaqessfvniGeRtnvaGskkfrklikaedyeealkiakqqveeGaqilDinvDevllDgead 399 e+l i++++ + +GeRtnv Gs+kf+++i a+++eea +ia+ qv+ Gaq++Di+v + +D ad lcl|NCBI__GCF_000020385.1:WP_012470188.1 329 EPLFIEEDNRPILVGERTNVIGSRKFKNMIVAGQFEEAAEIARAQVKTGAQVIDICVANPDRDELAD 395 ******************************************************************* PP TIGR02082 400 mkkllsllasepdiakvPlmlDssefevleaGLkviqGkaivnsislkdGeerFlekaklikeyGaa 466 m+++l l ++ + + P+m+Ds++ +v e L+++qGk+++nsi+l+dGeerF++ a l+++yG a lcl|NCBI__GCF_000020385.1:WP_012470188.1 396 MEQMLVYLPKK--V-RAPIMIDSTDAKVTELALQRLQGKCVINSINLEDGEERFAKVAPLLRHYGGA 459 ****9999888..5.9*************************************************** PP TIGR02082 467 vvvmafDee...GqartadkkieiakRayklltekvgfppediifDpniltiatGieehdryaidfi 530 vvv +De+ G+a t+++k+e+a+R y+ll++k+g+ ped+ifDp++++++ G ++ a+++i lcl|NCBI__GCF_000020385.1:WP_012470188.1 460 VVVGCIDEDpanGMAVTRQRKLEVAQRSYDLLVNKYGLRPEDLIFDPLVFPVGSGDDNYIGSAVETI 526 ********6666******************************************************* PP TIGR02082 531 eaireikeelPdakisgGvsnvsFslrgndavRealhsvFLyeaikaGlDmgivnagklavyddidk 597 e +r i e++P++ + +G+snvsF+l a Re+l++vFLy+++kaGl +ivn+ kl+ y+ i++ lcl|NCBI__GCF_000020385.1:WP_012470188.1 527 EGVRLITETFPQCSTILGISNVSFGLP--AAGREVLNAVFLYHCVKAGLGYAIVNTEKLERYASIPE 591 ***************************..999*********************************** PP TIGR02082 598 elrevvedlildrrreatekLlelaelykgtkeksskeaqeaewrnlpveeRLeralvkGeregiee 664 r ++edli+ r +++ +a ++++k +a ++lp++eRL +++G ++g+ + lcl|NCBI__GCF_000020385.1:WP_012470188.1 592 AERRLAEDLIFWR---GADPVAAFAAAFRDKKPV-----SHAPAAELPLDERLPLYIIEGSKDGLTA 650 ***********98...57899********99887.....457899********************** PP TIGR02082 665 dleearkklkapleiiegpLldGmkvvGdLFGsGkmfLPqvvksarvmkkavayLePylekekeedk 731 dl++a+ + +pleii+gpL++Gm vG LF + ++ + v++sa+ mk+av++LeP+l k++ ++ lcl|NCBI__GCF_000020385.1:WP_012470188.1 651 DLDAALTRGDKPLEIINGPLMAGMAEVGRLFNDNQLIVAEVLQSAEAMKAAVSHLEPHLSKDE--TS 715 ***********************************************************9887..99 PP TIGR02082 732 skGkivlatvkGDvhDiGknivdvvLscngyevvdlGvkvPvekileaakkkkaDviglsGLivksl 798 skGk++latvkGDvhDiGkn+v+++Ls+ng+ev++lG+kv e+++ aakk+++D iglsGL+vks lcl|NCBI__GCF_000020385.1:WP_012470188.1 716 SKGKMLLATVKGDVHDIGKNLVEIILSNNGFEVINLGIKVGPEELIAAAKKENPDFIGLSGLLVKSA 782 ******************************************************************* PP TIGR02082 799 demvevaeemerrgvkiPlllGGaalskahvavkiaekYkgevvyvkdaseavkvvdkllsekkkae 865 +m+ a +++ g++ Pll+GGaals+a+ ++ i ++Y+g+v+y+kda+++++++++l++ +++ lcl|NCBI__GCF_000020385.1:WP_012470188.1 783 LQMIVTAADLKAAGIDAPLLVGGAALSRAFADTRITPEYNGPVLYAKDAMAGLELANQLVDPALRQQ 849 *************************************************************998877 PP TIGR02082 866 elekikeeyee.irekfgekkeklialsekaarkevfaldrsedlevpapkflGtkvleas.ieell 930 + + +++e ++ + + + ++ + + s+++ + ++++l +++l+++ + +++ lcl|NCBI__GCF_000020385.1:WP_012470188.1 850 LMLDLARQQEAsAKIAAAK-AAGQSTV-----ATGSTHSAISSNAPILVAPDLEQHILRDIpVGQII 910 7655544433303333322.2222222.....2233344555578899999**************** PP TIGR02082 931 kyiDwkalFv.qWelrgkypkilkdeleglearklfkdakelldklsaekllrargvvGlfPaqsvg 996 +y++ ++l++ + +l g+ +k l ++++a+kl+ +++++l+++ +e+l+++++++ ++ a+ g lcl|NCBI__GCF_000020385.1:WP_012470188.1 911 PYLNRQMLYTkHLGLTGSVDKLLAG--QDEKATKLHLTVEAMLERVLQEGLIKPQAIYRFYQANGDG 975 **********899*******88875..67899*********************************** PP TIGR02082 997 ddieiytdetvsqetkpiatvrekleqlrqqsdrylclaDfiaskesGikDylgallvtaglgaeel 1063 +d+ ++++++ a +++ +rq+s+++lc+aDf+++ + +kD +++++vt g+g++el lcl|NCBI__GCF_000020385.1:WP_012470188.1 976 NDLILFNQDGS------EA---TRFTLPRQKSGEQLCVADFVRPLSGTEKDTMALFAVTCGQGVREL 1033 ******96665......33...35899*****************9999******************* PP TIGR02082 1064 akkleakeddydsilvkaladrlaealaellhervRkelwgyaeeenldkedllkerYrGirpafGY 1130 +++++a+ d+ +s l++ala +aea ae+lh+r+R wg +++ l +++l+ +Y+Gir +fGY lcl|NCBI__GCF_000020385.1:WP_012470188.1 1034 SEQWKAEGDYLNSHLLQALALEMAEATAEYLHKRIRTS-WGIVDDQTLTMKQLFNAEYQGIRVSFGY 1099 ************************************87.666699********************** PP TIGR02082 1131 pacPdhtekatlleLleaeriGlklteslalaPeasvsglyfahpeakYfa 1181 pacP+++++++l++Ll++e+iG++lte ++ Peasvs+l+f+hpe kYf lcl|NCBI__GCF_000020385.1:WP_012470188.1 1100 PACPNLDDQKKLFSLLKPEQIGINLTEGDMMDPEASVSALVFHHPEGKYFD 1150 **************************************************6 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (1182 nodes) Target sequences: 1 (1169 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.07u 0.03s 00:00:00.10 Elapsed: 00:00:00.10 # Mc/sec: 13.57 // [ok]
This GapMind analysis is from Apr 10 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