Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	213,341	0	A	G	50.0%	56.5 / 20.3	53	G337G (GGA→GGG)	tilS	tRNA(Ile)‑lysidine synthetase
Reads supporting (aligned to +/- strand):  ref base A (1/25);  new base G (26/0);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 1.09e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTGCC  >  NC_000913/213200‑213490                                                                                                                                               |                                                                                                                                                       aCCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACa                                                                                                                                                <  1:917863/150‑1 (MQ=255)   ccTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGGCGGGGAGGGTaaa                                                                                                                                                >  2:475555/1‑146 (MQ=255)    cTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGGTGGGGAGGGGACAGcct                                                                                                                                            >  1:930267/1‑149 (MQ=255)      gTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCtt                                                                                                                                            <  1:872334/149‑1 (MQ=255)        ttACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGGAGGGGACacactaa                                                                                                                                          >  1:880274/1‑143 (MQ=255)         tACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGTGCTCGGGAGGGTACAGc                                                                                                                                              >  1:634368/1‑144 (MQ=255)         tACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGAAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTTCCGGC‑GGGGCTGGGAAGTGTACAGCTTAATg                                                                                                                                        <  2:293310/150‑1 (MQ=255)           cGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGGTGGGGAgtgggaagctttaagcg                                                                                                                                      >  2:341662/1‑136 (MQ=255)           cGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGAGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGGTGGGGAGGGAGCAGGTTAATgtg                                                                                                                                      >  2:850487/1‑148 (MQ=255)               tGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAgg                                                                                                                                 >  1:202422/1‑151 (MQ=255)               tGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCGGGGGAGGGGACAGCTTAATGCgggggg                                                                                                                                 >  1:128730/1‑151 (MQ=255)                gggCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAgg                                                                                                                                 <  1:85821/150‑1 (MQ=255)                      gTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATAtt                                                                                                                          <  2:938713/151‑1 (MQ=255)                       tttGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGGCGGGGAGTGTGCCGCTTAAATCGGGGGGCGGTAtt                                                                                                                          >  2:395003/1‑150 (MQ=255)                         tGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGGAGTGGAAAACTTAAAGCGGGAGGGGAAAttt                                                                                                                         >  1:623582/1‑148 (MQ=255)                                     cTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGGAGGGGACAGCTTTATTCGGGGGGGGATATTCGcccccc                                                                                                                  >  1:550651/1‑144 (MQ=255)                                      tATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCGGGGGAGGGTAAAGCTTTATGCGGGAGGGGGGAATCGcccccccc                                                                                                                >  1:830497/1‑143 (MQ=255)                                         cAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGGTGGGGAGGGGACAGCTTAATGCGGGAGGGGGTAATCGCCCCCCCCGTGGAg                                                                                                          >  2:844809/1‑148 (MQ=255)                                          aGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGGAGGGGACAGCTTAATTCGGGGGGGGATATTTGCCCTCCCCCTGCAGACg                                                                                                       >  1:885496/1‑150 (MQ=255)                                           gTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGGGGGGGAGGGGAAAACTTAAAGCGGGGGGGGGAATTCGCCCCCcgcgg                                                                                                              >  2:483697/1‑141 (MQ=255)                                              gCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGCGGGGCAG‑GTAAAACTTAATACCGGGGGGGATATTTCCCCCCCCCGTGctggacca                                                                                                     >  2:147500/1‑141 (MQ=255)                                                aaCTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGt                                                                                                <  2:425558/151‑1 (MQ=255)                                                        tggATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGGTGGGGAGGGTACAGCTTTATGCGGGGGGCGGTATTCGCCCCCCCCGTGCAGACGAAGCGGGCAggggg                                                                                         >  1:33995/1‑146 (MQ=255)                                                          gATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAgcgtgc‑gt                                                                                     <  1:681271/151‑1 (MQ=255)                                                             tAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑Gttt                                                                                   <  1:757736/150‑1 (MQ=255)                                                                aTCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGGAGGGTACAGCTTTATGCGGGGGGCGATATTTGCCCCCCCCGTGCCGACGAAACGGGCAGCGGGC‑GGTTTaa                                                                                >  1:145582/1‑150 (MQ=255)                                                                  ccGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCGGGGGAGTGGACAGACTAATTAGGGGGGGGATATTCGGCCTCCGCGGGTCGGCGAAGCGGTCacggcgcggtgtaa                                                                                 >  2:403717/1‑135 (MQ=255)                                                                        ccGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGGAGGGTACAGCTTAATGCGGGGGGGGGTATTCGCCCTCCCCCTGGAGGCGAAGCGGGCCGCGGGC‑GGTTTaaa                                                                               >  2:205123/1‑143 (MQ=255)                                                                         cGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGGAGAGGAAAGCCTAATGCGGGGGGGGAATTTTGCCCCCCCCGTGCAGACGAAGCGGGCAGCGGGCGGTTTCAAA‑CGCCAgg                                                                       >  1:698984/1‑150 (MQ=255)                                                                           ggCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGAtt                                                                    <  1:463553/151‑1 (MQ=255)                                                                           ggCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGAtt                                                                    <  2:267286/151‑1 (MQ=255)                                                                             cAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTGCAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGAtt                                                                    <  2:610939/149‑1 (MQ=255)                                                                                    aaaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCatat                                                           <  1:562941/151‑1 (MQ=255)                                                                                      aaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTg                                                         <  2:4438/151‑1 (MQ=255)                                                                                            gtgCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGGTGGGGAGGGGACAGCTTTATGCGGGGGGGGGTATTCGCCCCCCCCGGGGAGACGGAGCGGGCAGCGGGC‑GGTTCAAAACGCCAGGGTTGCTGCATATTGTTggg                                                    >  2:462050/1‑150 (MQ=255)                                                                                             tgCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGGTGGGGGGGGTACAACCTAAAGCGGGGGGCGGTATTTGCCCCCCGCGTGCCGACGGAAGGGGCAGCGGGC‑GGTTCAAAACGCCAGGATTGCTGCATATTTTCggggg                                                  >  2:12425/1‑149 (MQ=255)                                                                                             tgCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGGAGTGGACAGCTTAATGCGGGAGGGGGTATTTGCCCTCCGCGGGCAGACGGAGCGGGGAGCGGGG‑GTTTTAAAGCGCCAGGATTGCTGCATAATGTTGGGCg                                                  >  1:650865/1‑151 (MQ=255)                                                                                              gCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGt                                                 <  2:930142/151‑1 (MQ=255)                                                                                              gCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGAGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCg                                                  <  1:441160/150‑1 (MQ=255)                                                                                                     cAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGGAGTGTACAGCTTAATTCGGGGGGCGATATTTGCCCTCCCCGTGCAGACGAAGCGGGCAGCGTGG‑GGTTCCAAGCGCCAGGATTGGTGCCTATTTTCGGGGGTAAcgggg                                          >  1:25495/1‑149 (MQ=255)                                                                                                      aGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAAcggcgg                                         <  2:572238/151‑1 (MQ=255)                                                                                                        aCGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGAc                                       <  1:169177/151‑1 (MQ=255)                                                                                                         cGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGGAGTGTACAGCTTTAAGCGGGGGGGGGTATTTGCCCCCCGCGGGCAGACGGAGCGGGCAGCGGGC‑GGTTTAAAACGCCCGGGTTGCTTCCTATTTTTGGGGGGAACggcgggcg                                      >  2:364631/1‑151 (MQ=255)                                                                                                         cGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACg                                      <  2:686951/151‑1 (MQ=255)                                                                                                             gCTTCAACCGCTGGAATTACCGGC‑GGGGGTTGGGAGGGTACAGCTTTATGCGGGGGGGGATATTTGCCCCCCGCGGGGAGACGGAGCGGGCCGCGGGG‑GGTTCAAAACGCCAGGACTGCTGCATATTTTTGGGGGGAAcggggggaggaa                                   >  2:401164/1‑144 (MQ=255)                                                                                                                         ggAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAAt                       <  2:1008449/150‑1 (MQ=255)                                                                                                                          gAATTACCGGC‑GGGGCTGGGGAGGGGGCAGCTTAAAGCGGGAGGGGATATTTGCCCTCCCCGGGGAGACGGAGCGGGACGCGTGC‑GGTTCAAAAAGCCAAGGGTGCTGGCTATTTTTTGGGGTAACGGCGGAGGGAAGAGaaaaaaaaaa                      >  2:983297/1‑148 (MQ=255)                                                                                                                              tACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGc                  <  2:289383/150‑1 (MQ=255)                                                                                                                                 cggc‑ggGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAaga              <  1:316157/151‑1 (MQ=255)                                                                                                                                  ggc‑ggGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAagag             <  2:541009/151‑1 (MQ=255)                                                                                                                                       gggCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTgg         <  2:815663/151‑1 (MQ=255)                                                                                                                                         gCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGc       <  2:984235/151‑1 (MQ=255)                                                                                                                                            gggAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCg      <  1:1007232/149‑1 (MQ=255)                                                                                                                                               aaGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGCAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTgcc  >  2:993884/1‑150 (MQ=255)                                                                                                                                               |                                                                                                                                                       ACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTGCC  >  NC_000913/213200‑213490

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