Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,158,7790TG64.5% 26.4 / 31.9 31L97W (TTG→TGG) oxyRoxidative and nitrosative stress transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (6/5);  new base G (11/9);  total (17/14)
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.76e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGA  >  NC_000913_3_bme_pgi/4158630‑4158914
                                                                                                                                                     |                                                                                                                                       
gCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGttt                                                                                                                                         <  2:196761/150‑1 (MQ=255)
  tGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGa                                                                                                                                       >  1:372951/1‑150 (MQ=255)
        ggACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCca                                                                                                                                 >  1:372062/1‑150 (MQ=255)
                 gTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGAc                                                                                                                        >  1:464516/1‑150 (MQ=255)
                      gtgtTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTAc                                                                                                                   <  2:372951/150‑1 (MQ=255)
                           gttCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCt                                                                                                              >  1:341005/1‑150 (MQ=255)
                                   aGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCata                                                                                                      <  2:372062/150‑1 (MQ=255)
                                            tgctgcTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTa                                                                                             <  2:464516/150‑1 (MQ=255)
                                               tgctgGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGtt                                                                                                                            <  2:97718/116‑1 (MQ=255)
                                               tgctgGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGtt                                                                                                                            >  1:97718/1‑116 (MQ=255)
                                               tgctgGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTAtgc                                                                                          >  2:351366/1‑150 (MQ=255)
                                                      ggATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCa                                                                                                        >  2:234160/1‑129 (MQ=255)
                                                      ggATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCa                                                                                                        <  1:234160/129‑1 (MQ=255)
                                                                        gctgcGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGa                                                                 >  2:290640/1‑150 (MQ=255)
                                                                          tgcGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGaaa                                                               <  2:341005/150‑1 (MQ=255)
                                                                             gTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGaa                                                                >  4:141884/1‑146 (MQ=255)
                                                                             gTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGaa                                                                <  3:141884/146‑1 (MQ=255)
                                                                                     aaaGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCAt                                                    >  1:462692/1‑150 (MQ=255)
                                                                                              aaaGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCAtattat                                                                                                  <  1:384078/95‑1 (MQ=255)
                                                                                              aaaGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCAtattat                                                                                                  >  2:384078/1‑95 (MQ=255)
                                                                                                        cAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTc                                                                           >  4:215017/1‑108 (MQ=255)
                                                                                                        cAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTc                                                                           <  3:215017/108‑1 (MQ=255)
                                                                                                                 agGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCAtatt                                                                                                    <  1:358778/74‑1 (MQ=255)
                                                                                                                 agGGCGAGACGATGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCAtatt                                                                                                    >  2:358778/1‑74 (MQ=255)
                                                                                                                           gATGTCCGGACCGCTGAACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAg              >  2:356936/1‑150 (MQ=255)
                                                                                                                            aTGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTcc                                                                          >  2:458707/1‑89 (MQ=255)
                                                                                                                            aTGTCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTcc                                                                          <  1:458707/89‑1 (MQ=255)
                                                                                                                               tCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCCCCAGTTCCTGGCGCAACTGGACAGCGGc          <  2:462692/150‑1 (MQ=255)
                                                                                                                                   gACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAAc      >  3:93686/1‑150 (MQ=255)
                                                                                                                                       gCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGa  <  1:351366/150‑1 (MQ=255)
                                                                                                                                        cTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAAt                                                              <  2:40652/89‑1 (MQ=255)
                                                                                                                                        cTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAAt                                                              >  1:40652/1‑89 (MQ=255)
                                                                                                                                                     |                                                                                                                                       
GCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGA  >  NC_000913_3_bme_pgi/4158630‑4158914

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: ‑