Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,158,8850AC55.9% 39.2 / 71.4 59L132F (TTA→TTC) oxyRoxidative and nitrosative stress transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (15/11);  new base C (20/13);  total (35/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATG  >  NC_000913_3_bme_pgi/4158737‑4159028
                                                                                                                                                    |                                                                                                                                               
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)
                          ccGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCg                                                                                                                     <  1:334686/151‑1 (MQ=255)
                             cTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGAtt                                                                                                                  >  1:397399/1‑151 (MQ=255)
                                      ggTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATc                                                                                                          <  4:245520/150‑1 (MQ=255)
                                           gATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTcgcg                                                                                                    <  2:41599/151‑1 (MQ=255)
                                               cccACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTgg                                                                                                >  1:76264/1‑151 (MQ=255)
                                                ccACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGt                                                                                               >  2:624136/1‑151 (MQ=255)
                                                  acaGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTg                                                                                              >  4:12817/1‑150 (MQ=255)
                                                   caGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGaa                                                                                            <  1:169240/151‑1 (MQ=255)
                                                           ccGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCg                                                                                    >  2:955674/1‑151 (MQ=255)
                                                             gTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGaa                                                                                  >  2:232729/1‑151 (MQ=255)
                                                              tACCTGCTACCGCATATTATCCGTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGaa                                                                                  <  3:498435/150‑1 (MQ=255)
                                                                   gCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGcattc                                                                            <  1:232089/151‑1 (MQ=255)
                                                                       ccGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGgc                                                                                                                                            <  1:845597/83‑1 (MQ=255)
                                                                       ccGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGgc                                                                                                                                            >  2:845597/1‑83 (MQ=255)
                                                                                atCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGcat                                                                              <  3:15753/136‑1 (MQ=255)
                                                                                atCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGcat                                                                              >  4:15753/1‑136 (MQ=255)
                                                                                    cTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGTGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCgttgtt                                                           >  2:59874/1‑151 (MQ=255)
                                                                                       tgctgcACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTtga                                                        >  1:254816/1‑151 (MQ=255)
                                                                                       tgctgcACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTtga                                                        >  2:160489/1‑151 (MQ=255)
                                                                                         ctgcACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTtgatg                                                      <  2:397399/151‑1 (MQ=255)
                                                                                          tgcACCAGACCTTACCAAAGCTGGAAATAAATCTGCATGAAGAACAAACCCACCAGTTACTGGCGCTACGGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGATGAAAGAGAGCGAAGCATTCATGGAAGTGCCGGTGTTtgatg                                                      >  4:134926/1‑150 (MQ=255)
                                                                                              ccAGACCTTTGCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCa                                                 >  1:650771/1‑151 (MQ=255)
                                                                                              ccAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGcc                                                  >  3:313238/1‑150 (MQ=255)
                                                                                               cAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTg                                                                                              >  1:223574/1‑105 (MQ=255)
                                                                                               cAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTg                                                                                              <  2:223574/105‑1 (MQ=255)
                                                                                                  aCCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGt                                             >  2:204023/1‑151 (MQ=255)
                                                                                                    cTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCgg                                                                                                                             <  2:578431/69‑1 (MQ=255)
                                                                                                    cTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCgg                                                                                                                             >  1:578431/1‑69 (MQ=255)
                                                                                                       tCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTg                                        <  1:315606/151‑1 (MQ=255)
                                                                                                        ccAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTg                                        <  3:580415/150‑1 (MQ=255)
                                                                                                         cAAAGCTGGAAATGTATCTGCATGAAGCACCGACCCACCAGTTACTGGCGCAACTGGACAGCGGCACACTCGATTGCGTGATCCTCGCGCTGGGGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGc                                      >  2:721181/1‑151 (MQ=255)
                                                                                                          aacgCTGGAAATGTATCTGAATGAAGCAGAGACCCACAAGTTCCTGGCGCAACTAGACAGCGGCAAACTAGATTACGTGATCCTCGCGCTGGTGGAAGAGAGCGAGGCATGCATTGTAGTGCCGTTGTTTGATGAGCTAATGTTGCTGGc                                      >  4:264613/4‑150 (MQ=255)
                                                                                                             gCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGctatc                                  >  2:450207/1‑151 (MQ=255)
                                                                                                               tGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGctatcta                                >  2:800859/1‑151 (MQ=255)
                                                                                                               tGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGctatcta                                <  1:888489/151‑1 (MQ=255)
                                                                                                                ggAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGctatcta                                >  4:91522/1‑150 (MQ=255)
                                                                                                                     tGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGa                          >  1:392933/1‑151 (MQ=255)
                                                                                                                           tGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCAccc                    <  1:624136/151‑1 (MQ=255)
                                                                                                                              aTGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGt                  >  3:573899/1‑150 (MQ=255)
                                                                                                                              aTGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGt                  >  3:323969/1‑150 (MQ=255)
                                                                                                                              aTGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTg                 >  2:854380/1‑151 (MQ=255)
                                                                                                                               tGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTtga                                                        >  3:189268/1‑111 (MQ=255)
                                                                                                                               tGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTtga                                                        <  4:189268/111‑1 (MQ=255)
                                                                                                                                gAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTggg               <  2:254816/151‑1 (MQ=255)
                                                                                                                                     acaGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAAc          <  2:274824/151‑1 (MQ=255)
                                                                                                                                         aCCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACcgc       >  3:592087/1‑150 (MQ=255)
                                                                                                                                            caccaGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGa            >  2:252316/1‑142 (MQ=255)
                                                                                                                                            caccaGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGa            <  1:252316/142‑1 (MQ=255)
                                                                                                                                              ccaGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATg  >  3:410470/1‑150 (MQ=255)
                                                                                                                                              ccaGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGtt                                            >  3:241011/1‑108 (MQ=255)
                                                                                                                                              ccaGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGtt                                            <  4:241011/108‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATG  >  NC_000913_3_bme_pgi/4158737‑4159028

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