Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,927,6690GC52.6% 18.2 / 28.9 19A398G (GCG→GGG) eddphosphogluconate dehydratase
Reads supporting (aligned to +/- strand):  ref base G (8/1);  new base C (0/10);  total (8/11)
Fisher's exact test for biased strand distribution p-value = 1.19e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-05
Rejected as consensus: Frequency below/above cutoff threshold.

ACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGT  >  NZ_CP009273/1927587‑1927748
                                                                                  |                                                                               
cccGTTTAACTTTTTTCCCCCCCCATGAGGAGAAAAGGTTTGTTCGAAGGAAGCGATCCCTTTGCTGTGGAGTGTTTTTTCCCCCCCtcc                                                                          <  1:423992/89‑3 (MQ=255)
 ccGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTcc                                                                         >  1:281377/1‑90 (MQ=255)
 ccGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTcc                                                                         >  1:372869/1‑90 (MQ=255)
 ccGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTcc                                                                         >  1:484831/1‑90 (MQ=255)
 ccGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTcc                                                                         >  1:57901/1‑90 (MQ=255)
     ttAACCTTTTTGCCCCCCCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACTTTGCTGTCGAGTGTTTTTCCCCCCCCTTCCCGcc                                                                     <  1:199595/90‑1 (MQ=255)
     taaaCTCTTTTGCCCCCCCAGGAGGAGAAAAAGGTTGTTCAAAGGAAGCGATCACATTGCTGTCGAGTGTTTTTCCCCCCCCTTCCCGcc                                                                     <  1:169622/88‑1 (MQ=255)
      tAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGcca                                                                    >  1:275884/1‑90 (MQ=255)
               tGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTc                                                           <  2:372869/90‑1 (MQ=255)
                    ccccaTGAGGAGAAAAAGTTTGTCCGAAGGAAGCGACCACTTTGCTGTGGTGTTTTTTTCCCCCCCCTCCCCGCCAGTCCAGTTCACCat                                                      <  1:435516/88‑1 (MQ=255)
                    caccaTGATGAGAAAAAGGTTTTTCGAAGGAAGCGATCACTTTGCTGTGGAGTGATTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCat                                                      <  2:281377/90‑1 (MQ=255)
                         tgatgaGAAAAGGGTTGTTCGAAGGAAGCGATCACTTTGCTGTCGAGTGTTTTTTCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTc                                                 <  1:142433/90‑1 (MQ=255)
                         tgaggaGAGAAGGGTTGTTCGAAGGAAGCGATCCCTTTCCTGTGGAGTGTTTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTc                                                 <  1:426964/90‑1 (MQ=255)
                               gagaAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCAt                                           >  1:481729/1‑90 (MQ=255)
                                          ttCGAAGGAAGCGATCCCTTTGTTGTCGAGTTTTTTTCCCCCCCCTCCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGCTCAAGGGt                                <  2:484831/90‑1 (MQ=255)
                                                             ttGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCaaaa             >  1:427202/1‑90 (MQ=255)
                                                               gCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAAcc           >  1:213627/1‑90 (MQ=255)
                                                                     gAGTGTTTTTCCCCCCCCTCCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCAc     <  1:461919/90‑1 (MQ=255)
                                                                        tatttttcCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGt  <  2:229272/87‑1 (MQ=255)
                                                                                  |                                                                               
ACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGT  >  NZ_CP009273/1927587‑1927748

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