Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,634,3410CT56.2% ‑3.3 / 13.3 16F40F (TTC→TTTdtpBdipeptide/tripeptide permease DtpB
Reads supporting (aligned to +/- strand):  ref base C (2/5);  new base T (9/0);  total (11/5)
Fisher's exact test for biased strand distribution p-value = 4.81e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTC  >  NZ_CP009273/3634252‑3634419
                                                                                         |                                                                              
cagcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTttc                                                                                <  1:285792/90‑1 (MQ=255)
cagcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTttc                                                                                <  1:42431/90‑1 (MQ=255)
  gcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTtttttt                                                                              >  1:466578/1‑90 (MQ=255)
                                              aGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTCTTTTAAAAACTTTGGTTTTCGCCAGAGCAGGCGTTTTTCACtt                                  >  2:140259/1‑90 (MQ=255)
                                               gCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTAAAAAGCTTTGGTTTTCGCAAAAGCAGGGGTTTTTTACttt                                 >  1:502561/1‑90 (MQ=255)
                                                cGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTTAACAGCTTGGATTTTCGCAAGAGCAGGCTTTTTTCACcttt                                >  1:170258/1‑90 (MQ=255)
                                                cGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTTGTTTAAAAACTTTGGTTTTCCCCAAAGAAGGCGTTTTTTACtttt                                >  2:275650/1‑90 (MQ=255)
                                                      ggCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCt                          <  1:183334/90‑1 (MQ=255)
                                                          actacGGCGTGCAGGGCGTACTGGCGGGTTTTTTTGTTTAAAAGCTTTGGTTTTCCCCAGAGCCGGGTTTTTTTTCTTTTTGGGGttttt                      >  2:128189/1‑89 (MQ=255)
                                                           ctacGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTtgc                     <  2:502561/90‑1 (MQ=255)
                                                                gCGTGCAGGGCGTACTGGCGGGTTTTTTTCTTTAACAGCTTGGGTTTTCCCAAAAGCAGGGTTTTTTTCCTTTTTGGGCTTTTTCTgcgc                >  2:67835/1‑90 (MQ=255)
                                                                 cGTGCAGGGCGTACTGGCGGGTTTTTTCGTTAAACAACTTGGATTCTCGCCAGAGCAGGCTTTTGTCACTTTTTGTGCTTTTTCTGCGCt               >  1:460090/1‑90 (MQ=255)
                                                                        ggCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTAt        <  1:408179/90‑1 (MQ=255)
                                                                         gCGTACTGGCGGTTTTTTTCGTTTAACAGCTTGGATTTTCGCAAGAGCAGGCTTTTTTCACTTTTTGTGCTTTTTCTTCGCTGGGCTATg       >  2:12455/1‑90 (MQ=255)
                                                                            tACTGGCGGTTTTTTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTTTTACTTTTTGTGCTTTTTCTTCGCTGGTTTTTGGcc    >  2:121381/1‑90 (MQ=255)
                                                                            tACTGGCGGTTTTTTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTTTTACTTTTTGGGCTTTTTCTGCGCTGGGCTATGGcc    >  2:142188/1‑90 (MQ=255)
                                                                            tACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTTTCACTTTTGGGGCTTTTTCTGCGCTGGTCTATGGcc    >  1:260093/1‑90 (MQ=255)
                                                                            tACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGcc    <  2:466578/90‑1 (MQ=255)
                                                                              cTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTc  <  1:477179/90‑1 (MQ=255)
                                                                                         |                                                                              
CAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTC  >  NZ_CP009273/3634252‑3634419

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