Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,269,3990TG53.3% 12.6 / 18.1 15G331G (GGT→GGGghxPguanine/hypoxanthine transporter GhxP
Reads supporting (aligned to +/- strand):  ref base T (1/6);  new base G (8/0);  total (9/6)
Fisher's exact test for biased strand distribution p-value = 1.40e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-03
Rejected as consensus: Frequency below/above cutoff threshold.

TGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCCTGC  >  NZ_CP009273/4269273‑4269520
                                                                                                                                |                                                                                                                           
tGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGACGGCGGGGGGGGGGAAAAACGGgttt                                                                                                                 >  1:1751/1‑138 (MQ=255)
       atcatcGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGACGGCGGCGGGGGGGGAAAACGGTGTGACGGcga                                                                                                          >  1:93426/1‑137 (MQ=255)
                  tGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGGCGGCGGGGGGCGGGGAAAAACGGTTGGCGGCCATCACCGGTGGc                                                                                             >  2:115731/1‑139 (MQ=255)
                                              agcagcGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGTGACGGCGGCGGGGGGGAAAACACGGTTGTCGGCGATTACCGCTTGGGTGCTTTTGCTCCTGATTCTGTTCctct                                                                 >  1:357327/1‑139 (MQ=255)
                                                  gcGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGACGGCGGCGGGGGGGGAAACCCGGTTGTCGGCTATCAACCGTTGGGGGGTGGTTCTCCTTGTTTTTTTTCTCTCTcc                                                             >  2:412049/1‑139 (MQ=255)
                                                   cGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGGGGGGAAAACCGGTTTGAGGCGAACACCGTTGGTGGGCTGGTTCTCCCGGTTTTTTTTCTCTCTccc                                                            >  1:57669/1‑138 (MQ=255)
                                                   cGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCg                                                            <  2:99164/139‑1 (MQ=255)
                                                    gTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGACGGCGGCGGGGGGGTAAAACGGGTTTGAGGGCATTACCGGTGGGGGGGGTGTTCTCCCTGTTTTTTTTCTCTCTCCGc                                                           >  2:425442/1‑139 (MQ=255)
                                                             ggCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGGGGGAAAACCGGGTTGACGGCTCTCACCCTTGGGGGGCTGTTCCTCCTGATTCTGTTTCTCTCTCCGCTCTCTTAcc                                                  >  1:389218/1‑139 (MQ=255)
                                                             ggCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGACGGCGGCGGGCGGGGAAAACGGGTTGGAGGCCGTCCCCCGTGGGGTGCGGTTTCTCCTGGTTCTGTTTCTCTCTCCCCTCTCTTAcc                                                  >  2:149149/1‑139 (MQ=255)
                                                                      ggtgCGGCTCCGGCAGCGTTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCg                                          <  1:15124/138‑1 (MQ=255)
                                                                                     gCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATGCAACGGCTc                          <  1:241351/139‑1 (MQ=255)
                                                                                               tCGAGTCTGCGGCGGGTACGGCGGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATGCAACGGCTCCGGCGCTGAt                <  1:115731/139‑1 (MQ=255)
                                                                                                             ggTACGGCGGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCctgc  <  2:437965/139‑1 (MQ=255)
                                                                                                             ggTACGGCGGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCctgc  <  2:54533/139‑1 (MQ=255)
                                                                                                                                |                                                                                                                           
TGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCCTGC  >  NZ_CP009273/4269273‑4269520

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