Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,380,3410GA55.0% 1.5 / 24.3 20intergenic (+31/‑72)yfbL/yfbMM28 family metallopeptidase/YfbM family protein
Reads supporting (aligned to +/- strand):  ref base G (1/8);  new base A (11/0);  total (12/8)
Fisher's exact test for biased strand distribution p-value = 7.14e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.80e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGC  >  NZ_CP009273/2380258‑2380427
                                                                                   |                                                                                       
aaaTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACAt                                                                                   <  2:1115577/90‑1 (MQ=255)
  aTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGc                                                                                 <  2:1133158/90‑1 (MQ=255)
  aTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGc                                                                                 <  2:948416/90‑1 (MQ=255)
              gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCACaatttttttt                                                                     >  1:1505872/1‑90 (MQ=255)
               tGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACatttttttttt                                                                    >  1:641412/1‑89 (MQ=255)
                 acaggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTccc                                                                  <  1:1249905/87‑1 (MQ=255)
                    ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACATTTTTTTTTTCCTcc                                                               >  1:465528/1‑89 (MQ=255)
                      tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                             <  1:1088514/90‑1 (MQ=255)
                      tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                             <  2:572957/90‑1 (MQ=255)
                               tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGGAAAAACGGCATGCACATTTTATTTTTCCTCCTCTTGTTTTag                                                    >  1:93444/1‑89 (MQ=255)
                               tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACATTTTTTTTTTTCTCATCtttttttta                                                   >  1:161136/1‑87 (MQ=255)
                               tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACAATTTATTTTTCCTCCTCTGtttttaa                                                   >  1:182097/1‑87 (MQ=255)
                               tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACAATTTTTTTTTCCTCCTCtttttttaa                                                   >  2:255106/1‑87 (MQ=255)
                               tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACATTTTATTTTTCCTCCTCTGtttttaa                                                   >  1:562353/1‑87 (MQ=255)
                               tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGAAAAGAAAAATTTTTTTTTTCTCATCTGttttttt                                                   >  2:723996/1‑87 (MQ=255)
                                          acaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCACATTTTTTTTTTCCTCATCTGTTTTTAAATGCGaaaaaa                                        >  1:485571/1‑87 (MQ=255)
                                             gTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCCTCTTTTTTTAAATGCGAAAAATGaa                                     >  1:1208035/1‑90 (MQ=255)
                                                    ataaAAATTTCTCGAGGCGGATCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAAc                              >  1:615948/1‑90 (MQ=255)
                                                                              aaaaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTg    <  1:629494/90‑1 (MQ=255)
                                                                                aaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGc  <  1:723996/90‑1 (MQ=255)
                                                                                   |                                                                                       
AAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGC  >  NZ_CP009273/2380258‑2380427

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