Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,042,9930CT60.0% ‑5.3 / 17.9 26F193F (TTC→TTTsrkACpx stress response Thr/Ser protein kinase; MazF antagonist protein
Reads supporting (aligned to +/- strand):  ref base C (8/2);  new base T (15/0);  total (24/2)
Fisher's exact test for biased strand distribution p-value = 1.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-01

CCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGAT  >  NC_000913/4042847‑4043142
                                                                                                                                                  |                                                                                                                                                     
ccGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGATTTTAc                                                                                                                                                     >  2:715644/1‑149 (MQ=255)
    ccATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTTAcccgtct                                                                                                                                               >  2:689414/1‑146 (MQ=255)
    ccATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAAATTTACCCGtt                                                                                                                                                 >  2:706389/1‑149 (MQ=255)
          gTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCCTTTTGCgg                                                                                                                                           >  2:30122/1‑149 (MQ=255)
           tttGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCgg                                                                                                                                           >  2:125730/1‑148 (MQ=255)
                  gAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTTCCGCTACa                                                                                                                                      >  1:721939/1‑146 (MQ=255)
                  gAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCTACATgg                                                                                                                                   >  2:67409/1‑149 (MQ=255)
                     tatCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAGTTTACCCCTCTGCGGCCACATgggg                                                                                                                                 >  2:163785/1‑146 (MQ=255)
                      atCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGATATCCCCCTTTTGCGGCATCCTGgagaa                                                                                                                               >  2:669469/1‑148 (MQ=255)
                      atCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTTTCGATCTGCGGCTACAAggggcccg                                                                                                                             >  1:817034/1‑145 (MQ=255)
                        ctcATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCCTTCTGCGGCTACATGGAGACTTcc                                                                                                                           >  2:300856/1‑151 (MQ=255)
                             tGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTTACCGGTTTGCGGCTATACGGCGACCGc                                                                                                                            >  1:281103/1‑145 (MQ=255)
                                            gTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGGATATtttc                                                                                                          >  1:660206/1‑148 (MQ=255)
                                            gTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTg                                                                                                        <  2:754461/150‑1 (MQ=255)
                                                        tACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGCTTTTACCGTTCTGCGGGTACATTGAGAATGCCACGCCGGGGATATTTTCTTGcgcggggggac                                                                                            >  1:473240/1‑143 (MQ=255)
                                                                      ccGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTAACGTTCCTCGGGTACATGGAGACTGCCCCACCCGGGAAATTTTCTTGCGCGATGGTCCCATGTTTGTTGAttt                                                                             >  1:327473/1‑149 (MQ=255)
                                                                            tGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGGTTTTACCGGTCTTCGGCTACAAGGGGACCGCCACGCCGGGAATATTTTCTGGGGCGGTGGGCCAATGGTTGTTTATCTGggtg                                                                        >  2:182361/1‑150 (MQ=255)
                                                                             gAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCCTTTTGCGGCCACACGGAGACTGCCACGCCGGGGATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGggtg                                                                        >  2:530645/1‑149 (MQ=255)
                                                                              aaaGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCTACATGGGGACTGCCCCCCCGGGGATATTCTCTGGCGCGGTGGTCCAATGTTTTTTGATCTGggtgat                                                                      >  2:64537/1‑150 (MQ=255)
                                                                                                              gATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTtgt                                       <  1:531744/149‑1 (MQ=255)
                                                                                                                      ccgTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGGATATTTTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGGTGATGCACGGAATGGGCCAGCCGTTCAGGGTTTTTGGATGTTg                               >  2:225488/1‑149 (MQ=255)
                                                                                                                       cgTTACAGCACACTGGCGGGGAGAATTTACCGTTCTGCGGCTACATTGAGGCCGCCACGCCGGGGATATTTTCCGGCGCGGTGGTCCACTGTTTGTTGTTTTGGGTTAAGCCCGCAATGGTCCCGCCGTTCTGGATTTTTGGGTGTTtctc                            >  1:221133/1‑151 (MQ=255)
                                                                                                                       cgTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGGGACCGCCCCGCCGGGGATATTCTCTGGCGCGAAGGGCCAATGTTTGTTTGTCTGGGTGATGCACGTAATGGTACAGCCCGTCAGGGTTTTTGGATGTTGCTc                            >  2:212283/1‑151 (MQ=255)
                                                                                                                                         ggAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCg            >  1:251828/1‑149 (MQ=255)
                                                                                                                                             gATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGagcagc      >  2:245973/1‑151 (MQ=255)
                                                                                                                                                 tCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGAt  >  1:677166/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGAT  >  NC_000913/4042847‑4043142

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: