Align methionine synthase; EC 2.1.1.13 (characterized)
to candidate WP_052368278.1 U743_RS15540 methionine synthase
Query= CharProtDB::CH_090726 (1227 letters) >NCBI__GCF_000733765.1:WP_052368278.1 Length = 1255 Score = 1499 bits (3880), Expect = 0.0 Identities = 755/1246 (60%), Positives = 937/1246 (75%), Gaps = 30/1246 (2%) Query: 6 EQLRAQLNERILVLDGGMGTMIQSYRLNEADFRGERFADWPCDLKGNNDLLVLSKPEVIA 65 E++R L+ERILVLDG MGT IQSY+L EAD+RGERF DW DL GNNDLLVL++PE+I Sbjct: 15 EEIRKHLSERILVLDGAMGTRIQSYKLGEADYRGERFKDWESDLAGNNDLLVLTQPELIK 74 Query: 66 AIHNAYFEAGADIIETNTFNSTTIAMADYQMESLSAEINFAAAKLARACADEWTARTPEK 125 +H AY +AGAD+I TNTFN I+MADY ME L+ E+N AAKLA AD + T E+ Sbjct: 75 EVHGAYLDAGADVIGTNTFNCQFISMADYHMEELAYELNVEAAKLACEAAD--SRATAER 132 Query: 126 PRYVAGVLGPTNRTASISPDVNDPAFRNITFDGLVAAYRESTKALVEGGADLILIETVFD 185 PR+VAG +GPTNRTASISPDVNDP RN +FD LV AY E L++GG DLILIET+FD Sbjct: 133 PRFVAGAVGPTNRTASISPDVNDPGKRNTSFDELVRAYSEQIHGLIDGGVDLILIETIFD 192 Query: 186 TLNAKAAVFAVKTEFEALGVELPIMISGTITDASGRTLSGQTTEAFYNSLRHAEALTFGL 245 TLNAKAA+FA T F ++LPI+ISGTITDASGRTLSGQ TEAF+NS+RHAE L GL Sbjct: 193 TLNAKAAIFAAHTVFAEKEIKLPIIISGTITDASGRTLSGQVTEAFWNSIRHAEPLAVGL 252 Query: 246 NCALGPDELRQYVQELSRIAECYVTAHPNAGLPNAFGEYDLDADTMAKQIREWAQAGFLN 305 NCALG ++R Y++ELS++A+C V+ +PNAGLPN FGEYD + + A+ + E+AQAGF+N Sbjct: 253 NCALGGRDMRPYIEELSKLADCNVSCYPNAGLPNTFGEYDEEPEETAEIVAEFAQAGFIN 312 Query: 306 IVGGCCGTTPQHIAAMSRAVEGLAPRKLPEIPVACRLSGLEPLNIGEDSLFVNVGERTNV 365 VGGCCGT P HI A++ AV + PR P RLSGLEPLN+ +D F+N+GERTN+ Sbjct: 313 FVGGCCGTQPSHIEAIANAVREIKPRVPPAHEPITRLSGLEPLNVTKDKGFINIGERTNI 372 Query: 366 TGSAKFKRLIKEEKYSEALDVARQQVENGAQIIDINMDEGMLDAEAAMVRFLNLIAGEPD 425 TGS KF+ LIK+ Y+ ALDVAR QVENGA I+DINMDEGM+D AM F L+A EPD Sbjct: 373 TGSTKFRNLIKDGDYATALDVARSQVENGAAIMDINMDEGMIDGVEAMRTFCQLVASEPD 432 Query: 426 IARVPIMIDSSKWDVIEKGLKCIQGKGIVNSISMKEGVDAFIHHAKLLRRYGAAVVVMAF 485 I+ +P+M+DSSKW+VIE GLKC+QGK +VNSISMKEG + F+ A L R YGAAVVVMAF Sbjct: 433 ISCIPMMVDSSKWEVIETGLKCLQGKCVVNSISMKEGEEKFLREANLCRHYGAAVVVMAF 492 Query: 486 DEQGQADTRARKIEICRRAYKILTEEVGFPPEDIIFDPNIFAVATGIEEHNNYAQDFIGA 545 DE GQAD R+ EIC RAYK+LTE++ FPPEDIIFDPNIFAVATGIEEH Y DFI A Sbjct: 493 DEDGQADNLERRKEICGRAYKLLTEKLDFPPEDIIFDPNIFAVATGIEEHARYGLDFIEA 552 Query: 546 CEDIKRELPHALISGGVSNVSFSFRGNDPVREAIHAVFLYYAIRNGMDMGIVNAGQLAIY 605 C IK PH ISGGVSN+SFSFRGN+ VREAIH+VFLY+A + GMDMGIVNAGQL +Y Sbjct: 553 CAWIKENCPHVHISGGVSNISFSFRGNNKVREAIHSVFLYHATKAGMDMGIVNAGQLTVY 612 Query: 606 DDLPAELRDAVEDVILNRRDD---GTERLLELAEKYRGSKTDDTANAQQAE---WRSWEV 659 +D+ ELRDA+EDVILNR++D TE+LLELAEKY+G D + A Q + WR + Sbjct: 613 EDVEQELRDAIEDVILNRKEDPQEATEKLLELAEKYKG----DASEADQGDRLAWRENDA 668 Query: 660 NKRLEYSLVKGITEFIEQDTEEARQQAT----RPIEVIEGPLMDGMNVVGDLFGEGKMFL 715 +R++++L+KGI ++ E+DTEE R + + +PIEVIEGPLMDGM VVGDLFG GKMFL Sbjct: 669 RERIKHALIKGIDQYAEEDTEELRAEISANGGKPIEVIEGPLMDGMGVVGDLFGAGKMFL 728 Query: 716 PQVVKSARVMKQAVAYLEPFIEASKEQGK---TNGKMVIATVKGDVHDIGKNIVGVVLQC 772 PQVVKSARVMK+AVAYL PFIEA KE+G+ NG +++ATVKGDVHDIGKNIVGVVLQC Sbjct: 729 PQVVKSARVMKKAVAYLIPFIEAEKEEGEDPSDNGTVIMATVKGDVHDIGKNIVGVVLQC 788 Query: 773 NNYEIVDLGVMVPAEKILRTAKEVNADLIGLSGLITPSLDEMVNVAKEMERQGFTIPLLI 832 NNYE++DLGVM PA+KIL AKE +AD+IGLSGLITPSLDEMV+ AKEMERQG IPLLI Sbjct: 789 NNYEVIDLGVMCPAQKILDAAKEHDADIIGLSGLITPSLDEMVHFAKEMERQGLDIPLLI 848 Query: 833 GGATTSKAHTAVKIEQNYSGPTVYVQNASRTVGVVAALLSDTQRDDFVARTRKEYETVRI 892 GGATTS+AHTAVKIE Y GP V+V++ASR+V VA LL +R + + +YE +R Sbjct: 849 GGATTSRAHTAVKIEPKYHGPVVWVKDASRSVPTVAKLLDKNERGAMLEELKADYEKLRE 908 Query: 893 QHGRKKPRTPPVTLEAARDNDFAFDWQAYTPPVAHRLGVQEVEAS---------IETLRN 943 +H K+ + ++ AR N DW+ Y PPV + +A + LR Sbjct: 909 RHAGKEKADKLLPIDKARANRDVPDWEGYEPPVPQLISEAGDKAEYVKVFKDYPLTELRE 968 Query: 944 YIDWTPFFMTWSLAGKYPRILEDEVVGVEAQRLFKDANDMLDKLSAEKTLNPRGVVGLFP 1003 IDW PFF W + G++P +L + G A++L+++AN MLD++ EK L G++GLFP Sbjct: 969 CIDWGPFFNAWEIKGRFPDVLNNPASGETARKLYEEANAMLDQIIEEKWLTANGIIGLFP 1028 Query: 1004 ANRV-GDDIEIYRDETRTHVINVSHHLRQQTE-KTGFANYCLADFVAPKLSGKADYIGAF 1061 AN+V +DI +Y DE+R+ + H LRQQ+E + G N L+D+VAP +G DY+GAF Sbjct: 1029 ANQVDNEDIAVYTDESRSDELLRLHCLRQQSEHREGVPNRALSDYVAPAETGLKDYVGAF 1088 Query: 1062 AVTGGLEEDALADAFEAQHDDYNKIMVKALADRLAEAFAEYLHERVRKVYWGYAPNENLS 1121 AVT GL + F +DDYN I++++LADRLAEAFAE LH+RVRK +WGY P+E L Sbjct: 1089 AVTAGLGSSEKCEEFRKANDDYNAILLESLADRLAEAFAERLHQRVRKEFWGYVPDEKLD 1148 Query: 1122 NEELIRENYQGIRPAPGYPACPEHTEKATIWELLEVEKHTGMKLTESFAMWPGASVSGWY 1181 N+ LI+E YQGIRPAPGYPACPEHTEK IW+LL+VE +TGM+LTES AMWPGASVSG + Sbjct: 1149 NDALIQEKYQGIRPAPGYPACPEHTEKQAIWDLLDVENNTGMQLTESMAMWPGASVSGLF 1208 Query: 1182 FSHPDSKYYAVAQIQRDQVEDYARRKGMSVTEVERWLAPNLGYDAD 1227 F+HPDS+Y+ + ++ RDQ+ DYA+RKG E ERWL+PNL Y+ + Sbjct: 1209 FAHPDSRYFVLGRVNRDQIADYAKRKGWGQREAERWLSPNLAYEPE 1254 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: 3818 Number of extensions: 163 Number of successful extensions: 10 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: 1255 Length adjustment: 48 Effective length of query: 1179 Effective length of database: 1207 Effective search space: 1423053 Effective search space used: 1423053 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: 59 (27.3 bits)
Align candidate WP_052368278.1 U743_RS15540 (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.362841.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 1721.2 0.0 0 1721.0 0.0 1.0 1 NCBI__GCF_000733765.1:WP_052368278.1 Domain annotation for each sequence (and alignments): >> NCBI__GCF_000733765.1:WP_052368278.1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 1721.0 0.0 0 0 2 1180 .. 22 1217 .. 21 1219 .. 0.97 Alignments for each domain: == domain 1 score: 1721.0 bits; conditional E-value: 0 TIGR02082 2 nkrilvlDGamGtqlqsanLteadFrge.eadlarelkGnndlLnltkPeviaaihrayfeaGaDivetnt 71 ++rilvlDGamGt +qs++L ead+rge ++d++++l GnndlL+lt+Pe+i+++h ay++aGaD++ tnt NCBI__GCF_000733765.1:WP_052368278.1 22 SERILVLDGAMGTRIQSYKLGEADYRGErFKDWESDLAGNNDLLVLTQPELIKEVHGAYLDAGADVIGTNT 92 79********************************************************************* PP TIGR02082 72 FnsteialadYdledkayelnkkaaklarevadeftltpekkRfvaGslGPtnklatlspdverpefrnvt 142 Fn+ i++adY++e++ayeln +aakla e+ad t+e++RfvaG++GPtn++a++spdv++p+ rn+ NCBI__GCF_000733765.1:WP_052368278.1 93 FNCQFISMADYHMEELAYELNVEAAKLACEAADSRA-TAERPRFVAGAVGPTNRTASISPDVNDPGKRNTS 162 *********************************999.********************************** PP TIGR02082 143 ydelvdaYkeqvkglldGGvDllLietvfDtlnakaalfaveevfeekgrelPilisgvivdksGrtLsGq 213 +delv aY+eq++gl+dGGvDl+Liet+fDtlnakaa+fa+++vf+ek+++lPi+isg+i+d+sGrtLsGq NCBI__GCF_000733765.1:WP_052368278.1 163 FDELVRAYSEQIHGLIDGGVDLILIETIFDTLNAKAAIFAAHTVFAEKEIKLPIIISGTITDASGRTLSGQ 233 *********************************************************************** PP TIGR02082 214 tleaflaslehaeililGLnCalGadelrefvkelsetaealvsviPnaGLPnalgeYdltpeelakalke 284 +eaf +s++hae l++GLnCalG ++r++++els+ a++ vs++PnaGLPn++geYd++pee+a++++e NCBI__GCF_000733765.1:WP_052368278.1 234 VTEAFWNSIRHAEPLAVGLNCALGGRDMRPYIEELSKLADCNVSCYPNAGLPNTFGEYDEEPEETAEIVAE 304 *********************************************************************** PP TIGR02082 285 faeegllnivGGCCGttPehiraiaeavkdikprkrqeleeksvlsglealkiaqessfvniGeRtnvaGs 355 fa++g++n vGGCCGt+P hi aia+av++ikpr +++e+ ++lsgle+l++++++ f+niGeRtn++Gs NCBI__GCF_000733765.1:WP_052368278.1 305 FAQAGFINFVGGCCGTQPSHIEAIANAVREIKPRVPPAHEPITRLSGLEPLNVTKDKGFINIGERTNITGS 375 *********************************************************************** PP TIGR02082 356 kkfrklikaedyeealkiakqqveeGaqilDinvDevllDgeadmkkllsllasepdiakvPlmlDssefe 426 kfr+lik++dy al++a++qve+Ga i+Din+De++ Dg+++m+++++l+asepdi+ +P+m+Dss++e NCBI__GCF_000733765.1:WP_052368278.1 376 TKFRNLIKDGDYATALDVARSQVENGAAIMDINMDEGMIDGVEAMRTFCQLVASEPDISCIPMMVDSSKWE 446 *********************************************************************** PP TIGR02082 427 vleaGLkviqGkaivnsislkdGeerFlekaklikeyGaavvvmafDeeGqartadkkieiakRaykllte 497 v+e+GLk++qGk++vnsis+k+Gee+Fl++a+l ++yGaavvvmafDe+Gqa++ +++ ei+ Raykllte NCBI__GCF_000733765.1:WP_052368278.1 447 VIETGLKCLQGKCVVNSISMKEGEEKFLREANLCRHYGAAVVVMAFDEDGQADNLERRKEICGRAYKLLTE 517 *********************************************************************** PP TIGR02082 498 kvgfppediifDpniltiatGieehdryaidfieaireikeelPdakisgGvsnvsFslrgndavRealhs 568 k++fppediifDpni+++atGieeh ry++dfiea+ +ike+ P+++isgGvsn+sFs+rgn++vRea+hs NCBI__GCF_000733765.1:WP_052368278.1 518 KLDFPPEDIIFDPNIFAVATGIEEHARYGLDFIEACAWIKENCPHVHISGGVSNISFSFRGNNKVREAIHS 588 *********************************************************************** PP TIGR02082 569 vFLyeaikaGlDmgivnagklavyddidkelrevvedlildrr...reatekLlelaelykgtkeksskea 636 vFLy+a kaG+Dmgivnag+l+vy+d+++elr+++ed+il+r+ +eatekLlelae+ykg ++ ++ NCBI__GCF_000733765.1:WP_052368278.1 589 VFLYHATKAGMDMGIVNAGQLTVYEDVEQELRDAIEDVILNRKedpQEATEKLLELAEKYKGDASEADQ-G 658 ******************************************988889***************999777.5 PP TIGR02082 637 qeaewrnlpveeRLeralvkGeregieedleear....kklkapleiiegpLldGmkvvGdLFGsGkmfLP 703 + +wr+ + +eR+++al+kG+ ++ eed+ee r + +p+e+iegpL+dGm vvGdLFG+GkmfLP NCBI__GCF_000733765.1:WP_052368278.1 659 DRLAWRENDARERIKHALIKGIDQYAEEDTEELRaeisANGGKPIEVIEGPLMDGMGVVGDLFGAGKMFLP 729 5668****************************99444445568**************************** PP TIGR02082 704 qvvksarvmkkavayLePylekekeed...kskGkivlatvkGDvhDiGknivdvvLscngyevvdlGvkv 771 qvvksarvmkkavayL+P++e+ekee + +G++++atvkGDvhDiGkniv+vvL+cn+yev+dlGv+ NCBI__GCF_000733765.1:WP_052368278.1 730 QVVKSARVMKKAVAYLIPFIEAEKEEGedpSDNGTVIMATVKGDVHDIGKNIVGVVLQCNNYEVIDLGVMC 800 ************************98878899*************************************** PP TIGR02082 772 PvekileaakkkkaDviglsGLivksldemvevaeemerrgvkiPlllGGaalskahvavkiaekYkgevv 842 P++kil+aak+++aD+iglsGLi++sldemv+ a+emer+g++iPll+GGa++s+ah+avki++kY+g+vv NCBI__GCF_000733765.1:WP_052368278.1 801 PAQKILDAAKEHDADIIGLSGLITPSLDEMVHFAKEMERQGLDIPLLIGGATTSRAHTAVKIEPKYHGPVV 871 *********************************************************************** PP TIGR02082 843 yvkdaseavkvvdkllsekkkaeelekikeeyeeirekfgekkeklialsekaarkevfaldrsedlevpa 913 +vkdas++v v+kll+++++ ++le++k++ye++re++ k++ + l++++ar ++ ++d++ +e+p+ NCBI__GCF_000733765.1:WP_052368278.1 872 WVKDASRSVPTVAKLLDKNERGAMLEELKADYEKLRERHAGKEKADKLLPIDKARANRDVPDWE-GYEPPV 941 ****************************************************************.999999 PP TIGR02082 914 pkfl...G.....tkvleas.ieellkyiDwkalFvqWelrgkypkilkdeleglearklfkdakelldkl 975 p+ + G +kv++++ + el++ iDw ++F +We++g++p +l+++ g+ arkl+++a+++ld++ NCBI__GCF_000733765.1:WP_052368278.1 942 PQLIseaGdkaeyVKVFKDYpLTELRECIDWGPFFNAWEIKGRFPDVLNNPASGETARKLYEEANAMLDQI 1012 99773333444445899999*************************************************** PP TIGR02082 976 saekllrargvvGlfPaqsv.gddieiytdetvsqetkpiatvrekleqlrqqsdr.....ylclaDfias 1040 ++ek l+a+g++GlfPa++v ++di++ytde++s e + + rqqs++ + +l+D++a+ NCBI__GCF_000733765.1:WP_052368278.1 1013 IEEKWLTANGIIGLFPANQVdNEDIAVYTDESRSDELLRLHCL-------RQQSEHregvpNRALSDYVAP 1076 *****************9762579******9996665555555.......444444444569********* PP TIGR02082 1041 kesGikDylgallvtaglgaeelakkleakeddydsilvkaladrlaealaellhervRkelwgyaeeenl 1111 e+G kDy+ga++vtaglg +e +++ + +ddy++il+++ladrlaea+ae lh+rvRke+wgy +e+l NCBI__GCF_000733765.1:WP_052368278.1 1077 AETGLKDYVGAFAVTAGLGSSEKCEEFRKANDDYNAILLESLADRLAEAFAERLHQRVRKEFWGYVPDEKL 1147 *********************************************************************** PP TIGR02082 1112 dkedllkerYrGirpafGYpacPdhtekatlleLleaer.iGlklteslalaPeasvsglyfahpeakYf 1180 d++ l++e+Y+Girpa+GYpacP+htek+ +++Ll++e+ +G++ltes+a++P asvsgl+fahp+++Yf NCBI__GCF_000733765.1:WP_052368278.1 1148 DNDALIQEKYQGIRPAPGYPACPEHTEKQAIWDLLDVENnTGMQLTESMAMWPGASVSGLFFAHPDSRYF 1217 ***************************************9****************************** PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (1182 nodes) Target sequences: 1 (1255 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.02u 0.01s 00:00:00.03 Elapsed: 00:00:00.02 # Mc/sec: 53.29 // [ok]
This GapMind analysis is from Jul 25 2024. The underlying query database was built on Jul 25 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