Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,889,0080AC55.6% 10.9 / 18.5 18V44G (GTG→GGG) ygcQelectron transfer flavoprotein subunit alpha/FixB family protein
Reads supporting (aligned to +/- strand):  ref base A (3/5);  new base C (0/10);  total (3/15)
Fisher's exact test for biased strand distribution p-value = 6.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-04
Rejected as consensus: Frequency below/above cutoff threshold.

ACAGTTCATCGCCAAACGTTCCCGGCGGAAAGAGAACAACGTCCGCAGGCGTTCGCTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTT  >  NZ_CP009273/2888925‑2889091
                                                                                   |                                                                                   
ccAGTTCATCCCCAACCGTCCCCGGGGGAAAGAAAACACCGCCCGCGGGCGTTCGCTGCCACTGTTCCACCAACGCACCCAGCCCCTGtt                                                                               <  2:321815/89‑1 (MQ=255)
  aGTTCACCGCCAACCGTCCCCGGCGGAAAGAAAACAACGCCCGCGGGCGTTCGCTGCCCCTGTTCCACCAACGCACCCAGCCCCTGTTcc                                                                             <  2:725981/90‑1 (MQ=255)
        tCGCCAAACGTTCCCGGCGGAAAGAGAACAACGTCCGCAGGCGTTCGCTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACt                                                                       >  2:331474/1‑90 (MQ=255)
          gCCAAACGTTCCCGGCGGAAAGAGAACAACGTCCGCAGGCGTTCGCTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTAc                                                                     <  1:261532/90‑1 (MQ=255)
             aaaCGTTCCCGGCGGAAAGAGAACAACGTCCGCAGGCGTTCGCTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAgg                                                                  >  1:319268/1‑90 (MQ=255)
                cGTCCCCGGCGGAAAGAAAACACCGCCCCCGGCCGTCCGCTGCCACTGTTCCACCAACGCATCCAGCCCCTGTTCCGCCACTACAGGCTg                                                               <  2:730373/90‑1 (MQ=255)
                    cccGGCGGAAAGAGAACAACGTCCGCAGGCGTTCGCTGCCACTGTTCCACCAACGCATCCAGCCCCTGTTCCGCCACTACAGGCTGAGGt                                                           <  1:56313/90‑1 (MQ=255)
                    cccGGCGGAAAGAGAACAACGTCCGCAGGCGTTCGCTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGt                                                           <  1:1597/90‑1 (MQ=255)
                    cccGGCGGAAAGAAAACACCGCCCGCGGGCGTTCGCTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGt                                                           <  2:319268/90‑1 (MQ=255)
                        gcggAAAGAGAACAACGTCCGCAGGCGTTCGCTGCCACTGTTCCCCCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCaa                                                       <  1:175790/90‑1 (MQ=255)
                                      ccgcccGCGGGCGTCCGCTGCCACTGTTCCCCCAACGCACCCAGCCCCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGcc                                         <  2:655792/86‑1 (MQ=255)
                                      acctCCGCGGGCTTTCGCTGCCACTGTTCCACCAACGCACCCAGCCCCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGcc                                         <  2:7429/87‑1 (MQ=255)
                                      acGTCCGCGGGCGTTCGCTGCCACTGTTCCCCCAACGCATCCAGCCCCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGcc                                         <  1:5686/90‑1 (MQ=255)
                                      acGTCCGCAGGCGTTCGCTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGcc                                         <  1:331474/90‑1 (MQ=255)
                                                 cGTTCGCTGCCCCTGTTCCACCAACGCATCCAGCCCCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAAcc                              <  1:112865/90‑1 (MQ=255)
                                                      gCTGCCACTGTCCCACCAACCCCTCCAGCCCCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGa                         <  1:712625/90‑1 (MQ=255)
                                                                       aaCGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCgccagc        >  1:540717/1‑90 (MQ=255)
                                                                             tCCAGCCCCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCtt  <  1:485624/90‑1 (MQ=255)
                                                                                   |                                                                                   
ACAGTTCATCGCCAAACGTTCCCGGCGGAAAGAGAACAACGTCCGCAGGCGTTCGCTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTT  >  NZ_CP009273/2888925‑2889091

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