Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,594,1670CT70.0% 4.6 / 10.8 20S64F (TCT→TTT) bglJDNA‑binding transcriptional activator BglJ
Reads supporting (aligned to +/- strand):  ref base C (0/6);  new base T (14/0);  total (14/6)
Fisher's exact test for biased strand distribution p-value = 2.58e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.25e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGG  >  NZ_CP009273/4594083‑4594256
                                                                                    |                                                                                         
ataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCa                                                                                      >  1:316908/1‑90 (MQ=255)
     cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCCTgaga                                                                                 >  1:261124/1‑90 (MQ=255)
       ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAg                                                                               >  1:124768/1‑90 (MQ=255)
              tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTGCATTAGAAGTGAgggg                                                                        >  2:27187/1‑87 (MQ=255)
              tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcgc                                                                        <  2:408807/90‑1 (MQ=255)
               ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTg                                                                                         >  1:90017/1‑72 (MQ=255)
               ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTg                                                                                         <  2:90017/72‑1 (MQ=255)
               ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTAATTTTTTTTTTTTTTGTCATGAGAAGGGAgcggc                                                                       >  1:297378/1‑88 (MQ=255)
                    aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTTCCATTAGAAGGGAGGGCagagag                                                                  >  1:239331/1‑90 (MQ=255)
                    aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTGCCTGAGAAGTGAGGGCcgagag                                                                  >  1:137377/1‑90 (MQ=255)
                    aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCagagag                                                                  >  1:83007/1‑90 (MQ=255)
                    aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGcc                                                                                       >  1:334089/1‑69 (MQ=255)
                    aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGcc                                                                                       <  2:334089/69‑1 (MQ=255)
                      tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCCTGGGAAGGGAGCGCagggaggg                                                                >  2:270925/1‑90 (MQ=255)
                       cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAgggg                                                               >  2:384773/1‑89 (MQ=255)
                       cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGa                                                               <  1:133578/90‑1 (MQ=255)
                                     ggATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGAGAGGGGTTATCTTTCctgac                                                 >  2:390861/1‑90 (MQ=255)
                                           cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACt                                           <  1:392975/90‑1 (MQ=255)
                                                        tcGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGAGAGGGGTTTTCTTTCCTGACTGAACTGGCGGTTTAGttt                              >  1:242419/1‑90 (MQ=255)
                                                                   ggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCg                   <  2:83007/90‑1 (MQ=255)
                                                                            ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGt          <  1:202899/90‑1 (MQ=255)
                                                                                ttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTAtt      <  2:137377/90‑1 (MQ=255)
                                                                                    ctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCgg  <  1:270925/90‑1 (MQ=255)
                                                                                    ctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCgg  <  2:124768/90‑1 (MQ=255)
                                                                                    |                                                                                         
ATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGG  >  NZ_CP009273/4594083‑4594256

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