Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,224,226	(A)5→4	intergenic (+319/+53)	ycgI → / ← minE	protein YcgI/Z‑ring positioning protein MinE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,224,222	0	A	.	100.0%	94.4 / NA	23	intergenic (+315/+57)	ycgI/minE	protein YcgI/Z‑ring positioning protein MinE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (8/15);  total (8/15)

AGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGC  >  NC_000913/1224145‑1224296                                                                              |                                                                           aGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGaaaa                                                                         >  2:218440/1‑81 (MQ=255)     ggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGaaaactc                                                                      >  1:218409/1‑77 (MQ=255)      gAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCaa                                                                   <  1:218410/81‑1 (MQ=255)              aaTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCa                                                               >  1:218434/1‑77 (MQ=255)                tAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGAt                                                          >  1:218429/1‑80 (MQ=255)                  gCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAAc                                                       >  1:218433/1‑81 (MQ=255)                     aTGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAActct                                                    <  1:218411/81‑1 (MQ=255)                          gagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCtt                                               <  2:218412/81‑1 (MQ=255)                               tCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGa                                          <  1:218413/81‑1 (MQ=255)                               tCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGa                                          <  2:218414/81‑1 (MQ=255)                                cAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGa                                          >  2:218436/1‑80 (MQ=255)                                  gCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATaa                                       <  2:218415/81‑1 (MQ=255)                                   cAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATaaa                                      >  2:218437/1‑81 (MQ=255)                                     aGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATg                                    >  2:218441/1‑81 (MQ=255)                                       aaTACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgcg                                  <  1:218416/81‑1 (MQ=255)                                               ccgcTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACa                          <  1:218418/81‑1 (MQ=255)                                               ccgcTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACa                          <  2:218417/81‑1 (MQ=255)                                                cgcTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAg                         >  1:218435/1‑81 (MQ=255)                                                    ttACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCggg                     >  2:218443/1‑81 (MQ=255)                                                      aCCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGc                    <  1:218419/80‑1 (MQ=255)                                                                     agagGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttct    <  1:218423/81‑1 (MQ=255)                                                                     agagGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttct    <  1:218422/81‑1 (MQ=255)                                                                     agagGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttct    <  1:218421/81‑1 (MQ=255)                                                                     agagGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttc     <  2:218420/80‑1 (MQ=255)                                                                       agGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGc  <  1:218424/81‑1 (MQ=255)                                                                              |                                                                           AGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGC  >  NC_000913/1224145‑1224296

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