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%	112.9 / NA	33	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/16);  total (17/16)

TATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATACATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGATAATAATGAATATAAGCATGACATGAGAATAAGGTT  >  NC_000913/1212122‑1212364                                                                                                                |                                                                                                                                    tataAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAAt                                                                                                          <  2:1502474/139‑1 (MQ=255)    aaaataaaatTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTaa                                                                                                       <  2:1496508/139‑1 (MQ=255)      aataaaatTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACNTCCNGTANTNaaan                                                                                                     >  1:1496508/1‑138 (MQ=255)      aataaaatTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTaaaa                                                                                                     >  1:1363678/1‑139 (MQ=255)      aataaaatTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGATTAGCTATAATTTTATGCATTCAAATGACTGAACATCCTATAATTaaaa                                                                                                     >  2:1030820/1‑139 (MQ=255)              ttAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCt                                                                                             <  1:5688/139‑1 (MQ=255)                       ggAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTAc                                                                                    >  1:1112367/1‑139 (MQ=255)                                   gATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCt                                                                                                               >  1:454608/1‑100 (MQ=255)                                   gATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCt                                                                                                               <  2:454608/100‑1 (MQ=255)                                   gATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAAtt                                                                        >  2:1085782/1‑139 (MQ=255)                                     tGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCt                                                                                             >  1:1127134/1‑116 (MQ=255)                                     tGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCt                                                                                             <  2:1127134/116‑1 (MQ=255)                                      gTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAAtt                                                                                                         <  2:1176462/103‑1 (MQ=255)                                      gTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAAtt                                                                                                         >  1:1176462/1‑103 (MQ=255)                                      gTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGc                                                                                               >  2:875010/1‑113 (MQ=255)                                      gTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGc                                                                                               <  1:875010/113‑1 (MQ=255)                                         aaaCAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTaa                                                                  >  1:185437/1‑139 (MQ=255)                                                 tGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGATCATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATg                                                          >  1:971953/1‑139 (MQ=255)                                                        ttGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCt                                                   <  1:1085782/139‑1 (MQ=255)                                                          gtgAAAGCAATAAATTATTGACTTCGTTTTAGATTAGTTTAGCTATAATGTTATGCATTCAAACGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCttt                                                 <  2:971953/139‑1 (MQ=255)                                                                           ttGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACGTAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGAtaataa                                >  1:316110/1‑139 (MQ=255)                                                                           ttGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCt                                                                    >  1:203645/1‑103 (MQ=255)                                                                           ttGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCt                                                                    <  2:203645/103‑1 (MQ=255)                                                                            tGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGATaataat                               <  1:1030820/139‑1 (MQ=255)                                                                                ttCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGt                                                                                                             <  2:1564362/57‑1 (MQ=255)                                                                                ttCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGt                                                                                                             >  1:1564362/1‑57 (MQ=255)                                                                                ttCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTACACATTATGGTTGCCTTTTTATATAACGataata                                 >  1:1325997/1‑133 (MQ=255)                                                                                ttCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGataata                                 <  2:1325997/133‑1 (MQ=255)                                                                                        aGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGATAATAATGAATATAAGcat                   <  2:185437/139‑1 (MQ=255)                                                                                             tGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACGTAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGATAATAATGAATATAAGcatgacat              <  2:316110/139‑1 (MQ=255)                                                                                             tGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGATAATAATGAATATAAGcatgaca               >  1:397053/1‑138 (MQ=255)                                                                                             tGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGATAATAATGAATATAAGcatgaca               <  2:397053/138‑1 (MQ=255)                                                                                                         aTGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGATAATAATGAATATAAGCATGACATGAGAATAAGGtt  >  1:914462/1‑139 (MQ=255)                                                                                                                |                                                                                                                                    TATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATACATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGATAATAATGAATATAAGCATGACATGAGAATAAGGTT  >  NC_000913/1212122‑1212364

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