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

CCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCAC‑GCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGAT  >  NC_000913/4042846‑4043142
                                                                                                                                                   |                                                                                                                                                      
cccGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTAcc                                                                                                                                                     >  2:640016/1‑151 (MQ=255)
 ccGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGATTTTAcc                                                                                                                                                     >  1:925380/1‑150 (MQ=255)
  cGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTTACCg                                                                                                                                                    >  1:445780/1‑150 (MQ=255)
   taCCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAGTTTACCGtt                                                                                                                                                  >  1:400032/2‑151 (MQ=255)
    aCCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTTACCgct                                                                                                                                                  >  1:503457/1‑148 (MQ=255)
     ccATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCt                                                                                                                                                <  2:876704/151‑1 (MQ=255)
      cATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCCTTCt                                                                                                                                                >  1:868328/1‑150 (MQ=255)
      cATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTCACCGtt                                                                                                                                                  >  2:935707/1‑148 (MQ=255)
           gTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTTAACGTTCttggg                                                                                                                                            >  2:936953/1‑145 (MQ=255)
            tttGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCgg                                                                                                                                            >  1:897859/1‑148 (MQ=255)
            tttGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTTCGGCt                                                                                                                                          >  2:634509/1‑150 (MQ=255)
               gAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACAt                                                                                                                                      <  2:57893/151‑1 (MQ=255)
                      tatCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCCACATgggg                                                                                                                                  >  2:299498/1‑146 (MQ=255)
                            aTTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCAc                                                                                                                          <  1:1016245/150‑1 (MQ=255)
                                 gCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGATTTTACCGTTTTGTGGGTTCATGGGGACATCCCA‑GCCggg                                                                                                                   >  1:736499/1‑151 (MQ=255)
                                     cGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCAC‑GCCGGGAat                                                                                                                <  2:861303/150‑1 (MQ=255)
                                        aaGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGGTCTGCGGCTACATGGAGAATGCCAC‑GCCGGGAAAAtt                                                                                                             >  2:584048/1‑150 (MQ=255)
                                        aaGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCAC‑GCCGGGAATATtc                                                                                                            <  2:114573/151‑1 (MQ=255)
                                            tGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGGCACATGGAGACTGCCAC‑GCCGGGGAAATTctct                                                                                                         >  2:970731/1‑150 (MQ=255)
                                             gTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTTCCGCTCTTCGGCTACATGGAGAATTCCCA‑GCCCGGAATATTTTCTgg                                                                                                       >  2:311904/1‑151 (MQ=255)
                                                tGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTTTGCGGCTACATGGGGACTGCCCC‑GCCGGGAATATTctct                                                                                                         >  2:523686/1‑146 (MQ=255)
                                                  aGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCAC‑GCCGGGGATATTCTCTGgcg                                                                                                     >  2:419351/1‑148 (MQ=255)
                                                    gACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCCACCTGGAGACTGCCCC‑CCCGGGGAAATTCTCTGGcgcggt                                                                                                 >  1:917570/1‑148 (MQ=255)
                                                       gCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCTACATGGAGACTGCCAC‑GCCGGGGATATTTTCTTGCGCGGTgg                                                                                               >  2:92181/1‑149 (MQ=255)
                                                          acacTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCAC‑GCCGGGAATATTCTCTGGCGCGATGGTCCaa                                                                                          <  2:976103/151‑1 (MQ=255)
                                                              tGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGAATTTACCGTTCTGGGGCTACCTGGGGACTGCCCC‑GCCGGGGATATTTTCTTGGGCGGTGGTCCCAtgtt                                                                                      >  1:365360/1‑151 (MQ=255)
                                                                 tACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCTACATGGGGACTGCCAC‑GCCGGGGATATTCTCTGGCGCGATGGTCCAAtgtttgt                                                                                   >  2:508603/1‑151 (MQ=255)
                                                                       ccGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCAC‑GCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCt                                                                             <  1:18124/151‑1 (MQ=255)
                                                                         gggTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCCC‑GCCGGGGATATTTTCTGGCGCGGTGGTCCAATGTTTGTTGGTCTg                                                                            >  1:405883/1‑150 (MQ=255)
                                                                             tGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGGTTTGCGGCTACATTGAGACTGCCCC‑GCCGGGGATATTTTCTGGCGCGATGGGCCAATGTTTGTTGATCTGgatg                                                                        >  2:679859/1‑150 (MQ=255)
                                                                                  gCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCCC‑GCCGGGGATATTTTCTGGCGCGGTGGGCCAAAGGTTGTTGATCTGggtg                                                                        >  1:827957/1‑145 (MQ=255)
                                                                                      cATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCAC‑GCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTaa                                                              <  1:935430/151‑1 (MQ=255)
                                                                                      cATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCAC‑GCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTaa                                                              <  1:508423/151‑1 (MQ=255)
                                                                                          ccTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCAC‑GCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATgg                                                           <  2:579470/150‑1 (MQ=255)
                                                                                             gAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCTACATGGAGACTGCCAC‑GCCGGGAATATTCTCTTGGGCGGTGGTCCAATGTTTGTTGGTCTGGGTGATGCACGTAATGGTcc                                                        >  2:706528/1‑150 (MQ=255)
                                                                                             gAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCAC‑GCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCa                                                       <  1:522237/151‑1 (MQ=255)
                                                                                              aaaGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCAC‑GCCGGGGATATTTTCTGGCGCGATGGTCCCATGTTTGTTGATCTGGATGATGCACGTAATGGTcc                                                        >  2:49667/1‑149 (MQ=255)
                                                                                                  cGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGGTCTGCGGCTACCTGGAGACTGCCAC‑GCCGGGGAAATTTTCTTGCGCGATGGTCCCATGTTTGTTGATCTGGGTGATGCACGTAATGGTCCAGcc                                                    >  1:128457/1‑149 (MQ=255)
                                                                                                   gACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCAC‑CCCGGGGATATTCTCTGGCGCGGTGGTCCAATGTTTGTTGATTTGGGTGATGCACGGAATGGGCCAGCCGtt                                                 >  2:282178/1‑151 (MQ=255)
                                                                                                    aCAGATGGGCGGATTGCCGCCGTTACAGCACCCTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCAC‑GCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGtt                                                 <  1:915132/150‑1 (MQ=255)
                                                                                                     cAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCAC‑GCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCa                                               <  2:500121/151‑1 (MQ=255)
                                                                                                     cAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTTACCGTTCTGGGGCCACATGGAGACTGCCAA‑GCCGGGGATATTCTCTTGCGCGAAGGTCCAATGTTTGTTGGTATGGATGATGCACGTAATGGTCCAGCCGtt                                                 >  1:996407/1‑149 (MQ=255)
                                                                                                      aGATGAGCTGCTTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCCCCCTTCTGCGGCTACCTGTAGGCCCCCACTGCCGGG‑ATATTTTCTCGCGCGCGGGTTCAATGTTTTTTGATCTTGATGATGTGCGGAATAGTCccgccggt                                                 >  1:183496/1‑146 (MQ=255)
                                                                                                          gAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCAC‑GCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGAt                                           <  1:448781/150‑1 (MQ=255)
                                                                                                             caaaTTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCAC‑GCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGAtt                                          <  2:697731/145‑1 (MQ=255)
                                                                                                                      gccgTTACAGCACACTGGCGGGAAGAATTTACCGTTCTGCGGCTACATGGAGACTGCCCC‑GCCGGGAATATTTTCTGGCGCGATGGGCCCAAGTTTGTTGATCTGGGTGATGCACGTAATTGTCCAGCCGTTCAGGGTTTTTGGGTGtt                                >  1:917822/1‑149 (MQ=255)
                                                                                                                       ccgTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCAC‑GCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCt                             <  2:886489/151‑1 (MQ=255)
                                                                                                                            aCAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCAC‑GCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATg                        <  1:205277/151‑1 (MQ=255)
                                                                                                                              aGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCAC‑GCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGc                      <  2:1005135/151‑1 (MQ=255)
                                                                                                                                    cTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCAC‑GCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGATGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATaa                 <  1:634306/150‑1 (MQ=255)
                                                                                                                                                  tCACCGTTCTGCGGCTACATGGAGACTGCCAC‑GCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGa   <  2:464402/150‑1 (MQ=255)
                                                                                                                                                  tCACCGTTCTGCGGCTACATGGAGACTGCCAC‑GCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGa   <  1:359058/150‑1 (MQ=255)
                                                                                                                                                   cACCGTTCTGCGGCTACATGGAGACTGCCAC‑GCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGAt  <  2:254731/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
CCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCAC‑GCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGAT  >  NC_000913/4042846‑4043142

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