Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,212,233	A→G	intergenic (‑230/+470)	iraM ← / ← ycgX	RpoS stabilzer during Mg starvation, anti‑RssB factor/DUF1398 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,212,233	0	A	G	100.0%	78.5 / NA	23	intergenic (‑230/+470)	iraM/ycgX	RpoS stabilzer during Mg starvation, anti‑RssB factor/DUF1398 family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/13);  total (10/13)

AATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATACATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGATAATAATGAAT  >  NC_000913/1212110‑1212339                                                                                                                            |                                                                                                           aaTGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGa                                                                                             <  2:1270112/139‑1 (MQ=255)      tgGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATc                                                                                        >  2:1656345/1‑139 (MQ=255)        gCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCt                                                                                      <  1:1656345/139‑1 (MQ=255)                                cAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGc                                                              <  2:1148509/139‑1 (MQ=255)                                   ggAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTAc                                                           >  2:1044600/1‑139 (MQ=255)                                             tGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCaa                                                 >  1:694325/1‑139 (MQ=255)                                                   ttAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCt                                           >  1:585240/1‑139 (MQ=255)                                                   ttAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCt                                           <  2:694325/139‑1 (MQ=255)                                                    tAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTa                                          <  1:870322/139‑1 (MQ=255)                                                        cAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACa                                      <  2:950575/139‑1 (MQ=255)                                                                            gCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTtatat                  <  2:585240/139‑1 (MQ=255)                                                                                 aaaTTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGa             <  2:687741/139‑1 (MQ=255)                                                                                      attGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATT‑AAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCttt                        >  1:484980/1‑122 (MQ=255)                                                                                      attGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATT‑AAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCttt                        <  2:484980/122‑1 (MQ=255)                                                                                       ttGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTg                             <  2:9754/117‑1 (MQ=255)                                                                                       ttGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTg                             >  1:9754/1‑117 (MQ=255)                                                                                            ttCGTTTTAGATTTGTGTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTATACTTTATGGTTGCCTTTTTATATAACGATAATAATGAat  >  1:330165/1‑139 (MQ=255)                                                                                              cGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCt                                                                    <  2:857879/71‑1 (MQ=255)                                                                                              cGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCt                                                                    >  1:857879/1‑71 (MQ=255)                                                                                              cGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTg                             <  2:3024/110‑1 (MQ=255)                                                                                              cGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTg                             >  1:3024/1‑110 (MQ=255)                                                                                              cGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGataata        <  1:1022459/131‑1 (MQ=255)                                                                                              cGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGataata        >  2:1022459/1‑131 (MQ=255)                                                                                                                            |                                                                                                           AATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATACATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGATAATAATGAAT  >  NC_000913/1212110‑1212339

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