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%	149.9 / NA	52	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 (22/30);  total (22/30)

CCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTG  >  NC_000913_3_bme_pgi/4233600‑4233888                                                                                                                                                   |                                                                                                                                               ccACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCTCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCCTCACATTTTTCTGTGACTGGCGCTACAGTCTTCCATAGTCACAATTCTCAAAATCAgaggag                                                                                                                                             >  1:523225/1‑150 (MQ=255) ccACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAgaggag                                                                                                                                             <  2:780149/150‑1 (MQ=255)  cACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGTGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGt                                                                                                                                            <  2:581822/150‑1 (MQ=255)   aCTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTa                                                                                                                                           >  1:276228/1‑150 (MQ=255)   aCTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTa                                                                                                                                           >  1:343836/1‑150 (MQ=255)      ccTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTg                                                                                                                                        <  1:156639/150‑1 (MQ=255)            gttCAGTAGTTCATCTAATAATTGCTTCACTGCGCTAAGGGTTTACACTGCTCATTTAGCTATCGTCACTAAAACCATCACATTTTTCTGTGACTGGCGCTCCAATCTTGCAAAGTCACATTTCTCAAAATCAGAGGAGTATTGCTAATg                                                                                                                                  <  2:801661/148‑1 (MQ=255)             ggcAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATgg                                                                                                                                 <  2:580042/148‑1 (MQ=255)             gTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATgg                                                                                                                                 <  2:276228/150‑1 (MQ=255)             gTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATgg                                                                                                                                 <  2:343836/150‑1 (MQ=255)                  tagtTCAGCTAATAAATTCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCTCATTTTTCTGGGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAaga                                                                                                                            <  1:524841/150‑1 (MQ=255)                    gtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAgaga                                                                                                                          <  1:665507/150‑1 (MQ=255)                          cTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGcgac                                                                                                                    >  1:420655/1‑150 (MQ=255)                            aataaATGCTTCACTGCGCTAAGGGTATACACTCAACATTACGCTAACGCCACTAAAACCATCACATTATTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTCCGAATGGCAACAGACGcgacga                                                                                                                  >  1:478657/1‑150 (MQ=255)                                   gcatCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGa                                                                                                           <  2:854903/147‑1 (MQ=255)                                      tCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGc                                                                                                        <  2:308065/150‑1 (MQ=255)                                         cTGCGCTAAGTGTTTACACTCAACATTACGCTAACGGAACTACAACCATCACATATTTCTGTGACTGGCGCTACACTCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGAAAGAGACGCGACGAAGCTTGAAGCAAc                                                                                                     <  1:559076/150‑1 (MQ=255)                                               tatgggTTTACACTCAACAATACGCTAACGGCACTAAAACCATCACATTTTTATGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGATTATATCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGc                                                                                               <  1:850824/147‑1 (MQ=255)                                               tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGc                                                                                               <  1:583232/150‑1 (MQ=255)                                               tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGc                                                                                               >  1:22879/1‑150 (MQ=255)                                               tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGc                                                                                               >  1:895721/1‑150 (MQ=255)                                                aatggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTAAAATCTTCCATAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCa                                                                                              <  1:247251/147‑1 (MQ=255)                                                  gggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCaaa                                                                                            >  2:491528/1‑150 (MQ=255)                                                          acTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaa                                                                                    >  1:770/1‑150 (MQ=255)                                                               aCATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCAt                                                                               >  2:374637/1‑150 (MQ=255)                                                               aCATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCAt                                                                               <  2:938045/150‑1 (MQ=255)                                                                 aTTACGCTTACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTCCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTg                                                                             >  1:162624/1‑150 (MQ=255)                                                                 aTTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTg                                                                             <  2:803938/150‑1 (MQ=255)                                                                  ttACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCTAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTgg                                                                            >  2:298318/1‑150 (MQ=255)                                                                       cTAACGTCAATAAAACAATCACATTTTTCTGTGTCTGGCGCTAAAAGATTACAATGTCACAATTCTCAAAATCAGAGTAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAATCTTAAAAAGCATTGGGCggt                                                                       <  2:660312/150‑2 (MQ=255)                                                                            ggCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCgg                                                                  >  2:759652/1‑150 (MQ=255)                                                                             gCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGc                                                                 >  2:356341/1‑150 (MQ=255)                                                                             gCACTAAAACCATCACATTTTTCTGTGACTCGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCCGAGGAGTATAGCTACTGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCTCAAAGCTGAAAAAGCATTGGGCGGAATCGGc                                                                 <  2:875185/150‑1 (MQ=255)                                                                               aCTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAAAGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGccc                                                               >  1:699316/1‑150 (MQ=255)                                                                                  aaaaCCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgc                                                            >  2:889224/1‑150 (MQ=255)                                                                                   aaaCCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGGCACAATTCGCAAAATCACAGGAGTATTGCTACTGGCAAGAGACGCGACGAAGCTGGAAGCAACGGCCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgcg                                                           <  2:609498/150‑1 (MQ=255)                                                                                   aaaCCATCACATTTTTCTGTGACTGGCGCTACAAACTTCCATAGTCACAATTCTCAAAATCAGAGGAGAATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACCGTCGCAAAGCTGAAAACGCATTGGGCGGAATCGGCCCcgcg                                                           >  2:915244/1‑150 (MQ=255)                                                                                    aaCCATCACATTTTTCTGAGACTGGCGCTACACTCTTCCAAAGTCACAATTCTCAAAATCAGAGGCGTATTGCTAATGGCATGAGACGCGACGAAGCTTGAAGCAACGGTCGAAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgcgc                                                          <  2:125033/150‑1 (MQ=255)                                                                                          cacaTTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATg                                                    >  2:260016/1‑150 (MQ=255)                                                                                                   cTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCa                                           >  2:477309/1‑150 (MQ=255)                                                                                                   cTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCa                                           >  2:808268/1‑150 (MQ=255)                                                                                                    tgtgACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCAt                                          <  2:172759/150‑1 (MQ=255)                                                                                                        aCTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAg                                      <  2:420655/150‑1 (MQ=255)                                                                                                          tGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAgcg                                    >  1:843346/1‑150 (MQ=255)                                                                                                            gcgcTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACCAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGcc                                  <  2:803504/150‑1 (MQ=255)                                                                                                                 aCAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGAccc                             <  1:253352/150‑1 (MQ=255)                                                                                                                 aCAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCAAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGAccc                             <  1:764563/150‑1 (MQ=255)                                                                                                                      cTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccggcc                        <  2:895721/150‑1 (MQ=255)                                                                                                                      cTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccggcc                        <  2:22879/150‑1 (MQ=255)                                                                                                                                cacaATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCg              >  1:806735/1‑150 (MQ=255)                                                                                                                                 acaAATCTCAATATCTCAGGAGCATTGCTAATGGCAAGAGACTCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAACCCAGTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCTGGCGa             <  1:530857/150‑1 (MQ=255)                                                                                                                                   aaTTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGata           <  1:614910/150‑1 (MQ=255)                                                                                                                                       ctcAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGGCGCAAAGCTGAAAAAGCATTGGGCGGAAGCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCg       <  1:465172/150‑1 (MQ=255)                                                                                                                                            aaTCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTg  >  2:150701/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                               CCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTG  >  NC_000913_3_bme_pgi/4233600‑4233888

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