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%	129.3 / NA	36	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 (17/19);  total (17/19)

GGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATACATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGATAATAATGAATATAAGCATGACATGAGAATA  >  NC_000913/1212116‑1212359                                                                                                                      |                                                                                                                               ggCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATcc                                                                                                           >  1:431525/1‑139 (MQ=255)       tataAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAAt                                                                                                     >  2:368106/1‑139 (MQ=255)        ataaaataaaATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAAtt                                                                                                    <  2:1118195/139‑1 (MQ=255)          aaaataaaatTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTaa                                                                                                  <  2:1093647/139‑1 (MQ=255)            aataaaatTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTaaaa                                                                                                >  2:148632/1‑139 (MQ=255)               aaaatTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACAt                                                                                             <  2:431525/139‑1 (MQ=255)               aaaatTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACAt                                                                                             <  1:980720/139‑1 (MQ=255)                 aatTTAAAGCAAGGAATGTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTAGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATg                                                                                                                      >  1:697677/1‑112 (MQ=255)                 aatTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATg                                                                                                                      <  2:697677/112‑1 (MQ=255)                  atTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGc                                                                                          <  1:560982/139‑1 (MQ=255)                                         gATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGc                                                                                          >  1:322879/1‑116 (MQ=255)                                         gATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGc                                                                                          <  2:322879/116‑1 (MQ=255)                                         gATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCt                                                                                        >  1:130502/1‑118 (MQ=255)                                         gATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCt                                                                                        <  2:130502/118‑1 (MQ=255)                                            gTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAAt                                                                    >  1:729045/1‑135 (MQ=255)                                            gTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAAt                                                                    <  2:729045/135‑1 (MQ=255)                                            gTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGc                                                                >  1:364049/1‑139 (MQ=255)                                             ttAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCt                                                               <  1:148632/139‑1 (MQ=255)                                             ttAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCt                                                               <  2:364049/139‑1 (MQ=255)                                                 aCAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGa                                                                                                                 <  2:722056/85‑1 (MQ=255)                                                 aCAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGa                                                                                                                 >  1:722056/1‑85 (MQ=255)                                                    aaaTGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCt                                                                                        <  2:414004/107‑1 (MQ=255)                                                    aaaTGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCt                                                                                        >  1:414004/1‑107 (MQ=255)                                                     aaTGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGc                                                                                          >  1:1218598/1‑104 (MQ=255)                                                     aaTGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGc                                                                                          <  2:1218598/104‑1 (MQ=255)                                                       tGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATg                                                     >  2:926431/1‑139 (MQ=255)                                                                   aaaGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTta                                         <  1:926431/139‑1 (MQ=255)                                                                                  tGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTtatat                                      <  1:81573/127‑1 (MQ=255)                                                                                  tGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTtatat                                      >  2:81573/1‑127 (MQ=255)                                                                                  tGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGATaataat                          >  2:1141648/1‑139 (MQ=255)                                                                                        cGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATg                                                                                                                      <  2:350501/41‑1 (MQ=255)                                                                                        cGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATg                                                                                                                      >  1:350501/1‑41 (MQ=255)                                                                                                   tGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGc                                                                                  <  2:881972/66‑1 (MQ=255)                                                                                                   tGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGc                                                                                  >  1:881972/1‑66 (MQ=255)                                                                                                         gCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGATAATAATGAATATAAGCATGACATGAGAAt   <  1:453047/139‑1 (MQ=255)                                                                                                          cTATAATGTTATGCATTCAAATGACTGAACATTCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGATAATAATGAATATAAGCATGACATGAGAATa  >  1:724109/1‑139 (MQ=255)                                                                                                                      |                                                                                                                               GGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATACATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGATAATAATGAATATAAGCATGACATGAGAATA  >  NC_000913/1212116‑1212359

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