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

TTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAA  >  NZ_CP009273/5074‑5239
                                                                                      |                                                                                
tttATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATa                                                                              >  2:192954/1‑90 (MQ=255)
    tGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaacaa                                                                          >  1:65714/1‑90 (MQ=255)
     gaagaaATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaacaac                                                                         >  2:111724/1‑90 (MQ=255)
                  agaAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATc                                                            >  2:471454/1‑90 (MQ=255)
                      aaaTGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGt                                                        >  1:512105/1‑90 (MQ=255)
                             agggaaaaagggaaaaTTTCCCAATAAATGGGTTATTTTAGAGATTGGTTTTGGGGAAAATAACACCCGCCGTTCTCATCGAGTAAtctc                                                  <  2:237573/90‑1 (MQ=255)
                                   aaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGatat                                           >  1:320075/1‑90 (MQ=255)
                                     agggAAAAATTCTCAAAAATTGTGGTTATTTGAAGATTGGTTTTGGGGAAACCACCAACCGCCGTTCTCATCGAGTAATCTCCGGATATc                                          <  2:33134/88‑1 (MQ=255)
                                                 tcAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACgg                              >  1:83391/1‑90 (MQ=255)
                                                 tcAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACgg                              >  1:58176/1‑90 (MQ=255)
                                                   aaTAAACGCGGTAACTAAGAGATTTGGATTGCGGAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGc                            <  2:166395/90‑1 (MQ=255)
                                                     tAATTGGGGTAATTTAGATATTGGGTTTGCGGAAAAAAAAACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa                          <  1:58688/90‑1 (MQ=255)
                                                     aaaaTGGGGTAACTAAGAAATTAGGTTGGGGGAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa                          <  2:512105/89‑1 (MQ=255)
                                                     aaaaTGCGGTATCTTAAAGTTTAGGATGGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa                          <  2:202002/89‑1 (MQ=255)
                                                            ggTAATTAAGAGATTAGGTTGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                   <  2:65714/90‑1 (MQ=255)
                                                                 tttAGAGTTTGGGTTGGCGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAg              <  2:527320/89‑1 (MQ=255)
                                                                   tAGAAATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTa            <  2:320075/90‑1 (MQ=255)
                                                                   aagagATTAGGATGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTa            <  2:123398/89‑1 (MQ=255)
                                                                          tAGGATGGGGGAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgtg     <  1:160255/90‑1 (MQ=255)
                                                                             gATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaa  <  2:97920/90‑1 (MQ=255)
                                                                                      |                                                                                
TTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAA  >  NZ_CP009273/5074‑5239

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