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

CACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTC‑GATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCAT  >  NC_000913/2326314‑2326579
                                                                                                                                        |                                                                                                                                  
cACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATg                                                                                                                                 <  1:270382/140‑1 (MQ=255)
     ttAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATg                                                                                                                                 >  2:529494/1‑135 (MQ=255)
     ttAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATg                                                                                                                                 <  1:529494/135‑1 (MQ=255)
     ttAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGaaaa                                                                                                                            >  1:5645/1‑140 (MQ=255)
             cGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGtt                                                                                                                    <  2:610020/140‑1 (MQ=255)
             cGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGtt                                                                                                                    <  1:311618/140‑1 (MQ=255)
             cGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAggt                                                                                                                    >  2:311618/1‑138 (MQ=255)
               ggCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGgttt                                                                                                                  >  2:306412/1‑139 (MQ=255)
                  tGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGcc                                                                                                               <  1:306412/140‑1 (MQ=255)
                   ggCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGccc                                                                                                              <  1:287419/140‑1 (MQ=255)
                    gCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGcccc                                                                                                             <  1:261085/140‑1 (MQ=255)
                             gtgtGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCtacttact                                                                                                    <  2:413303/140‑1 (MQ=255)
                             gtgtGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCtacttact                                                                                                    <  2:503221/140‑1 (MQ=255)
                                        cACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTC‑GATGCAAAAg                                                                                        >  1:442962/1‑140 (MQ=255)
                                               gacagaGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTC‑GATGCAAAAGCACGctc                                                                                 <  1:551558/135‑1 (MQ=255)
                                               aaTAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATGAGTTTAGCCCCCCACCTACTC‑GATGCCAAAGCACGctc                                                                                 >  2:41397/1‑140 (MQ=255)
                                                         gtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTC‑GATGCAAAAGCACGCTCTGGTTATCGt                                                                       >  1:663802/1‑140 (MQ=255)
                                                         gtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATTAGGTTAGTCCCCCACTTACTC‑GATGCAAAAAAACGCCCTGGTTATCGt                                                                       >  2:496588/1‑140 (MQ=255)
                                                           ggTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTC‑GATGCAAAAGCACGCTCTGGTTATCGTCt                                                                     <  1:496588/140‑1 (MQ=255)
                                                              tCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGTTTAGCCCCCTACTTACTC‑GATGCAAAAGCACGCTCTGGTTATCGTCTTgg                                                                  >  1:81643/1‑140 (MQ=255)
                                                               cGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTC‑GATGCAAAAGCACGCTCTGGTTATCGTCTTGga                                                                 <  2:5645/140‑1 (MQ=255)
                                                               cGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATGAGGTTTGCCCCCCACTTACTC‑CATGCAAAAACACGCCCTGGGTATCGTCCTgtc                                                                 >  2:551558/1‑138 (MQ=255)
                                                                  gTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAAATGGTTTATGCCCCCCCTTATTCGGATTC‑AAATAACGCGCTTGTTATTGTCTTGgggaa                                                              >  2:481327/1‑139 (MQ=255)
                                                                              tgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTC‑GATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTAt                                                  >  1:229422/1‑140 (MQ=255)
                                                                                         gccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTC‑GATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCt                                       <  2:442962/140‑1 (MQ=255)
                                                                                              ggCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTC‑GATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGcgcg                                  <  1:627265/140‑1 (MQ=255)
                                                                                               gCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTAGCCCCCTACTTACTC‑GATGCAAAAGCACGCTCTGGTTATCGTCTTTGGGACGGACAGGGTTTTGACGTAATCCTGCGCGa                                 >  2:605886/1‑140 (MQ=255)
                                                                                                   ttCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTC‑GATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGc                             <  2:663802/140‑1 (MQ=255)
                                                                                                              cAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTC‑GATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTgcgc                  <  1:481327/140‑1 (MQ=255)
                                                                                                                   gcAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTC‑GATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCAccc             <  1:41397/140‑1 (MQ=255)
                                                                                                                    cAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTC‑GATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCa            <  2:502147/140‑1 (MQ=255)
                                                                                                                              tgGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTC‑GATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCat  <  2:378106/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                  
CACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTC‑GATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCAT  >  NC_000913/2326314‑2326579

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: ‑