Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,007,6930TC53.3% 9.3 / 14.0 16P328P (CCT→CCCfliFflagellar basal‑body MS‑ring/collar protein FliF
Reads supporting (aligned to +/- strand):  ref base T (4/3);  new base C (0/8);  total (4/12)
Fisher's exact test for biased strand distribution p-value = 2.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAA  >  NZ_CP009273/2007610‑2007772
                                                                                   |                                                                               
aGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCa                                                                           <  2:845277/90‑1 (MQ=255)
            tATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAAc                                                               >  2:795700/1‑90 (MQ=255)
              tCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCg                                                             >  2:533099/1‑90 (MQ=255)
                  ggcggcGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCcag                                                         >  1:461064/1‑90 (MQ=255)
                   gcggcgTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCcagc                                                        >  1:90298/1‑90 (MQ=255)
                      gcgTCCGGGGGGCGTTGTCAAATAACCCGGCCCCTGCAAAAACCGCGCCAATCAGCACGCCCCCGGCAAATAAAAATAACCGCcagcagc                                                     <  1:885347/90‑1 (MQ=255)
                      gcgTACGGGGGGCGTTGTCGAAACACCCGGCACCTGCGAAAAACGCGCCAATCAGCACGCCCCCGGCAAATAAAAATAACCGCcagcagc                                                     <  1:725289/90‑1 (MQ=255)
                            cGGGGGCGTTGTCGAATCACCCGGCCCCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAAAAACCGCCAGCAGCAGGCGa                                               <  1:833940/90‑1 (MQ=255)
                              ggggCGTTGTCAAACCAACCGGCCCCTGCAAATAACGCGCCAATCAGCACCCCCCCGGCAAATAAAAATAACCGCCAGCAGCAGGCGAGc                                             <  1:533099/90‑1 (MQ=255)
                                  atgtGTCGAATCAACAGGCCGCTGCGAAAAGCGCCCAAATAAGCACGCTGCCGGCAAACAAAAAAAACCGCCAGCAGCAGGCGAGcacca                                         <  1:674050/87‑1 (MQ=255)
                                                   gccccTGCGAATAACGCGCCAATCAGCCCGCCCCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCcg                        <  1:752749/87‑1 (MQ=255)
                                                   gccccTGCGAATAACGCGCAAATCCGCACGCCCCCGGGAAACTAAAATATCCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCcg                        <  2:556703/87‑1 (MQ=255)
                                                           gAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGcac                <  2:72254/90‑1 (MQ=255)
                                                                aCGCGCCAACCACCACCCCCCCGGAAAATCAAAAAAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGc           <  1:165036/90‑1 (MQ=255)
                                                                    cccAATCAGCACCCCCCCGGAAAAAAAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGaa       <  1:589838/89‑1 (MQ=255)
                                                                         tCAGCACGCCCCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGaaa  <  2:778922/90‑1 (MQ=255)
                                                                                   |                                                                               
AGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAA  >  NZ_CP009273/2007610‑2007772

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