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%	93.9 / NA	21	intergenic (+315/+57)	ycgI/minE	protein YcgI/Z‑ring positioning protein MinE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (12/9);  total (12/9)

AGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTT  >  NC_000913/1224145‑1224298                                                                              |                                                                             aGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGaaaa                                                                           >  2:228020/1‑81 (MQ=255) aGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGaaaa                                                                           >  1:228025/1‑81 (MQ=255)    gggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGaaaact                                                                         >  2:228013/1‑78 (MQ=255)     ggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGaaaactc                                                                        >  1:228043/1‑77 (MQ=255)           ggAAATAGCAATGAGGAGTATCAGCAATAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCaatg                                                                   >  1:228037/1‑76 (MQ=255)              aaTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAgaga                                                             <  1:228000/81‑1 (MQ=255)               aTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGAt                                                            <  1:228001/81‑1 (MQ=255)                 aGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATaa                                                          >  1:228038/1‑81 (MQ=255)                 aGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATaa                                                          >  1:228027/1‑81 (MQ=255)                  gCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAAc                                                         <  2:228002/81‑1 (MQ=255)                     aTGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAActct                                                      >  2:228029/1‑81 (MQ=255)                       gaggagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGc                                                    >  2:228041/1‑81 (MQ=255)                        aggagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGcc                                                   <  1:228003/81‑1 (MQ=255)                              aTCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTg                                                >  1:228023/1‑78 (MQ=255)                              aTCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACACTGCCTTg                                                >  1:228028/1‑78 (MQ=255)                               tCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAg                                             <  1:228005/80‑1 (MQ=255)                               tCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAg                                             <  1:228004/80‑1 (MQ=255)                                    aaGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAAt                                       >  1:228042/1‑81 (MQ=255)                                       aaTACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgcg                                    >  1:228035/1‑81 (MQ=255)                                                 gcTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGc                          >  2:228030/1‑81 (MQ=255)                                                          caACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAtt                 <  2:228006/81‑1 (MQ=255)                                                           aaCGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAttt                >  2:228046/1‑81 (MQ=255)                                                                        gggATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCt   >  2:228032/1‑81 (MQ=255)                                                                         ggATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCtt  <  2:228008/81‑1 (MQ=255)                                                                         ggATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCtt  <  1:228007/81‑1 (MQ=255)                                                                              |                                                                             AGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTT  >  NC_000913/1224145‑1224298

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