Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,571,6470TC41.0% 63.1 / 10.8 39E147G (GAA→GGA) gntRgluconate operon transcriptional repressor GntR
Reads supporting (aligned to +/- strand):  ref base T (17/6);  new base C (0/16);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 2.06e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAAC  >  NZ_CP009273/3571509‑3571782
                                                                                                                                           |                                                                                                                                        
cGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGtt                                                                                                                                          <  2:54156/139‑1 (MQ=255)
 gTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTc                                                                                                                                         <  2:188974/139‑1 (MQ=255)
 gTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTc                                                                                                                                         <  2:390636/139‑1 (MQ=255)
     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:422557/1‑139 (MQ=255)
                                                                               ctTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTc                                                          >  1:708631/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                                     >  2:694968/1‑139 (MQ=255)
                                                                                                    cGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGa                                     >  1:926373/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)
                                                                                                                                           |                                                                                                                                        
CGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAAC  >  NZ_CP009273/3571509‑3571782

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