Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273970,3790CT58.3% 10.4 / 25.5 24S202F (TCC→TTC) mukFchromosome partition protein MukF
Reads supporting (aligned to +/- strand):  ref base C (1/9);  new base T (14/0);  total (15/9)
Fisher's exact test for biased strand distribution p-value = 7.65e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-05
Rejected as consensus: Frequency below/above cutoff threshold.

CTGACGCAACGTCTGATGGACGAACAGCAGCAGCAGGTGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTGC  >  NZ_CP009273/970291‑970462
                                                                                        |                                                                                    
cTGACGCAACGTCTGATGGACGAACAGCAGCAGCAGGTGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTcc                                                                                     <  1:727763/90‑1 (MQ=255)
        aCGTCTGATGGACGAACAGCAGCAGCAGGTGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGGGGTTTTCCGCtgtgg                                                                            >  2:1176595/1‑89 (MQ=255)
        aCGTCTGATGGACGAACAGCAGCAGCAGGTGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGGGGGTTTCCGCtgtgg                                                                            >  2:631390/1‑89 (MQ=255)
                   aCGAACAGCAGCAGCAGGTGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGGGGGGGTTTTCCGCTGGGGAATGTTGGttt                                                                 >  1:232690/1‑90 (MQ=255)
                    cGAACAGCAGCAGCAGGTGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGGGATTTTCAGCTGGGAAATGTTGCtttt                                                                >  2:615699/1‑89 (MQ=255)
                        cagcagcagcagGTGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGGATTGTTGCTTTCTGaa                                                            >  2:1028409/1‑90 (MQ=255)
                        cagcagcagcagGTGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGaa                                                            <  1:1019679/90‑1 (MQ=255)
                        cagcagcagcagGTGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGaa                                                            <  2:823224/90‑1 (MQ=255)
                            agcagcagGTGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGGGGGGGTTTTCAGCTGTGGATTTTTTGTTTTTGGaaatt                                                         >  1:318083/1‑88 (MQ=255)
                                gcagGTGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGGTTTTCCGCTGTGAATTGTTGCTTTTGGAAAATTcctg                                                    >  2:904673/1‑88 (MQ=255)
                                gcagGTGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTTCAGCTGTGAATTGTTGCTTTTGGAAAAtttccg                                                    >  2:153586/1‑86 (MQ=255)
                                       aaGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTg                                             <  2:318083/90‑1 (MQ=255)
                                        aGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGGTTTTCAGCTGTGGATTGTTGCTTTCGGAAACTTCCGGGAcgccgc                                            >  1:186637/1‑90 (MQ=255)
                                        aGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGGTTTTCAGCTGTGAATTGTTTCTTTTTGAAACTTCCCGGACGCTGc                                            >  2:858996/1‑90 (MQ=255)
                                             gATATCGCCCAGTTGCTGAACAAAGACTGGCGGGGGGGGGTTTTCAGCCGTGAATTTTTTGTTTTTGAAAATTCCCGAACGCcgcgggga                                       >  1:249355/1‑86 (MQ=255)
                                                 tCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGc                                   <  2:186637/90‑1 (MQ=255)
                                                   gCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTTCAGCTGTGAATTTTTGCTTTTGGAAACTTCCGGAAAGCTGCGGGAATTTCAg                                 >  2:1215184/1‑90 (MQ=255)
                                                      cAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGAt                              <  2:650333/90‑1 (MQ=255)
                                                        gTTGCTGAACAAAGACTGGCGGGGGGCGGTTTTCAGCTGTGAATTTTTGGTTTTTGAAACTTCCGGGACGCGGCGTGGATTGGAGGATaa                            >  1:903382/1‑89 (MQ=255)
                                                                  aaaGACTGGCGGGCGGCGATTTTCCGCTGTGAATTTTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTTCAGGATAAGCTGGGAGCg                  >  2:725203/1‑90 (MQ=255)
                                                                        tGGCGGGCGGCGATTTTCCGCTGGGGATTTTTTCTTTTTGAAAAATCCGGGACGCCGCCGGAATTTGAGGGTACGCCGGGAGGGGCgggg            >  1:556856/1‑89 (MQ=255)
                                                                             ggcggcGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACaa       <  1:416798/90‑1 (MQ=255)
                                                                             ggcggcGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACaa       <  1:1083935/90‑1 (MQ=255)
                                                                                  cgATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTGc  <  1:858996/90‑1 (MQ=255)
                                                                                  cgATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTGc  <  2:46080/90‑1 (MQ=255)
                                                                                        |                                                                                    
CTGACGCAACGTCTGATGGACGAACAGCAGCAGCAGGTGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTGC  >  NZ_CP009273/970291‑970462

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