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%	192.7 / NA	44	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 . (21/23);  total (21/23)

TCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGA  >  NC_000913_3_bme_pgi/4235259‑4235548                                                                                                                                                     |                                                                                                                                              tCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGa                                                                                                                                               <  2:111707/149‑1 (MQ=255)       ccTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCAtcg                                                                                                                                        >  2:409829/1‑149 (MQ=255)        cTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCACGCGCGACGGAAAGCTTG‑ATCAtcgt                                                                                                                                       >  2:519132/1‑149 (MQ=255)         ttCGATCAAGGGGGTGGTGAACGGGGTAAGGAACTGGCGACGGAATTTGTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGTGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCAtcgtc                                                                                                                                      <  2:38885/149‑1 (MQ=255)           ggATCAAGGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATAGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCTAAGGAAAGCTTG‑ATCATCTTCGa                                                                                                                                    <  2:722775/148‑1 (MQ=255)                    gggCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATAGGTAGGc                                                                                                                           >  1:561448/1‑149 (MQ=255)                      cttccgAACTGCATAAGGAACTGGCGACGGAACTTTTGCCGGTTCTATCGGGCAAGGAAGGGGCGCGTGGTCGCGACGCTTCTACGCAAGGGCTTGTCGCGCATCTGCATGCGCCAAGGAAAGCTTG‑ATCATCGTTGATCTGTAGGCCg                                                                                                                         <  2:623724/144‑1 (MQ=255)                              cttggTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACTCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCAACCTAAATCACG‑TTTATCGTCGATATGTAGGCCGGATAAGGc                                                                                                                 <  1:584260/146‑1 (MQ=255)                                      ggAACTGGGGAGGCAACTGTTGACGTTCGTATCGTGCAAGGAAGGGCCGAGTGGGCGCGACGCTTCCGCGCAAGGGCTGGTCGCGCAGCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgc                                                                                                         <  2:423902/149‑1 (MQ=255)                                                     aCTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCTCTACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCg                                                                                          >  1:615065/1‑149 (MQ=255)                                                       ttttGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAt                                                                                        >  1:630862/1‑149 (MQ=255)                                                           gCCGGTGGTAGCGGGCACGGAAGGGGCGAGGGGCCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCt                                                                                    <  2:584878/149‑1 (MQ=38)                                                              ggTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGAt                                                                                 <  1:405926/149‑1 (MQ=255)                                                                   tATCGGGCAAGGAAGGGGCGAGTGGTCGCGAAGCATCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAc                                                                            <  1:732776/149‑1 (MQ=255)                                                                   tATCGGGCAAGGAAGGGGAGAGTGGTCGCGACGGTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAc                                                                            <  1:373690/149‑1 (MQ=255)                                                                         gCAAGGAGGGGGCGAGGGGTCGCGACGCTCCAACGCAAGGGCTTGTCGCGCATCTGCACGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTc                                                                      <  2:55798/149‑1 (MQ=38)                                                                         gCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGAGTTAAGGCCGCATCGGGCAACCGATGCCTGATGCGACGCGGTc                                                                      >  1:415938/1‑149 (MQ=255)                                                                             ggaaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGCAAGCTTG‑ATCATCGTCGATATGTAGGCCGCATAAAGCGTTCACGCCGCATCCCGCAACCACAGCCTGCTGCGACGCGGTAGCGt                                                                  >  2:28676/1‑149 (MQ=255)                                                                             ggaaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCCTCCGGCACCCGATGCCGGATGTGACGCGGGCGCGt                                                                  >  1:309754/1‑149 (MQ=38)                                                                               aaggGGCGAGCGGTCGCGACGCTTCCACGCAAGGGCTGGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCt                                                                <  2:346690/149‑1 (MQ=255)                                                                                    gggAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                           <  2:577393/147‑1 (MQ=255)                                                                                     cGAGTGGTCGCGACGCTTCCAAGCAAGCGCTTGTTCTGTTTCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGGAAGCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                          <  1:19415/149‑1 (MQ=25)                                                                                      caGTCGTCGCGACGCTGCCACGCAAGCGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATAGGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGTGTCTTATCAga                                                         <  2:296538/148‑2 (MQ=37)                                                                                           gtagcgAGGCTTCCATGCAGGGGCTGTTCGGGCATCTGCGTGGGCGAAGGAAAGGTTG‑ATCAGTGTCGATATGTAGGCCGGATAAGGGGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAGGCGGTCGCTTCTTATCAGGCTTAc                                                    <  2:53143/146‑1 (MQ=25)                                                                                           gTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTAc                                                    <  2:698867/149‑1 (MQ=255)                                                                                            tCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGCCTTATCCGGCCTAcc                                                   >  2:340402/1‑148 (MQ=255)                                                                                                 aCGCTTCCACGCAAGGGCGTGTCGCGCAGCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCg                                              <  2:518700/149‑1 (MQ=255)                                                                                                   gCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGGCGCGTCTTATCAGGCCTACAGGTCGAt                                            >  1:261322/1‑149 (MQ=255)                                                                                                    cTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATg                                           >  1:42625/1‑149 (MQ=255)                                                                                                         aCGCAAGGGCTTGTCGCGCAACTGCATACGCGAAGGAAACCTTG‑ATCTGTTCCGCCAGCCGGGCCGAACAGGGCGTAGAGGCCGAATCCGGCAAACTATGCCTGAAGCGACGCGGTCGGGTCGTCCGATGACTACAGGCCGAAGCAGat                                      >  2:91248/1‑149 (MQ=255)                                                                                                           gCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCCTCTTTTCAGACCTACAGGTCGATGCCGatat                                    >  1:730312/1‑149 (MQ=255)                                                                                                            cAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATg                                   >  1:336668/1‑149 (MQ=255)                                                                                                            cAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATg                                   >  1:336673/1‑149 (MQ=255)                                                                                                               gggCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTAc                                <  2:42625/149‑1 (MQ=255)                                                                                                               gggCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTAc                                >  1:496701/1‑149 (MQ=255)                                                                                                                tttccttCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACa                               >  2:258029/7‑149 (MQ=255)                                                                                                                 gCTTGCCGCGCATCTGCATGCCCGAAGGAAAGCTTG‑ATCATCCTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACAt                              <  1:498933/149‑1 (MQ=255)                                                                                                                  cgtGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATc                             <  1:226617/147‑1 (MQ=255)                                                                                                                            aTCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCa                   >  1:652942/1‑149 (MQ=255)                                                                                                                             tCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCaa                  <  2:229275/149‑1 (MQ=255)                                                                                                                              cTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAAt                 <  1:128165/149‑1 (MQ=255)                                                                                                                                   tGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATa            >  1:336532/1‑149 (MQ=255)                                                                                                                                      gcgAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCCCGTCTTAGCACGCCTACACGTCGGTGCCGAAATGTACATCGCATTAGGCAATTGATACAt         >  1:558256/1‑149 (MQ=25)                                                                                                                                       cgAAGGAAAGCTTG‑ATCATCGTCGAGATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCGATTAATACATa        <  1:494712/149‑1 (MQ=255)                                                                                                                                             aaaGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGCTATGTACATCGTATTCGGCAATTAATACATGGCACga  >  2:627329/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                              TCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGA  >  NC_000913_3_bme_pgi/4235259‑4235548

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