Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,258,4430AG53.8% 8.8 / 16.3 13N368D (AAT→GAT) setBsugar efflux transporter SetB
Reads supporting (aligned to +/- strand):  ref base A (1/5);  new base G (7/0);  total (8/5)
Fisher's exact test for biased strand distribution p-value = 4.66e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CCACGCTCTATACCAACACTTCGCGCGTGGGCTGGATCATCGCAGGATCAGTGGCGGGCATCGTCGCCGAGATCTGGAATTATCACGCTGTGTTCTGGTTTGCGATGGTGATGATTATCGCCACTCTGTTTTGCTT  >  NZ_CP009273/2258366‑2258501
                                                                             |                                                          
ccacGCTCTATACCAACACTTCGCGCGTGGGCTGGATCATCGCAGGATCAGTGGCGGGCATCGTCGCCGAGATCTGGAATTATCACGCtg                                                <  1:427408/90‑1 (MQ=255)
              aaCACTTCGCGCGTGGGCTGGATCATCGCAGGATCAGTGGCGGGCATCGTCGCCCAGATCTGGGATTATTACGCTGTGTTCTGGTTTTCg                                  >  1:289077/1‑90 (MQ=255)
               acacTTCGCGCGTGGGCTGGATCATCGCAGGATCAGTGGCGGGCCTCGTCCCCCGGATCTGGGATTTTTACGCCGTGGTTTTGGTTGCgg                                 >  2:166403/1‑89 (MQ=255)
               acacTTCGCGCGTGGGCTGGATCATCGCAGGATCAGTGGCGGGCATCGTCGCCGAGATCTGGAATTATCACGCTGTGTTCTGGTTTGCGa                                 <  2:287197/90‑1 (MQ=255)
               acacTTCGCGCGTGGGCTGGATCATCGCAGGATCAGTGGCGGGCATCGTCGCCCGGATCTGGGATTATTACGCCGTGTTTTGGGTTgggg                                 >  2:29291/1‑87 (MQ=255)
                 acTTCGCGCGTGGGCTGGATCATCGCAGGATCAGTGGCGGGCATCGTCGCCGGGGTCTGGGATTATTACGCTGGGTTTTGGTTTGCGGTg                               >  1:297679/1‑90 (MQ=255)
                 acTTCGCGCGTGGGCTGGATCATCGCAGGATCAGTGGCGGGCAACGTCGCCGGGGTTTTGGATTATTAAGCTGTGGTTTGGGTTTCGGTg                               >  2:76448/1‑90 (MQ=255)
                  cTTCGCGCGTGGGCTGGATCATCGCAGGATCAGTGGCGGGCATCGTCGCCGAGATCTGGAATTATCACGCTGTGTTCTGGTTTGCGATgg                              <  1:150967/90‑1 (MQ=255)
                   ttCGCGCGTGGGCTGGATCATCGCAGGATCAGTGGCGGGCATCGTCGCCGGGATCTGGGATTATTACGCTGTGTTTTGGGTTTCGATggg                             >  2:224756/1‑89 (MQ=255)
                                       tcGCAGGATCAGTGGCGGGCATCGTCGCCGAGATCTGGAATTATCACGc                                                  >  1:170849/1‑49 (MQ=255)
                                       tcGCAGGATCAGTGGCGGGCATCGTCGCCGAGATCTGGAATTATCACGc                                                  <  2:170849/49‑1 (MQ=255)
                                            ggATCAGTGGCGGGCATCGGCGCCGGGATCTGGGATTATTACGCTGTGTTTTGGTTTTTGGTGGGGGTGATTTTTGCCCCCCTTTTTTgt    >  2:283766/1‑89 (MQ=255)
                                              aTCAGTGGCGGGCATCGTCGCCGAGATCTGGAATTATCACGCTGTGTTCTGGTTTGCGATGGTGATGATTATCGCCACTCTGTTTTGCtt  <  1:283766/90‑1 (MQ=255)
                                                                             |                                                          
CCACGCTCTATACCAACACTTCGCGCGTGGGCTGGATCATCGCAGGATCAGTGGCGGGCATCGTCGCCGAGATCTGGAATTATCACGCTGTGTTCTGGTTTGCGATGGTGATGATTATCGCCACTCTGTTTTGCTT  >  NZ_CP009273/2258366‑2258501

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