Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,571,6500AC45.7% 47.2 / 11.1 35V146G (GTG→GGG) gntRgluconate operon transcriptional repressor GntR
Reads supporting (aligned to +/- strand):  ref base A (15/4);  new base C (0/16);  total (15/20)
Fisher's exact test for biased strand distribution p-value = 1.90e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAAC  >  NZ_CP009273/3571514‑3571782
                                                                                                                                         |                                                                                                                                     
gtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTccac                                                                                                                                      <  1:185073/138‑1 (MQ=255)
 tcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTccacca                                                                                                                                    >  2:104646/1‑139 (MQ=255)
 tcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTccacca                                                                                                                                    >  1:483292/1‑139 (MQ=255)
      aCGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGa                                                                                                                               <  2:323605/139‑1 (MQ=255)
          gcgcCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATAc                                                                                                                          >  2:514173/1‑139 (MQ=255)
                       cAATGTGGCGATGCCCGCGAGAAATAATGGCAGTGGTCATCGGGCGGGCTGCTCAAAATTTATCAAAACCGCCGGCGATATAAGGGATGGGCGATTTGCTGTCCATCATTCCCCCCAGGGGAATACCCGCCACTTCAat                                                                                                             <  1:665045/139‑1 (MQ=255)
                       cAATGTGGCGATCCCCGCGAGAAATAATGGCAGTGGTCATCTGGGGTGCTGCTTAAAAGTTATCAAACCCGCGGGCGATATCAAGGCAGGGCGACTTGCTGTCCATCAGTCCCCCCACGGAAATCCCCGCCACTTCAat                                                                                                             <  1:159587/139‑1 (MQ=255)
                       cAATGGGGCGATGCCCGCGAGAAATAAGGGCAGTGGTCATTGGGCGTGCTGCTTAAAAGTTAAAAAAACCGACGGCGATATAAAGGCAGGGCGACTTGCTGCCCATCAGTCCCCCCAGGGGAATACCCGCCACTTCAat                                                                                                             <  1:100206/139‑1 (MQ=255)
                       cAATGGGGCGAGGCCCGCGAGAAATAATGGCAGTGGTCATCTGGCGTGCTGTTAAAAAGTTATAAAACCCGCCGGCGATATCAAGGCATGGCGACTTGCTGCCCATCAGTTCCCCCACGGAAATACCCGCCACTTCAat                                                                                                             <  1:514173/139‑1 (MQ=255)
                       cAAGGTGGCGAGGCCCGGGAGAAATAAGGGCAGGGGTCCTCTGGCGTGCTGCTTAAAAGTTAACAAAACCGACGGCAAAACCAAGGCATGGCGACTTCCTCTCCACAAGTCCCCCCACGGGAATACCCGCCACTTCAat                                                                                                             <  2:43518/139‑1 (MQ=255)
                       aattggggCGTTCCCCGCGAGAAATAGTGGCGGTGGTCATTGGGCGTGCTGTTTCAAATTTAAAAAAACCGACGGCGATATCAAGGCTTGGCGTTTTCCTGTCCATCGGTCCCCCCACGGAAATCCCCGCCACTTCAat                                                                                                             <  2:635537/133‑1 (MQ=255)
                       aagtgtgGCGATGCCCGCGAGAAATAAGGGCATGGTTCACCGGGCGTGCTTCTGAAAATTTACAAAACCCGACGGCGATATCAAGGCAGGGCGTCTTGCTGCCCATAAGTCCCACCAGGGGAATCCCCGCCACTTCAat                                                                                                             <  2:256443/136‑1 (MQ=255)
                               cGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCttt                                                                                                      >  2:740819/1‑138 (MQ=255)
                                        cgAGCAATAGGGGAAGTGGTCTTCTGGCGTGCTGCTTAAAATTTAACAAACCCGCCGGCGATACCAGGGCATGCCGACTTGCTGTCCATCATTCCCACCCCGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGcgcg                                                                                            <  2:36877/139‑1 (MQ=255)
                                        cgAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGcgcg                                                                                            >  2:860272/1‑139 (MQ=255)
                                        cgAGAAATAAGGGGGGTGGTCATCGGGCGTGCTGCTTAAAAGTTATCAAACCCGACGGCGATAAAAAGGCGTGGCGCCTTGCTGCCCATCAGTCCCCCCACGGAAATACCCGCCACTTCAATCATCTTTAAGGTGcgcg                                                                                            <  1:104646/139‑1 (MQ=255)
                                        cgAGAAATAAGGGCGGTGGTCTTCTGGCGTGTTGTTAAAAATTTAAAAAACCCGACGGCAATATCAGGGCATGGCGATTTGCTGCCCATCAGTCCCCCCACGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGcgcg                                                                                            <  1:215281/139‑1 (MQ=255)
                                        cgAACAATAAGGGGGGGGGTCATCTGGCGTGCTGTTAAAATGTTATAAAAACCGCGGGTGATAAAAAGGCATGGCGTTTTGTTGTCCAGCGTTCCCACCAGGGGAATACCCGCCACTTCAATCATCTTTAAGGTGcgcg                                                                                            <  2:937916/139‑1 (MQ=255)
                                              ataatGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCgtgt                                                                                      >  2:229738/1‑139 (MQ=255)
                                                     cAGTGGTCATCGGGCGTGCTGCTTAAAAGTTATAAAACCCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCATTCCCCCCGGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACg                                                                               <  1:511984/139‑1 (MQ=255)
                                                       gTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGtt                                                                             >  2:430194/1‑139 (MQ=255)
                                                           tCATTTGGCGTGCTGCTTAAAATTTACAAAACCCGACGGCGATATAAAGGCATGGCGTTTTGCTGTCCATCAGTCCCCCCACGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGt                                                                         <  1:155967/139‑1 (MQ=255)
                                                                     tgctgcTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAgg                                                               >  2:727107/1‑139 (MQ=255)
                                                                          ctTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTc                                                          >  1:708631/1‑139 (MQ=255)
                                                                          ctTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTc                                                          >  1:422557/1‑139 (MQ=255)
                                                                             aaaaGTTAAAAAACCCGACGGCGATACAAAGGCATGGCGATTTGCTGCCCACCAGTCCCCCCACGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGa                                                        <  1:860272/137‑1 (MQ=255)
                                                                                    aaaaaacCCGACGGGAATACCAGGGCAGGGCGCCTTGCTGCCCATCAGTCCCCCCAGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAg                                                <  2:926373/136‑1 (MQ=255)
                                                                                    aTCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAg                                                >  2:218358/1‑139 (MQ=255)
                                                                                     taaaacCCGGCGGCGATAAAAAGGCATGGCGACTTGCTGCCCATCAGTCCCCCCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAg                                                <  1:725871/136‑1 (MQ=255)
                                                                                            cGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCAt                                        >  2:253057/1‑139 (MQ=255)
                                                                                               cGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGa                                     >  1:926373/1‑139 (MQ=255)
                                                                                               cGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGa                                     >  2:694968/1‑139 (MQ=255)
                                                                                                                 ggCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCt                    >  1:324590/1‑139 (MQ=255)
                                                                                                                         gCTGTCCATCAGTTCCCCCACGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTc             <  1:210257/138‑1 (MQ=255)
                                                                                                                         gCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTcc            >  1:429552/1‑139 (MQ=255)
                                                                                                                           ggTCCAGCCGTCCCCCCGCGGGAACCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTc             <  1:576364/135‑1 (MQ=255)
                                                                                                                                  cAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATaa   >  1:718470/1‑139 (MQ=255)
                                                                                                                                    agtCCCCCCCGGGAAAACCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAAc  <  2:774955/137‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                     
GTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAAC  >  NZ_CP009273/3571514‑3571782

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