Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,594,1670CT57.1% ‑1.8 / 22.4 21S64F (TCT→TTT) bglJDNA‑binding transcriptional activator BglJ
Reads supporting (aligned to +/- strand):  ref base C (0/9);  new base T (12/0);  total (12/9)
Fisher's exact test for biased strand distribution p-value = 3.40e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.26e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATT  >  NZ_CP009273/4594088‑4594252
                                                                               |                                                                                     
cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATgaga                                                                             >  2:155078/1‑90 (MQ=255)
          ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCATTAGAAGTGAGCGCa                                                                   >  1:27089/1‑90 (MQ=255)
               aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCcgaggg                                                              >  2:137061/1‑88 (MQ=255)
                 tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGGGCAGAGAggg                                                            >  1:106894/1‑90 (MQ=255)
                 tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCAGAGAggg                                                            >  1:432233/1‑90 (MQ=255)
                 tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCAGAGAggg                                                            >  2:468936/1‑90 (MQ=255)
                 tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCCTGAGAAGATAGGGCagggaggg                                                            >  2:6298/1‑90 (MQ=255)
                  cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGa                                                           <  1:155078/90‑1 (MQ=255)
                     gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGGGGAGGGAGGGCAGAGggggggtt                                                        >  1:388949/1‑86 (MQ=255)
                     gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGGTTa                                                        >  2:38326/1‑90 (MQ=255)
                           tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGc                                                  <  1:137061/90‑1 (MQ=255)
                                   tGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGGGGGGTTATCTTGCctgactga                                          >  2:464008/1‑90 (MQ=255)
                                     ccATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGGGGCGGGGGGGGTTTTCTTTCCTGATTGaaa                                        >  1:82417/1‑89 (MQ=255)
                                      cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACt                                       <  1:245402/90‑1 (MQ=255)
                                      cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTCCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACt                                       <  1:353600/90‑1 (MQ=255)
                                                   tcGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGGGAAGGGGGCGCAGGGGGGGGGTTTTTTTCCTGACTGAACTGGGGGTTTAGttt                          >  1:236908/1‑90 (MQ=255)
                                                             cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCTATTAAGTTTCCGCGTAccc                <  1:407492/90‑1 (MQ=255)
                                                             cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCTATTAAGTTTCCGCGTAccc                <  1:386581/90‑1 (MQ=255)
                                                             cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAcac                <  2:388949/90‑3 (MQ=255)
                                                                tCATTTTTTCTTTTTCTGCCATGAGAAGTTAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCCCGTACCCGGc             <  1:374549/90‑1 (MQ=255)
                                                                tCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGc             <  2:82417/90‑1 (MQ=255)
                                                                           ttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTAtt  <  1:163697/90‑1 (MQ=255)
                                                                               |                                                                                     
CTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATT  >  NZ_CP009273/4594088‑4594252

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