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.

CGTCATCAGATTTTGCATAAACACCGGCGATCGCACTATTTGCTAAAATTT‑CATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATCCAGACGGGTTCAAAGAGA  >  NC_000913/2087346‑2087467
                                                             |                                                             
cGTCATCAGATTTTGCATAAACACCGGCGATCGCACTATTTGCTAAAATTT‑CATCCCGCCACCAGCCa                                                        >  1:2293312/1‑68 (MQ=255)
    atcaGTTTTTGCATAAACCCCGGCGATCGCATTATTTGCAAAAATTT‑CACCCCGCCCCCAGCCACTgg                                                    <  1:3960982/68‑1 (MQ=255)
         tttttGCATAACCCCCGGGGATCGCATTTTTTGCAAAATTTT‑CACCCCGCCCCCAGCCACTGGCAAGt                                               <  1:2479240/67‑1 (MQ=255)
         tttttGCATAACCCCCGGCGATCGCATTTTTTGCAAAATTTT‑CACCCCGCCCCCAGCCACTGGCAAGt                                               <  1:3723632/67‑1 (MQ=255)
         tttttGCAAAACCCCCGGCGATCGCATTATTTGCTAAAATTT‑CACCCCCCCCCCAGCCACTGGCAAGt                                               <  1:2862305/67‑1 (MQ=255)
         tttttGCAAAAACCCCGGCGAGCCCTTTTTTTGCAAAATTTT‑CACCCCGCCCCCAGCCACTGGCAAGt                                               <  1:2826606/67‑1 (MQ=255)
         tttttGAATAACCCCCGGCGATCGCATTTTTTGCTAAAATTT‑CACCCCGCCCCCAGCCACTNGCAAGt                                               <  1:2449169/67‑1 (MQ=255)
         aTTTTGCATAACCCCCGCCGACCCCACTTTTTGCAAAAATTT‑CCCCCCGCCCCCAGCCACTGGCAAGt                                               <  1:1751210/68‑1 (MQ=255)
         aTTTTCCATAAACCCGGGCGATCGCATTTTTTGCAAAAATTT‑CACCCCGCCCCCAGCCACTGGCAAGt                                               <  1:1123566/68‑1 (MQ=255)
           tttGCATAAACACCGGCGATCGCACTATTTGCTAAAATTT‑CATCCCGCCACCAGCCACTGGCAAGTcc                                             >  1:4040458/1‑68 (MQ=255)
            ttccATAACCCCCGGCGATCGCATTATTTGCAAAAATTT‑CACCCCGCCCCCAGCCACTGGCAAGTCCt                                            <  1:3847040/65‑1 (MQ=255)
                aaaaaCCCGGGCGATCGCCCGTTTTGCTAAAATTC‑CACCCCGCCACCAGCCACTGGCAAGTCCTGCgg                                        <  1:3606196/66‑1 (MQ=255)
                aTAAACCCCGGCGATCGCATTTTTGGCTAAATTTC‑CACCCCCCCCCCAGCCACTGGCAAGTCCTGCgg                                        <  1:3971523/68‑1 (MQ=255)
                aTAAACCCCGGCGATCGCATTATTTGCTAAAATTT‑CACCCCGCCCCCAGCCACTGGCAAGTCCTGCgg                                        <  1:1159804/68‑1 (MQ=255)
                    acacCGGCGATCGCTTTATTTGCAAAAATTT‑CACCCCGCCCCCAGCCACTGGCAAGTCCTGCGGTCTg                                    <  1:2550746/68‑1 (MQ=255)
                     cgccGGGGATCGCATTATTGGAAAAATTTC‑CACCCCGCCCCCAGCCACTGGCAAGTCCTGCGGTCTgt                                   <  1:2738904/66‑1 (MQ=255)
                     ccccGGCGCCCGCAGTATTTCCAAAAATTC‑CACCCCCCCCCCAGCCACTGGCAAGTCCTGCGGTCTgt                                   <  1:2364926/66‑1 (MQ=255)
                     ccccGGCAATCGCATTTTTTGCAAAATTTT‑CACCCCCCCCCCAGCCACTGGCAAGTCCTGCGGTCTgt                                   <  1:1394003/66‑1 (MQ=255)
                     cacCGGCCACCGATCTCTTTGCAAAATTTC‑CCTCCCGCCCCCAGCCACTGGCAAGTCCTGCGGTCTgt                                   <  1:782010/68‑1 (MQ=255)
                        cGGCGATCGCATTTTTTGCAAAATTTT‑CACCCCGCCCCCAGCCACTGGCAAGTCCTGCGGTCTgttcg                                <  1:3889566/68‑1 (MQ=255)
                         ggCGATCGCATTTTTTGCAAAAATTT‑CACCCCGCCCCCAGCCACTGGCAAGTCCTGCGGTCTgttcgt                               <  1:2959238/68‑1 (MQ=255)
                         ggCGATCGCATTTTTTGCAAAAATTT‑CACCCCGCCCCCAGCCACTGGCAAGTCCTGCGGTCTgttcgt                               <  1:2677149/68‑1 (MQ=255)
                         ggCGATCGCATTTTTGGAAAAATTTT‑CACCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTgttcgt                               <  1:2481381/68‑1 (MQ=255)
                         ggCGATCGCATTATTGTCTAAAAGTT‑CACCCCGCCCCCAGCCACTGGCAAGTCCTGCGGTCTgttcgt                               <  1:1281453/68‑1 (MQ=255)
                           cGATCGCACTATTTGCTAAAATTT‑CATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTgttcgttc                             >  1:3389556/1‑68 (MQ=255)
                                gCACTATTTGCTAAAATTT‑CATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGcc                        >  1:3391345/1‑68 (MQ=255)
                                  aCTATTTGCTAAAATTT‑CATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCaa                      >  1:567484/1‑68 (MQ=255)
                                    ttttttCTAAAATTTCCCCCCCCCCCCCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAAt                     <  1:1421462/66‑1 (MQ=255)
                                    ctttttCTAAAATTTCCCCCCCCCCCCCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAAt                     <  1:341032/66‑1 (MQ=255)
                                     aTTTGCTAAAATTT‑CATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATcc                   >  1:782981/1‑68 (MQ=255)
                                        tGCTAAAATTT‑CATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATccctg                >  1:3039403/1‑65 (MQ=255)
                                          cTAAAATTT‑CATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATCCAGACg              >  1:1018718/1‑68 (MQ=255)
                                          cTAAAATTT‑CATCCCGCCACCAGCCACTGGCAAGTCCTACGGTCTGTTCGTTCCAGCCAATCCAGACg              >  1:2484503/1‑68 (MQ=255)
                                            aaatTTC‑CACCCCCCCCCCAGCCACTGGCAAGGCCTGCGGTCTGTTCGTTCCAGCCAATCCAGACggg            <  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)
                                                             |                                                             
CGTCATCAGATTTTGCATAAACACCGGCGATCGCACTATTTGCTAAAATTT‑CATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATCCAGACGGGTTCAAAGAGA  >  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: ‑