Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913213,3410AG47.1% 33.1 / 10.4 34G337G (GGA→GGG) tilStRNA(Ile)‑lysidine synthetase
Reads supporting (aligned to +/- strand):  ref base A (5/13);  new base G (16/0);  total (21/13)
Fisher's exact test for biased strand distribution p-value = 9.84e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGG  >  NC_000913/213192‑213483
                                                                                                                                                     |                                                                                                                                               
gATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGaa                                                                                                                                                <  2:193469/151‑1 (MQ=255)
   gCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAggg                                                                                                                                             >  1:32864/1‑149 (MQ=255)
       cACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACa                                                                                                                                         <  1:693754/151‑1 (MQ=255)
         cccTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGc                                                                                                                                       <  1:552743/151‑1 (MQ=255)
          ccTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGGGTGTACa                                                                                                                                         >  1:694557/1‑148 (MQ=255)
             gTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCtt                                                                                                                                     <  1:198714/149‑1 (MQ=255)
                tACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCTTaa                                                                                                                                   >  1:407403/1‑148 (MQ=255)
                 aCGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGc                                                                                                                                <  2:346036/150‑1 (MQ=255)
                   gTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGGGGGTACAGCTTAATGCggg                                                                                                                             >  2:618261/1‑151 (MQ=255)
                          cgcgTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGCGGGGAGGGGACAGCTTAATGCgggggggggt                                                                                                                      >  2:593805/1‑147 (MQ=255)
                                      ccGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGAAAGAATAAAACGGGAGGCGGGATTCGCCCTccc                                                                                                          >  2:270769/1‑150 (MQ=255)
                                          cGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCTTAATGCGGGGGGGGGTATTTGCCCTCcgc                                                                                                         >  1:426753/1‑148 (MQ=255)
                                           gCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCTTAATGCGGGGGGGGGTATTTGCCCCCCCCGTg                                                                                                      >  2:610920/1‑150 (MQ=255)
                                              aTCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGa                                                                                                  <  1:487869/151‑1 (MQ=255)
                                               tCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGGGGGGAGGGTACAGCCTTATTCGGGAGGGGGTATTTCCCCTCCCCGTGccg                                                                                                   >  2:541218/1‑147 (MQ=255)
                                               tCAGTCGCAACTGTGGT‑GATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCCGGGAAGAGCACCGCATAAAGCGGGGGGGGGCATTCCCCCCCCCCgcgcccgcg                                                                                                >  1:667482/1‑146 (MQ=255)
                                                    cGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGT‑‑CAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCg                                                                                           <  2:436004/150‑1 (MQ=255)
                                                                                     aaaGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTCGGAGGGGTAATGCTTGATGTCGGGGGCGGTATTTCCCCTCCCCGGGGGGGCGAAGGGGTGCGCGAGC‑GGTTCCAAGCGCCAGGAtt                                                             >  2:14996/1‑148 (MQ=255)
                                                                                          gAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCata                                                     <  2:149826/151‑1 (MQ=255)
                                                                                           aaaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACACCTCTATTCGGGCGGGGATATTTTGCCGCCGGCGGGAGACGAAGCGGGCACCGTGCGGTTTGAAAGGGCCAGGAttttgggta                                                      >  2:411951/1‑143 (MQ=255)
                                                                                              aCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAATGCGGGAGGGGGTATTCGCCCTCCCCGTGCCGACGAAGCGGTCCGCGGGG‑GGTTCAAAGCGGCCGGAATGGTGCCTAtt                                                   >  2:425415/1‑149 (MQ=255)
                                                                                                aTTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCg                                               <  1:365754/151‑1 (MQ=255)
                                                                                                 ttGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCCCGGGCAGAAGAAGCGGGCAGCGTGC‑GGTTCAAAACGCCCGGATTGCTGCATATTGTCg                                               >  2:196018/1‑150 (MQ=255)
                                                                                                   gtgCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGGGGCGATATTTGCCCCCCGCGTGCAGACGAAGCGGGCAGCGTGC‑GGTTTAAAACGCCCGGATTGCTGCATATTGTCggg                                             >  1:762057/1‑150 (MQ=255)
                                                                                                       cGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAATGGAAAGCTTTATACGGGGGGGGGTATTTGCCCCCCGCGGGGCGACGGAGCGGGCCGCGTGG‑GGTTTAAAACGCCAGGGTTGCTGCATATTGTCggg                                             >  2:376873/1‑146 (MQ=255)
                                                                                                          ggCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAAcgg                                     <  1:699683/151‑1 (MQ=255)
                                                                                                            cAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGGGGCGATATTCGCCCCCCGCGTGCCGACGAAGCGGGCAGCGGGC‑GGTTTAAAACGCCAGGGTTGCTGCATATTTTTGGGCGGAAaggcg                                   >  1:185386/1‑151 (MQ=255)
                                                                                                              gACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGa                                 <  2:758415/151‑1 (MQ=255)
                                                                                                                     cTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGACAGCTTAATACGGGGGGCGATATTCGCCCTCCGCGGGCGGACGAAAGGGGCAGCGGGC‑GGTTTAAAACGCCCGGGTTGCTGCATATTTTCGGGCGGAACGGCGGaggcaa                            >  1:122529/1‑144 (MQ=255)
                                                                                                                       tCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGGACAGCTTAATGCGGGAGGGGATATTCGCCCTCCCCGTGCAGACGAAGCGGGCAGCGTGCGGTTTCAAA‑CGCCAGGATTGCTGCATATTGTCGGGGGGAAcggcgggaggaa                            >  1:481808/1‑141 (MQ=255)
                                                                                                                         aaCCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTaaa                      <  1:469464/151‑1 (MQ=255)
                                                                                                                          aCCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCCCCGCGTGCAGACGAAGCGGGCAGCGTGC‑GTTTCAAAGCGCCAGGGTTGCTGCATATTTTCGGGCGGAACGGCGGGCGTAAGCTaaa                      >  1:692385/1‑150 (MQ=255)
                                                                                                                                 gAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCt              <  1:618265/151‑1 (MQ=255)
                                                                                                                                         ggcggGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAagag      <  1:485270/151‑1 (MQ=255)
                                                                                                                                             gggCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTgg  >  2:59089/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGG  >  NC_000913/213192‑213483

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