Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,158,7790TG67.2% 81.6 / 54.3 58L97W (TTG→TGG) oxyRoxidative and nitrosative stress transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (10/9);  new base G (18/21);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 7.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATC  >  NC_000913_3_bme_pgi/4158641‑4158924
                                                                                                                                          |                                                                                                                                                 
ccAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAg                                                                                                                                        <  3:30338/150‑1 (MQ=255)
 cAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGtt                                                                                                                                      >  1:941983/1‑151 (MQ=255)
  aGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTg                                                                                                                                     <  1:63128/151‑1 (MQ=255)
   gCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGTAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTg                                                                                                                                     <  4:32528/150‑1 (MQ=255)
   gCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTgg                                                                                                                                    <  1:530724/151‑1 (MQ=255)
   gCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTgg                                                                                                                                    <  1:436106/151‑1 (MQ=255)
      gTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGAcc                                                                                                                                 >  1:330251/1‑151 (MQ=255)
       tAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCg                                                                                                                                >  2:297569/1‑151 (MQ=255)
            tgttgtTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCt                                                                                                                           >  1:18345/1‑151 (MQ=255)
                 ttCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTAc                                                                                                                      >  2:511345/1‑151 (MQ=255)
                       cAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCata                                                                                                                >  2:315606/1‑151 (MQ=255)
                        aGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCatat                                                                                                               >  2:658283/1‑151 (MQ=255)
                        aGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCatat                                                                                                               >  2:841865/1‑151 (MQ=255)
                         ggCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCcac                                                                                                                                          >  1:284127/1‑123 (MQ=255)
                         ggCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCcac                                                                                                                                          <  2:284127/123‑1 (MQ=255)
                               aaTGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATccc                                                                                                         >  3:498618/1‑150 (MQ=255)
                                  gctgctGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTAtg                                                                                                     <  2:516365/151‑1 (MQ=255)
                                   ctgctgGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTAtgc                                                                                                    >  2:440622/1‑151 (MQ=255)
                                      ctgGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGGAAGTCCTTAAAGATATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTAtgctgc                                                                                                 >  2:810392/1‑151 (MQ=255)
                                        ggtggATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCAc                                                                                               >  2:334686/1‑151 (MQ=255)
                                         gtggATCAGGCGCGTATCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCAc                                                                                               <  3:182375/150‑1 (MQ=255)
                                                        aCCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAg                                                                                <  3:138268/150‑1 (MQ=255)
                                                           gTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTgg                                                                            <  1:297569/151‑1 (MQ=255)
                                                            tgctgcGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTgg                                                                            >  3:218233/1‑150 (MQ=255)
                                                             gctgcGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGaa                                                                          <  1:511345/151‑1 (MQ=255)
                                                               tgcGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAAt                                                                        <  2:798437/151‑1 (MQ=255)
                                                                gcGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATg                                                                       <  2:18345/151‑1 (MQ=255)
                                                                gcGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATg                                                                       <  1:657262/151‑1 (MQ=255)
                                                                gcGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATg                                                                       >  1:41599/1‑151 (MQ=255)
                                                                          aaaGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCAtn                                                             <  2:941983/151‑2 (MQ=255)
                                                                           aaGTCCTTAAAGAGCTGGCAAGCCAGCAGGGCGAGACGATGTTCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGTATCGTACCTGTTACCGCATATTATTTCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCagg                                                             <  3:146599/150‑3 (MQ=255)
                                                                            aGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCatat                                                                                                               >  1:603147/1‑99 (MQ=255)
                                                                            aGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCatat                                                                                                               <  2:603147/99‑1 (MQ=255)
                                                                             gTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAg                                                          >  2:232089/1‑151 (MQ=255)
                                                                             gTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAg                                                          >  2:169240/1‑151 (MQ=255)
                                                                                    aaGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCatat                                                                                                               >  1:181407/1‑91 (MQ=255)
                                                                                    aaGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCatat                                                                                                               <  2:181407/91‑1 (MQ=255)
                                                                                      gagaTGGAAAGCCAGCAGGGCGAGCCGATGTCAGGACCGCTGCAAATTGGTTGGATTACCACAGTTCGACCGTCTCTGATTCCGCCTATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGAccc                                                 <  1:810392/151‑1 (MQ=255)
                                                                                          tGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACccacca                                             <  1:841865/151‑1 (MQ=255)
                                                                                          tGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACccacca                                             <  1:658283/151‑1 (MQ=255)
                                                                                            gCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAg                                            >  4:580415/1‑150 (MQ=255)
                                                                                               aGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCAcc                                                                                              <  3:154430/97‑1 (MQ=255)
                                                                                               aGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCAcc                                                                                              >  4:154430/1‑97 (MQ=255)
                                                                                                gCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCt                                       <  1:218081/151‑1 (MQ=255)
                                                                                                 ccAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTg                                      <  1:440622/151‑1 (MQ=255)
                                                                                                  cagcagGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTgg                                     <  1:600048/151‑1 (MQ=255)
                                                                                                    gcagGGCGAGACGATGTCCGGACCGCTGCACATTGGTCTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGACGCACAGACCCACCAGTTACTGgc                                    >  4:506143/1‑150 (MQ=255)
                                                                                                         gCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACt                              >  2:768923/1‑151 (MQ=255)
                                                                                                             gaCGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGa                           <  4:498618/150‑1 (MQ=255)
                                                                                                                   gTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGc                    >  2:888489/1‑151 (MQ=255)
                                                                                                                     ccGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCCAAGCTGGAAATGTATCCCCACGACAACCAGACCCACCCGTTCTAGGCGCGACGGGTCAgggaaaa                  >  2:82218/1‑145 (MQ=255)
                                                                                                                          ccGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCa                                                                                   >  1:649893/1‑81 (MQ=255)
                                                                                                                          ccGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCa                                                                                   <  2:649893/81‑1 (MQ=255)
                                                                                                                          ccGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCg             <  1:334686/151‑1 (MQ=255)
                                                                                                                             cTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGAtt          >  1:397399/1‑151 (MQ=255)
                                                                                                                                   aTTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAg                                                                                >  3:51106/1‑75 (MQ=255)
                                                                                                                                   aTTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAg                                                                                <  4:51106/75‑1 (MQ=255)
                                                                                                                                      ggTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATc  <  4:245520/150‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                 
CCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATC  >  NC_000913_3_bme_pgi/4158641‑4158924

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑