Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,043,0300AG44.7% 43.9 / 10.2 38N206D (AAT→GAT) srkACpx stress response Thr/Ser protein kinase; MazF antagonist protein
Reads supporting (aligned to +/- strand):  ref base A (12/9);  new base G (17/0);  total (29/9)
Fisher's exact test for biased strand distribution p-value = 1.95e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGA  >  NC_000913/4042884‑4043163
                                                                                                                                                  |                                                                                                                                     
gCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAGTTCACCGTTCTGCGGGTACATGGAGACTGCCACGCCGGGGat                                                                                                                                     >  2:479408/1‑149 (MQ=255)
  aaGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGATTTTAACCTTCTGCGGCTACACAGAGACTGCCACCCCGGGGAtttt                                                                                                                                  >  1:588620/1‑150 (MQ=255)
          tGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTCCCCGTTCTTCGGCTACATGGAGACTGCCACGCCGGGGATATTCTCTTgcgc                                                                                                                         >  2:366689/1‑151 (MQ=255)
                       cTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTTACCGTTCTGCGGCTACATGGAGACTGCCCCGCCGGGAATATTTTCTGGCGCGATGGTCCAAtg                                                                                                             >  2:13292/1‑150 (MQ=255)
                                   gggTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTgg                                                                                                <  1:939900/151‑1 (MQ=255)
                                    ggTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGga                                                                                               <  2:124557/151‑1 (MQ=255)
                                    ggTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGGTCTGCGGGTACATGGGGACTGCCCCGCCGGGGATATTCTCTTGCGCGATGGTCCAATGTTTGTTGATCTgg                                                                                                >  2:625519/1‑150 (MQ=255)
                                      ttGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGGAAATTTTCTGGCGCGATGGTCCAATGTTTGTTGTTCt                                                                                                  >  1:402894/1‑146 (MQ=255)
                                        gAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTTATCGGTCTGCCGATACATCGGGACTTGCACCCCCGGAATTTTTTCTGGCCCGGTGGGCCAAAGGTTGTTGTTCTgggt                                                                                              >  1:381195/1‑146 (MQ=255)
                                         aaaGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGATTTTACCGGTCTGCGGCCACATGGGGAACTCCACGCCGGGGAGATTTTCTGGCGCGAAGGGCCCATGTTTGTTGATCTTggtgat                                                                                           >  1:485463/1‑150 (MQ=255)
                                          aaGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTCACCGGTCTTCGGCTTCATTGAGACTGCCACGCCGGGGATATTCTCTTGCGCGGTGGGCCCATGTTTGTTGATCTgg                                                                                                >  2:655104/1‑144 (MQ=255)
                                                 aTTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCTACATGGGGACTGCCCCGCCGGGGATATTCTCTGGCGCGGTGGGCCAATGTTTTTTGATTTGGATGGTGGACGTaa                                                                                   >  2:142770/1‑150 (MQ=255)
                                                      tGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGATTTTACCGTTCTGCGGCTACATTGAGACTGCCACGCCGGGGAAATTTTCTTGCGCGGTGGGTCCATGTTTGTTGATCTGGGTGATGCACGTAATgg                                                                                >  1:656574/1‑148 (MQ=255)
                                                      tGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCTACCTTGAGACTGCCCCGCCGGGGATATTCTCTGGCGCGGTGGGCCCATGTTTGTTGGTTTGGATGATGCACGTAATgg                                                                                >  1:533766/1‑148 (MQ=255)
                                                       gAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGCAGATTTCACCGCTTTGCGGCCACTTAGAGACAGCCACGCCGGGGAAATTTTCTTGCGGGGTGGTCCAATGTTTGTTGATATGGGTGATGCACGGaaaag                                                                                >  1:465972/1‑144 (MQ=255)
                                                         aaGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGc                                                                          <  2:443884/151‑1 (MQ=255)
                                                          aGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGATTTTACCCTTTTTCGGCTACATGGGGACTGCCACCCCCGGGATATTTTCTGGCGCGGGGGGCCCATGTTTGTTTATCTGGGTGATGCACGTAATGGTCCAGcc                                                                         >  1:105038/1‑151 (MQ=255)
                                                           gCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCCCTTTTCGGCTACATGGAGACTGCCCCCCCGGGAAAATTTTCTGGCGCGGTGGGCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGcc                                                                         >  1:612763/1‑150 (MQ=255)
                                                             gACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAAGTTACCGTTCTGCGGCTTCATGGTGAAGGCCACGCCCGGGATATTTTCTGGCGCGATGGGCCCATGTTTTTTGTTCTGGAGGAAGCACGAAAAAGTCCAGCCg                                                                        >  1:99678/1‑149 (MQ=255)
                                                                     aGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTTTTCGGCTACATGGAGAATGCCACGCCGGGGAAATTTTCTTGCGCGGTGGGCCAATGTTTGTTGAATTGGATGATGGACGTTATGGGCCAGCCGTTCAGGAttt                                                              >  2:1047326/1‑151 (MQ=255)
                                                                          aTTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTAACGCTCTTCGGGTATATCGTGGGTGACCCGCCGGGGATATTTTTCTGCGCGGCCGTGCAACGATTGTTGTTCTGGAAGCAGGAAGTAAAGGTTCAGCTGGTCAGGATTtgttg                                                          >  2:677551/1‑148 (MQ=255)
                                                                               cgccgTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTg                                                    <  2:921993/151‑1 (MQ=255)
                                                                                  cgTTACAGCACACTGGCGGGGAGATTTTACCGTTCTGCGGGTACATGGAGACTGCCACGCCGGGGATATTTTCTGGCGCGATGGTCCCATGTTTTTTTATCTGGGTGATGCACGTAATGGTCCAGGCGTTCTGGGTTTGTGGGTGtt                                                     >  1:942971/1‑147 (MQ=255)
                                                                                  cgTTACAGCACACTGGCGGGAAGATTTCACCGGTTTGCGGCTACATGGAGACCGCCACGCCGGGGATATTTTCTGGCGCGATGGGCCCATGTTTGTTGGTCTGGGTGATGCACGTAATGGTCCAGCCGTTTAGGGTTTGTGGGTGtt                                                     >  2:472189/1‑147 (MQ=255)
                                                                                                    ggAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCg                                 >  1:144251/1‑149 (MQ=255)
                                                                                                      aaGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGagc                              >  1:463772/1‑150 (MQ=255)
                                                                                                           ttCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCgg                         <  1:840663/150‑1 (MQ=255)
                                                                                                             cACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATg                      <  1:711390/151‑1 (MQ=255)
                                                                                                              aCCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGc                     >  1:659689/1‑151 (MQ=255)
                                                                                                                   tCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGGCCAATGTTTGTTGATCTGGATTATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGGTCAATGGCGATAAAGCCGAGCAGCGGATGCAACt                 >  1:635839/1‑150 (MQ=255)
                                                                                                                      gCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGaa             <  2:1028107/151‑1 (MQ=255)
                                                                                                                      gCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa              >  1:218205/1‑150 (MQ=255)
                                                                                                                       cGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGaa             >  2:980196/1‑150 (MQ=255)
                                                                                                                        ggCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAAc           >  2:212560/1‑151 (MQ=255)
                                                                                                                             cATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACtatta      <  2:117189/151‑1 (MQ=255)
                                                                                                                                 gagaCTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGa  <  2:646076/151‑1 (MQ=255)
                                                                                                                                   gaCTGCCACGCCGGGACTATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTg   >  1:122693/1‑148 (MQ=255)
                                                                                                                                   gaCTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGa  >  2:980039/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                     
GCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGA  >  NC_000913/4042884‑4043163

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