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

CCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTG  >  NC_000913/213196‑213488                                                                                                                                                  |                                                                                                                                                    ccTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAgtgt                                                                                                                                                <  2:824882/151‑1 (MQ=255)   tCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAgtg                                                                                                                                                 >  2:904889/1‑148 (MQ=255)    cACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGTGGGGCTGGGAAGTGTAc                                                                                                                                              <  2:862197/150‑1 (MQ=255)     aCCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAggg                                                                                                                                                 >  1:556635/1‑144 (MQ=255)         tGTTTACGTTTGGGCGCGTTTGAAATCCGGCGCTATCAGTCGCAACTGTGGTGGATTAAGTCCGTCGGCGGGCAAAGCGAAAGCATTGTGCCGTGGGAGACGTGGCTTCAGCCGGTGGAATTACCGGCGGGGGTGGGAAgtgt                                                                                                                                                <  2:269511/143‑1 (MQ=255)            ttACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGCGGGGAGGGGACCGCttta                                                                                                                                       >  2:854756/1‑147 (MQ=255)             tACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGATGGGGAgtgggcaagttaa                                                                                                                                       >  1:781221/1‑138 (MQ=255)                   tGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAAAGCgggggg                                                                                                                              >  1:118539/1‑151 (MQ=255)                    gggCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAgg                                                                                                                              <  2:378949/150‑1 (MQ=255)                           tttGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAATGCGGGGGGGGATAtt                                                                                                                       >  1:347118/1‑150 (MQ=255)                               aaaTCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGTAACGCCTAATGCGGGGGGGGGGATTCGccc                                                                                                                  >  2:396106/1‑151 (MQ=255)                               aaaTCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCTTAATGCGGGGGGGGATATTTGccc                                                                                                                  >  1:553390/1‑151 (MQ=255)                                       cGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGGGTAAAGGTTTATTCGGGGGGTGATATTTGcccccccc                                                                                                             >  2:350079/1‑146 (MQ=255)                                       cGCTATCAGTCGCAACTGTGGTGGATTAAATCCGCCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGCATTACCGGCGGGGCTGGGGACGGTATAACCCTATCCGGGAGGTGCTATTCGCCCTc                                                                                                                >  2:786554/1‑145 (MQ=255)                                        gCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGGATGTACACCTTTACTAGGGAGGGGGGAGTCACCCTCCCCGTg                                                                                                          >  1:681303/1‑150 (MQ=255)                                        gCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGAGGGGAGGGTACAGCTTAATGCGGGAGGGGGGATTCGCCCTCcgcgcg                                                                                                          >  2:333849/1‑148 (MQ=255)                                               gTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGaa                                                                                                  <  2:764783/151‑1 (MQ=255)                                                 cGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAg                                                                                                 <  2:347028/150‑1 (MQ=255)                                                 cGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGc                                                                                                <  2:875780/151‑1 (MQ=255)                                                    aaCTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGAAAGCTCAAAGTGGGGGGGGATGCTCGGCCTCCGCGGGCAGAAAAAGc                                                                                                >  2:547308/1‑148 (MQ=255)                                                       tgtgGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCCGGGAGGGGGACAGCTTTAATCGGGAGGGGATATTCGCCCCCCCCGTGCCGACGAAGCgg                                                                                              >  2:521798/1‑147 (MQ=255)                                                         tggtggATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTTATGCGGGGGGGGAGATTTCCCCTCCGCGGGCAGACGAAGCGGGGAgcg                                                                                        >  2:102293/1‑151 (MQ=255)                                                             ggATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGTTGGGGAGTGGGAAACCTTATGCGGGGGGGGGTTATCGTCCCCCGCGTGCAGAAGAAGCggggagcccgg                                                                                     >  1:620530/1‑141 (MQ=255)                                                                      ccGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAg                                                                            <  2:69502/150‑1 (MQ=255)                                                                         tCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGCTTAATGCGGGAGGCGGTATTCGCCCCCCCCGTGCAGACGAAGCGGTCAGCGGGCGGTTCAAAACGcc                                                                        >  1:690936/1‑151 (MQ=255)                                                                               ggCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGAtt                                                                  <  2:232921/151‑1 (MQ=255)                                                                                 cAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAATACGGGGGGGGGTATTCGCCCCCCCCGGGCAGACGAAGCGGGCAGCGGGCGGTTTAAAACGCCAGGAtt                                                                  >  1:674323/1‑149 (MQ=255)                                                                                 cAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATtgc                                                                <  1:579759/151‑1 (MQ=255)                                                                                   aaGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCCCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAACGCCAGGATtgct                                                               >  1:859882/1‑150 (MQ=255)                                                                                      cGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGCCTAATTCGGGGGGGGGTATTCGCCCTCCCCGGGCAGGCGAAGCGGGGAGCGGGGGGTTTCAAGCGCCCggggttctggc                                                            >  2:795703/1‑141 (MQ=255)                                                                                             aTTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACCGCTTTAAGCCGGGGGCGCTATTTGCCCCCCGCGTGCCGACGAAGCGGTTAGCGTGCGTTTTAAAACGCCAGGATTGCTGCATATTGt                                                      >  1:646686/1‑149 (MQ=255)                                                                                                  gCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCg                                                <  1:358136/150‑1 (MQ=255)                                                                                                  gCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGt                                               <  1:542164/151‑1 (MQ=255)                                                                                                     gTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGGGGCGATATTCGCCCTCCCCGTGCAGACGAAACGGGCAGCGGGCGGTTTAAAACGCCAGGATTGGTGCATATTTTTggggggaa                                             >  1:168481/1‑145 (MQ=255)                                                                                                         cAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAAcggcg                                        <  1:2184/151‑1 (MQ=255)                                                                                                            aCGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGAc                                     <  1:659535/151‑1 (MQ=255)                                                                                                                    tCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGc                               <  2:781044/149‑1 (MQ=255)                                                                                                                     cAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGACAGGTTAATGCGGGAGGGGATATTTGCCCTCCGCGGGCAGACGAAGCGGGCAGCGGGCGGTTTAAAACGCCAGGGGTGCTGCATATTGTTGGGCGGAACGGGGGGCGGAAGcaa                             >  1:328020/1‑148 (MQ=255)                                                                                                                     cAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTa                             <  2:144943/150‑1 (MQ=255)                                                                                                                                     cggcggGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAaga            <  2:198734/151‑1 (MQ=255)                                                                                                                                     cggcggGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAaga            <  2:681155/151‑1 (MQ=255)                                                                                                                                      ggcggGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAagag           <  1:858838/151‑1 (MQ=255)                                                                                                                                      ggcggGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGGTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCa               >  2:714304/1‑147 (MQ=255)                                                                                                                                       gcggGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGc          <  2:791014/151‑1 (MQ=255)                                                                                                                                         ggggCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTg        >  2:144243/1‑151 (MQ=255)                                                                                                                                         ggggCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTg        <  1:907993/151‑1 (MQ=255)                                                                                                                                          gggCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTg        >  2:615350/1‑150 (MQ=255)                                                                                                                                           ggCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTggg      >  1:342205/1‑151 (MQ=255)                                                                                                                                             cTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCg    >  2:535650/1‑151 (MQ=255)                                                                                                                                              tGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCa                                                            >  2:615818/1‑94 (MQ=255)                                                                                                                                              tGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCa                                                            <  1:615955/94‑1 (MQ=255)                                                                                                                                               gggAAGTGTACAGCTTCTTGCGGGAGGCGATATTCGCCCTCCGCGGGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCCAAAGAAAATCTGGCAGGAGAAGGGCGTg  >  2:510677/1‑151 (MQ=255)                                                                                                                                               gggAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCg    <  2:14680/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTG  >  NC_000913/213196‑213488

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: ‑