Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,571,6470TC55.0% 8.2 / 24.1 20E147G (GAA→GGA) gntRgluconate operon transcriptional repressor GntR
Reads supporting (aligned to +/- strand):  ref base T (7/2);  new base C (0/11);  total (7/13)
Fisher's exact test for biased strand distribution p-value = 4.64e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCG  >  NZ_CP009273/3571514‑3571784
                                                                                                                                      |                                                                                                                                         
gtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAATTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTccac                                                                                                                                       <  1:151052/138‑1 (MQ=255)
gtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTccac                                                                                                                                       >  2:151052/1‑138 (MQ=255)
                     ggaaaTGGGGCGATGCCCGCGACAAATAATGGCGGGGGTCATCTGGCGTGCTGCTTAAAAGTTATCAAACCCGACGGCTATACCAGGGCATGGCGCTTTGCTGCCCATCATTCCCCCCACGGGAATCCCCGCCACTTCa                                                                                                                 <  1:47406/136‑1 (MQ=255)
                       cattgtgGCGAGGCCCGCGAGAAATAAGGGCGGTGGTCATCGGGCGTGCTGCTTAAAAGTTAAAAAACCCGACGGCGATATAAGGGCAGGGCGATTGGCTGTCCATCAGTCCCACCAGGGGAATACCCGCCACTTCAat                                                                                                               <  1:240579/136‑1 (MQ=255)
                          tgtgGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAatcat                                                                                                            >  2:240579/1‑139 (MQ=255)
                                gATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCttt                                                                                                        >  1:380716/1‑137 (MQ=255)
                                      cgcgAAAAATAAGGGAAGTGGTCATCTGGCGTGCTCCTTAAAAGTTACCAAACCCGACGGCAATATCAAGGCAGGGCGATTTGCTGTCCATCAGTCCCCCCACGGGAATCCCCGCCACTTCAATCATCTTTAAGGTgcg                                                                                                <  1:240611/139‑1 (MQ=255)
                                        ggagaAATAGGGGCATTGGTCATCTGGCGTGCGTCTAAAAATTAAAAAAACCCGACGGAGATACAAGGCCTGGGCGATTTGCTGTCCATCAGTCCCCCCGGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGcgcg                                                                                              <  2:399322/138‑1 (MQ=255)
                                        ggagaAATAATGGCGTTGGTTATCGGGCGTGCTGCTTAAAATTTACAAAAACCGACGGCGATATCAGGGCATGGCGTCTTGCTGTCCATCGTTCCCACCCGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGcgcg                                                                                              <  1:396868/138‑1 (MQ=255)
                                        ggagaAATAATGGCAGTGGTCCTCGGGCGTGCTGCTCAAAAGTTATAAAAACCGACGGCGATATCAGGGCAGGGCGCCTTGCTGCCCAGCGGTCCCACCGGGGGAATACCCGCCACTTCAATCATCTTTAAGGTGcgcg                                                                                              <  2:55848/138‑1 (MQ=255)
                                        caagaaataagGCCAGTGGTCATCTGGCTTGCTTCTACAAATTTAACAAACCCGACGGCGAAAACAGGCCTGGGCGACTTGCTGTCCATCATTCCCACCACGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGcgcg                                                                                              <  1:34404/134‑1 (MQ=255)
                                        caaGCAATAGTGCCGGTGGTCTTGTGGCGGTATGCTTAAATTTTAGCAAACCCGCCGGCAATATCAGGTCGTGGCGCCTTGCTGTCCAGCGGTCCCACCCCGGGAATACCCGCCACTTCAATCATCTTTAAGGTGcgcg                                                                                              <  1:63929/137‑1 (MQ=255)
                                          ataaataatGGCGTTGGTCGTCTGGCGTGCTGTTAAAAATTTAAAAAACCCGCCGGCGATAACAAGCCATGGCGATTTGCGTTCCATCTGTCCCCCCCCGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGcgcg                                                                                              <  2:413316/134‑1 (MQ=255)
                                                          gTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCgg                                                                            >  2:142641/1‑139 (MQ=255)
                                                                                  ttATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTc                                                     >  2:305527/1‑138 (MQ=255)
                                                                                                      atCAAGGCGTGGCGTCTTGCTGCCCATCAGTTCCCCCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAgg                                <  1:44737/139‑1 (MQ=255)
                                                                                                               tGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGc                       >  2:323710/1‑139 (MQ=255)
                                                                                                                      cTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCAtt                 >  2:32749/1‑139 (MQ=255)
                                                                                                                             tCCACCAGTCCCCCCAGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGtt          <  2:275384/139‑1 (MQ=255)
                                                                                                                               cATCGGTCCCCCCACGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGtt          <  2:100281/137‑1 (MQ=255)
                                                                                                                                     ttCCACCACGGGAATACCCGCCACTTCAATCATCTTTACGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGCGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCg  >  2:93796/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                         
GTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCG  >  NZ_CP009273/3571514‑3571784

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