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

GATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACA‑GCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTG  >  NC_000913/213192‑213488                                                                                                                                                      |                                                                                                                                                     gATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGaa                                                                                                                                                     <  2:730101/151‑1 (MQ=255)        cACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTa                                                                                                                                                <  1:384655/149‑1 (MQ=255)        cACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACa                                                                                                                                              <  1:554730/151‑1 (MQ=255)         aCCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACA‑g                                                                                                                                            <  2:161737/151‑1 (MQ=255)            cTGTTTCCGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACaa                                                                                                                                             >  1:393273/1‑147 (MQ=255)                 tACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACA‑GCTTTAt                                                                                                                                      >  1:365163/1‑149 (MQ=255)                       tGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACA‑GCTTAATGCGGGAg                                                                                                                               <  2:66097/150‑1 (MQ=255)                              gTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCACAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGCGGGGTGAGGAAC‑ACTTAAAGCGGGAGGAGatat                                                                                                                        >  1:590067/1‑150 (MQ=255)                                    aaTCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTGCA‑GCTTAAAAAGGGGGGGGATATTCGccc                                                                                                                  >  1:480743/1‑150 (MQ=255)                                    aaTCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACA‑GCTTAATGCGGGAGGCGATATTCGCCCt                                                                                                                 <  2:192195/151‑1 (MQ=255)                                         gACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACA‑GCTTAATGCGGGAGGCGATATTCGCCCTCcgcg                                                                                                            <  2:395231/151‑1 (MQ=255)                                          ccGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGGGTGAA‑GCTTAAAGCGGGGGGCGGTATTcccccccctc                                                                                                             >  1:142055/2‑147 (MQ=255)                                           cGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACA‑GCCTAATGCGGGGGGGGGTATTCGCCCTcc                                                                                                               >  2:20035/1‑146 (MQ=255)                                             cTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACA‑GCTTAATGCGGGGGGGGATATTTGCCCCCcgcggg                                                                                                          >  2:29656/1‑147 (MQ=255)                                                tCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACA‑GCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGa                                                                                                      <  1:700803/150‑1 (MQ=255)                                                tCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCGGGGGAGGGGAAAGGCTTAA‑ACGGGGGGGGGTATTCTCTCTCcgcggg                                                                                                          >  1:197710/1‑144 (MQ=255)                                                     cGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACA‑GCTTAATGCGGGAGGGGATATTCGCCCTCCGCGTGCAGACGaa                                                                                                  >  2:592890/1‑149 (MQ=255)                                                      gCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGAGGAAA‑GCTTTATGCGGGGGGCGGTATTTGCCCCCCGCGTGCCGACGaa                                                                                                  >  2:400008/1‑148 (MQ=255)                                                            gtggtgGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACA‑GCTTTATGCGGGAGGCGATATTCGCCCCCCGCGTGCAGACGAAGCGGTCa                                                                                           >  2:536208/1‑149 (MQ=255)                                                              ggtggATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACA‑GCTTAATGCGGGGGGGGATATTTGCCCTCCGCGGGCAGACggaacg                                                                                               >  2:523774/1‑138 (MQ=255)                                                              ggtggATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACA‑GCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAgcg                                                                                        <  2:134723/150‑1 (MQ=255)                                                                     tAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACA‑GCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGttt                                                                                 <  1:604611/150‑1 (MQ=255)                                                                       aaTCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGACTTACCGGCGGGGGTGGGGAGAGTACA‑GTCTAATGCGGGGGGGGGTAATACCCCCCCCCGTGGAGAGGAAACGGTGAGCGGGCGgtt                                                                                 >  1:406162/1‑148 (MQ=255)                                                                           cGTCACCGGGCAAAGCGAAAACATTGTGCCGTGTCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAAAATACA‑GCTTTATTGGGGAGGCGCTATTTTTCCTTCCCGTGCAGTCGAAACGGTCAGTGTGTGTTTTaata                                                                            >  1:377426/1‑147 (MQ=255)                                                                             tCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGACA‑GCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAACGGGCAGCGTGCGGTTCAAAGAGcc                                                                        >  2:128901/1‑151 (MQ=255)                                                                              cACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACA‑GCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCa                                                                       <  1:129573/151‑1 (MQ=255)                                                                                   ggCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACA‑GCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGAtt                                                                  <  2:165597/151‑1 (MQ=255)                                                                                    gCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACA‑GCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATtg                                                                 <  2:505799/151‑1 (MQ=255)                                                                                      aaaGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACA‑GCTTTATGCGGGGGGCGAAATTCCCCCCCCCCGTGCAGAAAAAGCGGGCAGCGTGCGGTTTAAAACGCCCGGAtt                                                                  >  1:34779/1‑148 (MQ=255)                                                                                           gAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAATGGACA‑GCCTAATGCGGGGGGGGAGATTTGCCCCCCCCGCGGAGACGAAGCGGGCAGCGGGGGGTTCTAAACGCCAGGAGtgct                                                               >  1:490089/1‑146 (MQ=255)                                                                                           gAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACA‑GCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGTATTGCTGCata                                                          <  1:525417/151‑1 (MQ=255)                                                                                             aaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGACA‑GCTTAAAGGGGGTGGGGATATTCGCCCTCCGCGTGCAGACGAAGCGGGCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCatat                                                         >  2:640422/1‑150 (MQ=255)                                                                                                cATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACA‑GCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTc                                                     <  2:681330/151‑1 (MQ=255)                                                                                                       ccGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACA‑GCTTAAATCGGGGGGGGGTATTCGCCCTCCCCGTGCAGAAGAAGCGGGCCGCGGGGGGTTTAAAACGCCAGGAGTGCTGCATTTTGTCGGGCg                                                >  2:270431/1‑149 (MQ=255)                                                                                                        cGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGACA‑GCTTAATGCGGGGGGGGGTATTCGCCCTCCGCGGGCAGACGAAGCGGGCAGCGTGCGGTTTAAAACGCCAGGGTTGCTGCATATTTTCggggg                                                >  2:413060/1‑146 (MQ=255)                                                                                                        cGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTGCA‑GCCTAATGAGGGGGGCGATATTCGCCCTCCGCGTGCGGACGAAGCGGTCAGCGTGCGGTTTAAAACGCCAGGATTGATGCCTATTGTCTGGCGt                                               >  2:771615/1‑149 (MQ=255)                                                                                                             cAGTCGTGGCCTCAACCGCTTGAATTACCTGCGGGGCTGGGAAGGGTACA‑GCTTAATGCGGGAGGCGTTATTCGCCCGCAGCCTGCAGACGACGCCCTCCCTCTGCATCTCAAAGCGCCAGCATTCCTGCATATTGTCGCGCATAAcggc                                         <  1:184637/150‑1 (MQ=255)                                                                                                             cAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACA‑GCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAAcggcg                                        <  2:142685/151‑1 (MQ=255)                                                                                                               gACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGCGGGGGGGGGTGA‑GGTTAATGGGGGGGGGGGTATTTGCCCCCCCCGTTCCGACGGAAGGGGGAGCGCGCGGTTTAAAGGGCACGGGGTGGTGGATATTGTCGGGGGTAAaggcgg                                       >  1:796317/1‑150 (MQ=255)                                                                                                                  gTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACA‑GCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGt                                   <  2:495932/151‑1 (MQ=255)                                                                                                                   tGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACA‑GCTTAATGCGGGGGGCGATATTTGCCCCCCCCGTGCAGACGAAGCGGGCAGCGTGCGGTTTAAAACGCCCGGGTTGCTGCATATTTTCGGGCGTAAcggcgg                                       >  2:290561/1‑146 (MQ=255)                                                                                                                    ggCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACA‑GCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTaa                                 <  2:363178/151‑1 (MQ=255)                                                                                                                     gCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACA‑GCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAg                                <  1:261538/151‑1 (MQ=255)                                                                                                                       ttCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGGACA‑GCTTTATTCGGGGGGGGGTATTCGCCCCCCCCGCGCGGACGAAGGGGGCAGCGGGGGGTTTAAAGCGCCCGGATTTCTGCATATTGTCGGGGGGAACGGGGGGGGGAAGCt                              >  2:353458/1‑151 (MQ=255)                                                                                                                        tCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACA‑GCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCt                              <  1:336937/150‑1 (MQ=255)                                                                                                                         cAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACA‑GCTTAATGCGGGAGGGGATATTTGCCCTCCGCGGGCAGACGAAGCGGGCCGCGTGCGTTTTAAAACGCCCGGGTTGCTGCCTATTGTTGGGCGTAACGGGGGGCGGAAg                                >  2:123321/1‑147 (MQ=255)                                                                                                                         cAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACA‑GCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTaa                            <  1:32793/151‑1 (MQ=255)                                                                                                                              gCTGGAATTACCGGCGGGGCGGGGGAGGGGACA‑GCCTAATGCGGGGGGGGGTATTTCCCCCCCCCGGGCGGACGAAACGGGCAGCGCGCGCTTCAAAGAGCCCGGGTTGGTGCCTATTTTCGGGCGGGACggcgggcggga                                 >  1:261023/1‑138 (MQ=255)                                                                                                                                 ggAATTACCGGCGGGGCTGGGAAGTGTACA‑GCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATc                    <  2:68982/151‑1 (MQ=255)                                                                                                                                  gAATTACCGGCGGGGCTGGGAAGTGTACA‑GCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCt                   <  1:151289/151‑1 (MQ=255)                                                                                                                                       aCCGGCGGGGCTGGGAAGTGTACA‑GCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCaa              <  1:415963/151‑1 (MQ=255)                                                                                                                                          ggcggGGCTGGGAAGTGTACA‑GCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAag             >  2:92101/1‑149 (MQ=255)                                                                                                                                           gcggGGCTGGGAAGTGTACA‑GCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGc          >  2:560780/1‑151 (MQ=255)                                                                                                                                                    ggAAGTGTACA‑GCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTg  >  2:660073/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACA‑GCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTG  >  NC_000913/213192‑213488

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