Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,042,9930CT44.0% 49.6 / 21.3 50F193F (TTC→TTT) srkACpx stress response Thr/Ser protein kinase; MazF antagonist protein
Reads supporting (aligned to +/- strand):  ref base C (10/18);  new base T (22/0);  total (32/18)
Fisher's exact test for biased strand distribution p-value = 7.68e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGA  >  NC_000913/4042846‑4043141
                                                                                                                                                    |                                                                                                                                                    
cccGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTTAcc                                                                                                                                                   >  2:1036667/1‑151 (MQ=255)
  cGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTTAcc                                                                                                                                                   >  2:986785/1‑149 (MQ=255)
   gACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTCAcc                                                                                                                                                   >  2:529412/1‑148 (MQ=255)
       aTCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAG‑ATTTCACCCTTCt                                                                                                                                              >  2:773640/1‑149 (MQ=255)
        tCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTCAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTCACCGTTCTGc                                                                                                                                            <  2:269433/150‑1 (MQ=255)
           gTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAG‑ATTTTACCGTTCTg                                                                                                                                             >  2:604328/1‑146 (MQ=255)
           gTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTTACCGctctac                                                                                                                                            >  2:957335/1‑145 (MQ=255)
               gAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTTACCGTTCTTCGGCCACAt                                                                                                                                    >  2:47773/1‑151 (MQ=255)
               gAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTCACCGTTCTGCGGCTACAt                                                                                                                                    <  1:551406/151‑1 (MQ=255)
                 aCGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGAATTTTACCG‑TCTGCCGCTTCaacg                                                                                                                                  >  1:286382/1‑148 (MQ=255)
                        tctcATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑AATTTACCCGTCTTGGGGCACATGGAGACt                                                                                                                             >  2:581558/1‑149 (MQ=255)
                          tcATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTTACCGTTTTGCGGCTACATTGAGACTGccc                                                                                                                         >  2:673384/1‑150 (MQ=255)
                           cATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTCACCGTTCTGCGGCTACATGGAGACTGCCAc                                                                                                                        <  2:150621/151‑1 (MQ=255)
                             ttGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTTACCGTTTTTCGGCTACATGGAGAATGCCAc                                                                                                                        >  2:574491/1‑149 (MQ=255)
                               gAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGcc                                                                                                                     <  2:833385/150‑1 (MQ=255)
                                    aCGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTTACCGTTCTTCGGCTACATGGAGACTGCCACGCCgggga                                                                                                                >  1:698962/1‑148 (MQ=255)
                                         aGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAG‑AATTCACCGTTCTGCGGGCACATTGAGACTGCCCCGCCCGGGATAtt                                                                                                            >  2:586157/1‑149 (MQ=255)
                                            tGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTTACCGTTTTGCGGCTCCATGGAGACTGCCCCCCCGGGGATAtttt                                                                                                          >  2:762635/1‑146 (MQ=255)
                                            tGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTTACCCTTCTGCGGCTACATGGAGACTGCCACGCCGGGGATATTttct                                                                                                        >  1:420649/1‑150 (MQ=255)
                                                       gCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTTACCGTTCTGCGGCTACATGGAGACTGCCCCGCCGGGGATATTCTCTGGCGCGGTgg                                                                                              >  2:225543/1‑149 (MQ=255)
                                                        cTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTTGCCGTTCTGCGGCCACATTGAGACCTCCACCCCGGGGATATTTTTTTGGGCggtggt                                                                                             >  1:533534/1‑149 (MQ=255)
                                                         tACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑AATTCACCGCTCTTCGGGTACATGGGGAATGGCACGCCGGGGAAATTCTCTTGGGCGGTGGGcc                                                                                           >  2:904068/1‑150 (MQ=255)
                                                           cacTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTTACCGTTCTGCGGCTACATGGGGAGTGCCACGCCGGGGATATTTTCTTGCGCGAGGGTccca                                                                                         >  1:222856/1‑148 (MQ=255)
                                                                    cTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTTACCGTTCTGCGGCTAAATGGAGACTGCCACGCCCGGGATATTCTCTGGGGCGGTGGTCCAATGTttgttg                                                                                >  1:45606/1‑150 (MQ=255)
                                                                      tCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTTACCGTTTTGCGGCTACATGGAGGCCGCCACGCCGGGGATATTTTCTGGCGCGATGGGCCAAtgtttgtt                                                                                 >  1:115244/1‑147 (MQ=255)
                                                                      tCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAA‑ATTTCACCGTTCTGCGGCTACATGGTGACTGCCACGCCCGGAATATTCTCTTGCGGGAAGGGCCAAtgcttttt                                                                                 >  1:464312/1‑147 (MQ=255)
                                                                       ccGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGAt                                                                              <  2:941380/149‑1 (MQ=255)
                                                                                         tCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTCACCGTTCTGCGGCTACATGGGGACTGCCCCGCCGGGGATATTCTCTGGCGCGGTGGGCCAATGTTTTTTGATTTTGATGAGGCACGGAATgg                                                          >  2:880022/1‑151 (MQ=255)
                                                                                         tCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGGATATTCTCTGGCGCGGTGGTCCAATGTTTGTTGATCTGGATGATGCACGTAAt                                                            >  1:981383/1‑149 (MQ=255)
                                                                                          ccTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGt                                                         <  2:579047/151‑1 (MQ=255)
                                                                                              aaaGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAg                                                     <  1:662663/151‑1 (MQ=255)
                                                                                                    aCAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTc                                               <  1:146543/151‑1 (MQ=255)
                                                                                                         tGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTTACCGTTCTTCGGCGACCTGGAGACTGCCACGCCGGGGATATTCTCTGGCGCGGTGGGCCCATGTTTGTTGGTCTGGGTGATGCACGTAATGGGCCAGCCGtttcgg                                            >  2:1023035/1‑145 (MQ=255)
                                                                                                          gAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGAtt                                         <  2:485006/151‑1 (MQ=255)
                                                                                                          gAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGAtt                                         <  2:185510/151‑1 (MQ=255)
                                                                                                          gAGCTGATGGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGa                                           <  2:196504/149‑1 (MQ=255)
                                                                                                             cTGATTGCCGCCGTTACAGCACACTGGCGGGAAA‑ATTTTACCCGTCTGCGGCCACATGGAGACTGCCACGCCGGGAATATTCTTTTGCGCGATGGGCCCATGTTTGTTGGTCTGGGGGATGCACGTAATGGTCCCGCCGTTCAGGGTTtgt                                      >  2:851651/1‑151 (MQ=255)
                                                                                                                     cgccgTTACAGCACACTGGCGGGGAG‑ATTTTCCCGTTTTGCGGCTACATGGAGAGTGCCACGCCGGGGATATTTTTTGGCGCGATGGGCCCATTTTTGTTGGTCTGGATGATGCACGTAATGGGCCAGCCGTTCCGGGTTTTTGGGTGtt                               >  1:11327/1‑150 (MQ=255)
                                                                                                                     cgccgTTACAGCACACTGGCGGGAAG‑ATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTg                              <  2:259901/151‑1 (MQ=255)
                                                                                                                       ccgTTACAGCACACTGGCGGGAAA‑ATTTTACCCTTCTGCGGGTACATTGGGAGCGGCACGCCCGGAATATTTTCTGGCGCGGGGGTCCCATGTTTGGTGATTTTGGTGATGCACGGAAAGGTCCAGCCGTGCGGGATTTGTGGGTGtt                               >  2:184587/1‑148 (MQ=255)
                                                                                                                             cAGCACACTGGCGGGGAG‑ATTTTACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGGATATTTTCTGGCGCGATGGGCCAATGTTTGTTGTTCTGGGTGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGGTCAATgg                      >  2:1071034/1‑151 (MQ=255)
                                                                                                                                cacaCTGGCGGGAAG‑ATTTCACCCCTCTGCGGCTACATGGAGACTGCCCCGCCGGGGATATTCTCTGGCGCGATGGGCCCATGTTTGTTGGTCTGGATGATGCACGGAATGGTCCCGCCGTTCCGGATTTTTGGGTGTTGCTCAATgg                      >  2:544453/1‑148 (MQ=255)
                                                                                                                                  cacTGGCGGGAAG‑ATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATa                 <  1:773791/151‑1 (MQ=255)
                                                                                                                                       gCGGGAAG‑ATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGcc            <  2:76943/151‑1 (MQ=255)
                                                                                                                                           gAAG‑ATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGagc        <  1:561403/151‑1 (MQ=255)
                                                                                                                                           gAAG‑ATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGagc        >  2:755584/1‑151 (MQ=255)
                                                                                                                                            aaG‑ATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGagca       <  1:574622/151‑1 (MQ=255)
                                                                                                                                            aaG‑ATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGagc        >  1:568574/1‑150 (MQ=255)
                                                                                                                                              g‑ATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGagcagc     <  2:583569/151‑1 (MQ=255)
                                                                                                                                                   tCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGa  <  1:1023219/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGA  >  NC_000913/4042846‑4043141

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