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→TTTsrkACpx 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

CCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGAT  >  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)
                                 gCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGATTTTACCGTTTTGTGGGTTCATGGGGACATCCCAGCCggg                                                                                                                   >  1:736499/1‑151 (MQ=255)
                                     cGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAat                                                                                                                <  2:861303/150‑1 (MQ=255)
                                        aaGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGGTCTGCGGCTACATGGAGAATGCCACGCCGGGAAAAtt                                                                                                             >  2:584048/1‑150 (MQ=255)
                                        aaGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATtc                                                                                                            <  2:114573/151‑1 (MQ=255)
                                            tGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGGCACATGGAGACTGCCACGCCGGGGAAATTctct                                                                                                         >  2:970731/1‑150 (MQ=255)
                                             gTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTTCCGCTCTTCGGCTACATGGAGAATTCCCAGCCCGGAATATTTTCTgg                                                                                                       >  2:311904/1‑151 (MQ=255)
                                                tGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTTTGCGGCTACATGGGGACTGCCCCGCCGGGAATATTctct                                                                                                         >  2:523686/1‑146 (MQ=255)
                                                  aGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGGATATTCTCTGgcg                                                                                                     >  2:419351/1‑148 (MQ=255)
                                                    gACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCCACCTGGAGACTGCCCCCCCGGGGAAATTCTCTGGcgcggt                                                                                                 >  1:917570/1‑148 (MQ=255)
                                                       gCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGGATATTTTCTTGCGCGGTgg                                                                                               >  2:92181/1‑149 (MQ=255)
                                                          acacTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCaa                                                                                          <  2:976103/151‑1 (MQ=255)
                                                              tGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGAATTTACCGTTCTGGGGCTACCTGGGGACTGCCCCGCCGGGGATATTTTCTTGGGCGGTGGTCCCAtgtt                                                                                      >  1:365360/1‑151 (MQ=255)
                                                                 tACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCTACATGGGGACTGCCACGCCGGGGATATTCTCTGGCGCGATGGTCCAAtgtttgt                                                                                   >  2:508603/1‑151 (MQ=255)
                                                                       ccGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCt                                                                             <  1:18124/151‑1 (MQ=255)
                                                                         gggTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCCCGCCGGGGATATTTTCTGGCGCGGTGGTCCAATGTTTGTTGGTCTg                                                                            >  1:405883/1‑150 (MQ=255)
                                                                             tGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGGTTTGCGGCTACATTGAGACTGCCCCGCCGGGGATATTTTCTGGCGCGATGGGCCAATGTTTGTTGATCTGgatg                                                                        >  2:679859/1‑150 (MQ=255)
                                                                                  gCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCCCGCCGGGGATATTTTCTGGCGCGGTGGGCCAAAGGTTGTTGATCTGggtg                                                                        >  1:827957/1‑145 (MQ=255)
                                                                                      cATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTaa                                                              <  1:935430/151‑1 (MQ=255)
                                                                                      cATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTaa                                                              <  1:508423/151‑1 (MQ=255)
                                                                                          ccTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATgg                                                           <  2:579470/150‑1 (MQ=255)
                                                                                             gAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTTGGGCGGTGGTCCAATGTTTGTTGGTCTGGGTGATGCACGTAATGGTcc                                                        >  2:706528/1‑150 (MQ=255)
                                                                                             gAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCa                                                       <  1:522237/151‑1 (MQ=255)
                                                                                              aaaGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGGATATTTTCTGGCGCGATGGTCCCATGTTTGTTGATCTGGATGATGCACGTAATGGTcc                                                        >  2:49667/1‑149 (MQ=255)
                                                                                                  cGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGGTCTGCGGCTACCTGGAGACTGCCACGCCGGGGAAATTTTCTTGCGCGATGGTCCCATGTTTGTTGATCTGGGTGATGCACGTAATGGTCCAGcc                                                    >  1:128457/1‑149 (MQ=255)
                                                                                                   gACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACCCCGGGGATATTCTCTGGCGCGGTGGTCCAATGTTTGTTGATTTGGGTGATGCACGGAATGGGCCAGCCGtt                                                 >  2:282178/1‑151 (MQ=255)
                                                                                                    aCAGATGGGCGGATTGCCGCCGTTACAGCACCCTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGtt                                                 <  1:915132/150‑1 (MQ=255)
                                                                                                     cAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCa                                               <  2:500121/151‑1 (MQ=255)
                                                                                                     cAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTTACCGTTCTGGGGCCACATGGAGACTGCCAAGCCGGGGATATTCTCTTGCGCGAAGGTCCAATGTTTGTTGGTATGGATGATGCACGTAATGGTCCAGCCGtt                                                 >  1:996407/1‑149 (MQ=255)
                                                                                                      aGATGAGCTGCTTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCCCCCTTCTGCGGCTACCTGTAGGCCCCCACTGCCGGGATATTTTCTCGCGCGCGGGTTCAATGTTTTTTGATCTTGATGATGTGCGGAATAGTCccgccggt                                                 >  1:183496/1‑146 (MQ=255)
                                                                                                          gAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGAt                                           <  1:448781/150‑1 (MQ=255)
                                                                                                             caaaTTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGAtt                                          <  2:697731/145‑1 (MQ=255)
                                                                                                                      gccgTTACAGCACACTGGCGGGAAGAATTTACCGTTCTGCGGCTACATGGAGACTGCCCCGCCGGGAATATTTTCTGGCGCGATGGGCCCAAGTTTGTTGATCTGGGTGATGCACGTAATTGTCCAGCCGTTCAGGGTTTTTGGGTGtt                                >  1:917822/1‑149 (MQ=255)
                                                                                                                       ccgTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCt                             <  2:886489/151‑1 (MQ=255)
                                                                                                                            aCAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATg                        <  1:205277/151‑1 (MQ=255)
                                                                                                                              aGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGc                      <  2:1005135/151‑1 (MQ=255)
                                                                                                                                    cTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGATGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATaa                 <  1:634306/150‑1 (MQ=255)
                                                                                                                                                  tCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGa   <  2:464402/150‑1 (MQ=255)
                                                                                                                                                  tCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGa   <  1:359058/150‑1 (MQ=255)
                                                                                                                                                   cACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGAt  <  2:254731/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
CCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGAT  >  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: