Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,398,0510GT46.2% ‑0.3 / 14.1 13Q258K (CAA→AAA) abgAp‑aminobenzoyl‑glutamate hydrolase subunit AbgA
Reads supporting (aligned to +/- strand):  ref base G (1/6);  new base T (6/0);  total (7/6)
Fisher's exact test for biased strand distribution p-value = 4.66e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCTTCCTGCCTGCATAACGCCCACGTTTACTCTGGAAGCTCCTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATT  >  NZ_CP009273/1397968‑1398134
                                                                                   |                                                                                   
gCTTCCTGCCTGCATAACGCCCACGTTTACTCTGGAAGCTCCTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTTtggtgtg                                                                              >  2:373973/1‑89 (MQ=255)
 cTTCCTGCCTGCATAACGCCCACGTTTACTCTGGAAGCTCCTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTgccg                                                                             <  2:103919/91‑1 (MQ=255)
     ctgcctgcATAACGCCCACGTTTACTCTGGAAGCTCCTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCcaa                                                                         <  1:4571/91‑1 (MQ=255)
           gcATAACGCCCACGTTTACTCTGGAAGCTCCTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTTTTCTGTCGCCCACAAGGc                                                                   >  2:18549/1‑91 (MQ=255)
                     cACGTTTACTCTGGAAGCTCCTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTTTGCTTCCGCCAACAAAGGATTTTGAcgg                                                         >  2:507128/1‑89 (MQ=255)
                            aCTCTGGAAGCTCCTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTg                                                  <  2:224333/91‑1 (MQ=255)
                              tctGGAAGCTCCTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTTTGCTGCCGCCCACAAGGCATTGTGACCGTCTTCTGGt                                                >  2:451520/1‑91 (MQ=255)
                                 ggAAGCTCCTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTTTTTTGCCGCCAACAAAGGACTTGGGGCGTCTTCTGGtttt                                             >  2:44028/1‑91 (MQ=255)
                                 ggAAGCTCCTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGtttt                                             >  1:660782/1‑91 (MQ=255)
                                         cTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTTTGTTGCCGCCCACAAAGCATTGTGACCGTCTTTTTGTTTTTCGCCTGc                                     >  2:418318/1‑91 (MQ=255)
                                           tCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTTTTCTGCCGCCCACAAGGCATTGTGGCCGTCTTCTGGTTTTTTGCCCGCGt                                   >  2:66057/1‑91 (MQ=255)
                                             gCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGa                                 <  1:66057/91‑1 (MQ=255)
                                                     gggCGATTGCATGCAGTGCAAGAGTGGCTTTTGTTGCCGCCCCCAAAGGATTTTGACCGGCTTCTTTTTTTTTGCCTGCGTGGGgggcggg                         >  2:145126/1‑90 (MQ=255)
                                                       gCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTAc                       <  1:834212/91‑1 (MQ=255)
                                                       gCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTAc                       <  1:18549/91‑1 (MQ=255)
                                                                    gTGCAAGAGTGGCTTTTTTTGCCGCCCACAAAGCATTGTGACCGGCTTCTTGTTTTGTGCCCGCGTGGGCGGCGGGACCCGGGGAGTgggg          >  2:21877/1‑88 (MQ=255)
                                                                            gTGGCTTTTTTTGCCGCCCACAAAGCATTGTGACCGTCTTCTTGTTTTTTGCCTGCGTGAGCGGCGGGACCGGTGAAGGGGGCGTCaaaat  >  2:176780/1‑89 (MQ=255)
                                                                                   |                                                                                   
GCTTCCTGCCTGCATAACGCCCACGTTTACTCTGGAAGCTCCTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATT  >  NZ_CP009273/1397968‑1398134

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