Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,158,8850AC63.6% 33.8 / 33.8 33L132F (TTA→TTC) oxyRoxidative and nitrosative stress transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (7/5);  new base C (10/11);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAAT  >  NC_000913_3_bme_pgi/4158753‑4159027
                                                                                                                                    |                                                                                                                                              
gATGTCCGGACCGCTGAACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAg                                                                                                                               >  2:356936/1‑150 (MQ=255)
    tCCGGACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCCCCAGTTCCTGGCGCAACTGGACAGCGGc                                                                                                                           <  2:462692/150‑1 (MQ=255)
        gACCGCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAAc                                                                                                                       >  3:93686/1‑150 (MQ=255)
            gCTGCACATTGGTTGGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGa                                                                                                                   <  1:351366/150‑1 (MQ=255)
                                ccACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTgg                                                                                               >  2:23375/1‑150 (MQ=255)
                                            cGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCg                                                                                   >  1:68842/1‑150 (MQ=255)
                                                   gCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGa                                                                                                           <  2:43835/119‑1 (MQ=255)
                                                   gCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGa                                                                                                           >  1:43835/1‑119 (MQ=255)
                                                     tACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGcattca                                                                          >  3:18708/1‑150 (MQ=255)
                                                      aCCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGcattcat                                                                         <  1:290640/150‑1 (MQ=255)
                                                       ccGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACt                                                                                                                      <  1:273077/104‑1 (MQ=255)
                                                       ccGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACt                                                                                                                      >  2:273077/1‑104 (MQ=255)
                                                               tatCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGcc                                                                <  1:356936/150‑1 (MQ=255)
                                                                          tgcACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTtgatg                                                     <  1:23375/150‑1 (MQ=255)
                                                                              ccAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCgttgtt                                                          >  1:422215/1‑141 (MQ=255)
                                                                              ccAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCgttgtt                                                          <  2:422215/141‑1 (MQ=255)
                                                                               cAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCa                                                >  3:260248/1‑150 (MQ=255)
                                                                                 gACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGAtt                                                                                                                 <  4:134246/83‑1 (MQ=255)
                                                                                 gACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGAtt                                                                                                                 >  3:134246/1‑83 (MQ=255)
                                                                                    cTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCgttgtt                                                          >  4:106961/1‑135 (MQ=255)
                                                                                    cTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCgttgtt                                                          <  3:106961/135‑1 (MQ=255)
                                                                                     tttCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTg                                          <  1:339111/150‑1 (MQ=255)
                                                                                                          cTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCAc                     <  2:418844/150‑1 (MQ=255)
                                                                                                          cTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCAc                     >  1:418844/1‑150 (MQ=255)
                                                                                                             cATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCg                  >  3:268116/1‑150 (MQ=255)
                                                                                                               tGAAGCACAGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTg                <  3:13922/150‑1 (MQ=255)
                                                                                                                   gCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCg            >  4:1284/1‑150 (MQ=255)
                                                                                                                    cacaGACCCACCAGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGa           <  4:93686/150‑1 (MQ=255)
                                                                                                                           ccaccaGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTg                                                                  >  1:21475/1‑88 (MQ=255)
                                                                                                                           ccaccaGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTg                                                                  <  2:21475/88‑1 (MQ=255)
                                                                                                                             accaGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCgtt                                                             >  2:423454/1‑91 (MQ=255)
                                                                                                                             accaGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCgtt                                                             <  1:423454/91‑1 (MQ=255)
                                                                                                                             accaGTTCCTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAAt  >  1:36681/1‑150 (MQ=255)
                                                                                                                                    |                                                                                                                                              
GATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAAT  >  NC_000913_3_bme_pgi/4158753‑4159027

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