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

GCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGC  >  NC_000913/213195‑213485
                                                                                                                                                  |                                                                                                                                                 
gCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAgt                                                                                                                                                <  2:230572/150‑1 (MQ=255)
gCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAgt                                                                                                                                                <  1:297832/150‑1 (MQ=255)
   tCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCCTGGGAGGGTAc                                                                                                                                            >  1:56096/1‑151 (MQ=255)
        cTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCtt                                                                                                                                       >  2:403921/1‑151 (MQ=255)
                   tGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAg                                                                                                                             <  2:9013/150‑1 (MQ=255)
                       cgcgTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGAAAGCTTAATGCGGGGGGCg                                                                                                                          >  2:128082/1‑149 (MQ=255)
                         cgTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCTTAATACGGGGGGCg                                                                                                                          >  1:542799/1‑147 (MQ=255)
                           tttGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGAGAGGTTTATGCCGGGGGCgga                                                                                                                        >  1:398865/1‑145 (MQ=255)
                             tGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCg                                                                                                                   <  1:131203/150‑1 (MQ=255)
                                           aTCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGa                                                                                                    <  2:367120/151‑1 (MQ=255)
                                                 cGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGAGAGCTTAATGCGGGGGGGGGGATTTGCCCCCTGCGTGCAGACgga                                                                                                >  1:16414/1‑147 (MQ=255)
                                                        gtggtgGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCa                                                                                         <  1:442676/149‑1 (MQ=255)
                                                               aTTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGAAATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGtt                                                                               <  2:398902/151‑1 (MQ=255)
                                                                  aaaTCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAAAGCGGGGGGGGGTATTTTCCCCCCCCGTGCCGACGAAGCGGGCAGC‑GGGCGgtt                                                                              >  1:666527/1‑149 (MQ=255)
                                                                    aTCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAAAGCGGGGGGGGGTATTTGCCCTCCGCGTGCAGACGAAGCGGGCAGC‑GTGGGGTTCa                                                                            >  2:20455/1‑149 (MQ=255)
                                                                      ccGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAgc                                                                        <  1:496507/151‑1 (MQ=255)
                                                                       cGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGAAAGCCTAATGCGGGGGGGGGGATTTGCCCTCCCCGTGCAGACGAAAGGGGCAgc‑ggggggtttaaaacg                                                                       >  1:263991/1‑137 (MQ=255)
                                                                            ccGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAg                                                                   <  2:316513/150‑1 (MQ=255)
                                                                             cGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGAAAGCCTAATTCGGGGGGGGGTATTTGCCCTCCCCGTGCCGACGAAGCGGGGAGCGGGGGGTTTCaaaagcc                                                                      >  1:122835/1‑143 (MQ=255)
                                                                               ggCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGAtt                                                               <  2:538883/151‑1 (MQ=255)
                                                                                gCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATtg                                                              <  2:345419/151‑1 (MQ=255)
                                                                                       gAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGGGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGGCAGC‑GTGCGGTTCAAAGCGCCCGGATTTCTGCa                                                         >  2:629274/1‑149 (MQ=255)
                                                                                       gAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGAAAGCCTAATGCGGGGGGGGATATTTGCCCCCCGCGGGCCGACGAAGGGGGCCGC‑GTGGGGTTCAAAACGCCCGGGTTGCTGCat                                                        >  1:103568/1‑150 (MQ=255)
                                                                                                   ccGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTa                                           <  1:273638/151‑1 (MQ=255)
                                                                                                     gTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAAc                                         <  2:522857/151‑1 (MQ=255)
                                                                                                       ggCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAAcg                                        <  2:103588/150‑1 (MQ=255)
                                                                                                        gCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAAcggc                                      <  1:44656/151‑1 (MQ=255)
                                                                                                          aGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGCCTAATGCGGGGGGGGGAATTTGCCCCCCCCGCGCAGACGAAGCGGGCAGC‑GGGCGGTTCAAAACGCCCGGATTGCTGCATATTTTTGGGCGGAAcggcgg                                    >  1:66854/1‑151 (MQ=255)
                                                                                                            aCGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGAc                                  <  1:240586/151‑1 (MQ=255)
                                                                                                              gTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGTACAGCTTAAATCGGGGGGCGATATTTGCCCTCCGCGTGCAGGCGAAGCGGGCCGC‑GTGGGTTTCAAAGCGCCCGGGATGCTGCCTATTTTCGGGCGTAACCGCGGACg                                 >  2:215137/1‑150 (MQ=255)
                                                                                                                 gCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTTATGCGGGGGGGGGTATTCGCCCTCCGCGTGCAGACGAAGCGGGCAGC‑GGGCGTTTCAAAACGCCCGGATTGCTGCATAATTTCGGGCGGAAcggcgg                                    >  1:101586/1‑144 (MQ=255)
                                                                                                                    tCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGGGGCGGTATTCCCCCCCCCCCGGCAGAAGAAGCGGGCAGC‑GGGGGGTTCAAAACGCCCGGATTGCTGCCTATTGTCGGGCGGAACGGCGGGCGGAAggt                           >  1:149270/1‑148 (MQ=255)
                                                                                                                    tCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCt                           <  1:33956/150‑1 (MQ=255)
                                                                                                                      aaCCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGAGGCGATATTTGCCCTCCGCGGGCAGACGAAGCGGGCAGC‑GGGCGGTTTAAAACGCCCGGGTTGCTGCATATTGTCGGGCGGAAAGGCGGGCGGAAGCTaa                         >  2:75508/1‑150 (MQ=255)
                                                                                                                             ggAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGaaaa                   <  1:364961/149‑1 (MQ=255)
                                                                                                                             ggAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATc                 <  1:5278/151‑1 (MQ=255)
                                                                                                                                            gCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGc  <  2:407897/151‑1 (MQ=255)
                                                                                                                                            gCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGc  >  2:673628/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGC  >  NC_000913/213195‑213485

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