Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273529,7710CT40.0% 43.3 / 10.1 30R134C (CGT→TGT) gclglyoxylate carboligase
Reads supporting (aligned to +/- strand):  ref base C (2/16);  new base T (12/0);  total (14/16)
Fisher's exact test for biased strand distribution p-value = 1.05e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTCCGCTTCTGCTGATTCCATTCCTATTCTGTGCATTACCGGCCAGGCACCGCGCGCCCGTCTGCATAAAGAAGATTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGGTTACAGTTCGTGAAGCGGCGCTGGTGCCTCGC‑GTGCTGCAACAGGCATTTCACC‑TGATGCGTTCTGGTCGTCCGGGTCCGGTACTGGTGGATTTACCGTTCGACGTTCAGGTTGCGGAAATCGAGTTTGATCCTGACATGTA  >  NZ_CP009273/529635‑529904
                                                                                                                                        |                                                                                                                                       
ttCCGCTTCTGCTGATTCCATTCCTATTCTGTGCATTACCGGCCAGGCACCGCGCGCCCGTCTGCATAAAGAAGATTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGGTTACAGTTCg                                                                                                                                        <  1:682329/138‑1 (MQ=255)
        ctgctgATTCCATTCCTATTCTGTGCATTACCGGCCAGGCACCGCGCGCCCGTCTGCATAAAGAAGATTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGGTTACAGTTCGTGAAgcgg                                                                                                                                <  1:177808/138‑1 (MQ=255)
        ctgctgATTCCATTCCTATTCTGTGCATTACCGGCCAGGCACCGCGCGCCCGTCTGCATAAAGAAGATTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGGTTACAGTTCGGGGAggggc                                                                                                                               >  2:554589/1‑139 (MQ=255)
                               tgCATTACCGGCCAGGCACCGCGCGCCCGTCTGCATAAAGAAGATTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGGTTTCAGTTTGTGAAAGGGCGCTGGGGCCCCGC‑GTGCTGCaa                                                                                                        >  1:424970/1‑138 (MQ=255)
                               tgCATTACCGGCCAGGCACCGCGCGCCCGTCTGCATAAAGAAGATTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGGTTACAGTTTGTGGAGCGGCGCTGGGGCCTCGC‑Gtgctgcca                                                                                                        >  1:163413/1‑136 (MQ=255)
                               tgCATTACCGGCCAGGCACCGCGCGCCCGTCTGCATAAAGAAGATTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGGTTAAAGTTTGTGAAACGGGGCTGGTGCCCCGC‑Gggctgcga                                                                                                        >  2:891463/1‑136 (MQ=255)
                                      ccGGCCAGGCACCGCGCGCCCGTCTGCATAAAGAAGATTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGGTTACAGTTCGTGAAGCGGCGCTGGTGCCTCGC‑GTGCTGCAACAGGCAt                                                                                                 <  1:608359/138‑1 (MQ=255)
                                         gCCAGGCACCGCGCGCCCGTCTGCATAAAGAAGATTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGGTTTCAGTTTGTGGAGCGGCGCCGGGGGCCCGCGGGGCGGCACAAGGCAttta                                                                                              >  1:96672/1‑138 (MQ=255)
                                         gCCAGGCACCGCGCGCCCGTCTGCATAAAGAAGATTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGGTTAAAATTTGGGGAGCGGGGCCGGGGCCTCGC‑GTGGTGCAACAAGGAATTCa                                                                                             >  2:369184/1‑139 (MQ=255)
                                                    cgcgcCCGTCTGCATAAAGAAGATTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGGTTAAAGTTTGTGAAGAGGGGCTGGCGCCCCGC‑GCGCTGCAAAAAGAATTTCACC‑TGAAGCggt                                                                                 >  2:672725/1‑137 (MQ=255)
                                                                cATAAAGAAGATTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGGTTACAGTTCGTGAAGCGGCGCTGGTGCCTCGC‑GTGCTGCAACAGGCATTTCACC‑TGATGCGTTCTGGTCgtcc                                                                       <  2:459813/137‑1 (MQ=255)
                                                                cATAAAGAAGATTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGGTTACAGTTCGTGAAGCGGCGCTGGTGCCTCGC‑GTGCTGCAACAGGCATTTCACC‑TGATGCGTTCTGGTCgtcc                                                                       <  1:73144/137‑1 (MQ=255)
                                                                    aagaagATTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAGTGGCGGTTACAGTTCGTGAAGCGGCGCTGGTGCCTCGC‑GTGCTGCAACAGGCATTTCACC‑TGATGCGTTCTGGTCgtccgggtc                                                                  <  1:913312/138‑1 (MQ=255)
                                                                       aagaTTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGGTTTAAATTTGTGAAGAGGGGCTGGGGGCCCCC‑GTGCTGCAAACGGCATTTTTCC‑TGATGCGGTTTGGGCgtcccgggccggg                                                              >  2:951601/1‑138 (MQ=255)
                                                                          aTTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGGTTTCAGCTTGTGGAGAGGCGCTGGGGCCGCGC‑GGGCTGGAGAAGGAATTTTACC‑TGATGCGCTCTGGGCgtccggggcccgggat                                                           >  2:857441/1‑133 (MQ=255)
                                                                          aTTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGGTTGCAGTTCGTGAAGCGGCGCTGGTGGCCCGC‑GGGCTGCCACAGGCATTTTACC‑CGCTGCGTTCTTGTCGTCCGGGGCCCGGACt                                                           >  2:666243/1‑139 (MQ=255)
                                                                                      tAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGGTTACAGTTCGTGAAGCGGCGCTGGTGCCTCGC‑GTGCTGCAACAGGCATTTCACC‑TGATGCGTTCTGGTCGTCCGGGTCCGGTACTGGTGGATTTAcc                                               <  1:554589/139‑1 (MQ=255)
                                                                                                   aaTTGCTAAACCGGTCAGCAAAATGGCGGTTTCAGTTTGTGAAACGGGGCTGGGGGCTCGC‑GTGCTGCAACAGGGAATTTACC‑TGATGCGGTCTGGTCGTCCGGGGCCGGGACTGGGGGGTTTTCCGTTCGGCGTTCAg                                  >  1:410186/1‑139 (MQ=255)
                                                                                                                gTCAGCAAAATGGCGGTTTCAGTTTGTGAAGCGGCGCTGGTGCCTCGC‑GTGCTGCAACAAGGATTTTTCC‑CGATGGGGTCTGGTCGTCCCGGGCCGGGACTGGGGGGTTTACCGTTCGACGTTCAGGGTGCGGGGATCg                     >  2:841257/1‑139 (MQ=255)
                                                                                                                gTCAGCAAAATGGCGGTTTCAGTTTGGGAAGCGGGGCTGGGGCCTCGC‑GTGCTGCCACAAGGATTTCACC‑TGATGCGGTCTGGGCGTCCGGGGCCGGGACTGGGGGGGTTTACGTTCGACGTTCCGGTTGCggggaacg                     >  2:569226/1‑133 (MQ=255)
                                                                                                                gTCAGCAAAATGGCGGTTACAGTTTTGGGAACGGCGCCGGGGCCTCGC‑GCGCTGTCACAAGGCTTTCACCTTGAT‑CGGTCTGGGCGTCCCGGGCCGGGGCCGGGGGGTTTTCACTTCTACACTCCGGTTGTGGAAAACg                     >  2:663571/1‑139 (MQ=255)
                                                                                                                    gCAAAATGGCGGTTACAGTTCGTGAAGCGGCGCTGGTGCCTCGC‑GTGCTGCAACAGGCATTTCACC‑TGATGCGTTCTGGTCGTCCGGGTCCGGTACTGGTGGATTTACCGTTCGACGTTCAGGTTGCGGAAATCGAGtt                 <  1:369184/139‑1 (MQ=255)
                                                                                                                    gCAAAATGGCGGTTACAGTTCGTGAAGCGGCGCTGGTGCCTCGC‑GTGCTGCAACAGGCATTTCACC‑TGATGCGTTCTGGTCGTCCGGGTCCGGTACTGGTGGATTTACCGTTCGACGTTCAGGTTGCGGAAATCGAGtt                 <  1:569226/139‑1 (MQ=255)
                                                                                                                    gCAAAATGGCGGTTACAGTTCGTGAAGCGGCGCTGGTGCCTCGC‑GTGCTGCAACAGGCATTTCACC‑TGATGCGTTCTGGTCGTCCGGGTCCGGTACTGGTGGATTTACCGTTCGACGTTCAGGTTGCGGAAATCGAGtt                 <  2:5128/139‑1 (MQ=255)
                                                                                                                    gCAAAATGGCGGTTACAGTTCGTGAAGCGGCGCTGGTGCCTCGC‑GTGCTGCAACAGGCATTTCACC‑TGATGCGTTCTGGTCGTCCGGGTCCGGTACTGGTGGATTTACCGTTCGACGTTCAGGTTGCGGAAATCGAGtt                 <  1:474962/139‑1 (MQ=255)
                                                                                                                         aTGGCGGTTACAGTTCGTGAAGCGGCGCTGGTGCCTCGC‑GTGCTGCAACAGGCATTTCACC‑TGATGCGTTCTGGTCGTCCGGGTCCGGTACTGGTGGATTTACCGTTCGACGTTCAGGTTGCGGAAATCGAGTTTGATc            <  2:424970/139‑1 (MQ=255)
                                                                                                                            gCGGTTACAGTTCGTGAAGCGGCGCTGGTGCCTCGC‑GTGCTGCAACAGGCATTTCACC‑TGATGCGTTCTGGTCGTCCGGGTCCGGTACTGGTGGATTTACCGTTCGACGTTCAGGTTGCGGAAATCGAGTTTGATCCt          <  2:163413/138‑1 (MQ=255)
                                                                                                                             cGGTTACAGTTCGTGAAGCGGCGCTGGTGCCTCGC‑GTGCTGCAACAGGCATTTCACC‑TGATGCGTTCTGGTCGTCCGGGTCCGGTACTGGTGGATTTACCGTTCGACGTTCAGGTTGCGGAAATCGAGTTTGATCCTGa        <  1:614634/139‑1 (MQ=255)
                                                                                                                              ggTTACAGTTCGTGAAGCGGCGCTGGTGCCTCGC‑GTGCTGCAACAGGCATTTCACC‑TGATGCGTTCTGGTCGTCCGGGTCCGGTACTGGTGGATTTACCGTTCGACGTTCAGGTTGCGGAAATCGAGTTTGATCCTGAc       <  2:644899/139‑1 (MQ=255)
                                                                                                                                   cAGTTCGTGAAGCGGCGCTGGTGCCTCGC‑GTGCTGCAACAGGCATTTCACC‑TGATGCGTTCTGGTCGTCCGGGTCCGGTACTGGTGGATTTACCGTTCGACGTTCAGGTTGCGGAAATCGAGTTTGATCCTGACATGTa  <  1:891463/139‑1 (MQ=255)
                                                                                                                                        cGTGAAGCGGCGCTGGTGCCTCGC‑GTGCTGCAACAGGCATTTCACC‑TGATGCGTTCTGGTCGTCCGGGTCCGGTACTGGTGGATTTACCGTTCGACGTTCAg                                  <  1:915249/102‑1 (MQ=255)
                                                                                                                                        cGTGAAGCGGCGCTGGTGCCTCGC‑GTGCTGCAACAGGCATTTCACC‑TGATGCGTTCTGGTCGTCCGGGTCCGGTACTGGTGGATTTACCGTTCGACGTTCAg                                  >  2:915249/1‑102 (MQ=255)
                                                                                                                                        |                                                                                                                                       
TTCCGCTTCTGCTGATTCCATTCCTATTCTGTGCATTACCGGCCAGGCACCGCGCGCCCGTCTGCATAAAGAAGATTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGGTTACAGTTCGTGAAGCGGCGCTGGTGCCTCGC‑GTGCTGCAACAGGCATTTCACC‑TGATGCGTTCTGGTCGTCCGGGTCCGGTACTGGTGGATTTACCGTTCGACGTTCAGGTTGCGGAAATCGAGTTTGATCCTGACATGTA  >  NZ_CP009273/529635‑529904

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