Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	213,341	0	A	G	43.9%	53.0 / 11.1	41	G337G (GGA→GGG)	tilS	tRNA(Ile)‑lysidine synthetase
Reads supporting (aligned to +/- strand):  ref base A (4/19);  new base G (18/0);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 5.49e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGC  >  NC_000913/213191‑213485                                                                                                                                                       |                                                                                                                                                   aGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGa                                                                                                                                                    <  2:658609/151‑1 (MQ=255)  gATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGaa                                                                                                                                                   <  1:644951/151‑1 (MQ=255)         cACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTa                                                                                                                                              <  2:181277/149‑1 (MQ=255)              tGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCtt                                                                                                                                        >  2:700613/1‑150 (MQ=255)               gTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTaa                                                                                                                                      <  1:254076/151‑1 (MQ=255)                   aCGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCg                                                                                                                                  <  2:588867/151‑1 (MQ=255)                      tttGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCTTTATGCggg                                                                                                                                >  1:349142/1‑150 (MQ=255)                       ttGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCCTAATACggggg                                                                                                                              >  2:51460/1‑149 (MQ=255)                         gggCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAATGCgggggg                                                                                                                             >  2:202272/1‑150 (MQ=255)                                  tGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCg                                                                                                                   <  2:301042/150‑1 (MQ=255)                                  tGAAATCCGACGCTATCAGTCGCAACTGTGGTGGAGTAAATCCGTCACCGGGCAAAGCGAAAACAGTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCGGGGAAGTGTACAGCTTAATGCGGGAGGCGat                                                                                                                         <  1:116205/145‑1 (MQ=255)                                     aaTCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCAGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGCTTTATGAGGGGGGGGGTA‑TTCTccc                                                                                                                >  2:478996/1‑150 (MQ=255)                                      aTCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGGGGGGGTA‑TTCGCCCTc                                                                                                              >  2:482983/1‑151 (MQ=255)                                              cTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCGGGGGAGGGCTCAGCTTAATGCGGGGGGCGATATTTCGcccccgccg                                                                                                          >  2:194572/1‑142 (MQ=255)                                                aTCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAg                                                                                                     <  1:684472/150‑1 (MQ=255)                                                 tCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCTTAAAGCGGGGGGGGGTA‑TTCGCCCCCCCCGGGCa                                                                                                      >  1:722991/1‑148 (MQ=255)                                                      cGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAAt                                                                                                                                     <  2:168782/113‑1 (MQ=255)                                                        cAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGCTTAATGCGGGGGGGGGTA‑TTTGCCCTCCGCGTGGAGACGAAGCgg                                                                                            >  2:407973/1‑151 (MQ=255)                                                                      tAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAAGGTATAGCCTAAAACGGGGGGGGAAA‑TTTGCCCCCCCCGCGCAGGAGAAGCggggagagtgtggt                                                                                >  1:741228/1‑137 (MQ=255)                                                                       aaaTCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGttt                                                                              <  1:176472/149‑1 (MQ=255)                                                                               cACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAAGGTACAGCTTAATGTGGGGGGCGGTA‑TTTCCCCTCCCCCTGGCGAAGAAGCGGGCAGC‑GGGGGGTTTAAAACGcc                                                                     >  1:112608/1‑150 (MQ=255)                                                                                 ccGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGAGTACAGCTTAATTCGGGGGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGGCAGC‑GTGCGTTTCAAAGCGCCAgg                                                                  >  2:230115/1‑151 (MQ=255)                                                                                    ggCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGAtt                                                               <  2:104763/151‑1 (MQ=255)                                                                                        aaGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATtgctg                                                           <  2:421790/151‑1 (MQ=255)                                                                                          gCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCCGGGGAGGGGACAGCCTTAAGCGGGGGGGGGTA‑TTTGCCCCCCCCGGGCAGACGAAAGGGGCCGC‑GGGGGGTTCAAAACGCCCGGGGtgct                                                            >  1:475936/1‑148 (MQ=255)                                                                                           cGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGACAGCTTAATGCGGGGGGGGATA‑TTTGCCCCCCGCCGGCAGACGAAGCGGGCCGCGGGGGGTTTCAAA‑CGCCCGGGTtggt                                                            >  1:45603/1‑145 (MQ=255)                                                                                                 cATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGACAGCTTTAATCGGGGGGGGGTA‑TTTCCCCTCCCCGGGGCGACGAAGCGGGCAGC‑GGGGGGTTTAAAACGCCCGGGGTGCTGCATAttt                                                    >  1:775306/1‑148 (MQ=255)                                                                                                    tgtgCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCCTAATGCGGGGGGCGGAA‑TTTGCCCCCCGCGGGCAGACGAAGCGGGCAGC‑GGGCGGTTCAAAGCGCCCGGATTGCTGCATATTGTCggg                                               >  1:362103/1‑151 (MQ=255)                                                                                                            ggCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAAcgg                                       <  1:114653/151‑1 (MQ=255)                                                                                                             gCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTGCAGCTTTATGCGGGGGGGGGTA‑TTTGCCCCCCGCGTGCCGGCGAAGCGGGCAGC‑GGGGGGTTCAAAAAGGCCGGATTGTTGCCTATTGTCGGGGGGAAcgg                                       >  1:176624/1‑150 (MQ=255)                                                                                                             gCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAAcggc                                      <  2:290620/151‑1 (MQ=255)                                                                                                                 aCGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGAc                                  <  1:202240/151‑1 (MQ=255)                                                                                                                         tCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCt                           <  1:511768/150‑1 (MQ=255)                                                                                                                               gCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGaaa                    <  2:171146/151‑1 (MQ=255)                                                                                                                                 tGGAATTACCGGCGGGGGTGGGGAGGGTACAGCTTTATGCGGGGGGCGGTA‑TTTGCCCCCCCCCGGCAGACGAAGCGGGCCGC‑GTGCGGTTCAAAACGCCCGGATTTCTGCATATTGTCGGGCGTAACGGCGGGCGGAAACTaaaaaaaa                   >  2:217600/1‑150 (MQ=255)                                                                                                                                   gAATTACCGGCGGGGCTGGGGAGTGGACAGCCTAAAGGGGGGGGCGGTA‑TTTCCCCCCCCCGTGCAGACGAAGCGGGCAGC‑GGGGGGTTCCAAACGCCCGGATTTCTGCATATTTTCGGGGGGAACGGCGGGAGTAAGCTaaaaaaaa                   >  1:615238/1‑148 (MQ=255)                                                                                                                                    aaTTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTg               <  1:260793/151‑1 (MQ=255)                                                                                                                                         ccGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGc             >  2:244820/1‑148 (MQ=255)                                                                                                                                            gcggGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGc       <  1:748225/151‑1 (MQ=255)                                                                                                                                             cggGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCt      >  1:745609/1‑151 (MQ=255)                                                                                                                                               gggCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTgg    <  1:150580/151‑1 (MQ=255)                                                                                                                                                ggCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTggg   >  1:708712/1‑151 (MQ=255)                                                                                                                                                 gCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGc  <  2:536720/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                   AGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGC  >  NC_000913/213191‑213485

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