Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	4,235,408	Δ1 bp	intergenic (‑2175/‑499)	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,235,407	0	A	.	100.0%	224.2 / NA	50	intergenic (‑2174/‑500)	lysC/yjbE	lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (27/23);  total (27/23)

ATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCAC  >  NC_000913_3_bme_pgi/4235262‑4235546                                                                                                                                                  |                                                                                                                                            aTGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGatca                                                                                                                                          >  1:177475/1‑146 (MQ=255)     cTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTAGTCTGCGGCAGTGAAGGGGCGAGCGTTCGCGACGCTTCCAAGCAAGGGCTGGGCGCGGATCTGCATGCGCGAAGGAAAGCTTG‑ATCAtcgt                                                                                                                                     >  1:187210/1‑149 (MQ=255)      ttCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCAtcgtc                                                                                                                                    >  1:401005/1‑149 (MQ=255)      ttCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTGGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCAtcgtc                                                                                                                                    <  1:624350/149‑1 (MQ=255)         gATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTGGCGAAGCTGGATGCGCGAAGGAAAGCTTG‑ATCCTCGTCGat                                                                                                                                 >  1:622212/1‑149 (MQ=255)          aTCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGCGTGGTCGCGACGCTTCCTCGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGata                                                                                                                                >  2:88001/1‑149 (MQ=255)           tCAATGGGGCGTGGAACGGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGAGAGTGGTCGCGACGCTTCCAAGCACGGGCTTGTCGCGCATCTGCATGCGCGGAGGGAAGCTTG‑ATCATCGTCGctat                                                                                                                               >  1:268320/1‑149 (MQ=255)               tGGGGCGTGGATCTGTGTAAGGAGCTGGCGACGGAACTTTTGCCTGTTGTATCGGGCAAGGAAGGGGCGGGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCGGCATGCGCGAAGGAAAGCTTG‑ATCATCATCGATATGTAg                                                                                                                           <  2:475586/149‑1 (MQ=255)                ggggCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAgg                                                                                                                          >  1:500197/1‑149 (MQ=255)                  ggCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGcc                                                                                                                        <  1:508652/149‑1 (MQ=255)                    cGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAGACTTG‑ATCATCGTCGATACGTAGGCTgg                                                                                                                      >  2:493294/1‑149 (MQ=255)                              ggTAGGGAACTGGCGAGGGAAGTTGTGCGGGGTGTATCGGGAAAGGAAGGGGCGAGTGGTCGCGTCGCTTCCACGCAAGGGTTTGTCGCGCATCTGCATGCGCGAAGGAGAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                                                            <  2:68541/149‑1 (MQ=255)                              ggTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                                                            >  1:801420/1‑149 (MQ=255)                              ggTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGACGCGCATCTGCATGAGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                                                            <  1:707137/149‑1 (MQ=255)                                tAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGACGGGGCGAGGGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAAGCCGGAAAAGGCGtgcc                                                                                                          >  2:29098/1‑146 (MQ=255)                                           cGACGGAACTTTTGCAGATTGTGTGGGGCATGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                               <  2:603009/149‑1 (MQ=39)                                               ggAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTCGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                                           >  1:591942/1‑149 (MQ=255)                                                  aCTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGACCAGGAAAGCTTG‑ATAATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACCCg                                                                                        >  1:589971/1‑149 (MQ=38)                                                        gCCGGTTGTATGGGGCAAGGACGGGGCGAGTGGTCGGGACGCTTCCACGCAAGGGCTTGGCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCt                                                                                  <  2:327519/149‑1 (MQ=39)                                                         ccGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTg                                                                                 >  2:156728/1‑149 (MQ=255)                                                          cGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCCCGCCGCATCCGGCAACCGATTCCTGa                                                                                >  1:30855/1‑149 (MQ=255)                                                           ggtggtATCGGGCAAGCAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGAt                                                                               <  2:744978/149‑1 (MQ=255)                                                                 aTCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcg                                                                         >  1:463500/1‑149 (MQ=255)                                                                 aTCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCAATGCCTGATGCGAcg                                                                         >  1:428520/1‑149 (MQ=255)                                                                     ggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGt                                                                     >  2:72235/1‑149 (MQ=255)                                                                        atagaaggGGTGGGGGGGCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTTGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgc                                                                  <  2:561107/146‑1 (MQ=37)                                                                         aggaaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGACGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgcc                                                                 >  2:466304/1‑148 (MQ=255)                                                                          ggaaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGt                                                                <  1:156728/149‑1 (MQ=255)                                                                           gatggggCGGGTTGTCGCGACGCTTCCACGCAATGGCTTGTCGCGCATCTGCATTCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGTATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTc                                                               <  2:232548/146‑1 (MQ=38)                                                                             aggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCtt                                                             <  1:72235/149‑1 (MQ=255)                                                                               gggCGAGAAGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTAt                                                           >  2:496475/1‑149 (MQ=255)                                                                                    aGTGGCCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc                                                      <  1:159471/149‑1 (MQ=255)                                                                                        gTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTAc                                                  >  2:318493/1‑149 (MQ=255)                                                                                        gTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTAc                                                  >  1:47142/1‑149 (MQ=255)                                                                                         tCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGCCGCGGTCGCGTCTTATCAGGCCTACa                                                 >  1:692258/1‑149 (MQ=255)                                                                                         tCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACa                                                 <  1:669106/149‑1 (MQ=255)                                                                                          ggcgACGCTTCCACGGAAGGGCTTGGCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATACGGCGTTCTCGTCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAg                                                <  1:425420/148‑1 (MQ=37)                                                                                          cgcgACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAg                                                >  2:228454/1‑149 (MQ=255)                                                                                            cgACGCTTCCACGGAAGGGCTTGTCTCGCAGCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                              <  2:622212/149‑1 (MQ=255)                                                                                            cgACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                              <  1:389644/149‑1 (MQ=255)                                                                                              aCGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGGTTG‑ATCATCGTCGACATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCg                                            >  1:706776/1‑149 (MQ=255)                                                                                                  ttCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGATTCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGc                                        <  1:556972/149‑1 (MQ=255)                                                                                                  ttCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGc                                        >  1:697924/1‑149 (MQ=255)                                                                                                   tCCACGAAAGGGCTTGCCCCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGcc                                       <  1:325125/149‑1 (MQ=255)                                                                                                            gggCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGCTTCCGATATGTAc                              >  2:430695/1‑149 (MQ=255)                                                                                                                    cgcgCATCTGCATGCGCGAAGGAAAGCTTG‑AGCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTAtt                      <  2:504938/149‑1 (MQ=255)                                                                                                                       gcAGCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGGAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCgg                   <  1:494979/149‑1 (MQ=255)                                                                                                                              cctGCGCGAAGGGAAGCTTG‑ATCCTCGTGGATACGTAGGCCGGATAAGGCGTTCCCGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTaa            <  1:234409/147‑1 (MQ=25)                                                                                                                               aggcggGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCCGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAAt           <  2:401005/147‑1 (MQ=255)                                                                                                                                       aGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCa   >  1:418992/1‑149 (MQ=255)                                                                                                                                        ggAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGAACATCGTATTCGGCAATTAATACaatacac  >  2:612041/1‑143 (MQ=255)                                                                                                                                                  |                                                                                                                                            ATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCAC  >  NC_000913_3_bme_pgi/4235262‑4235546

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