Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273369,0930TG53.3% 16.1 / 19.1 15G239G (GGT→GGGmhpFacetaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (3/4);  new base G (8/0);  total (11/4)
Fisher's exact test for biased strand distribution p-value = 2.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-04
Rejected as consensus: Frequency below/above cutoff threshold.

ATTGAGCGACGAAGCTTCACAAGATGATATCGAAGCCTCAATCAATGAAATGGCTGAGGCGGTGCAGGCTTACGTACCGGGTTATCGCCTGAAACAGCGCGTGCAGTTTGAAGTTATCCCGCAGGATAAACCGGTCAATTTACCGGGCGTGGGGCAATTCTCCGGACTG  >  NZ_CP009273/369012‑369180
                                                                                 |                                                                                       
atTGAGCGACGAAGCTTCACAAGATGATATCGAAGCCTCAATCAATGAAATGGCTGAGGCGGTGCAGGCTTACGTACCGGGGTATCGCCt                                                                                 >  2:405099/1‑90 (MQ=255)
  tGAGCGACGAAGCTTCACAAGATGATATCGAAGCCTCAATCAATGAAATGGCTGAGGCGGTGCAGGCTTACGTACCGGGTTAt                                                                                      <  1:338081/83‑1 (MQ=255)
  tGAGCGACGAAGCTTCACAAGATGATATCGAAGCCTCAATCAATGAAATGGCTGAGGCGGTGCAGGCTTACGTACCGGGTTAt                                                                                      >  2:338081/1‑83 (MQ=255)
  tGAGCGACGAAGCTTCACAAGATGATATCGAAGCCTCAATCAATGAAATGGCTGAGGCGGGGGAGGGTTTCGTTCCGGGGTATCGCCTGa                                                                               >  1:441065/1‑90 (MQ=255)
  tGAGCGACGAAGCTTCACAAGATGATATCGAAGCCTCAATCAATGAAATGGCTGAGGCGGGGCAGGCTTACGTACCGGGGTATCGCCTGa                                                                               >  1:393996/1‑90 (MQ=255)
  tGAGCGACGAAGCTTCACAAGATGATATCGAAGCCTCAATCAACGAAATGGCTGAGGCGGTGCAGGCTTACGTACCGGGGTATCGCCTGa                                                                               >  1:4299/1‑90 (MQ=255)
   gAGCGACGAAGCTTCACAAGATGATATCGAAGCCTCAATCAATGAAATGGCTGAGGCGGTGCAGGGTTACGTACCGGGGTATCGCCTGaa                                                                              >  1:258023/1‑90 (MQ=255)
              cTTCACAAGATGATATCGAAGCCTCAATCAATGAAATGGCTGAGGCGGTGCAGGCTTACGTACCGGGTTATCGCCTGAAACAGCGCGTGc                                                                   <  1:417417/90‑1 (MQ=255)
                   caAGATGATATCGAAGCCTCAATCAATGAAATGGCTGAGGCGGTGCAGGCTTACGTACCGGGTTATCGCCTGAAACAGCGCGTGCAGttt                                                              <  2:343583/90‑1 (MQ=255)
                               gAAGCCTCAATCAATGAAATGGCTGAGGCGGTGCAGGCTTACGTACCGGGTTATCGCCTGAAACAGCGCGTGCAGTTTGAAGTTATCCCg                                                  >  2:18111/1‑90 (MQ=255)
                                            atGAAATGGCTGAGGCGGTGCAGGCTTACGTACCGGGGTATCGCCTGAAAAAGCGCGTGCAGTTTGAAGTTATCCCGCAGGATAAACCgg                                     >  1:169836/1‑90 (MQ=255)
                                                    gCTGAGGCGGTGCAGGCTTACGTACCGGGGTATCGGCTGAAACAACGCGTGCAGTTTGGAGATATCCCCCAGGATAAACCGGGCAATTTa                             >  1:469932/1‑90 (MQ=255)
                                                                  ggCTTACGTACCGGGGTATCGCCTGGAAAAGCGCGGGCAGGTTGAAGTTATCCCGCCGGATAAACCGGGCAATTTACCggggggggggga               >  1:393794/1‑88 (MQ=255)
                                                                              gggTTATCGCCTGAAACAGCGCGTGCAGTTTGAAGTTATCCCGCAGGATAAACCGGTCAATTTACCGGGCGTGGGGCAATTCTCCGGACt   >  2:145534/1‑90 (MQ=255)
                                                                               ggTTATCGCCTGAAACAGCGCGTGCAGTTTGAAGTTATCCCGCAGGATAAACCGGTCAATTTACCGGGCGTGGGGCAATTCTCCGGACTg  <  1:145534/90‑1 (MQ=255)
                                                                                 |                                                                                       
ATTGAGCGACGAAGCTTCACAAGATGATATCGAAGCCTCAATCAATGAAATGGCTGAGGCGGTGCAGGCTTACGTACCGGGTTATCGCCTGAAACAGCGCGTGCAGTTTGAAGTTATCCCGCAGGATAAACCGGTCAATTTACCGGGCGTGGGGCAATTCTCCGGACTG  >  NZ_CP009273/369012‑369180

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: