Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092735,1620TA45.5% 19.9 / 10.0 23intergenic (+142/‑72)thrC/yaaXthreonine synthase/DUF2502 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (10/2);  new base A (0/10);  total (10/13)
Fisher's exact test for biased strand distribution p-value = 1.04e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCA  >  NZ_CP009273/5027‑5296
                                                                                                                                         |                                                                                                                                      
ttCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaa                                                                                                                                      >  2:215488/1‑138 (MQ=255)
 tCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaac                                                                                                                                     >  1:5356/1‑138 (MQ=255)
             tCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTtct                                                                                                                         >  2:354637/1‑138 (MQ=255)
                     ggTTTGTTTTTAAGCAGCCGGGTTTTTTTATGAAAAAATGAGGAAAAAAAATGAGAGGAAAAAAGGAAAAATTCTCAATAAATGCGGTAACTGAGAGATTAGGATGGCGGAAAATAAAAACCGCCGTTCTCATCGAGTa                                                                                                                <  2:169900/139‑1 (MQ=255)
                     ggTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGt                                                                                                                 >  2:32531/1‑138 (MQ=255)
                                  gaaGCCCGCTTTTTTTAGGAAAAAATACGGGAAAAAACGGCCGGGAAAAAAGGAAAATTCCCAAACAAAGGCATTAAATTAGAAATCGGATTTGCGGAAAACAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATc                                                                                                   <  1:44651/137‑1 (MQ=255)
                                     gCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAc                                                                                                >  2:44651/1‑139 (MQ=255)
                                      cccGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                                                                               >  1:52098/1‑139 (MQ=255)
                                             tttttATGAAAAAATTATGGAAAAAAAGGAAAGGAAAAAAGGAAAAATTCTAAAAAAAGGCGTAAACTTAGAGATTAGGATTGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAAc                                                                                         <  1:469348/138‑1 (MQ=255)
                                             tttttATAAAAAAATTATGGAAAAAAATGAGAGGAAAAAAGGAAAAATTCAAAAAAAAGGCGGTAATTTAGAGATTAGGATGGGGGAGAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAAc                                                                                         <  1:398365/138‑1 (MQ=255)
                                                  ataaaaaaaTAATGGAAAAAAAGGAAAGAAAAAAAGGAAAAATTCAAAAAAAAGGTGGTAATTAAGAAATTAGGAGGGCGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCa                                                                                    <  2:611794/132‑1 (MQ=255)
                                                  aggaaaaaATTAGGGAAAAAAATGAAAGGAAAAAGGGAAAAATTCCAAAAAAGGGGGGAAACTAAGAATTGGGTATGGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCa                                                                                    <  1:462423/136‑1 (MQ=255)
                                                  aaaaaaaaaTTAGGAAAAAAAAAGACGGGAAAAAAGGAGAAACTCTAAAAAATGGCTGAAACATAGAAATTGGGGTTGGGGAGAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCa                                                                                    <  1:600859/132‑1 (MQ=255)
                                                                             ggaaaaaaGGAAAAATTCTAAAAAAATGCTGAAACTAAGAGTTTAGGATGGGGGAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAg                                                        <  2:52098/136‑1 (MQ=255)
                                                                                        gaAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt                                             >  1:290809/1‑139 (MQ=255)
                                                                                        aaaaTTCTAAAAAAAGGCGGTAACTAAGAGATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt                                             <  1:354637/138‑1 (MQ=255)
                                                                                                       aTGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACttt                               >  2:108572/1‑139 (MQ=255)
                                                                                                                    aataTTAGGATTGCGTAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGt                  <  1:215488/136‑1 (MQ=255)
                                                                                                                     agaTTAGGATGGGGGAGAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTc                 <  1:632121/139‑1 (MQ=255)
                                                                                                                     agaTAAGGTTTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTc                 <  1:253238/139‑1 (MQ=255)
                                                                                                                     aaaTTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTc                 <  2:5356/137‑1 (MQ=255)
                                                                                                                                  cGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATgg       >  2:194235/1‑137 (MQ=255)
                                                                                                                                       aaTAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagca  >  2:501435/1‑137 (MQ=255)
                                                                                                                                         |                                                                                                                                      
TTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCA  >  NZ_CP009273/5027‑5296

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