Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,326,4500TG35.7% 25.9 / 21.3 28G114G (GGT→GGG) atoBacetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/18);  new base G (10/0);  total (10/18)
Fisher's exact test for biased strand distribution p-value = 7.62e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTA  >  NC_000913/2326327‑2326575
                                                                                                                           |                                                                                                                             
cGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGtt                                                                                                               <  1:79951/140‑1 (MQ=255)
cGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATGAGtt                                                                                                               >  2:391231/1‑140 (MQ=255)
cGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGtt                                                                                                               >  2:81185/1‑140 (MQ=255)
  ggCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTa                                                                                                             <  2:32358/140‑1 (MQ=255)
  ggCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTa                                                                                                             <  1:327050/140‑1 (MQ=255)
  ggCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATTAGtttt                                                                                                             >  2:147488/1‑139 (MQ=255)
  ggCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTa                                                                                                             >  1:32358/1‑140 (MQ=255)
     tGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGcc                                                                                                          <  1:74779/140‑1 (MQ=255)
      ggCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGccc                                                                                                         <  1:147488/140‑1 (MQ=255)
       gCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGcccc                                                                                                        <  2:168837/140‑1 (MQ=255)
       gCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGcccc                                                                                                        <  1:391231/140‑1 (MQ=255)
                   tgCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGa                                                                                            >  1:55110/1‑140 (MQ=255)
                          tCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATGAGGTTAACCCCCCACTTACTCGATGCaaaa                                                                                     >  2:425050/1‑140 (MQ=255)
                           cACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAg                                                                                    <  1:294050/140‑1 (MQ=255)
                             cGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCa                                                                                  <  2:653487/140‑1 (MQ=255)
                             cGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGTTGGAAAAAATGAGTTTTGCCCCCCTCCTCCTCGATGCaaaaaaa                                                                                  >  2:291053/1‑137 (MQ=255)
                                  aaTAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATGAGTTTAGCCCCCCACTTACTCGATGCAAAAACACGccc                                                                             >  2:529839/1‑138 (MQ=255)
                                   aTAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGctct                                                                            >  1:101189/1‑140 (MQ=255)
                                            gtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTATCCCCCCACTTACTCGATGCAAAAACACGCCCCGGTTATCgg                                                                   >  2:269600/1‑139 (MQ=255)
                                            gtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGTTTAGCCCCCCACTTACTCGATGCAAAAACACGCCCCGGGTATTGt                                                                   >  2:197633/1‑140 (MQ=255)
                                              ggTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCt                                                                 <  1:529839/140‑1 (MQ=255)
                                              ggTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCt                                                                 >  1:200566/1‑140 (MQ=255)
                                                                gtgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTa                                               <  2:101189/140‑1 (MQ=255)
                                                                 tgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTAt                                              >  1:401599/1‑140 (MQ=255)
                                                                        tGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTaa                                       <  2:274841/140‑1 (MQ=255)
                                                                                 ggCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGcgcg                              <  2:401599/140‑1 (MQ=255)
                                                                                       tCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGGTTAGCCCCCCACTTTCTCGATGCAAAAACACGCTCTGGGTATCGTCTTGGAGACGGACAGGGTTATGACGTAATCCTGCGCGGTGGcc                        >  2:18706/1‑140 (MQ=255)
                                                                                                   ggCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGTTTAGCCCCCCACTTTCTCGATGCAAAAGCACGCCCTGGTTATTGTCTTGGAGACGGACAGGGTTTTGACGTAATCCTGCGCGGTGGGCTGATGTGCGccc            >  2:308966/1‑139 (MQ=255)
                                                                                                     cgcAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCAcc          <  1:291053/140‑1 (MQ=255)
                                                                                                       cAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGTGCCACCCa        <  1:197633/140‑1 (MQ=255)
                                                                                                         gagCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATg      <  2:200566/140‑1 (MQ=255)
                                                                                                            cATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGtt   <  2:55110/140‑1 (MQ=255)
                                                                                                             aTTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTa  <  1:269600/140‑1 (MQ=255)
                                                                                                                           |                                                                                                                             
CGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTA  >  NC_000913/2326327‑2326575

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑