Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	4,233,746	A→G	intergenic (‑513/‑2161)	lysC ← / → yjbE	lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	4,233,746	0	A	G	100.0%	159.5 / NA	51	intergenic (‑513/‑2161)	lysC/yjbE	lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/17);  total (34/17)

ACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCT  >  NC_000913_3_bme_pgi/4233602‑4233887                                                                                                                                                 |                                                                                                                                               aCTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTGCACTGCGCTAAGGGTTTACGCTCAGCATTACGCGAGCGGCCCTAAAACCATCACAATGTTCTGTCACTGGCGCTACAATCTTGCAAACTCACAACTTCCAGAATCAGAcgagc                                                                                                                                            >  2:645794/1‑148 (MQ=255)     cTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGCAGTATTg                                                                                                                                        >  2:624394/1‑149 (MQ=255)           gTCAGTAGTTCAGCTAATAAATGCTTCACTGCGTTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATg                                                                                                                                 <  2:620568/149‑1 (MQ=255)                   ttCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAcgaaa                                                                                                                         >  2:663650/1‑147 (MQ=255)                   ttCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAgaga                                                                                                                         <  2:615399/149‑1 (MQ=255)                        cTAATCAATGCTTCACTGCGGTAACGGTTTACTCTCAGCATGACGCTAACGGGACTAAAACCCTCACAGTTTTCTGTGACTGGCGCTACAATCTTCCGAGGTCTCCATTCTCGATTTCAGAGGAGTATTGCTA‑ATGGCAAGAGACGcga                                                                                                                    <  2:639040/149‑1 (MQ=255)                              aaggCTTCACGGCGGCAAGGGTGTACGCTCAACGTTACGCTAACGGGCCGAGAACGATCACAGTTGTCTGTGACTGGCGCTACCATTTTCCACACGCTCAATTCTCCAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGc                                                                                                              <  2:245238/146‑1 (MQ=255)                              aaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCTTCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAga                                                                                                              >  2:331804/1‑148 (MQ=255)                              aaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCCCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCTGAGGAATATTGCTA‑ATGGCAAGAGACGCGACGAAGc                                                                                                              >  1:51230/1‑149 (MQ=255)                              aaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGc                                                                                                              >  1:604812/1‑149 (MQ=255)                              aaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATATTCACGAGATACTATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAGGAGACGCGACGAAGc                                                                                                              >  1:586618/1‑149 (MQ=255)                                 gCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACGATGTTTAATTGACTCAATTCTCAAAATCGGAGGAGTATTGCTA‑ATGGCAAGAGATGCGACGAAGCTTg                                                                                                           >  1:588129/1‑149 (MQ=255)                                 gCTTCACTGCGCTAAAGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAATATTACTA‑ATGGCAAGAGACGCGACGAAGCTTg                                                                                                           >  2:294848/1‑149 (MQ=255)                                   ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCCCTCAAACCATCACATTTTTCTTTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGaa                                                                                                         >  1:39847/1‑149 (MQ=255)                                    tCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAg                                                                                                        >  1:101154/1‑149 (MQ=255)                                      aCTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAAGAGAGGTA‑AGGGCAAGAGACGCGAGGAAGCTTGAAGCa                                                                                                      >  1:750437/1‑149 (MQ=255)                                               aGGGTTTACACTCAACATTACGCTAACGGCACTAAGACCATCACATTTTTTTGTGTTTTGTGTTACAATCTTCAGAAGTCACAGTTCTCAAGATGGGGGGGGTGTTGCTA‑TTGGCAGTAGACGCGGGGAGGCTTGTGGCAACAGttcga                                                                                             >  2:334571/1‑145 (MQ=255)                                               aGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCa                                                                                             >  1:520021/1‑149 (MQ=255)                                                  gTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAg                                                                                          <  2:659943/149‑1 (MQ=255)                                                   tttACACTCAACATTACGCAAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAACCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGc                                                                                         <  2:642665/149‑1 (MQ=255)                                                    ttACACTCCACATTACGCTAACGGCACTAAAACCATCACCTTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCt                                                                                        <  1:604395/149‑1 (MQ=255)                                                    ttACACTCAACAGCACCCTAACGGAACTAAAACCATCACATGTTTATTGAAATGGTGGAACAATCTTCCAAAGTCACAATTCTTTTAATGAGAGGAGTTTGGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCt                                                                                        <  2:45928/149‑1 (MQ=255)                                                         cTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGAAGAAGCTTGAAGCAACAGTCGCAAGACTGaaaa                                                                                   >  1:204969/1‑149 (MQ=255)                                                          tCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCCCAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCAACGAAGCTTGAAGCAACGGGCACAAAGCGCaaaaa                                                                                  >  2:235369/1‑149 (MQ=255)                                                          tCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGATGCAACGGTCGCAAAGCTGaaaaa                                                                                  >  1:427932/1‑149 (MQ=255)                                                          tCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaaa                                                                                  >  2:595281/1‑149 (MQ=255)                                                          tCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaaa                                                                                  >  1:797696/1‑149 (MQ=255)                                                          tCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTACCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCGTGAAGCAACGGTCGCAAAGCTGaaaaa                                                                                  >  1:396298/1‑149 (MQ=255)                                                              cATTACGCTACCGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTAAACGCAACGTTCGAGCAGATGaaaatcgat                                                                              >  1:783793/1‑144 (MQ=255)                                                                     cTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAGCACCGGTCGCAAAGCTGAAAAAGCATTGGGCgg                                                                       >  2:519923/1‑149 (MQ=255)                                                                     cTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCgg                                                                       >  1:744288/1‑149 (MQ=255)                                                                       aagtgCACTATAAGCATCGCATTTTTCTGTTTCTCGCGCTACATTCGTCCAAAGTCACAAGTGTCAAAATCAGAGTAGTATTGCTA‑ATGGCTACGGACGCGATGAAGGTTGAAGCCACGGTCGCAAAGCGGAAAAAGCATTGGGCGGaa                                                                     <  2:5807/145‑1 (MQ=255)                                                                        aCGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAAc                                                                                                    >  2:267028/1‑117 (MQ=255)                                                                        aCGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAAc                                                                                                    <  1:267028/117‑1 (MQ=255)                                                                               tAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGcccc                                                             >  1:464443/1‑149 (MQ=255)                                                                                        ctcctttttCTGTGTAGGGCGCTACAAGCTTCCAAAGTCACAATTCTCAGAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTGGAAGGAACGGTCGGAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGatat                                                    <  2:428744/145‑1 (MQ=255)                                                                                                    tgACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGGCGTGTCTGCCCCGCGCGATATGCGCGCGGCAtt                                        <  1:595281/149‑1 (MQ=255)                                                                                                        tGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAgc                                    <  2:464443/149‑1 (MQ=255)                                                                                                                  aTCTTCCAAAGTCACAATTTTCATGATTGGAGGAGTATTGCTA‑GTGGCAAGAGACGCGACGAGGCTTGAAGCAACGGGCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATTTGCGCGTGGCATTCAGCGCCGACCccg                          <  2:300605/149‑1 (MQ=255)                                                                                                                    cTTCCAAAGTCACAAATCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAAACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTAGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccggc                        >  2:660319/1‑149 (MQ=255)                                                                                                                        cAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGt                    >  1:627506/1‑149 (MQ=255)                                                                                                                             tCACAATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCggg               <  1:444466/149‑1 (MQ=255)                                                                                                                               acaATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCg             >  1:801719/1‑149 (MQ=255)                                                                                                                               acaATTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTCGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGATTCGGGCg             >  1:711599/1‑149 (MQ=255)                                                                                                                               acaATTCTCAAAATCAGAGGAGTATTGCAAGACAGACAG‑GATGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCGCGGCCGGTTCGGGCg             >  1:794499/1‑149 (MQ=255)                                                                                                                                 aaTTCTCAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGat           >  1:534712/1‑149 (MQ=255)                                                                                                                                      tcAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGTCCCCGGCCGGTTCGGGCGATATTCg      >  2:623774/1‑149 (MQ=255)                                                                                                                                      tcAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCg      >  1:741983/1‑149 (MQ=255)                                                                                                                                       cAAAATCAGAGGAGTATTGCTA‑ATGGCAAGAGGCGCGACGAGGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATACGCGCGGGCTTCAGCACCGACACCGGCCGGTTCGGGCGATATTCGc     >  2:467687/1‑149 (MQ=255)                                                                                                                                         aaaTCAGAGTAGTATCGATA‑ATGGCAAGAGACGCAACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGctc   <  2:493776/149‑1 (MQ=255)                                                                                                                                          aaTCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCACCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGctct  >  1:768419/1‑149 (MQ=255)                                                                                                                                          aaTCAGAGGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCAGAAACGACTGCTGTCGATATGCGCGCGGCATTCAGTGCCGACCCCGGCCGGTTCGGGCGATATTCGctct  <  1:151749/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                               ACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCT  >  NC_000913_3_bme_pgi/4233602‑4233887

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