Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,255,3050TG57.9% 17.4 / 21.4 20A180A (GCA→GCCfruK1‑phosphofructokinase
Reads supporting (aligned to +/- strand):  ref base T (1/7);  new base G (11/0);  total (13/7)
Fisher's exact test for biased strand distribution p-value = 1.59e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-04
Rejected as consensus: Frequency below/above cutoff threshold.

CAGCTTCAATCACATCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTG  >  NZ_CP009273/2255220‑2255367
                                                                                     |                                                              
cAGCTTCAATCACATCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAGCCACGGTGCCg                                                            <  1:544312/90‑1 (MQ=255)
    ttCAATCACATCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGGTAGGGTTTACCCACCCCGGGGCCGCttt                                                        >  2:639821/1‑90 (MQ=255)
            caTCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAGCCACGGGGCCGCTTTCAAACCTg                                                >  1:686964/1‑90 (MQ=255)
                tttcatttcaGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCCGCCACGGGGCCGCTTTTAAAACTGctac                                            >  2:151093/1‑90 (MQ=255)
                      ttcaGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTTTGTTTCCCCCGCCCCGGGGCCCCCTTTAAACCTGCTACTaaacc                                      >  2:466374/1‑87 (MQ=255)
                          ggCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCa                                  <  1:154556/90‑1 (MQ=255)
                               tttACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGct                             <  1:466374/90‑1 (MQ=255)
                                 tACGGCCTGCCCAGATTTCCAGCTCGCGGCTGTTAGGTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGctac                           <  2:90368/90‑1 (MQ=255)
                                 tACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGctac                           <  2:232521/90‑1 (MQ=255)
                                 tACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGctac                           <  1:151093/90‑1 (MQ=255)
                                        tGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTTACCAGCCCCGGGGCCCCTTTCAAACCTGCTACTAACGCTTCACGGCCACTATCaaa                    >  2:308505/1‑90 (MQ=255)
                                        tGCCCAGATTTCCAGCTCGCGGCGGTTAGGGTTTACCAGCCCCGGGGCCGCCTTCAAAACCGCCACTAACGCCTCCCGGCTACTATCaaa                    >  2:628448/1‑90 (MQ=255)
                                        tGCCCAGATTTCCAGCTCGCGGCGGTTAGGGTTCACCAGCCACGGGGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCaaa                    >  1:507087/1‑90 (MQ=255)
                                        tGCCCAGATTTCCAGCTCGCGGCGGGTAGGTTTTACCAGCCCCGGGGCCCCCTTCAAACCCGCTCCTAACGCTTCCCGGGTACTATCaaa                    >  2:467202/1‑90 (MQ=255)
                                         gCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAGCCACGGGGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAg                   >  1:1491/1‑90 (MQ=255)
                                              gATTTCCAGCTCGCGGCGGTTAGGTTTCACCAGCCACGGGGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAAt              >  2:718386/1‑90 (MQ=255)
                                              gATTTCCAGCTCGCGGCGGTTAGGGTTCACCAGCCACGGCGCCGCTTTCAAACCTGCTGCTAACGCTTCACGGCTCCTATCAAAaataat              >  1:462062/1‑90 (MQ=255)
                                                        tCGCGGCGGTTAGGTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGacac    <  2:450795/90‑1 (MQ=255)
                                                        tCGCGGCGGTTAGGTTTCACCAGCCACGGTGCCCCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGacac    >  1:618426/1‑90 (MQ=255)
                                                          gcggcgGTTAGGGTTCACCAGCCACGGGGCCGCTTTCAAAACTGCTACTAACGCTTTACGGCTACTATCAAAGATAATGCAAAGacaccg  >  1:486804/1‑88 (MQ=255)
                                                                                     |                                                              
CAGCTTCAATCACATCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTG  >  NZ_CP009273/2255220‑2255367

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