Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,976,0040AG57.1% 2.2 / 14.0 15intergenic (‑136/+31)otsB/araHtrehalose‑phosphatase/arabinose ABC transporter permease AraH
Reads supporting (aligned to +/- strand):  ref base A (2/4);  new base G (8/0);  total (10/5)
Fisher's exact test for biased strand distribution p-value = 1.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGG  >  NZ_CP009273/1975869‑1976140
                                                                                                                                       |                                                                                                                                         
ccGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGaaaa                                                                                                                                        <  1:24475/139‑1 (MQ=255)
 cGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGTATGGTGGGgaaa                                                                                                                                        >  2:487585/1‑138 (MQ=255)
              tctcTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTGTGGGGGGGAAAAAAGCCCAattatt                                                                                                                          >  2:147279/1‑139 (MQ=255)
                               ttGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCGGTGGGGGGGAAAAAAGACAAATTGTTGTAGAAAAAAacaaca                                                                                                         >  2:94643/1‑138 (MQ=255)
                               ttGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGTTTTGGGGGGAAAAAAACACAAATGTTGGAGcaaaaaaaaa                                                                                                            >  2:254660/1‑132 (MQ=255)
                                                tACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCttt                                                                                       <  1:173036/138‑1 (MQ=255)
                                                        tGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGTCTGTGGTTGGGAAAAAAGCCAAAATTTTGTAGAAAAAAAGAACAGACAGGGCGGTTCGCttttttttttt                                                                               >  2:193851/1‑133 (MQ=255)
                                                                       tgtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTTGCGATGTTGGGGAAAAAAGCGAAAATGTTGTAGAAAAAAACAACAAAACAAGGGGTTTGCTTTTTTCTTGTTAAGGTCGCAGAACa                                                                >  1:189094/1‑139 (MQ=255)
                                                                       tgtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGGGGGGAAAAAAGCTAAAATGTTGCGGAAAAAAACATAAGACAGTGGGGTTAGCTTTTTTGTTTTAAAGGTCGGAGAta                                                                 >  1:192091/1‑137 (MQ=255)
                                                                       tgtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGTATGGGGGGGAAAAAAGCAAAATTGTTGCAGAAAAAAAAAACAGACAGTGCGGTTCGCGTTTTTTTTGTAAAGGGCGAAGATCa                                                                >  1:20651/1‑139 (MQ=255)
                                                                         tgACGGGGGTCGCCATTTGCAACGAAGACCAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCAcc                                                              <  2:127015/139‑1 (MQ=255)
                                                                                 gTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTTTGGTGGGGAAAAAAGCTAAAATTTTGCAGAAAAAAACAACAGACAGGGCGTTTCGCTTTTTTTTTGTTAAGGGCGAAGAACACCGCGGcccg                                                      >  1:232402/1‑137 (MQ=255)
                                                                                               cGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCcgc                                        <  1:254660/139‑1 (MQ=255)
                                                                                                                   tGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGaa                     <  1:487585/138‑1 (MQ=255)
                                                                                                                                  gtggGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGc     <  2:191174/139‑1 (MQ=255)
                                                                                                                                     ggAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGACATATTAAGCAgg  >  2:297920/1‑139 (MQ=255)
                                                                                                                                     ggAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAgg  >  2:528832/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                         
CCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGG  >  NZ_CP009273/1975869‑1976140

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