GapMind for Amino acid biosynthesis

 

Alignments for a candidate for metH in Algiphilus aromaticivorans DG1253

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.

Links

Downloads

Related tools

About GapMind

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:

where "other" refers to the best ublast hit to a sequence that is not annotated as performing this step (and is not "ignored").

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