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%	108.6 / NA	26	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 . (9/17);  total (9/17)

TCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAG  >  NC_000913_3_bme_pgi/4235259‑4235543                                                                                                                                                     |                                                                                                                                         tCAATGGCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCCAGGGCTTGTCCCGCATCTGCATGCGCGAAGGAAAGCTTGa                                                                                                                                          <  1:268165/149‑1 (MQ=255) tCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTACGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGCCGCTTCCACGCAAGGGATGGTGGCGCGTCTGCGTGCCGGAAGAGAAGCTTGa                                                                                                                                          >  1:20916/1‑149 (MQ=255)            gATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGat                                                                                                                              <  2:281464/149‑1 (MQ=255)              tCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGGGGTCGCGACGCTACCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCAtcgtcgcact                                                                                                                            >  1:518591/1‑145 (MQ=255)                         tGGAACTGTGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGTCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGAt                                                                                                                 <  1:41314/149‑1 (MQ=255)                              cTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGc                                                                                                            <  1:156892/149‑1 (MQ=255)                                     aGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcg                                                                                                     >  1:146302/1‑149 (MQ=255)                                      gggaCTGGCGAGGGAGCTTGTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGATGCTTCCACGCGAGGGCTTGTCGCGTATCTGCATGCGAGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgc                                                                                                    <  2:137937/146‑1 (MQ=255)                                           ctgCGACGGATCTTGGGCCGGTGGTAGCGGGCAAGGAAGGGTGGGGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTTCATGCGCGAAGGAACGCTTG‑ATCATCGTCGATATGTAGGCCCGATAAGGCGTTCACGCCGCAt                                                                                               <  2:480340/147‑1 (MQ=255)                                               gACGGAACTTTTGCCGGTTGTATCGGGCAAGGAATGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                           <  2:423582/149‑1 (MQ=255)                                               gACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                           >  1:7098/1‑149 (MQ=255)                                                       ttttGCCGGGTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTTCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAt                                                                                   <  2:247605/149‑1 (MQ=255)                                                                ttGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCAGCTACAGGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCGGCATCCGGCACCCGATccccgcggc                                                                          >  2:219909/1‑143 (MQ=38)                                                                         gCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTc                                                                 >  1:65355/1‑149 (MQ=255)                                                                                     cgggTGGGCGCGACGTTTCCCGGGAAGGGCTTGTCGCGGATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCCCGGCGCAGCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                     <  2:105399/146‑1 (MQ=25)                                                                                     agcggggTCGCGCTGCTTGCAGGGAGTGGCGTGGCGCGCTGGTGCATGTGCGAAGTATAGCTTG‑ATCATCGTCGGTATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTTGCGTCTTATCAg                                                     <  2:106285/144‑1 (MQ=12)                                                                                       aGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc                                                   <  1:294050/149‑1 (MQ=255)                                                                                                     ttCAACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGCTATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGc                                     <  1:422231/149‑1 (MQ=255)                                                                                                             aaGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGGCGCGGTCGCGTCTTATCAGGCCTGCAGGTGGATGCCGATATGt                             >  1:237058/1‑149 (MQ=255)                                                                                                                ggCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACGGATTCCTGATGCGCCGCCGGCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACa                          <  1:487515/149‑1 (MQ=37)                                                                                                                 ccTTGTCGCGCATCGGCGTGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACAt                         <  2:458517/148‑1 (MQ=255)                                                                                                                                 cATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTaa         >  1:189273/1‑149 (MQ=255)                                                                                                                                     agcgAAGGAACGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACa     <  1:91482/148‑1 (MQ=255)                                                                                                                                        gAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAg  <  1:455669/149‑1 (MQ=255)                                                                                                                                          aGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATAAGGCCTACAGGTCGATGCCGATATg                              <  2:36036/119‑1 (MQ=255)                                                                                                                                          aGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAACGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATg                              >  1:36036/1‑119 (MQ=255)                                                                                                                                           ggAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACa                                              >  2:154603/1‑102 (MQ=255)                                                                                                                                           ggAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACa                                              <  1:154603/102‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                         TCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAG  >  NC_000913_3_bme_pgi/4235259‑4235543

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