Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,380,3410GA59.1% ‑2.0 / 18.1 22intergenic (+31/‑72)yfbL/yfbMM28 family metallopeptidase/YfbM family protein
Reads supporting (aligned to +/- strand):  ref base G (2/7);  new base A (13/0);  total (15/7)
Fisher's exact test for biased strand distribution p-value = 2.11e-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.

ATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAAC  >  NZ_CP009273/2380260‑2380399
                                                                                 |                                                           
aTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGc                                                     <  2:512053/90‑1 (MQ=255)
  ggCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGACGCGGGTCCAGAAAAAGCGACATGcac                                                   >  2:467439/1‑90 (MQ=255)
     tCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACatt                                                <  1:175409/90‑1 (MQ=255)
      cAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACaatt                                               >  2:74905/1‑90 (MQ=255)
                    tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                 <  1:459670/90‑1 (MQ=255)
                    tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                 <  2:161267/90‑1 (MQ=255)
                    tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                 <  2:331464/90‑1 (MQ=255)
                    ggtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                 <  1:63877/89‑1 (MQ=255)
                      ttATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAAAACATTTTTTTTTTTCTCATCtt                               >  1:356170/1‑89 (MQ=255)
                             tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAGCACATTTTTTTTTTCCTCCTCTGtttttaa                       >  2:81657/1‑87 (MQ=255)
                             tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACATTTTTTTTTTCCTCCTCtttttttaa                       >  2:562516/1‑87 (MQ=255)
                             tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACAATTTTTTTTTTCTCCTCtttttttaa                       >  2:157004/1‑87 (MQ=255)
                             tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAGCACAATTTATTTTTCCTCCTCTGtttttaa                       >  1:594200/1‑87 (MQ=255)
                                   ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGa                 <  2:828572/90‑1 (MQ=255)
                                        acaGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACAAGGAAATTTTATTTTTTCTCATCCGCTTTTAGATGAGaaaaaa            >  1:384091/1‑87 (MQ=255)
                                        acaGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACAATTTTTTTTTCCTCATCTTTTTTTAAATGCGgaaaaa            >  2:242256/1‑87 (MQ=255)
                                          aGTAAATAATAAAAATTTCTCGAGGCGGGGGCAGAAAAAAAGACAAGCACAATTTTTTTTTTCTCCTCTTTTTTTAGATGGGAAACaaaa          >  1:637946/1‑87 (MQ=255)
                                          aGTAAATAATAAAAATTTCTCGAGGCGGGGGCAGAAAAAAAAACAAGCAAATTTTTTTTTTTCCCCTCTGTTTTTAGATGCGGAAAATGa          >  1:332376/1‑90 (MQ=255)
                                          aGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACAATTTTTTTTTCCTCATTTGTTTTTAAAATGGGaaaaaaa           >  2:652876/1‑86 (MQ=255)
                                                taataaAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACAAGCACATTTTATTTTTCCTCATCTGTTTTAAGATGCGAAACATGAAATTCa    >  1:57638/1‑90 (MQ=255)
                                                  ataaAAATTTCTCGAGGCGGGGCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCCTCTGTTTTTAGATGCGAAAAATGAAAATCaaa  >  2:372841/1‑89 (MQ=255)
                                                  ataaAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACATTTTTTTTTTCCTCCTCTTTTTTTAGATGCGGAAAATGAAAATTAAc  >  2:692851/1‑90 (MQ=255)
                                                                                 |                                                           
ATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAAC  >  NZ_CP009273/2380260‑2380399

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