Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,679,0060AC57.1% 9.0 / 16.1 14D164A (GAC→GCC) tusDNA replication terminus site‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (6/0);  new base C (0/8);  total (6/8)
Fisher's exact test for biased strand distribution p-value = 3.33e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCC  >  NZ_CP009273/1678917‑1679095
                                                                                         |                                                                                         
ccccACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGa                                                                                           >  2:412917/1‑90 (MQ=255)
                  ggAAGGGGTGCATCGTCTTTTGCGGGGGCTGACCCCCTTTATTCCTTCCCCCCCGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCtt                                                                         <  2:769245/89‑1 (MQ=255)
                        ggTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTttggttg                                                                   >  2:692131/1‑90 (MQ=255)
                         gTGCATCGTCTTTTGCGGGGGCTGACCCCCTTTAATGCTTCCCGCCCGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTttggttgg                                                                  <  2:789294/90‑1 (MQ=255)
                         gTGCATCGTCTTTTGCCGGGGCTGACCCCCCTTAATGCTTCCCGCACGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTttggttgg                                                                  <  2:88421/90‑1 (MQ=255)
                                gTCTTTTGCGGGGGCTGACCCCCCTTAATGCTTCCCGCACGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCtaata                                                           <  2:32186/90‑1 (MQ=255)
                                             gCTGACCCCCTTTATTGTTTCCCGCCCGCCCACCGTTCTCCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTaa                                              <  1:23417/90‑1 (MQ=255)
                                                               ttACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtgatga                            >  2:318813/1‑90 (MQ=255)
                                                               ttACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtgatga                            >  2:486872/1‑90 (MQ=255)
                                                               ttACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtgatga                            >  2:534444/1‑90 (MQ=255)
                                                                ccccGCACGCTCCCCGTTCTGCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGaa                           <  2:502053/88‑1 (MQ=255)
                                                                      aCGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTg                     >  1:257440/1‑90 (MQ=255)
                                                                       cGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTgg                    >  2:794768/1‑90 (MQ=255)
                                                                        ccTCACCTTTCTCCACCCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGc                   <  1:297491/89‑1 (MQ=255)
                                                                           ccccGTTCGCCACCCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGcaca                <  1:1864/88‑1 (MQ=255)
                                                                                         aCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGcc  >  1:304868/1‑90 (MQ=255)
                                                                                         |                                                                                         
CCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCC  >  NZ_CP009273/1678917‑1679095

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