Align cobalamin-dependent methionine synthase (EC 2.1.1.13) (characterized)
to candidate WP_013074393.1 BTUS_RS01655 methionine synthase
Query= metacyc::G18NG-11090-MONOMER (1221 letters) >NCBI__GCF_000092905.1:WP_013074393.1 Length = 1159 Score = 756 bits (1953), Expect = 0.0 Identities = 489/1214 (40%), Positives = 668/1214 (55%), Gaps = 99/1214 (8%) Query: 22 LANHVLIGDGAMGTQLQGFDLDVEKDFLD--LEGCNEILNDTRPDVLRQIHRAYFEAGAD 79 L +L+ DGAMGT +Q L DF EGCNE+LN TRPDV+ IH AY EAGAD Sbjct: 24 LKEAILVLDGAMGTMIQQQGLTAA-DFGGEAQEGCNEVLNLTRPDVITSIHEAYLEAGAD 82 Query: 80 LVETNTFGCNLPNLADYDIADRCRELAYKGTAVAREVADEMGPGRNGMRRFVVGSLGPGT 139 ++ETNTFG LA+Y + + + +AR AD G RFV G++GP T Sbjct: 83 IIETNTFGATRVVLAEYGLEGKAYAINRTAAELARRAAD--GYSTPDRPRFVAGAMGPTT 140 Query: 140 KLPSL-GHAPYADLRGHYKEAALGIIDGGGDAFLIETAQDLLQVKAAVHGVQDAMAELDT 198 K S+ G + L Y E A G++DGG D L+ETAQD L VKAA G+ A AEL Sbjct: 141 KTLSVTGGIDFDTLAAAYAEQAKGLLDGGVDLLLVETAQDTLNVKAAGVGIHRAFAELGF 200 Query: 199 FLPIICHVTVETTGTMLMGSEIGAALTALQPLGIDMIGLNCATGPDEMSEHLRYLSKHAD 258 +P++ T+E GT L G I A +++ L +GLNCATGP+ M++HLR LS A Sbjct: 201 EVPVMISGTIEPMGTTLAGQGIDAFYVSVEHLHPVAVGLNCATGPEFMTDHLRTLSGLAK 260 Query: 259 IPVSVMPNAGLPVLGKNGAEYPLEAEDLAQALAGFVSEYGLSMVGGCCGTTPEHIRAVRD 318 VS PNAGLP Y + LA+ +AGF + +++VGGCCGTTPEHIRA+ Sbjct: 261 TAVSCYPNAGLP---DEEGHYHESPQSLAEKMAGFARQGWVNVVGGCCGTTPEHIRALAQ 317 Query: 319 AVVGVPEQETSTLTKIPAGPVEQASREVEKEDSVASLYTSVPLSQETGISMIGERTNSNG 378 AV G+P + PA A VE +V L + ++GERTN G Sbjct: 318 AVRGLPP-------RTPAADHPPAVSGVE----------AVFLEGDNRPLLVGERTNVIG 360 Query: 379 SKAFREAMLSGDWEKCVDIAKQQTRDGAHMLDLCVDYVGRDGTADMATLAALLATSSTLP 438 S+ FRE + SG +E+ +IA+ Q R GA ++D+C+ RD ADM + +P Sbjct: 361 SRKFRELIQSGRFEEAAEIARAQVRGGAQIIDVCLADPDRDEYADMDQFLKYVTRMVKVP 420 Query: 439 IMIDSTEPEVIRTGLEHLGGRSIVNSVNFEDGDGPESRYQRIMKLVKQHGAAVVALTIDE 498 +MIDST+P V+ L+ G++I+NS+N EDG E R+ ++ L +++GAAVV TIDE Sbjct: 421 LMIDSTDPRVVELALQRSQGKAIINSINLEDG---EKRFAQVAPLARRYGAAVVVGTIDE 477 Query: 499 EGQARTAEHKVRIAKRLIDDITGSYGLDIKDIVVDCLTFPISTGQEETRRDGIETIEAIR 558 EG A T E K+ +A+R +T YG+ +D++ D L FP+ TG E ET+ A+R Sbjct: 478 EGMAVTRERKLAVAQRSYQLLTEEYGIPPEDLIFDPLVFPVGTGDENYVGAAAETVAALR 537 Query: 559 ELKKLYPEIHTTLGLSNISFGLNPAARQVLNSVFLNECIEAGLDSAIAHSSKILPMNRID 618 LK+ +P T LG+SN+SFGL PA R+VLN+VFL C +AGLD AI + K+ I Sbjct: 538 MLKERFPRSATILGISNVSFGLPPAGREVLNAVFLYHCTKAGLDYAIVNPEKLQRYASIP 597 Query: 619 DRQREVALDMVYDRRTEDYDPLQEFMQLFEGVSAADAKDARAEQLAAMPLFERLAQRIID 678 + +R +A ++++ E L F F+ A R E L+ L ERLA+ +++ Sbjct: 598 EEERRLAEALLFETGQE---TLARFHDRFKDRPA--GPPVRRETLS---LNERLARCVVE 649 Query: 679 GDKNGLEDDLEAGMKEKSPIAIINEDLLNGMKTVGELFGSGQMQLPFVLQSAETMKTAVA 738 G K GL DDL ++ SPI IIN L+ GM VG LF Q+ + VLQSAE MK AVA Sbjct: 650 GSKEGLLDDLSEALQSASPIEIINGPLMAGMDEVGRLFNDNQLIVAEVLQSAEVMKAAVA 709 Query: 739 YLEPFMEEEAEATGSAQAEGKGKIVVATVKGDVHDIGKNLVDIILSNNGYDVVNLGIKQP 798 YLEP ME+ A+ K KI++ATVKGDVHDIGKNLV+IILSNNG++VVNLGIK Sbjct: 710 YLEPHMEK-------AETAKKAKILLATVKGDVHDIGKNLVEIILSNNGFEVVNLGIKVA 762 Query: 799 LSAMLEAAEEHKADVIGMSGLLVKSTVVMKENLEEMNNAGASNYPVILGGAALTRTYVEN 858 ++EA + D IG+SGLLVKS M +++ AG S P+++GGAALT + Sbjct: 763 PEVLIEAVRREEPDAIGLSGLLVKSAQQMVVTAQDLRAAGVS-LPLLVGGAALTAKFAYT 821 Query: 859 DLNEVYTGEVYYARDAFEGLRLMDEVMAEKRGEGLDPNSPEAIEQAKKKAERKARNERSR 918 + Y G V YA+DA EGL + + + +DP + + + ERK R Sbjct: 822 RIAPEYEGLVLYAKDAMEGLEIANRL--------VDP-------EGRPELERKTAELR-- 864 Query: 919 KIAAERKANAAPVIV----PERSDVSTDTPTAAPPFWGTRIVKGLPLAEFLGNLDERALF 974 A+ AAP + P RS VS P A PP + + P ++ + L Sbjct: 865 --ASLTDRPAAPAVSTGSRPRRSSVSRTVPVAVPPDCRRHLWRDYPFDHVQPYINLQMLL 922 Query: 975 MGQWGLKSTRGNEGPSYEDLVETEGRPRLRYWLDRL-------KSEGILDHVALVYGYFP 1027 GL RGN D + E P+ W++RL K EG+L + VY +FP Sbjct: 923 GRHMGL---RGN-----VDKLLAERDPKAVEWVERLNALIAEAKGEGLL-RLHAVYQFFP 973 Query: 1028 AVAEGDDVVILESPDPHAAERMRFSFPRQQRGRFLCIADFIRPREQAVKDGQVDVMPFQL 1087 + +EGDD++I + P RF+FPRQ R +LC++D++RP V+ G++D + F Sbjct: 974 SQSEGDDILIYDPEQPGRVIE-RFTFPRQGREPYLCLSDYVRP----VESGEMDYVAFFA 1028 Query: 1088 VTMGNPIADFANELFAANEYREYLEVHGIGVQLTEALAEYWHSRVRSELKLNDGGSVADF 1147 VT G I D A E A EY + + + ++ EA AE H +R Sbjct: 1029 VTAGAGIRDKAEEWKARGEYVKSHFLQALALETAEAFAERLHGMIRDLWGFP-------- 1080 Query: 1148 DPEDKT--KFFDLDYRGARFSFGYGSCPDLEDRAKLVELLEPGRIGVELSEELQLHPEQS 1205 DP D T + F YRG R SFGY +CP+LED+AKL LL P IG+ L++ + PE S Sbjct: 1081 DPPDMTMRERFQAKYRGIRVSFGYPACPNLEDQAKLFRLLRPEEIGIRLTDGFMMDPEAS 1140 Query: 1206 TDAFVLYHPEAKYF 1219 A V HPEA YF Sbjct: 1141 VSALVFSHPEAVYF 1154 Lambda K H 0.316 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: 3147 Number of extensions: 153 Number of successful extensions: 14 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: 1221 Length of database: 1159 Length adjustment: 47 Effective length of query: 1174 Effective length of database: 1112 Effective search space: 1305488 Effective search space used: 1305488 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.6 bits) S2: 59 (27.3 bits)
Align candidate WP_013074393.1 BTUS_RS01655 (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.28638.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 1262.9 0.0 0 1262.7 0.0 1.0 1 lcl|NCBI__GCF_000092905.1:WP_013074393.1 BTUS_RS01655 methionine synthase Domain annotation for each sequence (and alignments): >> lcl|NCBI__GCF_000092905.1:WP_013074393.1 BTUS_RS01655 methionine synthase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 1262.7 0.0 0 0 2 1181 .. 25 1155 .. 24 1156 .. 0.95 Alignments for each domain: == domain 1 score: 1262.7 bits; conditional E-value: 0 TIGR02082 2 nkrilvlDGamGtqlqsanLteadFrgeeadlarelkGnndlLnltkPeviaaihrayfeaGaDive 68 ++ ilvlDGamGt++q++ Lt+adF ge ++G+n++Lnlt+P+vi++ih+ay+eaGaDi+e lcl|NCBI__GCF_000092905.1:WP_013074393.1 25 KEAILVLDGAMGTMIQQQGLTAADFGGE------AQEGCNEVLNLTRPDVITSIHEAYLEAGADIIE 85 688************************7......89******************************* PP TIGR02082 69 tntFnsteialadYdledkayelnkkaaklarevadeftltpekkRfvaGslGPtnklatlspdver 135 tntF++t ++la+Y+le kay++n++aa+lar++ad ++ tp+++RfvaG++GPt+k+ +++ lcl|NCBI__GCF_000092905.1:WP_013074393.1 86 TNTFGATRVVLAEYGLEGKAYAINRTAAELARRAADGYS-TPDRPRFVAGAMGPTTKTLSVT----- 146 ***************************************.*********************9..... PP TIGR02082 136 pefrnvtydelvdaYkeqvkglldGGvDllLietvfDtlnakaalfaveevfeekgrelPilisgvi 202 +++d+l +aY eq+kglldGGvDllL+et +Dtln+kaa ++++++f+e g e+P++isg+i lcl|NCBI__GCF_000092905.1:WP_013074393.1 147 ---GGIDFDTLAAAYAEQAKGLLDGGVDLLLVETAQDTLNVKAAGVGIHRAFAELGFEVPVMISGTI 210 ...78***********************************************************875 PP TIGR02082 203 vdksGrtLsGqtleaflaslehaeililGLnCalGadelrefvkelsetaealvsviPnaGLPnalg 269 + +G+tL Gq ++af++s+eh + +++GLnCa+G++ + +++++ls a++ vs++PnaGLP++ g lcl|NCBI__GCF_000092905.1:WP_013074393.1 211 -EPMGTTLAGQGIDAFYVSVEHLHPVAVGLNCATGPEFMTDHLRTLSGLAKTAVSCYPNAGLPDEEG 276 .99**************************************************************** PP TIGR02082 270 eYdltpeelakalkefaeegllnivGGCCGttPehiraiaeavkdikprkrqeleeksvlsglealk 336 +Y++ p++la++++ fa++g +n+vGGCCGttPehira+a+av+++ pr +++ ++sg+ea+ lcl|NCBI__GCF_000092905.1:WP_013074393.1 277 HYHESPQSLAEKMAGFARQGWVNVVGGCCGTTPEHIRALAQAVRGLPPRTPA-ADHPPAVSGVEAVF 342 **********************************************999765.5666788******* PP TIGR02082 337 iaqessfvniGeRtnvaGskkfrklikaedyeealkiakqqveeGaqilDinvDevllDgeadmkkl 403 ++ ++ + +GeRtnv Gs+kfr+li+ + +eea +ia+ qv+ Gaqi+D+++ + +D adm ++ lcl|NCBI__GCF_000092905.1:WP_013074393.1 343 LEGDNRPLLVGERTNVIGSRKFRELIQSGRFEEAAEIARAQVRGGAQIIDVCLADPDRDEYADMDQF 409 ******************************************************************* PP TIGR02082 404 lsllasepdiakvPlmlDssefevleaGLkviqGkaivnsislkdGeerFlekaklikeyGaavvvm 470 l+ ++ + kvPlm+Ds++ v+e L++ qGkai+nsi+l+dGe+rF++ a l+++yGaavvv lcl|NCBI__GCF_000092905.1:WP_013074393.1 410 LKYVTRM--V-KVPLMIDSTDPRVVELALQRSQGKAIINSINLEDGEKRFAQVAPLARRYGAAVVVG 473 **99988..5.******************************************************** PP TIGR02082 471 afDeeGqartadkkieiakRayklltekvgfppediifDpniltiatGieehdryaidfieaireik 537 ++DeeG+a t+++k+ +a+R y+llte++g+pped+ifDp++++++tG e+ a +++ a+r +k lcl|NCBI__GCF_000092905.1:WP_013074393.1 474 TIDEEGMAVTRERKLAVAQRSYQLLTEEYGIPPEDLIFDPLVFPVGTGDENYVGAAAETVAALRMLK 540 ******************************************************************* PP TIGR02082 538 eelPdakisgGvsnvsFslrgndavRealhsvFLyeaikaGlDmgivnagklavyddidkelrevve 604 e++P + +G+snvsF+l +a Re+l++vFLy++ kaGlD +ivn+ kl+ y+ i++e r ++e lcl|NCBI__GCF_000092905.1:WP_013074393.1 541 ERFPRSATILGISNVSFGLP--PAGREVLNAVFLYHCTKAGLDYAIVNPEKLQRYASIPEEERRLAE 605 ********************..********************************************* PP TIGR02082 605 dlildrrreatekLlelaelykgtkeksskeaqeaewrnlpveeRLeralvkGeregieedleeark 671 l+++ +e L ++ +++k+ + l+++eRL+r +v+G +eg+ +dl ea+ lcl|NCBI__GCF_000092905.1:WP_013074393.1 606 ALLFET---GQETLARFHDRFKDRPAGPP-----VRRETLSLNERLARCVVEGSKEGLLDDLSEAL- 663 ****85...67889*********999833.....3456799*************************. PP TIGR02082 672 klkapleiiegpLldGmkvvGdLFGsGkmfLPqvvksarvmkkavayLePylekekeedkskGkivl 738 + +p+eii+gpL++Gm+ vG LF + ++ + v++sa+vmk+avayLeP++ek + + +k ki+l lcl|NCBI__GCF_000092905.1:WP_013074393.1 664 QSASPIEIINGPLMAGMDEVGRLFNDNQLIVAEVLQSAEVMKAAVAYLEPHMEKAE--TAKKAKILL 728 899*************************************************9876..788999*** PP TIGR02082 739 atvkGDvhDiGknivdvvLscngyevvdlGvkvPvekileaakkkkaDviglsGLivksldemveva 805 atvkGDvhDiGkn+v+++Ls+ng+evv+lG+kv e ++ea +++++D iglsGL+vks ++mv a lcl|NCBI__GCF_000092905.1:WP_013074393.1 729 ATVKGDVHDIGKNLVEIILSNNGFEVVNLGIKVAPEVLIEAVRREEPDAIGLSGLLVKSAQQMVVTA 795 ******************************************************************* PP TIGR02082 806 eemerrgvkiPlllGGaalskahvavkiaekYkgevvyvkdaseavkvvdkllsekkkaeelekike 872 ++++ gv++Pll+GGaal+ ++ + ia++Y+g v+y+kda+e+++++++l++ + +e +++ lcl|NCBI__GCF_000092905.1:WP_013074393.1 796 QDLRAAGVSLPLLVGGAALTAKFAYTRIAPEYEGLVLYAKDAMEGLEIANRLVDPEG-RPE---LER 858 *****************************************************9995.444...344 PP TIGR02082 873 eyeeirekfgekkeklialsekaarkevfaldrsedlevpapkflGtkvleas.ieellkyiDwkal 938 + e+r + + + a+s+ r+ + ++ r+ + v++p++ ++ +++ +++ +yi++++l lcl|NCBI__GCF_000092905.1:WP_013074393.1 859 KTAELRASLTD-RPAAPAVST-GSRPRRSSVSRT--VPVAVPPDCRRHLWRDYpFDHVQPYINLQML 921 44556655554.444555554.445566666666..88999************************** PP TIGR02082 939 Fv.qWelrgkypkilkdeleglearklfkdakelldklsaekllrargvvGlfPaqsvgddieiytd 1004 + +lrg+ +k l + ++ +a + +++++l+ + + e+llr ++v+ +fP+qs+gddi iy++ lcl|NCBI__GCF_000092905.1:WP_013074393.1 922 LGrHMGLRGNVDKLLAE--RDPKAVEWVERLNALIAEAKGEGLLRLHAVYQFFPSQSEGDDILIYDP 986 99899*******98887..567788999**************************************6 PP TIGR02082 1005 etvsqetkpiatvrekleqlrqqsdrylclaDfiaskesGikDylgallvtaglgaeelakkleake 1071 e+ v+e++ ++rq ++ ylcl+D++++ esG+ Dy+++++vtag g+++ a++++a+ lcl|NCBI__GCF_000092905.1:WP_013074393.1 987 EQP-------GRVIERFTFPRQGREPYLCLSDYVRPVESGEMDYVAFFAVTAGAGIRDKAEEWKARG 1046 654.......6677899************************************************** PP TIGR02082 1072 ddydsilvkaladrlaealaellhervRkelwgyaeeenldkedllkerYrGirpafGYpacPdhte 1138 ++ +s ++ala aea+ae lh +R lwg + + +++ ++ +YrGir +fGYpacP++++ lcl|NCBI__GCF_000092905.1:WP_013074393.1 1047 EYVKSHFLQALALETAEAFAERLHGMIRD-LWGFPDPPDMTMRERFQAKYRGIRVSFGYPACPNLED 1112 ***************************96.6777688888999************************ PP TIGR02082 1139 katlleLleaeriGlklteslalaPeasvsglyfahpeakYfa 1181 +a+l++Ll +e+iG++lt+ +++ Peasvs+l+f+hpea Yf lcl|NCBI__GCF_000092905.1:WP_013074393.1 1113 QAKLFRLLRPEEIGIRLTDGFMMDPEASVSALVFSHPEAVYFR 1155 ******************************************6 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (1182 nodes) Target sequences: 1 (1159 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.04s 00:00:00.12 Elapsed: 00:00:00.12 # Mc/sec: 11.04 // [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