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

GCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTGCC  >  NC_000913/213195‑213490                                                                                                                                                   |                                                                                                                                                      gCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAgt                                                                                                                                                    <  2:551158/150‑1 (MQ=255) gCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAggg                                                                                                                                                   >  1:1399/1‑149 (MQ=255)          tGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGAAGTGGACAGCtt                                                                                                                                           >  1:747480/1‑150 (MQ=255)             ttACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTaa                                                                                                                                         >  1:194264/1‑149 (MQ=255)              tACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATg                                                                                                                                       <  2:716878/150‑1 (MQ=255)               aCGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTAAAGCTTAATGc                                                                                                                                      >  1:620886/1‑150 (MQ=255)               aCGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGc                                                                                                                                      <  1:84947/150‑1 (MQ=255)                         gcgTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGata                                                                                                                           <  1:652432/151‑1 (MQ=255)                           gTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGCTTAATACGGGGGGGGGTAtt                                                                                                                         >  2:31412/1‑151 (MQ=255)                           gTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAAGTTAATGGGGGGGGCGGTAtt                                                                                                                         >  1:228/1‑151 (MQ=255)                           gTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGatat                                                                                                                          <  1:160647/150‑1 (MQ=255)                            tttGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGCTTAATGGGGGAGGGGGTAtt                                                                                                                         >  1:129219/1‑150 (MQ=255)                              tGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGGTTAATGCGGGAgggggga                                                                                                                           >  2:98949/1‑141 (MQ=255)                               gAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGc                                                                                                                      <  1:336618/150‑1 (MQ=255)                                 aaTCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGAGGGGATATTCGccc                                                                                                                    >  2:663522/1‑150 (MQ=255)                                        cGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAACCTAATGCGGGAGGCGGTATTTGCCCCCCGCCTg                                                                                                            >  1:367141/1‑151 (MQ=255)                                         gCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTgt                                                                                                           <  1:131425/151‑2 (MQ=255)                                            aTCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAg                                                                                                         <  1:403330/150‑1 (MQ=255)                                             tCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTTCAGCCTAATGCGGGGGGGGGTATTCGCCCTCCCCGGGCAg                                                                                                         >  2:263280/1‑149 (MQ=255)                                                 tCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCCTAATGCGGGGGGGGGTATTTGCCCTCCGCCTGCCGa                                                                                                        >  2:17936/1‑146 (MQ=255)                                                         gtggtgGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCa                                                                                             <  2:229274/149‑1 (MQ=255)                                                         gtggtgGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAgc                                                                                           <  1:686725/151‑1 (MQ=255)                                                          tggtggATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGAGTACAGCTTAATGCGGGGGGCGATATTCTCCCCCCCCGCGCCGACGAAACgggcgagg                                                                                          >  1:195622/1‑145 (MQ=255)                                                               gATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGCTTAATGCGGGAGGCGATATTTGCCCTCCGCGCGCCGAAGAAACGGGCAgcgggcg                                                                                      >  1:495158/1‑150 (MQ=255)                                                                     aTCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCaaa                                                                               <  2:317067/151‑1 (MQ=255)                                                                       ccGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAgc                                                                             <  1:812396/151‑1 (MQ=255)                                                                           cACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGCGGGAAGGGTACAGCCTCATTCCGGCGGCCATATTCTGCCTCCTCGGGCAGCCCACGAGGTCAGCGCGCGGTTCCAAGCGccc                                                                         >  1:700158/1‑150 (MQ=255)                                                                           cACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCa                                                                         <  2:605208/151‑1 (MQ=255)                                                                                 gCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATtg                                                                   <  1:225184/151‑1 (MQ=255)                                                                                  cAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGAAGGGTACAGCTTAATGCGGGGGGCGATATTTGCCCCCCCCGTGCAGACGAAGCGGGCAGCGGGCGTTTCAAAGCGCCAGGAtt                                                                    >  1:228390/1‑149 (MQ=255)                                                                                    aaGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGCGGGGAGGGGAAAGCTTAATTCGGGGGGCGAAATTCGCCCCCCCCGTGCAGGCGAAACGGGCAGCGGGGGGTTCACAACGCCAGGGtt                                                                    >  2:188118/1‑147 (MQ=255)                                                                                    aaGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATtgctg                                                                <  1:64729/151‑1 (MQ=255)                                                                                        gAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGTACAGCTTTAAGCGGGGGGCGGGATTTGCCCCCCGCGGGCAGACGAAGCGGGGAGCGGGCGGTTCAAAGAGCCAGGGGtgctgcct                                                             >  2:357278/1‑148 (MQ=255)                                                                                         aaaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCTTTAAACGGGGGGGGGGATTTGCCCTCCCCGTGCAGACGAAGCGGGCAGCGCGGGGTTCAAAACGCCCGGGGtgctgggt                                                             >  2:642692/1‑146 (MQ=255)                                                                                          aaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTGCAGCTTAAATCGGGGGGCGATAATTGCCCCCCGCGGGCAGACGAAACGGGCAGGGGGCGGTTCTAAGCGCCAGGGATGCCGCCTAtt                                                          >  2:38613/1‑151 (MQ=255)                                                                                          aaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGCCTAATGCGGGGGGGGATATTTGCCCTCCGCGTGCAGAAGAAGCGGGCAGCGGGCGTTTCAAAGCGCCAGGATTGCTGCat                                                             >  1:643583/1‑148 (MQ=255)                                                                                                 gtgCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGc                                                   <  1:140019/151‑1 (MQ=255)                                                                                                      gTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGTACAGCTTAATGCGGGGGGGGATATTCGCCCCCCCCGTGCAGACGAAGCGGGCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTaa                                               >  1:914455/1‑150 (MQ=255)                                                                                                         gCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAAcggc                                           <  1:445236/151‑1 (MQ=255)                                                                                                          cAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGACAGCTTAATGCGGGGGGGGATATTCGCCCTCCGCGGGCAGAAGAAGCGGGCAGCGTGGGTTTCAAAACGCCAGGATTGGTGCATATTGTCGGGGGTAAcggcg                                          >  2:341392/1‑151 (MQ=255)                                                                                                           aGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAAcggcgg                                         <  1:444100/151‑1 (MQ=255)                                                                                                             aCGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGAc                                       <  1:188668/151‑1 (MQ=255)                                                                                                               gTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGGGGGAGGCGATATTTCCCCTCCGCGTGCAGACGAAGCGGGGCGCGGGCGTTTCTAAGCGCCAGGATTGCTGCATATTGTCGGGGGTAACggcgggcgg                                     >  1:807187/1‑150 (MQ=255)                                                                                                                tGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAATGAACAACTTAAAGCGGGAGGGCATATTTGCCCCCCTCGTGCAGAAAAAAAGGGCAGCGGGCGGGTTAAAGCAGCAGGTGTTCGGCTTATTGTTggggg                                                  >  2:348988/1‑135 (MQ=255)                                                                                                                   cTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGACAGCTTAATGCGGGGGGGGGTATTTGCCCTCCGCGTGGAGACGAAGCGGGCAGCGGGGGGTTCAAAGCGCCAGGATTGCTGCCTATTGTCGGGGGGAAAGGCGGACGTaa                                   >  1:328088/1‑149 (MQ=255)                                                                                                                     tCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGc                                 <  1:733914/149‑1 (MQ=255)                                                                                                                     tCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGc                                 <  1:426569/149‑1 (MQ=255)                                                                                                                         ccGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGa                           <  2:371702/151‑1 (MQ=255)                                                                                                                                aaTTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTg                    <  1:375181/151‑1 (MQ=255)                                                                                                                                    aCCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCaa                <  1:820565/151‑1 (MQ=255)                                                                                                                                    aCCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCaa                <  1:121261/151‑1 (MQ=255)                                                                                                                                                 ggAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGCCGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTg    >  2:204126/1‑150 (MQ=255)                                                                                                                                                  gAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTgcc  <  2:40476/151‑1 (MQ=255)                                                                                                                                                  gAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCAAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGCAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTgcc  >  2:730155/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTGCC  >  NC_000913/213195‑213490

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