Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913213,3410AG50.0% 36.6 / 14.6 38G337G (GGA→GGG) tilStRNA(Ile)‑lysidine synthetase
Reads supporting (aligned to +/- strand):  ref base A (1/18);  new base G (19/0);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 1.13e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCG‑CCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGC  >  NC_000913/213202‑213485
                                                                                                                                           |                                                                                                                                                 
ccTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAggggac                                                                                                                                            >  2:353391/1‑142 (MQ=255)
     ttACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAAt                                                                                                                                    <  1:335331/150‑1 (MQ=255)
          tttGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATTCggg                                                                                                                               >  2:94387/1‑150 (MQ=255)
            tGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAg                                                                                                                             <  2:176928/150‑1 (MQ=255)
             gggCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGc                                                                                                                           <  2:100750/151‑1 (MQ=255)
                                 gCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGGGGGGAGGGTACAGCTTAAAACGGGGGGGGGGATTTGCCCCCCCCGTg                                                                                                        >  1:567837/1‑150 (MQ=255)
                                   tATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTAAAGCTTAATACGGGGGGGGATATTCGCCCCCCCCGTGcc                                                                                                      >  1:65179/1‑149 (MQ=255)
                                    aTCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAg                                                                                                     <  1:116097/150‑1 (MQ=255)
                                      cAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGGACAGCTTAATGCGGGGGGGGGTATTTGCCCTCCGCGTGCAg                                                                                                     >  2:664565/1‑148 (MQ=255)
                                      cAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACg                                                                                                  <  2:523024/151‑1 (MQ=255)
                                         tCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGACAGCTTAATTCGGGGGGGGGTATTTGCCCTCCGCGGGCAGACGaa                                                                                                >  1:386520/1‑150 (MQ=255)
                                          cGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGaa                                                                                                >  2:107202/1‑149 (MQ=255)
                                          cGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCGCTCCGCGTGCAGACGAAGc                                                                                              <  1:352570/151‑1 (MQ=255)
                                          cGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAg                                                                                               <  1:477893/150‑1 (MQ=255)
                                          cGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAg                                                                                               <  1:488154/150‑1 (MQ=255)
                                              aCTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGGGGGGAGATTTGCCCCCCCCGGGCAGACGAAGCgg                                                                                            >  2:157137/1‑149 (MQ=255)
                                                 gtggtgGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAg                                                                                        <  2:364309/150‑1 (MQ=255)
                                                    gtggATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAgcgtg                                                                                    <  2:661460/151‑1 (MQ=255)
                                                               ccGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAgc                                                                         <  2:75029/151‑1 (MQ=255)
                                                               ccGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAg                                                                          <  1:236062/150‑1 (MQ=255)
                                                                   cACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATTCGGGAGGCGGTATTCGCCCTCCGCGTGCAGACGAAGCGGGCAGCGTGCGGTTCAAAACG‑cc                                                                     >  2:391442/1‑150 (MQ=255)
                                                                    aCCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGACAGCTTAATGCGGGGGGCGGTAATTGCCCCCCCCGTGCAGACGAAGCGGGGAGCGGGGGGTTCAAAACG‑cc                                                                     >  1:190/1‑149 (MQ=255)
                                                                      cGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCG‑CCAGGa                                                                 <  1:581458/151‑1 (MQ=255)
                                                                       gggCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCG‑CCAGGAt                                                                <  1:343394/151‑1 (MQ=255)
                                                                                 aaaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAAGGGACAGCTTAAAACGGGAGGCGATATTCGCCCTCCGCCTGCAGACGAAACGGGCAGCGGGCGTTTTCAAGCG‑CCAGGATTGCTGCata                                                       >  2:241970/1‑150 (MQ=255)
                                                                                  aaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCACAGCG‑CCAGGATTGCTGCATAtt                                                     <  1:444505/151‑1 (MQ=255)
                                                                                   aaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGAAAACTTTATGCGGGGGGGGGAAATTGCCCCCCGCGGGCCGGCGAAGAGGGCAGCGGGCGGTTCCAAACG‑CCCGGGTtgctggc                                                         >  1:506879/1‑144 (MQ=255)
                                                                                    aCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAACTTAATGCGGGAGGCGAAATTCTCCCCCCCCGTGGAGAAGAAACGGGCCGCGTGCGGTTCAAAGAG‑CCCGGGTTGCTGCATATTGt                                                   >  1:531683/1‑151 (MQ=255)
                                                                                        tgtgCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGCGGGGAATGGACAGCCTTATGCGGGGGGAGATATTACCCCTCACCGTGAGGACGAAGCGGTCCGCGCGCGCGTCAAACCG‑CCCGGATTGCTGCCTATTGTCggg                                               >  2:423059/1‑151 (MQ=255)
                                                                                               tGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGAAAGCTTTATTCGGGGGGGGATATTTGCCCCCCCCGGGCCGGCGAAACGGGCAGCGGGCGTTTCAAAACG‑CCCGGGGTGCTGCCTATTTTCggggggaa                                          >  1:4791/1‑144 (MQ=255)
                                                                                                ggCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGAGGCGATATTTGCCCTCCCCGTGCAGACGAAACGGGCAGCGGGCGGTTTAAAGCG‑CCAGGGTTGCTGCATATTTTCGGGCGGAAcgg                                       >  2:145667/1‑151 (MQ=255)
                                                                                                 gCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCG‑CCAGGATTGCTGCATATTGTCGGGCGTAAcggc                                      <  2:140811/151‑1 (MQ=255)
                                                                                                   aGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTTATGCGGGGGGGGATATTCGCCCTCCCCGCGCCGACGAAGCGGTTAGCGGGGGTTTTCAAGCG‑CCCGGATTGCTTCCTATTTTCGGGGGGAAcggcgg                                    >  1:328839/1‑151 (MQ=255)
                                                                                                       gTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGAGTTCAGCTTAATGCGGGGGGGGATATTTGCCCCCCGCGTGCAGACGAAGCGGTCAGCGGGCGTTTCAAAGCG‑CCAGGATTGCTGCATATTTTCGGGCGGAAcggcgg                                    >  2:459972/1‑147 (MQ=255)
                                                                                                             tCAACCGCTGGAATTACCGGCGGGGGTGGGGAGAGTACAGCTTAAAACGGGGGGCGATATTTGCCCTCCCCGCGCAGACAAAACGGGCAGCGGGCGTTTTAAAAAGCCCAGGA‑TGCTGGATATTTTTGGGGGGAAcggcggggggaagtt                           >  1:360594/1‑141 (MQ=255)
                                                                                                                 ccGCTGGAATTACCGGCGGGGGTGGGAAGGGGACAGCTTAATGCGGGAGGGGGAATTCGCCCTCCGCGTGCAGACGAAACGGGCAGCGGGCGGTTCAAAGCG‑CCAGGGTTGGTTCCTATTTTCGGGCGTAACGGCGGGCGGAAACTaa                         >  1:470980/1‑148 (MQ=255)
                                                                                                                                     gCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCG‑CCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGc  <  2:48825/151‑1 (MQ=255)
                                                                                                                                     gCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCG‑CCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGc  <  1:107199/151‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
CCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCG‑CCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGC  >  NC_000913/213202‑213485

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