Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,087,4060AC52.6% 4.8 / 28.3 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 = 5.61e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TTTCCGTCATCAGATTTTGCATAAACACCGGCGATCGCACTATTTGCTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATCCAGACGGGTTCAAAGAGA  >  NC_000913/2087342‑2087467
                                                                 |                                                             
tttCCGTCATCAGATTTTGCATAAACACCGGCGATCGCACTATTTGCTAAAATTTCATCCCGcca                                                               >  1:1882336/1‑65 (MQ=255)
 ttCCGTCATCAGATTTTGCATAAACACCGGCGATCGCACTATTTGCTAAAATTTCATCCCGccac                                                              >  1:3325464/1‑65 (MQ=255)
     ttcatcaGATTTTGCATAAACCCGGGCGATCGCACTTTTTGCAAAAATTTCACCCCGCCCCCAGc                                                          <  1:838195/64‑1 (MQ=255)
         tcaGATTTTGCAAAACCACCGGCGATCGCACTTTTGGCAAAATTTTCACCCCGCCCCCAGCCACt                                                      <  1:1997584/65‑1 (MQ=255)
           aGATTTTGCATAACCCCCGGCGACCGCATTTTTTGCTAAAATTTCACCCCGCCCCCAGCCACTgg                                                    <  1:4310974/65‑1 (MQ=255)
                ttccATAACCNCCGGCGATCGCGTTATTTGCTAAATTTTCACCCCCCCCCCAGCCACTGGCAAGt                                               <  1:1318933/62‑1 (MQ=255)
                ttccATAACCCCCGGAGATCGCTTTATTTGCAAAATTTTCACCCCGCCCCCAGCCACTGGCAAGt                                               <  1:330872/62‑1 (MQ=255)
                ttccATAAACCCCGGCGATCGCTTTATTTGCAAAAATTTCACCCCGCCCCCAGCCACTGGCAAGt                                               <  1:1237976/62‑1 (MQ=255)
                ttGCATAAACCCCGGCGATCGCATTTTTTGCAAAAATTTCACCCCGCCCCCAGCCACTGGCAAGt                                               <  1:4403376/65‑1 (MQ=255)
                ttGCAAAACCCCCGGCGACCGCTTTATTTGCAAAAATTTCACCCCGCCCCCAGCCACTGGCAAGt                                               <  1:2441847/65‑1 (MQ=255)
                ttGCAAAAACCCCGGCGATCGCATTATTGGCTAAATTTCCACCCCGCCCCCAGCCACTGGCAAGt                                               <  1:3201582/65‑1 (MQ=255)
                tggCATAAACCCCGGCGATCGCATTTTTTGCTAAATTTTCACCCCGCCCCCAGCCACTGGCAAGt                                               <  1:1587693/63‑1 (MQ=255)
                tggCATAAACACCGGCGACCGCCTTTTTTGCTAAATTTTCACCCCGCCCCCAGCCACTGGCAAGt                                               <  1:2995026/63‑1 (MQ=255)
                    aTAAACACCGGCGATCGCACTATTTGCTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTg                                           >  1:1827541/1‑65 (MQ=255)
                    aTAAACACCGGCGATCGCACTATTTGCTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTg                                           >  1:799823/1‑65 (MQ=255)
                       aaCCCCGGCGATCGCATTATTTCCAAAAATTTCACCCCGCCCCCAGCCACTGGCAAGTCCTGCgg                                        <  1:2246767/65‑1 (MQ=255)
                            cGGCGATCGCATTTTTGGCTAAATTTTCCCCCCGCCCCCAGCCACTGGCAAGTCCTGCGGTCTgt                                   <  1:3138813/65‑1 (MQ=255)
                            cGGCCACCGCATTTTTTGCAAAAATTTCGTCCCGCCCCCAGCCACTGGCAAGTCCTGCGGTCTgt                                   <  1:1905312/65‑1 (MQ=255)
                               cGATCGCACTATTTGCTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTgttcg                                >  1:2561024/1‑65 (MQ=255)
                                gctCGCATTTTTTGCTAAAATTTCACCCCGCCCCCAGCCACTGGCAAGTCCTGCGGTCTgttcgt                               <  1:3631490/63‑1 (MQ=255)
                                      aCTATTTGCTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGc                         >  1:1792871/1‑65 (MQ=255)
                                      aCTATTTGCTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGc                         >  1:2934697/1‑65 (MQ=255)
                                        tATTTGCTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCa                       >  1:760358/1‑65 (MQ=255)
                                         aTTTGCTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCaa                      >  1:1565038/1‑65 (MQ=255)
                                          tttGCAAAAATTTCCCCCCCCCCCCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAAt                     <  1:28995/65‑1 (MQ=255)
                                          tttGCAAAAATTCCCCCCCCCCCCCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAAt                     <  1:1277233/65‑1 (MQ=255)
                                           tttcTAAAATTTCCCCCCCCCCCCCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAAt                     <  1:1439185/62‑1 (MQ=255)
                                              cTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATCCAg                 >  1:1093101/1‑65 (MQ=255)
                                               aaaatTTTCCCCCCCCCCCCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATCCAGa                <  1:1629967/64‑1 (MQ=255)
                                                         accccACCCCCAGCCACTGGCGAGTCCTGCGGTCTGTTCGTTCCAGCCAATCCAGACGGGTTCaa       <  1:1142999/63‑1 (MQ=255)
                                                          cccccccccCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATCCAGACGGGTTCaaa      <  1:671086/64‑1 (MQ=255)
                                                              gCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATCCAGACGGGTTCAAagaga  >  1:2122590/1‑65 (MQ=255)
                                                                 |                                                             
TTTCCGTCATCAGATTTTGCATAAACACCGGCGATCGCACTATTTGCTAAAATTTCATCCCGCCACCAGCCACTGGCAAGTCCTGCGGTCTGTTCGTTCCAGCCAATCCAGACGGGTTCAAAGAGA  >  NC_000913/2087342‑2087467

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