Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,326,4500TG52.4% 18.7 / 25.8 21G114G (GGT→GGGatoBacetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/9);  new base G (10/1);  total (11/10)
Fisher's exact test for biased strand distribution p-value = 3.46e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.82e-05
Rejected as consensus: Frequency below/above cutoff threshold.

ATTGAACGTGCAAAAATCGATTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTT  >  NC_000913/2326220‑2326695
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
aTTGAACGTGCAAAAATCGATTCACAACACGTTGATGAAGCGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGGGttttgc                                                                                                                                                                                                                                    >  2:322990/1‑248 (MQ=255)
     aCGTGCAAAAATCGATTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCt                                                                                                                                                                                                                               <  2:324546/251‑1 (MQ=255)
              aaTCGATTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGTTTAACCCCCTCCTTATTCg                                                                                                                                                                                                                      >  1:38778/1‑251 (MQ=255)
                             cGTTGAGGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGct                                                                                                                                                                                                       <  1:206124/251‑1 (MQ=255)
                                            tATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCAGGCGCGACAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATGAGATTCGCCCCCTACTTACTCGATGCaaaaacaacccctggttatagtgttg                                                                                                                                                                                        >  2:309341/1‑229 (MQ=255)
                                                     cGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATTAGTTTTGCCCCCTCCTTTCTCGATGCAAAAACACGCTCTGGTTATCGTCTTGGAGACggtc                                                                                                                                                                               >  1:142219/1‑249 (MQ=255)
                                                       tgtTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTAAGGCAATGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTCAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGAGGGACAATAAGGGGTTAACCCCCTACTTAATCGATGGAAAAGCACGCAATGGTTATAGTCTTAGAGACgggc                                                                                                                                                                               >  1:219543/1‑247 (MQ=255)
                                                                    gCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTaa                                                                                                                                                                <  2:219543/251‑1 (MQ=255)
                                                                       ggggCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATc                                                                                                                                                              <  1:109144/250‑1 (MQ=255)
                                                                                                                                             caTAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGGGGTTAGCCCCCTACTTACTCGATGAAAAAGCAAGCCTTGGGTATCGTCTTGGAGACGGGCAGGGTTATTAAGTAATCCTGCGCGATGGCCTGCTGTGCGCCCCCCATGGTTATTCTTTGGGGGTTACCGCCCaaaaaggtgctaa                                                                                         >  2:248536/2‑239 (MQ=255)
                                                                                                                                                       gtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAAt                                                                             <  1:309341/250‑1 (MQ=255)
                                                                                                                                                          gTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAGAAGGGTTTAGCCCCCTACCTACTCGATGGCAAAGCACGACCTGGTTATCGTCTTGGAGACGGGCAGGCTTATAACGTAATTCTGCGAGATGGCCTGGTGTGGGCCACCCATGGCTATTATATGGGGGTTAACGCCGGAAAAGTGGCTCAagagaaagggaatacc                                                                         >  2:322770/1‑239 (MQ=255)
                                                                                                                                                                       aaGTGTGGCGCATGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGCTGGAAAAAATGCGTTTTGCCCCCCACTTACTCGGTGCAAAAACACGCTCTGGTTATCGGCTTTGGGACGGGCAGGTTTATGACGTAAACCTGCGCGATGGCCTGATGTGGGCCACCCATGGGTTTCATATGGGGATTACCGCCGCAAAAGCGGGTAAAAAGTACGGAATTACCCGGTAAATGCAgg                                                            >  2:187941/1‑251 (MQ=255)
                                                                                                                                                                           gtgGCGCTTGACGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGTTGGAAAATATGAGTTTAGGCCCCTACTTACTCGACGCAAAAGAACGCACCGGTTATCGTCTTGGAGGAGGACAGGGTTTTGACGGAATCCTGCGCTATGGCCTCATGTGGGCCCCCCCAGGGTATCATTTGGGGGTTACTCCCCGAAAACTGGCGAAAGAGGACGGAATTACCCGTGTAATGAAgggtgg                                                        >  2:32863/1‑247 (MQ=255)
                                                                                                                                                                                   tGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGgc                                                  <  1:322990/249‑1 (MQ=255)
                                                                                                                                                                                   cgccgccCATACCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGCCTTGGAGACGGACAAGTTTATGCCGTAATCCTGCGCGATGGCCTGAGGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAAcagacg                                                 <  2:8407/249‑6 (MQ=255)
                                                                                                                                                                                                                    gagCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCCAAAGCACGCCCTGGTTATCGTCTTTGAGACGGACCGGGTTATGACGCAATCCTTCGCGATGGCCTGATGTGCGCCACCCCTGGTTATCATATTGGGGTTATCGCCGAAAACGTGGCTAAAGAGCACGGAATTACCCCTGAAATGCAGGATGAACTGGCGCTACATTCACAACGTAAAGCGGCAGCCGCAAttt               >  1:310098/1‑250 (MQ=255)
                                                                                                                                                                                                                         gtgtgGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTcc          <  1:204181/250‑1 (MQ=255)
                                                                                                                                                                                                                           ttgggggggggTATGGAAAAATTGAGTTTAGCCCCCCACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTTATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCATCCGCAATTGAGTCCg         <  2:333055/245‑1 (MQ=255)
                                                                                                                                                                                                                             ggCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTg      <  1:143434/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gggggTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCttt  >  2:281372/1‑251 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
ATTGAACGTGCAAAAATCGATTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTT  >  NC_000913/2326220‑2326695

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: