Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,290,2430TG57.7% 21.4 / 30.1 26L87F (TTA→TTC) ccmBheme exporter protein CcmB
Reads supporting (aligned to +/- strand):  ref base T (2/9);  new base G (15/0);  total (17/9)
Fisher's exact test for biased strand distribution p-value = 1.76e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.78e-05
Rejected as consensus: Frequency below/above cutoff threshold.

TACCAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAAGCAGCGC  >  NZ_CP009273/2290159‑2290332
                                                                                    |                                                                                         
tACCAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGGAAGGGGAACgg                                                                                      >  1:215820/1‑90 (MQ=255)
  ccAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGca                                                                                    <  2:249437/90‑1 (MQ=255)
        gCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGGAAGGGGAAAGGgcacaaaa                                                                              >  1:571459/1‑88 (MQ=255)
            aaGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGGAAGGGGAACGGcaacaacaacaa                                                                          >  1:1267984/1‑90 (MQ=255)
             aGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAAt                                                                         <  2:1280841/90‑1 (MQ=255)
                aTGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGt                                                                      <  2:971676/90‑1 (MQ=255)
                 tGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGGCTGGGGAAAAGCAACAACATCAAATGtt                                                                     >  2:887587/1‑90 (MQ=255)
                          gCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGc                                                            <  2:978202/90‑1 (MQ=255)
                            aGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACAACAATTGTTCAAAAccgccg                                                          >  1:317742/1‑90 (MQ=255)
                            aGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAAAAACAACAATTTTTCAAGACTGCCg                                                          >  1:445375/1‑90 (MQ=255)
                             gACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGt                                                         <  1:638396/90‑1 (MQ=255)
                             gACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCCATTGTTCAAGACTTCCgg                                                         >  2:504560/1‑89 (MQ=255)
                             gACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAAAGGCAACAACAACAATTTTTTAAAACTGCCgc                                                         >  1:324347/1‑89 (MQ=255)
                                  gTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGGAAGGGGAACGGCCACAACATCAATTGTTCAAGACTGCCGGCCCGc                                                    >  1:865941/1‑90 (MQ=255)
                                  gTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGc                                                    >  2:570654/1‑90 (MQ=255)
                                        aTCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCGACAACATCCATTTTTCCAGAACGCCGCCCTGCAAAACg                                              >  1:590234/1‑90 (MQ=255)
                                           cAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCa                                           >  2:860828/1‑90 (MQ=255)
                                                 gcCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAAACGTCACGGAAc                                     >  1:424703/1‑90 (MQ=255)
                                                  ccATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACa                                    <  2:180556/90‑1 (MQ=255)
                                                  ccATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACa                                    <  2:399986/90‑1 (MQ=255)
                                                          ttCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTcc                            >  2:136469/1‑90 (MQ=255)
                                                                     aaCGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCaa                 <  1:1094195/90‑1 (MQ=255)
                                                                      aCGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAg                >  1:720/1‑90 (MQ=255)
                                                                           gggCCAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGGCCCGCAAAACGGCACGGGAAAGTCATTCCAGCGCCAACCAGGGTGa           >  1:879003/1‑90 (MQ=255)
                                                                                  ggTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAagcagc    <  2:1372738/90‑1 (MQ=255)
                                                                                   gTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAAgcagcg   >  2:1005399/1‑90 (MQ=255)
                                                                                    tAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAAgcagcgc  <  1:90608/90‑1 (MQ=255)
                                                                                    |                                                                                         
TACCAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAAGCAGCGC  >  NZ_CP009273/2290159‑2290332

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