Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	89,890	T→A	M86K (ATG→AAG)	mraZ →	DNA‑binding transcriptional repressor MraZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	89,890	0	T	A	100.0%	138.8 / NA	38	M86K (ATG→AAG)	mraZ	DNA‑binding transcriptional repressor MraZ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (14/24);  total (14/24)

ATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTG  >  NC_000913/89813‑89967                                                                              |                                                                              aTCGCGTCTGTCGAGCATGAACCCGGTTGAACGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgaagg                                                                             <  2:20100/80‑1 (MQ=255)    gcgTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGaaggatgg                                                                         >  2:20151/1‑81 (MQ=255)     cgTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGAAGGATGgc                                                                        <  2:20101/81‑1 (MQ=255)        cTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGAAGGATGgcg                                                                       >  1:20146/1‑79 (MQ=255)             gAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTc                                                                 >  1:20143/1‑80 (MQ=255)             gAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTc                                                                 <  2:20103/80‑1 (MQ=255)             gAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTc                                                                 <  1:20102/80‑1 (MQ=255)                      cccGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAAt                                                       <  2:20104/81‑1 (MQ=255)                          gTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATcgcg                                                   <  2:20105/81‑1 (MQ=255)                           ttGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCgcgc                                                  <  2:20106/81‑1 (MQ=255)                             gAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCgcgc                                                  >  2:20145/1‑79 (MQ=255)                             gAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGcc                                                 <  2:20107/80‑1 (MQ=255)                                 gctgtgtGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTAc                                            <  2:20108/78‑1 (MQ=255)                                 gcCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTAc                                            >  1:20147/1‑81 (MQ=255)                                  ccGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACt                                           <  2:20109/81‑1 (MQ=255)                                       tgCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCgg                                       <  2:20110/80‑1 (MQ=255)                                        gCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGc                                      >  2:20157/1‑80 (MQ=255)                                         cAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCa                                     <  1:20111/80‑1 (MQ=255)                                          aGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCaa                                    >  2:20152/1‑80 (MQ=255)                                          aGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAAc                                   <  2:20114/81‑1 (MQ=255)                                          aGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAAc                                   <  2:20113/81‑1 (MQ=255)                                          aGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAAc                                   <  2:20112/81‑1 (MQ=255)                                            cgcCTACTGTTAGGTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACAt                                 <  1:20115/81‑1 (MQ=255)                                            cgcCTACTGTTAGGTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACAt                                 <  1:20116/81‑1 (MQ=255)                                             gcCTACTGTTAGGTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATg                                <  2:20117/81‑1 (MQ=255)                                               cTACTGTTAGGTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGcc                              >  2:20153/1‑81 (MQ=255)                                                  cTGTTAGGTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCggg                           <  1:20119/81‑1 (MQ=255)                                                  cTGTTAGGTCACGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCg                             <  2:20118/79‑1 (MQ=255)                                                    gTTAGGTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCt                         >  2:20150/1‑81 (MQ=255)                                                    gTTAGGTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCt                         >  1:20162/1‑81 (MQ=255)                                                        ggTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGAc                      <  2:20122/80‑1 (MQ=255)                                                        ggTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACa                     <  1:20120/81‑1 (MQ=255)                                                        ggTCATGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACa                     <  1:20121/81‑1 (MQ=255)                                                            aTGCCAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACaa                    >  1:20154/1‑78 (MQ=255)                                                                cAGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGt             >  1:20168/1‑81 (MQ=255)                                                                 aGCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGt             >  1:20148/1‑80 (MQ=255)                                                                  gCGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGa           >  2:20149/1‑81 (MQ=255)                                                                   cGAATGTCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGAt          <  2:20123/81‑1 (MQ=255)                                                                         tCAGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATGCTgg     <  2:20124/80‑1 (MQ=255)                                                                           aGAAGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTg  >  1:20161/1‑81 (MQ=255)                                                                              |                                                                              ATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTG  >  NC_000913/89813‑89967

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