Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,087,4060AC52.6% 8.8 / 29.8 19W285G (TGG→GGG) yeeYputative DNA‑binding transcriptional regulator YeeY
Reads supporting (aligned to +/- strand):  ref base A (9/0);  new base C (0/10);  total (9/10)
Fisher's exact test for biased strand distribution p-value = 1.08e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-05
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGTCATCAGATTTTGCATAAACACCGGCGATCGCACTATTTGCTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATCCAGACGGGTTCAAAGAGA  >  NC_000913/2087346‑2087467
                                                             |                                                             
cGTCATCAGATTTTGCATAAACACCGGCGATCGCACTATTTGCTAAAATTTCATCCCGCCACCAGCCa                                                        >  1:2293312/1‑68 (MQ=255)
    atcaGTTTTTGCATAAACCCCGGCGATCGCATTATTTGCAAAAATTTCACCCCGCCCCCAGCCACTgg                                                    <  1:3960982/68‑1 (MQ=255)
         tttttGCATAACCCCCGGGGATCGCATTTTTTGCAAAATTTTCACCCCGCCCCCAGCCACTGGCAAGt                                               <  1:2479240/67‑1 (MQ=255)
         tttttGCATAACCCCCGGCGATCGCATTTTTTGCAAAATTTTCACCCCGCCCCCAGCCACTGGCAAGt                                               <  1:3723632/67‑1 (MQ=255)
         tttttGCAAAACCCCCGGCGATCGCATTATTTGCTAAAATTTCACCCCCCCCCCAGCCACTGGCAAGt                                               <  1:2862305/67‑1 (MQ=255)
         tttttGCAAAAACCCCGGCGAGCCCTTTTTTTGCAAAATTTTCACCCCGCCCCCAGCCACTGGCAAGt                                               <  1:2826606/67‑1 (MQ=255)
         tttttGAATAACCCCCGGCGATCGCATTTTTTGCTAAAATTTCACCCCGCCCCCAGCCACTNGCAAGt                                               <  1:2449169/67‑1 (MQ=255)
         aTTTTGCATAACCCCCGCCGACCCCACTTTTTGCAAAAATTTCCCCCCGCCCCCAGCCACTGGCAAGt                                               <  1:1751210/68‑1 (MQ=255)
         aTTTTCCATAAACCCGGGCGATCGCATTTTTTGCAAAAATTTCACCCCGCCCCCAGCCACTGGCAAGt                                               <  1:1123566/68‑1 (MQ=255)
           tttGCATAAACACCGGCGATCGCACTATTTGCTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTcc                                             >  1:4040458/1‑68 (MQ=255)
            ttccATAACCCCCGGCGATCGCATTATTTGCAAAAATTTCACCCCGCCCCCAGCCACTGGCAAGTCCt                                            <  1:3847040/65‑1 (MQ=255)
                aaaaaCCCGGGCGATCGCCCGTTTTGCTAAAATTCCACCCCGCCACCAGCCACTGGCAAGTCCTGCgg                                        <  1:3606196/66‑1 (MQ=255)
                aTAAACCCCGGCGATCGCATTTTTGGCTAAATTTCCACCCCCCCCCCAGCCACTGGCAAGTCCTGCgg                                        <  1:3971523/68‑1 (MQ=255)
                aTAAACCCCGGCGATCGCATTATTTGCTAAAATTTCACCCCGCCCCCAGCCACTGGCAAGTCCTGCgg                                        <  1:1159804/68‑1 (MQ=255)
                    acacCGGCGATCGCTTTATTTGCAAAAATTTCACCCCGCCCCCAGCCACTGGCAAGTCCTGCGGTCTg                                    <  1:2550746/68‑1 (MQ=255)
                     cgccGGGGATCGCATTATTGGAAAAATTTCCACCCCGCCCCCAGCCACTGGCAAGTCCTGCGGTCTgt                                   <  1:2738904/66‑1 (MQ=255)
                     ccccGGCGCCCGCAGTATTTCCAAAAATTCCACCCCCCCCCCAGCCACTGGCAAGTCCTGCGGTCTgt                                   <  1:2364926/66‑1 (MQ=255)
                     ccccGGCAATCGCATTTTTTGCAAAATTTTCACCCCCCCCCCAGCCACTGGCAAGTCCTGCGGTCTgt                                   <  1:1394003/66‑1 (MQ=255)
                     cacCGGCCACCGATCTCTTTGCAAAATTTCCCTCCCGCCCCCAGCCACTGGCAAGTCCTGCGGTCTgt                                   <  1:782010/68‑1 (MQ=255)
                        cGGCGATCGCATTTTTTGCAAAATTTTCACCCCGCCCCCAGCCACTGGCAAGTCCTGCGGTCTgttcg                                <  1:3889566/68‑1 (MQ=255)
                         ggCGATCGCATTTTTTGCAAAAATTTCACCCCGCCCCCAGCCACTGGCAAGTCCTGCGGTCTgttcgt                               <  1:2959238/68‑1 (MQ=255)
                         ggCGATCGCATTTTTTGCAAAAATTTCACCCCGCCCCCAGCCACTGGCAAGTCCTGCGGTCTgttcgt                               <  1:2677149/68‑1 (MQ=255)
                         ggCGATCGCATTTTTGGAAAAATTTTCACCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTgttcgt                               <  1:2481381/68‑1 (MQ=255)
                         ggCGATCGCATTATTGTCTAAAAGTTCACCCCGCCCCCAGCCACTGGCAAGTCCTGCGGTCTgttcgt                               <  1:1281453/68‑1 (MQ=255)
                           cGATCGCACTATTTGCTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTgttcgttc                             >  1:3389556/1‑68 (MQ=255)
                                gCACTATTTGCTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGcc                        >  1:3391345/1‑68 (MQ=255)
                                  aCTATTTGCTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCaa                      >  1:567484/1‑68 (MQ=255)
                                    ttttttCTAAAATTTCCCCCCCCCCCCCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAAt                     <  1:1421462/66‑1 (MQ=255)
                                    ctttttCTAAAATTTCCCCCCCCCCCCCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAAt                     <  1:341032/66‑1 (MQ=255)
                                     aTTTGCTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATcc                   >  1:782981/1‑68 (MQ=255)
                                        tGCTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATccctg                >  1:3039403/1‑65 (MQ=255)
                                          cTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATCCAGACg              >  1:1018718/1‑68 (MQ=255)
                                          cTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTACGGTCTGTTCGTTCCAGCCAATCCAGACg              >  1:2484503/1‑68 (MQ=255)
                                            aaatTTCCACCCCCCCCCCAGCCACTGGCAAGGCCTGCGGTCTGTTCGTTCCAGCCAATCCAGACggg            <  1:814701/68‑1 (MQ=255)
                                                       ccccccccCAGCGACTGGCAAGTCCGGCGGTCTCTTCGTTCCAGCCAATCCAGACGGGTTCAAagaga  <  1:1876468/68‑1 (MQ=255)
                                                       ccccccccCAGCCCCTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATCCAGACGGGTTCAAagaga  <  1:3140724/68‑1 (MQ=255)
                                                       ccccccccCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATCCAGACGGGTTCAAagaga  <  1:3107747/68‑1 (MQ=255)
                                                       ccccccccCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATCCAGACGGGTTCAAagaga  <  1:2360452/68‑1 (MQ=255)
                                                             |                                                             
CGTCATCAGATTTTGCATAAACACCGGCGATCGCACTATTTGCTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATCCAGACGGGTTCAAAGAGA  >  NC_000913/2087346‑2087467

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: