Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913213,3410AG42.6% 72.1 / 15.6 54G337G (GGA→GGGtilStRNA(Ile)‑lysidine synthetase
Reads supporting (aligned to +/- strand):  ref base A (10/21);  new base G (23/0);  total (33/21)
Fisher's exact test for biased strand distribution p-value = 9.33e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTGC  >  NC_000913/213206‑213489
                                                                                                                                        |                                                                                                                                                       
tttACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTaa                                                                                                                                         >  2:615666/1‑150 (MQ=255)
   aCGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGc                                                                                                                                      <  2:943503/150‑1 (MQ=255)
    cGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGAGGCTTCAACCGCTGGAATTACCGGCGGGGCCGGGAAATGGACAAGTTAAt                                                                                                                                        >  2:649631/1‑147 (MQ=255)
     gTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGCTTAATGCggg                                                                                                                                   >  1:506617/1‑151 (MQ=255)
      tttGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTCGGAAGGGGACacacttattccggg                                                                                                                                   >  1:468955/1‑138 (MQ=255)
            cgcgTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCTTAAAGCGGGAgggg                                                                                                                              >  2:793859/1‑147 (MQ=255)
             gcgTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGTACAGCTTAATgggggggggg                                                                                                                              >  2:520424/1‑146 (MQ=255)
                  tGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCg                                                                                                                      <  1:654386/150‑1 (MQ=255)
                     aaTCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAATGCGGGAGGGGATATTCGccc                                                                                                                   >  1:559139/1‑150 (MQ=255)
                          gACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGGGGCGATATTCGcccccccc                                                                                                              >  2:210704/1‑148 (MQ=255)
                             gCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGAGGCGGTATTCGCCCCCCGCGTg                                                                                                           >  1:383619/1‑150 (MQ=255)
                                aTCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCa                                                                                                         <  1:498748/149‑1 (MQ=255)
                                aTCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCa                                                                                                         <  1:592341/149‑1 (MQ=255)
                                    gTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGaa                                                                                                   <  2:1084464/151‑1 (MQ=255)
                                      cGCAACTGTGGAGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGCGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGAGGTGGGACGGGTCACGCTTAATAAGCACGGCGCTATTTGCCCACCGCCGGCAgcc                                                                                                      >  1:413975/1‑144 (MQ=255)
                                         aaCTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAAGCTAAAACGGGAGGCGGTATTTGCCCTCCCCGATCAGACGAAGAggg                                                                                              >  2:102693/1‑150 (MQ=255)
                                         aaCTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGt                                                                                              <  2:130647/151‑1 (MQ=255)
                                         aaCTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGCGTCGGGAAGGGTATAGACTAAATGGGGGGGCGATTTTCTCCCTCCGCTTGGGGACggaa                                                                                                  >  1:389126/1‑144 (MQ=255)
                                          aCTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTAAAGCTTAAAGCGGGGGGGGGTAATTGGCCCCCGCGGGCAg                                                                                                        >  2:517874/1‑140 (MQ=255)
                                          aCTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTc                                                                                             <  1:192406/151‑1 (MQ=255)
                                             gtggtgGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAgc                                                                                          <  2:1026943/151‑1 (MQ=255)
                                                      tAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGCTTTATTCGGGAGGCGATATTTGCCCTCCCCGTGCAGACGAAGCGGGCAGCGTGCGgtt                                                                                  >  2:311280/1‑150 (MQ=255)
                                                       aaaTCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAATGCGGGGGGCGGTATTTGCCCCCCGCGTGCAGACGAAACGGGCAGCGTGCGgtt                                                                                  >  2:634532/1‑149 (MQ=255)
                                                        aaTCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGGACAGCTTAAATCGGGGGGGGATATTCGCCCCCCGCGGGCAGAAGAAGCGGGCAGCGCGCGgtt                                                                                  >  2:1109753/1‑148 (MQ=255)
                                                        aaTCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGAGGGGATATTTTCCCCCCGCCTGCCGAAGAAACGGGCAGCGTGCGttt                                                                                  >  2:107660/1‑148 (MQ=255)
                                                           ccGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTTGGGAGGGTTACACCTTATCCGGGGGGCGATAATTTCCCCTCCCGTGCCGGAGAAGCGGGCTGCGGGCGGTTCaaa                                                                              >  1:636462/1‑149 (MQ=255)
                                                           ccGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAgc                                                                            <  2:1099810/151‑1 (MQ=255)
                                                             gTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAACTTAAAGCGGGGGGCGATAATTGCCCTCCCCGTGCAGACGAAGCGGGCAGCGTGCGGTTCcaa                                                                              >  1:424075/1‑147 (MQ=255)
                                                                       aaaGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATtgct                                                                <  2:133586/151‑1 (MQ=255)
                                                                        aaGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATtgctg                                                               <  1:514894/151‑1 (MQ=255)
                                                                            gAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGCCTAATTCGGGGGGGGGTATTTGCCCCCCCCGGGCCGGCGGAGGGGGCAGGGTGGGGTTTAAAACGCCCGGGAtgctgggt                                                            >  1:358717/1‑147 (MQ=255)
                                                                              aaaCATTGTGCCGTGTCAGACGTGGCTTCATCCGCTGGAATTACCGGGGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATAtt                                                         <  1:882022/151‑1 (MQ=255)
                                                                               aaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGCTTAATTCGGGAGGGGATATTTGCCCCCCCCGTGCAGACGAAGCGGGCAGCGGGCGTTTTAAAAGGCCAGGATTTCTGCATAtt                                                         >  2:376877/1‑150 (MQ=255)
                                                                                    tgtgCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGGACAGCTTAAATCGGGGGGGGATATTTGCCCCCCCCGTGCCGAAGAAGGGGGTAGGGCGGGGTTTAAAACGCCAGGAATGATGCATATTTTCggg                                                   >  2:633839/1‑151 (MQ=255)
                                                                                     gtgCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCggg                                                   <  1:179788/150‑1 (MQ=255)
                                                                                     gtgCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGc                                                  <  1:158157/151‑1 (MQ=255)
                                                                                      tgCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGAAGGGGACAGCTTAATGCGGGAGGCGAAATTTGCCTTCCGCGGGCAGACGAAGCGGGCAGCGGGCGGTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCg                                                 >  2:962817/1‑151 (MQ=255)
                                                                                             gCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGAGGGGGTATTTGCCCTCCGCGTGCAGACGAAGCGGGCAGCGTGCGTTTTAAAGCGCCCGGATTGGTGCATATTGTTGGGGGTAAcgg                                           >  2:375998/1‑150 (MQ=255)
                                                                                             gCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAAcggc                                          <  1:662524/151‑1 (MQ=255)
                                                                                             gCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAAcggc                                          <  2:528624/151‑1 (MQ=255)
                                                                                               aGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGGGGGGGTATTCGCCCCCCCCGTGCAGACGAAGAGGGCAGCGTGCGTTTCAAAGCGCCACGGTTTCTGCCTATTGTCGGGCGTAAacgcgg                                        >  1:846211/1‑151 (MQ=255)
                                                                                                 aCGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGAc                                      <  1:615574/151‑1 (MQ=255)
                                                                                                    tGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACCGCTTAATGCGGGAGGGGGTATTCGCCCTCCGCGTGGCGACGAAGCGGTCAGCGTGGGTTTCAAAGCGCCCGGAGTTCTGCATATTTTCGGGCGTAAcggggggagt                                    >  2:899771/1‑146 (MQ=255)
                                                                                                     ggCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTa                                   <  2:65173/150‑1 (MQ=255)
                                                                                                         tCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGc                                <  1:772568/149‑1 (MQ=255)
                                                                                                          cAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAACTTTAAGCGGGGGGCGGTATTTGCCCTCCGCGTGCAGACGAAGCGGGCAGCGGGCGTTTCAAAACGCCCGGATTGCTGCATATTTTTGGGGGGGAcggcgg                                        >  1:29421/1‑140 (MQ=255)
                                                                                                           aaCCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCTTAATGCGGGGGGCGATATTCGCCCCCCCCGTGCAGACGAAACGGGCAGCGGGGGGTTTAAAACGCCCGGGTTGCTGCATATTGTCGGGGGGAAAGGCGGGCGGaaagca                              >  2:483692/1‑145 (MQ=255)
                                                                                                           aaCCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTaaa                            <  1:643985/151‑1 (MQ=255)
                                                                                                            aCCGCTGGAATTACCGGCGGGGCTGGGGAGGGTTCCGGCTTAACCGGGGGGCGAATCTTGCCCCCCGCCGGCAGAAGAAACGGGGAGCGGGCGCTTCAAACCGGCACGAGTGGTGCAGAATGATGGGCGTGTCGGCGGGCGTTAg                                 >  1:911700/1‑145 (MQ=255)
                                                                                                                cTGGAATTACCGGCGGGGCTGGGAAAGGTACAGCTTAATGCGGGGGGCGAGATTTGCCCTCCCCGGGCAGACGAAACGGGGCGCCGGCGTTTCAAAAAGCCAGGGTTGCTGGATATTGTCGGGCGTAAAGGGGGACGGAAACTCaaaaaa                        >  1:381930/1‑150 (MQ=255)
                                                                                                                          cggcggGGCTGGGAGGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAaga             <  1:669889/151‑1 (MQ=255)
                                                                                                                             cggGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCt          >  2:802178/1‑151 (MQ=255)
                                                                                                                                    tGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCg     >  2:480535/1‑150 (MQ=255)
                                                                                                                                      ggAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTgc  >  2:189180/1‑151 (MQ=255)
                                                                                                                                        |                                                                                                                                                       
TTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTGC  >  NC_000913/213206‑213489

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