Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	374,326	G→A	K153K (AAG→AAA)	mhpE →	4‑hydroxy‑2‑oxovalerate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	374,326	0	G	A	100.0%	168.6 / NA	50	K153K (AAG→AAA)	mhpE	4‑hydroxy‑2‑oxovalerate aldolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/30);  total (20/30)

CGTGTGGCAACGCACTGTACCGAAGCTGATGTTTCCGCCCAGCATATTCAGTATGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAGCTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCGATGAACATGAGCGATATCCGTGACCGTTTCCGCGCCCTGAAAGCAGAGCTGAAACCAGAAACGCAAACTGGCATGCACGCT  >  NC_000913/374189‑374461                                                                                                                                          |                                                                                                                                        cGTGTGGCAACGCACTGTACCGAAGCTGATGTTTCCGCCCAGCATATTCAGTATGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACt                                                                                                                                       >  2:111138/1‑140 (MQ=255)     tgGCAACGCACTGTACCGAAGCTGATGTTTCCGCCCAGCATATTCAGTATGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATg                                                                                                                                   >  1:202901/1‑140 (MQ=255)        cAACGCACTGTACCGAAGCTGATGTTTCCGCCCAGCATATTCAGTATGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGaa                                                                                                                                >  2:757540/1‑140 (MQ=255)                  tACCGAAGCTGATGTTTCCGCCCAGCATATTCAGTATGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGa                                                                                                                                 >  2:715588/1‑129 (MQ=255)                  tACCGAAGCTGATGTTTCCGCCCAGCATATTCAGTATGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGa                                                                                                                                 <  1:715588/129‑1 (MQ=255)                  tACCGAAGCTGATGTTTCCGCCCAGCATATTCAGTATGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTg                                                                                                                      <  2:817864/140‑1 (MQ=255)                  tACCGAAGCTGATGTTTCCGCCCAGCATATTCAGTATGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTg                                                                                                                      <  1:547988/140‑1 (MQ=255)                    ccGAAGCTGATGTTTCCGCCCAGCATATTCAGTATGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCg                                                                                                                    <  2:512843/140‑1 (MQ=255)                     cGAAGCTGATGTTTCCGCCCAGGATATTCAGTATGCCCGCGAGCTCGGTATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGa                                                                                                                   <  2:268300/140‑1 (MQ=255)                      gAAGCTGATGTTTCCGCCCAGCTTATTCAGTATGCACGCTAGCTCGGACTGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGAc                                                                                                                  <  2:350704/140‑1 (MQ=255)                      gAAGCTGATGTTTCCGCCCAGCATATTCAGTATGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGAc                                                                                                                  <  1:120017/140‑1 (MQ=255)                      gAAGCTGATGTTTCCGCCCAGCATATTCAGTATGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGAc                                                                                                                  <  2:479166/140‑1 (MQ=255)                      gAAGCTGATGTTTCCGCCCAGCATATTCAGTATGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGAc                                                                                                                  <  1:807272/140‑1 (MQ=255)                      gAAGCTGATGTTTCCGCCCAGCATATTCAGTATGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGAc                                                                                                                  <  1:785672/140‑1 (MQ=255)                       aaGCTGATGTTTCCGCCCAGCATATTCAGTATGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGAcc                                                                                                                 <  1:695022/140‑1 (MQ=255)                       aaGCTGATGTTTCCGCCCAGCATATTCAGTATGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGAcc                                                                                                                 <  2:551647/140‑1 (MQ=255)                         gCTGATGTTTCCGCCCAGCATACTCGGTAGGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTg                                                                                                               <  1:329714/140‑1 (MQ=255)                            gATGTTTCCGCCCAGCATATTCAGTATGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCATGCATGTAGAACTGATGGAAGGGTACGGTGCGACCTGTAt                                                                                                            >  2:367801/1‑140 (MQ=255)                            gATGTTTCCGCCCAGCATATTCAGTATGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTAt                                                                                                            >  2:207538/1‑140 (MQ=255)                            gATGTTTCCGCCCAGCATATTCAGTATGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGAAGGCAAAACTGCTGGAAGGCTCCGGTGCGACCTGTAt                                                                                                            >  1:322225/1‑140 (MQ=255)                            gATGTTTCCGCCCAGCATATTCAGTATGCCCGCGAGCTCGGAATGGACACCATTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTAt                                                                                                            >  1:753391/1‑140 (MQ=255)                              tGTTTCCGCCCAGCATATTCAGTATGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTAttt                                                                                                          <  1:673596/140‑1 (MQ=255)                                tttCCGCCCAGCATATTCAGTATGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTATTTAt                                                                                                        <  2:202901/140‑1 (MQ=255)                                      cccAGCATATTCAGTATGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTATTTATgtggtg                                                                                                  <  1:25991/140‑1 (MQ=255)                                                 cAGTATGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCgg                                                                                       <  2:774502/140‑1 (MQ=255)                                                   gTATGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTg                                                                                     <  1:569143/140‑1 (MQ=255)                                                     aTGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGTTGCGACCGGTAATTATGTGGTTGATTCTGGCGGTGCg                                                                                   >  1:240160/1‑140 (MQ=255)                                                     aTGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCg                                                                                   <  1:644710/140‑1 (MQ=255)                                                     aTGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATAGAAGGCTAGGGTCCGAACTGTATTTATGTGGTGGATTCTGGCGGTGCg                                                                                   >  2:241224/1‑140 (MQ=255)                                                       gCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCg                                                                                   <  1:769539/138‑1 (MQ=255)                                                       gCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCg                                                                                   >  2:769539/1‑138 (MQ=255)                                                       gCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCGAt                                                                                 >  1:664494/1‑140 (MQ=255)                                                         ccGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCGATGa                                                                               <  1:757540/140‑1 (MQ=255)                                                         ccGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTATTTAGGTGGTGGATCCTGGCGGTGCGATGa                                                                               >  1:268886/1‑140 (MQ=255)                                                            cgAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCGATGAACa                                                                            <  2:810253/140‑1 (MQ=255)                                                            cgAGCTCGGAATGGACACCGTTGGTTATCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCGATGAACa                                                                            >  2:432894/1‑140 (MQ=255)                                                                  cGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCGATGAACATGAGCg                                                                      >  1:709382/1‑140 (MQ=255)                                                                  cGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCGATGAACATGAGCg                                                                      <  1:798146/140‑1 (MQ=255)                                                                                     tttCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCGATGAACATGAGCGATATCCGTGACCGTTTCcg                                                   <  1:89643/140‑1 (MQ=255)                                                                                        cTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCGATGAACATGAGCGATATCCGTGACCGTTTCCgcgc                                                >  1:182603/1‑140 (MQ=255)                                                                                           atgatgaGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCGATGAACATGAGCGATATCCGTGACCGTTTCCGCGCCCt                                             <  2:182603/140‑1 (MQ=255)                                                                                                     atatGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCGATGAACATGAGCGATATCCGTGACCGTTTCCGCGCCCt                                             <  1:801722/130‑1 (MQ=255)                                                                                                     atatGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCGATGAACATGAGCGATATCCGTGACCGTTTCCGCGCCCt                                             >  2:801722/1‑130 (MQ=255)                                                                                                             aCGCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCGATGAACATGAGCGATATCCGTGACCGTTTCCGCGCCCTGAAAGCAGAGCTGAAAcc                           >  2:826602/1‑140 (MQ=255)                                                                                                               gCCGGAGAATCTCGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCGATGAACATGAGCGATATCCGTGACCGTTTCCGCGCCCTGAAAGCAGAGCTGAAACCAg                         <  1:207538/140‑1 (MQ=255)                                                                                                                        tctcGCCAAGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCGATGAACATGAGCGATATCCGTGACCGTTTCCGCGCCCTGAAAGCAGAGCTGAAACCAGAAACGCaaa                <  1:826602/140‑1 (MQ=255)                                                                                                                                aGCAGGCAAAACTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCGATGAACATGAGCGATATCCGTGACCGTTTCCGCGCCCTGAAAGCAGAGCTGAGACCAGATACGCAAACTGGCATg        >  1:287937/1‑140 (MQ=255)                                                                                                                                    ggCAAAACTGATGGAAGGCTACCGTACGAGCTGTATTTATGTGGTGGATTCTGGCGGTGCGATGAACATGAGCGATATCCGTGACCGTTTCCGCGCCCTGAAAGCAGAGCTGAAACCAGAAACGCAAACTGGCATGCACg    <  1:351093/140‑1 (MQ=255)                                                                                                                                     gCAAAACTGATGGAAGGCTACGATGCGACCTGTATGTATGTGGTGGATTCTGGCGGTGCGATGAACATGAGCGATATCCGTGACCGTTTCCGCGCCCTGAAAGCAGAGCTGAAACCAGAAACGCAAACTGGCATGCACGc   <  2:322225/140‑1 (MQ=255)                                                                                                                                      cAAAACTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCGATGAACATGAGCGATATCCGTGACCGTTTCCGCGCCCTGAAAGCAGAGCTGAAACCAGAAACGCAAACTGGCATGCACGCt  >  2:480202/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                        CGTGTGGCAACGCACTGTACCGAAGCTGATGTTTCCGCCCAGCATATTCAGTATGCCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAGCTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCGATGAACATGAGCGATATCCGTGACCGTTTCCGCGCCCTGAAAGCAGAGCTGAAACCAGAAACGCAAACTGGCATGCACGCT  >  NC_000913/374189‑374461

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑