Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,642,689	0	T	G	48.1%	74.0 / 11.6	52	A52A (GCA→GCC)	rodZ	MreB assembly cytoskeletal protein
Reads supporting (aligned to +/- strand):  ref base T (5/22);  new base G (25/0);  total (30/22)
Fisher's exact test for biased strand distribution p-value = 3.07e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCC‑TTATCTTCTTCAATGTCGCGTACCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGT  >  NC_000913/2642540‑2642838                                                                                                                                                      |                                                                                                                                                        aGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTg                                                                                                                                                        <  1:148247/151‑1 (MQ=255)               cATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCC‑TTATCTTCTTCa                                                                                                                                         <  2:569318/151‑1 (MQ=255)               cATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCC‑TTATCTTCTTCa                                                                                                                                         <  1:42981/151‑1 (MQ=255)               cATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGC‑TTATcttctt                                                                                                                                           >  2:134158/1‑149 (MQ=255)                  cGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGC‑TTTTcttctttac                                                                                                                                        >  2:239579/1‑146 (MQ=255)                    gcgcAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCC‑TTATCTTCTTCAATGTc                                                                                                                                    <  2:705020/151‑1 (MQ=255)                      gcAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCC‑TTATCTTCTTCAATGTcgc                                                                                                                                  <  1:803287/151‑1 (MQ=255)                      gcAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGC‑TTATCTTCTTCAATGTcg                                                                                                                                   >  1:169234/1‑150 (MQ=255)                        aaCTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGC‑TTATcttcttttatttcg                                                                                                                                   >  1:884504/1‑140 (MQ=255)                         aCTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGCC‑CTTTCTTTTTCAATGt                                                                                                                                     >  1:371489/1‑145 (MQ=255)                          cTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGCC‑TTTTCTTTTTTAATTTCGCGTa                                                                                                                               >  1:672070/1‑150 (MQ=255)                           ttttGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGC‑TTTTCTTCTTCAATGTTGCGGAcc                                                                                                                             >  1:996406/1‑151 (MQ=255)                             ttGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCC‑TTATCTTCTTCAATGTCGCGTACCg                                                                                                                            <  1:320230/150‑1 (MQ=255)                             ttGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGCAATGTTGAAGCAAGATCGGCGGGGGCC‑TTATCTTTTTTAATTTCGCGGAcc                                                                                                                             >  1:360735/1‑149 (MQ=255)                                cAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGC‑TTTTTTTTTTCAATGTTGCGTACCgg‑gg                                                                                                                        >  1:899532/1‑150 (MQ=255)                                  gCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGC‑TTTTCTTTTTTAATGTCGCGGACCCT‑ggga                                                                                                                      >  1:647720/1‑148 (MQ=255)                                     cGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGCC‑TTTTTTTCTTCAATGTCGCGGACCGG‑GGaaaac                                                                                                                   >  1:709891/1‑148 (MQ=255)                                      gAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCC‑TTATCTTCTTCAATGTCGCGTACCGT‑GGAAACCtt                                                                                                                 <  1:210403/151‑1 (MQ=255)                                           ggAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCC‑TTATCTTCTTCAATGTCGCGTACCGT‑GGAAACCTTCAgg                                                                                                             <  2:688879/150‑1 (MQ=255)                                              gCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGCC‑TTATTTTCTTCAATGTTGCGTACCCTGGGAAACCtc                                                                                                                 >  2:180316/1‑143 (MQ=255)                                                  gCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCC‑TTATCTTCTTCAATGTCGCGTACCGT‑GGAAACCTTCAGGCAAAGTCg                                                                                                     <  1:178346/151‑1 (MQ=255)                                                         cAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGC‑TTATCTTCTTCAATGTCGCGGACCCT‑GGGAACCTTCCGGCAAAAACGCCCgg                                                                                                >  2:772838/1‑149 (MQ=255)                                                          aGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGC‑CTCTCTTCTTCAATAGCGAGTAGCGC‑GGAAACCTTCACGAAAAGCCGcgctg                                                                                                >  2:732719/1‑144 (MQ=255)                                                               tGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGCC‑TTATCTTTTTCAATGTCTCGTACCGT‑GGAAAACTTTAGGCAAAAACGCTCGGCCAAGGcc                                                                                        >  1:178919/1‑151 (MQ=255)                                                                ggCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCC‑TTATCTTCTTCAATGTCGCGTACCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCct                                                                                       <  2:697068/151‑1 (MQ=255)                                                                  cagcagTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGCC‑TTATCTTCTTTTATGTCGCGTACCGT‑GGGAACCTTCAGGCAAAAACGCTCCGCaaacggc                                                                                        >  2:433159/1‑143 (MQ=255)                                                                    gcagTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGCC‑TTATCTTCTTTAATTTTTCGTACCGT‑GGGAAACCTCAGGCAAAAACGCTCGGCAAAGGCccgc                                                                                     >  2:955190/1‑149 (MQ=255)                                                                           ttcttctGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGGC‑CTATTTTCTTTAATGTGGCGGACCGT‑GGGAAACTTCACGGAAAAACGCTCGGGAACGGGccgctgggta                                                                               >  1:566845/1‑144 (MQ=255)                                                                            tcttctGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGCC‑TTTTCTTCTTCAATTTTGCGGACCGT‑GGGAAACTTTAGGCAAAAACGCTCGGCAACGGCccgctggcca                                                                               >  1:232121/1‑143 (MQ=255)                                                                               tctGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCC‑TTATCTTCTTCAATGTCGCGTACCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTcc                                                                          <  2:169414/149‑1 (MQ=255)                                                                                  ggAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCC‑TTATCTTCTTCAATGTCGCGTACCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTa                                                                        <  2:557788/148‑1 (MQ=255)                                                                                  ggAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCC‑TTATCTTCTTCAATGTCGCGTACCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAgtt                                                                     <  2:976584/151‑1 (MQ=255)                                                                                   gAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCC‑TTATCTTCTTCAATGTCGCGTACCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAgttg                                                                    <  2:647648/151‑1 (MQ=255)                                                                                       atGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCC‑TTATCTTCTTCAATGTCGCGTACCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTcg                                                                <  1:898060/151‑1 (MQ=255)                                                                                         gTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCC‑TTATCTTCTTCAATGTCGCGTACCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTcgc                                                               <  2:273677/150‑1 (MQ=255)                                                                                         gTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGC‑TTATCTTTTTTCATGTTGCGGACCGG‑GGGAAACTTCAGGCAAAAACGCTCGGCAACGGCCCGCTGACTAAAACCTACTTGTTtgcg                                                              >  2:594681/1‑151 (MQ=255)                                                                                          tACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGC‑TTTTCTTCTTCCAAGTCGCGGACCGG‑GGGAAACTTCTGGCAAAAACGCCCCGCAAAGGCCCGCCGGCTCAAGCCTAGTTGTTtgcg                                                              >  2:624239/1‑150 (MQ=255)                                                                                          tACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGCC‑TTATCTTTTTTAATTTCGCGGACCCG‑GGGAAACTTCCGGCAAAAACGCCCGGCAACGGGCCGCTGACTAAATCCCCgttgtt                                                                  >  1:136655/1‑146 (MQ=255)                                                                                                aCGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCC‑TTATCTTCTTCAATGTCGCGTACCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTa                                                       <  1:674914/151‑1 (MQ=255)                                                                                                     cATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGC‑TTATCTTTTTCAATGTCGCGGACCGG‑GGGAAACCTCAGGCAAAAACGCTCCGCAAACGCCCGCTGGCTAAGTCCTAGTTTTTTGCGAGCATTACGc                                                    >  2:771480/1‑149 (MQ=255)                                                                                                        aaGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCC‑TTATCTTCTTCAATGTCGCGTACCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGc                                                <  2:191332/150‑1 (MQ=255)                                                                                                         agagCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGC‑TTATTTTCTTCAATTTCGCGGACCGT‑GGGAACCTTCAGGCAAAAGCGCCCCGCAACCGCCCGCTGACTAAGTCCTAGTTGTTTGCGAGGATTACGCAgg                                                 >  1:25240/1‑148 (MQ=255)                                                                                                                      ccGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCC‑TTATCTTCTTCAATGTCGCGTACCGT‑GGGAAACTTCAGGCAAAAACGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCCGGCGAGCGCCGggagga                                  >  1:154053/1‑148 (MQ=255)                                                                                                                      ccGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCC‑TTATCTTCTTCAATGTCGCGTACCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGgtagta                                  <  1:974919/150‑1 (MQ=255)                                                                                                                           cAGGAATGTTGAAGCAAGATCGGCGGGGGCCTTTTTTTTTTTCAAT‑TTGCGTACCGG‑GGGAAACCTCCGGCAAAAACGCCCCGCAACCGCCCGCTGGCTAAAGCCTAGTTTTTTTCGGGCATTTCGCCGGCGGGCGcccggaggaaagggt                            >  1:611958/1‑138 (MQ=255)                                                                                                                                aTGTTGAAGCAAGATCGGCGGGGGCC‑TTTTTTTCTTTAATATTGCGGACCGT‑GGAAACCTTCAGGCAAAAACGCCCGGCAAACGCCCGCTGGCTAAGGCCCCGTTTTTTGCGCGCATTACGCAGGCGGGCGGCGGGAGGAAGTGCTTCAtt                       >  2:794910/1‑151 (MQ=255)                                                                                                                                        gCAAGATCGGCGGGTGCC‑TTATCTTCTTCAATGTCGCGTACCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTcgt               <  1:650167/151‑1 (MQ=255)                                                                                                                                                ggcggGTGCC‑TTATCTTCTTCAATGTCGCGTACCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGc        >  2:980562/1‑150 (MQ=255)                                                                                                                                                 gcggGTGCC‑TTATCTTCTTCAATGTCGCGTACCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCtt      <  1:557058/151‑1 (MQ=255)                                                                                                                                                  cggGTGCC‑TTATCTTCTTCAATGTCGCGTACCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTc     <  1:134169/151‑1 (MQ=255)                                                                                                                                                   gggTGCC‑TTATCTTCTTCAATGTCGCGTACCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTc     >  2:790731/1‑150 (MQ=255)                                                                                                                                                     gTGCC‑TTATCTTCTTCAATGTCGCGTACCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGt  >  1:854115/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                        AGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCC‑TTATCTTCTTCAATGTCGCGTACCGT‑GGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGT  >  NC_000913/2642540‑2642838

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