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

TTTCCGTCATCAGATTTTGCATAAACACCGGCGATCGCACTATTTGCTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATCCAGACGGGTTCAAAGAGAG  >  NC_000913/2087342‑2087468
                                                                  |                                                              
tttCCGTCATCAGATTTTGCATAAACACCGGCGATCGCACTATTTGCTAAAATTTCATCCCGcca                                                                >  1:2869132/1‑65 (MQ=255)
 ttCCGTCATCAGATTTTGCATAAACACCGGCGATCGCACTATTTGCTAAAATTTCATCCCGccac                                                               >  1:1742209/1‑65 (MQ=255)
           acaTTTTGCATAACCCCCGGCGCCCGCATTTTTTGCAAAAATTTCTTCCCCCCCCCAGCCACTgg                                                     <  1:2453655/63‑1 (MQ=255)
            atatttGCAAAAACCACCGGCGATCCCCCTATTTGAAAAAATTCCACCCCGCCCCCAGCCACTgg                                                     <  1:4426894/62‑1 (MQ=255)
                ttGCATAACCCCCGGCGATCCCATTTTTTGCAAAATTTTCACCCCGCCCCCAGCCACTGGCAAGt                                                <  1:4484916/65‑1 (MQ=255)
                ttGCAAAACCCCCGGCGATCGCATTTTTTGCAAAATTTTCACCCCGCCCCCAGCCACTGGCAAGt                                                <  1:1070597/65‑1 (MQ=255)
                 tGCAAAAACCCCGGCCATCCCCTTTTTTTCTAAAATTTCCCCCCCGCCCCCAGCCACTGGCAAGt                                                <  1:3370723/65‑1 (MQ=255)
                  gCATAAACCCCGGCGATGGCGTTTTTTGCTAAAATTCCCTCCCGCCCCCAGCCACTGGCAAGTcc                                              <  1:640468/65‑1 (MQ=255)
                  gCAAAAACACCGGCGATCGCATTTTTTGCAAAATTTTCACCCCCCCCCCAGCCACTGGCAAGTcc                                              <  1:4034654/65‑1 (MQ=255)
                  ccATAAACACCGGCGATCGCATTTTTTGCAAAAATTTCACCCCGCCCCCAGCCACTGGCAAGTcc                                              <  1:662191/64‑1 (MQ=255)
                       aaCCCCGGCGATCGCCCTATTTGCAAAAATTTCACCCCGCCCCCAGCCACTGGTAAGTCCTGCgg                                         <  1:2463794/65‑1 (MQ=255)
                            cgggCGATCGCCTTTTTTGCTAAATTTTCACCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTg                                     <  1:4395477/63‑1 (MQ=255)
                             cGGCGATCGCCTTTTTTGCTAAATTTTCACCCCGCCCCCAGCCACTGGCAAGTCCTGCGGTCTgt                                    <  1:3971563/65‑1 (MQ=255)
                             cGGCGATCGCACTATTTGCTAAAATTTCTCCCCGCCCCCAGCCACTGGCAAGTCCTGCGGTCTgt                                    <  1:390225/65‑1 (MQ=255)
                               gCGATCGCACTATTTGCTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTgttc                                  >  1:2197555/1‑65 (MQ=255)
                                 gaccGCACTTTTGTCTAAATTTTCATCCCGCCCCCAGCCACTGGCAAGTCCTGCGGTCTgttcgt                                <  1:1845329/62‑1 (MQ=255)
                                 gATCCCACTATTTGCAAAATTTTCCCCCCGCCCCCAGCCACTGGCAAGTCCTGCGGTCTgttcgt                                <  1:4388237/65‑1 (MQ=255)
                                   tCGCACTATTTGCTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTgttcgttc                              >  1:587963/1‑65 (MQ=255)
                                     gCACTATTTGCTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCa                            >  1:1948216/1‑65 (MQ=255)
                                       aCTATTTGCTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGc                          >  1:252984/1‑65 (MQ=255)
                                         tATTTGCTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCa                        >  1:1319103/1‑65 (MQ=255)
                                           tttGCAAAATTTTCCCCCCCCCCCCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAAt                      <  1:4568966/65‑1 (MQ=255)
                                            ttGCTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATc                     >  1:4361370/1‑65 (MQ=255)
                                               cTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATCCAg                  >  1:326562/1‑65 (MQ=255)
                                                    aTTCCACCCCCCCCCCAGCCACGGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATCCAGACggg             <  1:362050/65‑1 (MQ=255)
                                                      ttCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATCCAGACGGGtt           >  1:2358/1‑65 (MQ=255)
                                                           ccccGCCCCCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATCCAGACGGGTTCaaa       <  1:806500/64‑1 (MQ=255)
                                                                ccaccaGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATCCAGACGGGTTCAAagagag  >  1:38739/1‑65 (MQ=255)
                                                                  |                                                              
TTTCCGTCATCAGATTTTGCATAAACACCGGCGATCGCACTATTTGCTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATCCAGACGGGTTCAAAGAGAG  >  NC_000913/2087342‑2087468

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