Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092735,1590GA50.0% 5.4 / 10.8 18intergenic (+139/‑75)thrC/yaaXthreonine synthase/DUF2502 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (9/0);  new base A (0/9);  total (9/9)
Fisher's exact test for biased strand distribution p-value = 4.11e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGA  >  NZ_CP009273/5079‑5243
                                                                                |                                                                                    
gaagaaATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTTTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaacaac                                                                             >  1:283386/1‑90 (MQ=255)
gaagaaATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaacaac                                                                             >  2:703101/1‑90 (MQ=255)
   aaaaTTAGGGAAAAAAATGAGGGGAAAAAAGGAAAAATTCTCAATAAATGGGGTAACTGAAAGATTAGGTTGGCGGAAAATAACAAccgc                                                                          <  1:555388/89‑1 (MQ=255)
       ttAGGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTCAATAAAGGCGGTAACTTAGAGATTAGGATTGCGGAAAATAACAACCGCCGtt                                                                      <  1:2277/90‑1 (MQ=255)
              gaAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCg                                                               >  1:43229/1‑90 (MQ=255)
                  aaTGACAGGAAAAAAGGAAAAATTCTAAAAAAAGGGGGTAACTTAGAGATTAGGATTGGGGAAAATAACAACCGCCGTTCTCATCGAGTa                                                           <  1:915001/90‑1 (MQ=255)
                           gAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGa                                                  >  1:393322/1‑90 (MQ=255)
                              aaaGAAAAATTTCTAAAAAAAGGGGGAAATTAAGAGATTAGGATGGGGGAAAAAACCAACCGCCGTTCTCATCGAGTAATCTCCGGatat                                               <  1:703101/90‑1 (MQ=255)
                               aaGGAGAAATTCTAAATAAATGCGGTAACTTAGAGATTAGGTTTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATc                                              <  1:407647/90‑1 (MQ=255)
                               aaGGAAAAATCCTAAAAAAATGGGTTAACTTAGAAATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATc                                              <  1:454894/90‑1 (MQ=255)
                                        ttCTAAAAAAAGGCGTTAATTAAGAAATTAGGATGGCGAAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATaa                                     <  2:14036/90‑1 (MQ=255)
                                           tcAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACgg                                  >  1:27149/1‑90 (MQ=255)
                                               aaaaTGCGGTAATATAGAAATTAGGATGGGGAAAAAAAAAACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa                              <  2:723395/89‑1 (MQ=255)
                                                      ggTAATTAAGAGATAAGGATGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                       <  2:27149/90‑1 (MQ=255)
                                                                   ttAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgt          >  1:847074/1‑90 (MQ=255)
                                                                   ttAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgt          >  1:696334/1‑90 (MQ=255)
                                                                     aGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa        >  1:555850/1‑90 (MQ=255)
                                                                           gCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGa  >  1:221927/1‑90 (MQ=255)
                                                                                |                                                                                    
GAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGA  >  NZ_CP009273/5079‑5243

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