Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,976,0000TG64.3% 5.3 / 11.3 14intergenic (‑132/+35)otsB/araHtrehalose‑phosphatase/arabinose ABC transporter permease AraH
Reads supporting (aligned to +/- strand):  ref base T (0/5);  new base G (9/0);  total (9/5)
Fisher's exact test for biased strand distribution p-value = 5.00e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.61e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGC  >  NZ_CP009273/1975924‑1976089
                                                                              |                                                                                         
aTGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTTGGGTTTGGGGGGaaaaaagct                                                                               >  1:492642/1‑87 (MQ=255)
         taGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTTTGGGGGGGAAAAAAAGTAAAatgttt                                                                     >  1:736424/1‑89 (MQ=255)
         taGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGTTTGGGGGGGGAAAAAACTAAAatgttt                                                                     >  1:352771/1‑89 (MQ=255)
               gtgtGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGTGATGGGGGGGGAAAAAGACAAATTTGTGCAGcaa                                                               >  1:676222/1‑90 (MQ=255)
               gtgtGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGTATTGGGGGGAAAAAAACCAAAATTTTGTAGaaa                                                               >  1:284871/1‑90 (MQ=255)
               gtgtGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGTATTGGGGGGAAAAAAACTAAAATTTTGTAGaaa                                                               >  1:215928/1‑90 (MQ=255)
               gtgtGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGGATGGGGGGGAAAAAAGCTAAATTGTTTTAGaaa                                                               >  1:339471/1‑90 (MQ=255)
               gtgtGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGGGGGGTTGGGGGGGAAAAAAACAAAAATTTTGTAGaaa                                                               >  1:775367/1‑90 (MQ=255)
                tgtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGGGGGGAAAAAAGCCAAAATTTTGGAGaaaa                                                              >  1:133629/1‑90 (MQ=255)
                 gtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGaaaaa                                                             <  2:696000/90‑1 (MQ=255)
                 gtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGaaaaa                                                             <  2:628680/90‑1 (MQ=255)
                 gtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGaaaaa                                                             <  2:47602/90‑1 (MQ=255)
                 gtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGaaaaa                                                             <  1:368074/90‑1 (MQ=255)
                                                     aCCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAAc                         <  2:489096/90‑1 (MQ=255)
                                                                              tggGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCcagc  <  2:352771/90‑1 (MQ=255)
                                                                              tggGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCcagc  <  2:588806/90‑1 (MQ=255)
                                                                              tggGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCcagc  <  1:179402/90‑1 (MQ=255)
                                                                              |                                                                                         
ATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGC  >  NZ_CP009273/1975924‑1976089

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