Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,374,970	G→A	R12W (CGG→TGG)	menE ←	O‑succinylbenzoate‑CoA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,374,970	0	G	A	100.0%	147.6 / NA	41	R12W (CGG→TGG)	menE	O‑succinylbenzoate‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/25);  total (16/25)

CCACGCCCGCAACATCACGCCGCTGCCCTCAACCACCCCCTGCACCGCAAATCCGGAGGCTAATTCATCGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCGCCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCAACTTCCACGACAGGCAGCGTGCTACCCGGCCAGCGACGTACCTGCTGCGCCTGCATCAGATCCAGCGTGTCCAGCCCTGGAATGGTGTCCGG  >  NC_000913/2374824‑2375116                                                                                                                                                   |                                                                                                                                                   ccACGCCCGCAACATCACGCCGCTGCCCTCAACCACCCCCTGCACCGCAAATCCGGAGGCTAATTCATCGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGccacca                                                                                                                                                 <  1:648316/150‑1 (MQ=255)      cccGCAACATCACGCCGCTGCCCTCAACCACCCCCTGCACCGCAAATCCGGAGGCTAATTCATCGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACg                                                                                                                                           >  2:556731/1‑151 (MQ=255)          cAATATCACGCCGCTGCCCTCAACCACCCCCTGCACCGCAAATCCGGAGGCTAATTCATCGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCAc                                                                                                                                       >  1:277045/1‑151 (MQ=255)          cAACATCACGCCGCTGCCCTCAACCACCCCCTGCACCGCAAATCCGGAGGCTAATTCATCGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCAc                                                                                                                                       >  1:171839/1‑151 (MQ=255)          cAACATCACGCCGCTGCCCTCAACCACCCCCTGCACCGCAAATCCGGAGGCTAATTCATCGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCAc                                                                                                                                       <  1:270728/151‑1 (MQ=255)           aaCATCACGCCGCTGCCCTCAACCACCCCCTGCACCGCAAATCCGGAGGCTAATTCATCGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACg                                                                                                                                      <  1:244306/151‑1 (MQ=255)           aaCATCACGCCGCTGCCCTCAACCACCCCCTGCACCGCAAATCCGGAGGCTAATTCATCGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACg                                                                                                                                      <  1:75525/151‑1 (MQ=255)           aaCATCACGCCGCTGCCCTCAACCACCCCCTGCACCGCAAATCCGGAGGCTAATTCATCGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACg                                                                                                                                      <  1:394810/151‑1 (MQ=255)                 aCGCCGCTGCCCTCAACCACCCCCTGCACCGCAAATCCGGAGGCTAATTCATCGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGt                                                                                                                                <  1:46856/151‑1 (MQ=255)                   gccgcTGCCCTCAACCACCCCCTGCACCGCAAATCCGGAGGCTAATTCATCGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCa                                                                                                                              <  2:448477/151‑1 (MQ=255)                    ccgcTGCCCTCAACCACCCCCTGCACCGCAAATCCGGAGGCTAATTCATCGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCag                                                                                                                             <  2:759549/151‑1 (MQ=255)                     cgcTGCCCTCAACCACCCCCTGCACCGCAAATCCGGAGGCTAATTCATCGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTc                                                                                                                               >  1:659376/1‑148 (MQ=255)                         gCCCTCAACCACCCCCTGCACCGCAAATCCGGAGGCTAATTCATCGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGagaa                                                                                                                         <  2:227087/150‑1 (MQ=255)                            cTCAACCACCCCCTGCACCGCAAATCCGGAGGCTAATTCATCGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAgaga                                                                                                                          >  1:228385/1‑146 (MQ=255)                               aaCCACCCCCTGCACCGCAAATCCGGAGGCTAATTCATCGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATa                                                                                                                  <  2:283065/151‑1 (MQ=255)                                 ccaccCCCTGCACCGCAAATCCGGAGGCTAATTCATCGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATaac                                                                                                                <  1:271481/151‑1 (MQ=255)                                   accCCCTGCACCGCAAATCCGGAGGCTAATTCATCGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATaaca                                                                                                               >  1:683558/1‑150 (MQ=255)                                               gCAAATCCGGAGGCTAATTCATCGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGc                                                                                                  >  1:170718/1‑151 (MQ=255)                                                   aTCCGGAGGCTAATTCATCGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCa                                                                                              <  2:677170/151‑1 (MQ=255)                                                             tAATTCATCGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCAACTTCCACg                                                                                     >  1:690215/1‑150 (MQ=255)                                                              aaTTCATCGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCAACTTCCACGAc                                                                                   <  1:85059/151‑1 (MQ=255)                                                                   aTCGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCAACTTCCACGACAGGCa                                                                              <  1:479137/151‑1 (MQ=255)                                                                    tCGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCAACTTCCACGACAGGc                                                                               >  2:290433/1‑149 (MQ=255)                                                                     cGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCaa                                                                                             <  1:735429/134‑1 (MQ=255)                                                                     cGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCaa                                                                                             >  2:735503/1‑134 (MQ=255)                                                                     cGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCaa                                                                                             >  2:380888/1‑134 (MQ=255)                                                                     cGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCaa                                                                                             <  1:380840/134‑1 (MQ=255)                                                                      gACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCAACTTCCACGACAGGCAGCg                                                                           <  2:440906/151‑1 (MQ=255)                                                                                 aaaGCTCGCGCCAGTTGAGTTGCTCATCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCAACTTCCACGACAGGCAGCGTGCTACCCGGc                                                                <  1:704744/151‑1 (MQ=255)                                                                                  aaGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCAACTTCCACGACAGGCAGCGTGCTACCCGGcc                                                               >  2:214573/1‑151 (MQ=255)                                                                                         cgcCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCAACTTCCACGACAGGCAGCGTGCTACCCGGCCAGCGACg                                                        <  2:20303/151‑1 (MQ=255)                                                                                              gttgagttgCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCAACTTCCACGACAGGCAGCGTGCTACCCGGCCAGCGACGTACct                                                   >  1:330821/1‑151 (MQ=255)                                                                                                tgagttgCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCAACTTCCACGACAGGCAGCGTGCTACCCGGCCAGCGACGTACctg                                                  <  2:603192/150‑1 (MQ=255)                                                                                                         cGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCAACTTCCACGACAGGCAGCGTGCTACCCGGCCAGCGACGTACCTGCTGCGCCTGc                                        <  1:640357/151‑1 (MQ=255)                                                                                                                     gTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCAACTTCCACGACAGGCAGCGTGCTACCCGGCCAGCGACGTACCTGCTGCGCCTGCATCAGATCCAGc                            <  2:13379/151‑1 (MQ=255)                                                                                                                     gTAAGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCAACTTCCACGACAGGCAGCGTGCTACCCAGCCAGCGACGTACCTGCTGCGCCTGCATCAGATCCAGc                            <  2:219629/151‑1 (MQ=255)                                                                                                                        aGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCAACTTCCACGACAGGCAGCGTGCTACCCGGCCAGCGACGTACCTGCTGCGCCTGCATCAGATCCAGCgtg                         <  2:159971/151‑1 (MQ=255)                                                                                                                        aGGCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCAACTTCCACGACAGGCAGCGTGCTACCCGGCCAGCGACGTACCTGCTGCGCCTGCATCAGATCCAGCgtg                         >  1:541501/1‑151 (MQ=255)                                                                                                                         ggCGATGGTTTCTCCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCAACTTCCACGACAGGCAGCGTGCTACCCGGCCAGCGACGTACCTGCTGCGCCTGCATCAGATCCAGCgtgt                        >  2:698166/1‑151 (MQ=255)                                                                                                                                    ctcCCCGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCAACTTCCACGACAGGCAGCGTGCTACCCGGCCAGCGACGTACCTGCTGCGCCTGCATCAGATCCAGCGTGTCCAGCCCTgg              >  1:342301/1‑150 (MQ=255)                                                                                                                                      ccccGCACTTGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCAACTTCCACGACAGGCAGCGTGCTACCCGGCCAGCGACGTACCTGCTGCGCCTGCATCAGATCCAGCGTGTCCAGCCCTGGaa            <  1:79445/150‑1 (MQ=255)                                                                                                                                               tGCCACCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCAACTTCCACGACAGGCAGCGTGCTACCCGGCCAGCGACGTACCTGCTGCGCCTGCATCAGATCCAGCGTGTCCAGCCCTGGAATGGTGTCCg   <  2:638170/150‑1 (MQ=255)                                                                                                                                                gccaccAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCAACTTCCACGACAGGCAGCGTGCTACCCGGCCAGCGACGTACCTGCTGCGCCTGCATCAGATCCAGCGTGTCCAGCCCTGGAATGGTGTCCgg  <  1:204404/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CCACGCCCGCAACATCACGCCGCTGCCCTCAACCACCCCCTGCACCGCAAATCCGGAGGCTAATTCATCGACGCGAGCACAAAGCTCGCGCCAGTTGAGTTGCTCGTCATTAAGACGTAAGGCGATGGTTTCTCCCCGCACTTGCCGCCAGTGACGCCACGGCCAGTCAGAGAAGATCATAACAACCGCTCCAGTGCATCAACTTCCACGACAGGCAGCGTGCTACCCGGCCAGCGACGTACCTGCTGCGCCTGCATCAGATCCAGCGTGTCCAGCCCTGGAATGGTGTCCGG  >  NC_000913/2374824‑2375116

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