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

AAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCT  >  NC_000913/1224144‑1224294                                                                               |                                                                         aaGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGaa                                                                         >  2:222714/1‑80 (MQ=255)  aGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGaaaa                                                                       >  1:222711/1‑81 (MQ=255)   gggggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGaaaac                                                                      <  1:222682/81‑2 (MQ=255)    ggggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGaaaact                                                                     <  1:222683/81‑3 (MQ=255)      ggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGaaaactc                                                                    >  1:222715/1‑77 (MQ=255)       gAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCaa                                                                 <  2:222684/81‑1 (MQ=255)        aGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAg                                                                >  1:222707/1‑81 (MQ=255)             gAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCag                                                            <  2:222685/80‑1 (MQ=255)                     aaTGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAActc                                                   >  1:222718/1‑81 (MQ=255)                      aTGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAActct                                                  >  2:222699/1‑81 (MQ=255)                           gagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCtt                                             <  2:222686/81‑1 (MQ=255)                           gagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCtt                                             >  2:222712/1‑81 (MQ=255)                             gTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGa                                           <  1:222687/81‑1 (MQ=255)                                tCAGCAAGGATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAg                                         <  2:222688/80‑1 (MQ=255)                                     aaGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAAt                                   <  2:222689/81‑1 (MQ=255)                                              cgccgcTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTa                          >  1:222720/1‑81 (MQ=255)                                               gccgcTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTAc                         <  2:222690/81‑1 (MQ=255)                                                   cTTTACCACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCg                     <  1:222691/81‑1 (MQ=255)                                                        ccACAACGTGGATGAGAGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCt                 <  2:222692/80‑1 (MQ=255)                                                                     gagaGGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttc   >  2:222719/1‑81 (MQ=255)                                                                      agagGGATG‑AAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttct  <  1:222693/81‑1 (MQ=255)                                                                               |                                                                         AAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCT  >  NC_000913/1224144‑1224294

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