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

TGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATG  >  NZ_CP009273/5089‑5245
                                                                      |                                                                                      
tGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTAAGAGATTAGGATTGCGGAGAATAACAACCGCCGTtctc                                                                     >  1:671548/1‑90 (MQ=255)
 ggAAAAAAAAGACAGGAAAAAAGAAAAAACTCCAAAAAATTGCGGTAACTTAGAGTTTAGGATGGCGGAAAAAAACAACCGCCGTTCTCa                                                                    <  2:480450/90‑1 (MQ=255)
   aaaaaaaTGAAGGGAAAAAAGAAAAAATTCTAAAAAAATGCGGTAACTTAGAGATTAGGTTGGGGGAGAAAAACAACCGCCGTTCTCATc                                                                  <  2:58929/88‑1 (MQ=255)
    gaAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCg                                                                 >  1:759674/1‑90 (MQ=255)
    aaaaaaTGACAGGAAAAAAGGAAAAATTCTCAAAAAATGCGGTAACTTAGAGATTAGGTTTGGGGAAAATAACAACCGCCGTTCTCATCg                                                                 <  1:73253/89‑1 (MQ=255)
          tGAAAGGAAAAAAGAAAAAATCCTAAAAAAATGGGGTAACTTAGAAATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAAt                                                           <  2:671548/90‑1 (MQ=255)
            aCAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAAtct                                                         >  2:273388/1‑90 (MQ=255)
                     aaGGAAAAATTCAAAAAAAGGGGGGTAATTAAGATTTTAGGTTGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATc                                                <  2:266071/90‑1 (MQ=255)
                         agaAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                            >  2:67853/1‑90 (MQ=255)
                          gaAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcac                                           >  1:675777/1‑88 (MQ=255)
                               tctcAATAAATGCGGTAACTTAGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAAc                                      <  1:134399/90‑1 (MQ=255)
                               tctcAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAAc                                      >  1:693121/1‑90 (MQ=255)
                               tcccAAAAAACGCGTTAACTAAGAGTTTAGGCTGGGGGAAAAAAAACACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAAc                                      <  1:755444/87‑1 (MQ=255)
                                            ggTAACTTAGAGTTTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                         <  1:9889/90‑1 (MQ=255)
                                             gTAACTTAGAGATTAGGTTGGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaaa                        <  1:273388/90‑1 (MQ=255)
                                                     gagaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAc                >  1:307212/1‑90 (MQ=255)
                                                      agaTTAGGATGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc               <  1:527406/90‑1 (MQ=255)
                                                           aGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa          >  1:734662/1‑90 (MQ=255)
                                                                   ggAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg  <  1:652270/90‑1 (MQ=255)
                                                                      |                                                                                      
TGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATG  >  NZ_CP009273/5089‑5245

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