Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913213,3410AG58.6% 12.8 / 34.6 29G337G (GGA→GGG) tilStRNA(Ile)‑lysidine synthetase
Reads supporting (aligned to +/- strand):  ref base A (1/11);  new base G (17/0);  total (18/11)
Fisher's exact test for biased strand distribution p-value = 3.47e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-06

GATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGC  >  NC_000913/213192‑213485
                                                                                                                                                     |                                                                                                                                                
gATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGaa                                                                                                                                                 <  2:128475/151‑1 (MQ=255)
 aTGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAg                                                                                                                                                <  2:510185/151‑1 (MQ=255)
         cccTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTa                                                                                                                                            <  2:144246/147‑1 (MQ=255)
                                           gCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGAACACCTTAAATCGGGAGGCGATATTAGCCCTcccccgg                                                                                                       >  1:139706/1‑145 (MQ=255)
                                                cAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACg                                                                                                 <  1:509820/151‑1 (MQ=255)
                                                      cAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAATACGGGAGGGGATAATCGCCCCCCCCGGGCCGACGAAACgg                                                                                           >  2:540059/1‑151 (MQ=255)
                                                           gtggtgGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGGGGCGGTATTCGCCCTCCGCGGGCAGACGAAGCgg                                                                                           >  2:415700/1‑146 (MQ=255)
                                                            tggtggATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAACTTAAAGCGGGGGGGGGTATTTGCCCTCCGCGGGCAGACGAAGCGGGCAg                                                                                       >  2:414152/1‑149 (MQ=255)
                                                                 gATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAgcgtgcgt                                                                                <  2:90290/151‑1 (MQ=255)
                                                                         ccGTCACCGGGCAAAGCGAAAACATTGTGCCGCGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGCGGGGAGGGTACCGCTTTATGCGGGGGGGGGACTTTGCCCCCCGCCGGGCGGCGAAGCGGGCCGCGGGCGGTTCaaa                                                                          >  1:392759/1‑149 (MQ=255)
                                                                          cGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCGCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAgcg                                                                       <  1:415937/151‑1 (MQ=255)
                                                                              aCCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGAGGCGGTATTCGCCCTCCGCGTGCAGACGAAGCGGGCAGCGGGCGGTTTaaa                                                                          >  2:217847/1‑144 (MQ=255)
                                                                               ccGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGCTTAATTCGGGAGGGGGTATTTGCCCCCCGCGGGCAGACGAAGCGGGCAGCGGGCGTTTCAAAGCGCCAgg                                                                  >  2:199125/1‑151 (MQ=255)
                                                                                 gggCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGAGGGGATATTTGCCCTCCGCGTGCAGAAGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAg                                                                   >  2:48763/1‑148 (MQ=255)
                                                                                          gAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGTAAAGCTTTATGCGGGGGGGGGTATTTGCCCCCCCCGTGCAGACGAAGCGGGCAGAGGGGGGTTCAAAACGCCAGGAGTGCTGCa                                                         >  1:484870/1‑149 (MQ=255)
                                                                                          gAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGGACAGCTTAATGCGGGAGGCGGTATTTGCCCTCCGCGGGCCGACGAAGCGGGCCGCGGGCGGTTTAAAGCGCCAGGGTTTCTGCat                                                        >  2:144226/1‑150 (MQ=255)
                                                                                          gAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGCTTAATTCGGGAGGGGATATTCGCCCTCCCCGGGCAGACGAAGCGGTCAGCGGGGGTTTCAAAGAGCCAGGATtggt                                                            >  2:65042/1‑144 (MQ=255)
                                                                                          gAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGCGGGGAGCGTAACGCCTAATGCGGGGGGCGATATTCGCCCCCCGCCGGCCGACGAAGCGGGCAGCGGGCGGTTTAAAACGCCAGGATtgct                                                            >  2:453960/1‑146 (MQ=255)
                                                                                          gAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGAGGGGGTATTCGCCCCCCGCGGGCAGACGAAGCGGGCAGCGGGCGGTTTAAAAAGCCCGGATTGCTGCat                                                        >  2:195275/1‑150 (MQ=255)
                                                                                                aTTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGGGGGGGTATTTGCCCTCCGCGTGCAGAAGAAGCGGTTCGCGGGGGGTTTTAAGCGCCAGGATTGGCGCATATTTTCg                                                 >  2:52706/1‑151 (MQ=255)
                                                                                                aTTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCg                                                 <  1:86317/151‑1 (MQ=255)
                                                                                                               aCGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGAc                                  <  1:291348/151‑1 (MQ=255)
                                                                                                                   ggCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTa                               <  1:195211/150‑1 (MQ=255)
                                                                                                                     cTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGAACGCTTTATGCGGGGGGCGATATTCGCCCTCCGCGGGCCGACGAAACGGGCAGCGGGCGGTTCCAAACGCCCGGAGTGCCGCATATTTTTGGGCGGAAAGGGGGGCGGaa                              >  1:105048/1‑149 (MQ=255)
                                                                                                                     cTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAg                             <  2:18716/150‑1 (MQ=255)
                                                                                                                           ccGCTGGAATTACCGGCGGGGGTGGGGAGTGTACAGCTTAATACGGGGGGCGATATTCCCCCCCCGCGGGCAGACGAAGCGGGCAGCGGGCGTTTCCAAGGGCCCGGGTTGCTGGATATTGTCGGGGGGAACGGCGGGCGTAAGCTaaa                        >  1:201299/1‑149 (MQ=255)
                                                                                                                            cGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCCTAATTCGGGGGGCGGAATTCGCCCTCCGCCGGCCGACGAAGGGGGCAGCGGGCGGTTCCAAGCGCCAGGGTTGCTGCCTATTTTCGGGGGGAACGGGGGGCGGAAACGaaa                        >  1:176384/1‑148 (MQ=255)
                                                                                                                                         ggcggGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAag          >  2:207743/1‑149 (MQ=255)
                                                                                                                                          gcggGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGc       <  2:231867/151‑1 (MQ=255)
                                                                                                                                                cTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTggtg  <  1:388414/150‑3 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGC  >  NC_000913/213192‑213485

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