Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,642,689	0	T	G	51.1%	64.0 / 80.9	46	A52A (GCA→GCC)	rodZ	MreB assembly cytoskeletal protein
Reads supporting (aligned to +/- strand):  ref base T (3/19);  new base G (23/0);  total (27/19)
Fisher's exact test for biased strand distribution p-value = 6.32e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.26e-09

GAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATG‑TCGCGTA‑CCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTC  >  NC_000913/2642539‑2642835                                                                                                                                                       |                                                                                                                                                      gAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGt                                                                                                                                                       <  2:650124/151‑1 (MQ=255)      aaaaaCTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGccct                                                                                                                                                  >  2:532886/1‑149 (MQ=255)       aaaaCTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTa                                                                                                                                                 <  1:789107/151‑1 (MQ=255)         aaCTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGCCtttt                                                                                                                                                >  1:713462/1‑148 (MQ=255)               gCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGCTTATcttttt                                                                                                                                          >  1:450132/1‑150 (MQ=255)                cATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCa                                                                                                                                        <  1:758688/151‑1 (MQ=255)                cATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGCCTTATcttctt                                                                                                                                          >  1:897448/1‑149 (MQ=255)                  tCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCaa                                                                                                                                       <  1:75631/150‑1 (MQ=255)                   cGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGCTTAtttttt                                                                                                                                           >  2:490432/1‑143 (MQ=255)                    ggCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATG‑t                                                                                                                                   <  1:687847/151‑1 (MQ=255)                      cgcAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATG‑Tcg                                                                                                                                 <  1:29464/151‑1 (MQ=255)                          aCTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGCCTTATtttctt                                                                                                                                          >  1:635848/1‑139 (MQ=255)                                gCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATG‑TCGCGTA‑CCGT‑g                                                                                                                      <  2:739503/150‑1 (MQ=255)                                        aaGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGCTTTTTTTTTTTCTTG‑TTGCGTA‑CCGTGGGAAACCt                                                                                                               >  2:293098/1‑150 (MQ=255)                                            ggAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATG‑TCGCGTA‑CCGT‑GGAAACCTTCAgg                                                                                                          <  2:602198/150‑1 (MQ=255)                                                    cTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGGTTTTTTTCTTCTATT‑TCGCGGACCCGG‑GGAAACCTCCGGGaaaa                                                                                                      >  2:530355/1‑147 (MQ=255)                                                        tCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGCCTTATCTTCTTCAATT‑TCGCGTA‑CCGT‑GGGAAACTTCAGGCAAAATCGCTCgg                                                                                             >  1:597775/1‑151 (MQ=255)                                                          cAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATG‑TCGCGTA‑CCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCa                                                                                           <  1:211560/151‑1 (MQ=255)                                                          cAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGCCTTATCTTCTTCAATG‑TCGCGTA‑CCGT‑GGAAACCTTCataccaa                                                                                                      >  1:473432/1‑134 (MQ=255)                                                           aGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGCTTATTTTCTTCAATG‑TCGCGGA‑CCGT‑GGGAACCTTCAGGCAAAATCGCTCg                                                                                              >  2:306401/1‑147 (MQ=255)                                                               cTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGCTTTTTTTTTTTAATG‑TTGCGTACCCGG‑GGAAACC‑TCAGGCAAAAACGCCCGGCaa                                                                                          >  1:2063/1‑147 (MQ=255)                                                               cTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGCCTTTCTTCTTTTATG‑TCGCGGA‑CCGG‑GGAAACCTTCTGGCAAAGACGCTCGGCCAc                                                                                         >  2:888343/1‑148 (MQ=255)                                                                       agTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGCTTTTTTTCTTTCATT‑TTGGGTA‑CCCC‑GGGAACCTTCAAGCAAAATCGCCCGGCAAccgcccgc                                                                                  >  2:742539/1‑147 (MQ=255)                                                                       agTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGCCCTATTTTTTTTCATGTTCGGCAC‑CCAG‑GGAAACTTTCA‑GCAAAAACACTCGGGAACGGcc                                                                                     >  1:579492/1‑144 (MQ=255)                                                                         ttcttcttcTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGCCTTTCTTTTTTTATT‑TCGCGCA‑CCCC‑GGAAAACTTCACGCCAAAACGCCCGGCACCGGCccgctgt                                                                               >  2:186481/1‑147 (MQ=255)                                                                         ttcttcttcTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGACCTATTTTTTTCAAAG‑TCGCGTA‑CCCT‑GGGAACCTTCCAGCAAAAGCGCTCGGCCACCGCccgctgg                                                                               >  2:119233/1‑147 (MQ=255)                                                                                 ctGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGCCTTATCTTCTTCAATT‑TTGCGTA‑CCGT‑GGGAAACTTCAGGCAAAAACGCCCGGCCACGGCCCCCCGACTAAGTcc                                                                       >  1:760225/1‑148 (MQ=255)                                                                                   ggAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATG‑TCGCGTA‑CCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTa                                                                     <  1:583968/148‑1 (MQ=255)                                                                                   ggAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGCTTATCTTTTTCAATT‑TCGCGTA‑CCGC‑GGGAAACTTCATGCAAAAAAGCTCCGCAACGGCCTGCTGACTAAATCCTAatt                                                                  >  1:90292/1‑151 (MQ=255)                                                                                             ccAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGCCTTTCTTCTTTAATG‑TCGCGTA‑CCCC‑GGGAAACATCACGCGAAATCGCTCGGCCACGGCCCGCTGACTAAATCCTActtgtt                                                               >  2:367981/1‑144 (MQ=255)                                                                                             ccAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGCCTTTTTTTTTTTCATG‑TCGCGTA‑CCGC‑GGGAAACCTCCTGCAAAATCGCCCGGCAACGGCCCGCCGGCTAAGACCTAATTGTTtgcggg                                                         >  1:109461/1‑148 (MQ=255)                                                                                                       aTAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCGTCAATG‑TCGCGTA‑CCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAgg                                              <  1:385759/151‑1 (MQ=255)                                                                                                         aaGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATG‑TCGCGTA‑CCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGc                                             <  2:567195/150‑1 (MQ=255)                                                                                                          agagCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATG‑TCGCGTA‑CCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGa                                           <  1:7342/151‑1 (MQ=255)                                                                                                          agagCGGATATATCCGCGCAGGAATATTGAAGCAAGATCGGCGGGGGGCCTTTCTTTTTCAATG‑TCGCGGA‑CCGG‑GGGAAACTTCACGGGAAAACGCCCGGGAACGGGCCGCGGGCTAAATCCTTGTTTTTTGCGGGCAATAAGCAgg                                              >  1:503426/1‑148 (MQ=255)                                                                                                            agCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTTTTTTCTTTAATT‑TCGCGGA‑CCCT‑GGGAAACCTCCGGCAAAAACGCTCGGCAACCGCCCCCTGGATAAGTCCCAGTTTTTtgcggggcttacac                                                 >  1:849456/1‑133 (MQ=255)                                                                                                              cGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGCCTTTTCTTTTTTAAAG‑TCGCGTA‑CCGC‑GGAAAACTTCCGGCAAAAACGCTCCGCAAACGCCCGCTGACTAAAGCCTCGGTGTTTGCGAGCATTACGCCgg                                              >  1:403345/1‑144 (MQ=255)                                                                                                                              ggAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATG‑TCGCGTA‑CCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTc                       <  2:729364/151‑1 (MQ=255)                                                                                                                                     tGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATG‑TCGCGTA‑CCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTg                 <  2:850700/150‑1 (MQ=255)                                                                                                                                         gCAAGATCGGCGGGTGCCTTATCTTCTTCAATG‑TCGCGTA‑CCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTcgt            <  2:804173/151‑1 (MQ=255)                                                                                                                                          cAAGATCGGCGGGTGCCTTATCTTCTTCAATG‑TCGCGTA‑CCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTcgtg           <  1:122407/151‑1 (MQ=255)                                                                                                                                              aTCGGCGGGTGCCTTATCTTCTTCAATG‑TCGCGTA‑CCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTcgtgcgtg       <  2:238449/151‑1 (MQ=255)                                                                                                                                              aTCGGCGGGTGCCTTATCTTCTTCAATG‑TCGCGTA‑CCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTcgtgcgtg       <  1:865355/151‑1 (MQ=255)                                                                                                                                                cggcggGTGCCTTATCTTCTTCAATG‑TCGCGTA‑CCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGc     <  2:690957/151‑1 (MQ=255)                                                                                                                                                    gggTGCCTTATCTTCTTCAATG‑TCGCGTA‑CCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTc  >  1:770020/1‑150 (MQ=255)                                                                                                                                                    gggAGCCTTATCTTCTTCAATG‑TCGCGTA‑CCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGc     >  1:214078/1‑147 (MQ=255)                                                                                                                                                     ggTGCCTTATCTTCTTCAATG‑TCGCGTA‑CCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTc  >  1:680960/1‑149 (MQ=255)                                                                                                                                                       |                                                                                                                                                      GAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATG‑TCGCGTA‑CCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTC  >  NC_000913/2642539‑2642835

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