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%	189.6 / NA	42	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 . (22/20);  total (22/20)

TCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGAT  >  NC_000913_3_bme_pgi/4235259‑4235549                                                                                                                                                     |                                                                                                                                               tCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGAAAAGCTTGa                                                                                                                                                >  1:277651/1‑149 (MQ=255)     tGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGAGGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTGGGGGATTTGCATGCGCGAAGGAAAGCTTGatcat                                                                                                                                            >  1:209341/1‑145 (MQ=255)      gcattCGACTAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑atcatc                                                                                                                                          <  2:318189/146‑1 (MQ=255)      gCCTTCGATCAAGGGGGCGTGGAACTGGGTAAGGAAGTGGCCGCGGATCTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑atcatc                                                                                                                                          <  1:363597/149‑1 (MQ=255)        cTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCAtcgt                                                                                                                                        <  2:471665/149‑1 (MQ=255)         ctCGATCAATGGGGCGGGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCAtcgtc                                                                                                                                       <  2:125308/148‑1 (MQ=255)              tCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGa                                                                                                                                     <  1:398488/146‑1 (MQ=255)              tCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCGTG‑ATCATCGTCGa                                                                                                                                     >  2:398488/1‑146 (MQ=255)               cAATGGGGCGGGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATg                                                                                                                                 <  1:181647/149‑1 (MQ=255)                aaTGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGt                                                                                                                                >  1:382763/1‑149 (MQ=255)                 aTGGGGCGGGTAAGTGGGTAAGGAACGGGGGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTa                                                                                                                               <  2:384612/149‑1 (MQ=255)                          ggAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGaga                                                                                                                      >  1:1437/1‑147 (MQ=255)                            aaCTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGATGTATCGGGCAAGGAAGGGCCGAGTGGTCGCGACGCTTCCACGGAAGGACTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGCAGTGCACGCCGGATCAg                                                                                                                    >  2:231590/1‑149 (MQ=255)                                  ggaaGGAACGGGCGAGGGAACTTGTGCGGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTc                                                                                                              <  2:266034/147‑1 (MQ=255)                                  gTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGGGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTc                                                                                                              >  1:346557/1‑149 (MQ=255)                                              cGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                  <  1:315657/149‑1 (MQ=255)                                                 cGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                               <  1:350938/149‑2 (MQ=255)                                                   gAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                             <  1:297736/149‑1 (MQ=255)                                                                ttGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAc                                                                                                            <  2:98118/121‑1 (MQ=255)                                                                ttGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAc                                                                                                            >  1:98118/1‑121 (MQ=255)                                                                     tCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgc                                                                           >  2:83410/1‑149 (MQ=255)                                                                       gggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGAcgcga                                                                         >  1:338609/1‑148 (MQ=255)                                                                         gCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTc                                                                       <  1:351202/149‑1 (MQ=255)                                                                                aggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCg                                                                                                                          >  1:185099/1‑91 (MQ=255)                                                                                aggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCg                                                                                                                          <  2:185099/91‑1 (MQ=255)                                                                                  gggCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCTTCTTAt                                                              >  1:82533/1‑149 (MQ=255)                                                                                    gCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTcgcatca                                                            <  2:338609/149‑7 (MQ=38)                                                                                      gAGTGGTCGCGACGCTTCCACGCAAAGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGTTGCGACGCGGTCGCGTCTCATCAgg                                                          >  1:478071/1‑149 (MQ=255)                                                                                            tCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACa                                                    >  1:291454/1‑149 (MQ=255)                                                                                              gcgACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCTCGTCTAATCAGACCTAAAgg                                                  >  1:132844/1‑149 (MQ=37)                                                                                               cgACGCTCCCACGCAAGGGCTTGCCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                 <  1:68010/149‑1 (MQ=255)                                                                                                gACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTc                                                >  1:316290/1‑149 (MQ=255)                                                                                                  ggCTTCCACGCAAGGGCGGGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCCCGCCGCATCCGGCAACCGAGGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGa                                              <  1:294195/148‑1 (MQ=37)                                                                                                   gCTTCCACGGAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGCCGTTCACGCCGCATCCGGCAACCGAGGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGAt                                             <  2:435998/149‑1 (MQ=255)                                                                                                      tCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGcc                                          >  1:173480/1‑149 (MQ=255)                                                                                                        cACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATTGTCGTTATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGa                                        <  1:436400/149‑1 (MQ=255)                                                                                                           acAAGGGCTTGGCGCTCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGatat                                     <  2:209341/148‑1 (MQ=255)                                                                                                            cAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATACCGCATGCGAAGCAGTCGCGTCTTATCAGGCAGACAGGACGATGCCGgtagg                                    >  1:466647/1‑147 (MQ=37)                                                                                                                              cTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGGTAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACCGGTCGATGCCGATATGTACATCGTATTCGGCAAt                  >  1:209086/1‑149 (MQ=255)                                                                                                                              cTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACGGGTCGATGCCGATTTGTACATAGTATTCGGCAAt                  >  1:25860/1‑149 (MQ=255)                                                                                                                                      gcgAAGGAAAGGTTG‑ATCATCGTCCATATGGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACAt          <  2:27354/149‑1 (MQ=255)                                                                                                                                      gcgAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACAt          >  2:20619/1‑149 (MQ=255)                                                                                                                                        gAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGAAATTAATACaaaa        >  1:192199/1‑146 (MQ=255)                                                                                                                                            gAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGCTGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTGGATGCCGGTATGTACATCGTATTCGGCAGTTTATACGTAGCCCg    >  1:248484/1‑149 (MQ=25)                                                                                                                                            gAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACCTAGCACg    >  1:284846/1‑149 (MQ=255)                                                                                                                                              aaGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCTCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgat  <  1:264578/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                               TCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGAT  >  NC_000913_3_bme_pgi/4235259‑4235549

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: ‑