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

GTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTA  >  NZ_CP009273/2380272‑2380429
                                                                      |                                                                                        
gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACattttttttt                                                                      >  1:364779/1‑90 (MQ=255)
   gATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTccc                                                                   <  1:120045/90‑1 (MQ=255)
   gATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTccc                                                                   <  2:1047358/90‑1 (MQ=255)
     tGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACAAGCACAATTTTTTTTTCCTc                                                                 >  1:886643/1‑90 (MQ=255)
      ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACAAGCACATTTTTTTTTtccccc                                                                >  1:856204/1‑87 (MQ=255)
      ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACAAGCACATTTTTTTTTTCCTCa                                                                >  1:217111/1‑90 (MQ=255)
      ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAGCGACATGGACATTTTTTTTTTCCTcc                                                                >  1:369024/1‑89 (MQ=255)
      ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAAAAAAAAAGATAAACACAATTTTTTTTTCCTcc                                                                >  1:1267879/1‑89 (MQ=255)
       gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAAAAAAAAAGACAAGCAAATTTTTTTTTTTCTCAt                                                               >  1:692769/1‑90 (MQ=255)
         gtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGGCATGCACATTTTTTTTTTCCTCctct                                                             >  1:516704/1‑90 (MQ=255)
          ttATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACATGCACATTTTTTTTTTTCTCctctt                                                            >  1:119690/1‑89 (MQ=255)
                 tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCTGAAAAAAAGACAAGCAAAATTTTTTTTTTCTCCTCTGttttttt                                                     >  2:28964/1‑87 (MQ=255)
                 tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACAATTTTTTTTTTCTCCTCtttttttaa                                                     >  2:442801/1‑87 (MQ=255)
                 tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACATTTTTTTTTTTCTCATTTGTTTTTAg                                                     >  1:966558/1‑90 (MQ=255)
                 tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGAAAAGAACATTTTTTTTTTTCTCCTCTGtttttat                                                     >  2:881828/1‑87 (MQ=255)
                 tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGAAAAACAAATTTTTTTTTTTCTCCTCtttttttat                                                     >  2:893369/1‑87 (MQ=255)
                         acaacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGaaa                                             <  2:119690/90‑1 (MQ=255)
                              aGTAAATAATAAAAATTTCTCGAGGCGGGCTCGCGAAAAACGAAAAGAACATTTTATTTTTCCCCCTCTGTTTTTAGATGCGAAACACGa                                        >  2:356206/1‑90 (MQ=255)
                              aGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGa                                        <  1:601527/90‑1 (MQ=255)
                                 aaaTAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAAACAAGCACAATTTTTTTTTCCTCCTCTGTTTTTAGAAGGGAAAAAAGaaaa                                     >  2:297752/1‑89 (MQ=255)
                                      ataaAAATTTCTCGAGGCGGGTCTAGAAAAAACGACATGCACATTTTATTTTTCCTCATCTGTTTTTAGATGCGAAAAATGAAATTTaaa                                >  1:509939/1‑89 (MQ=255)
                                            aTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACggagga                          <  2:1117443/90‑1 (MQ=255)
                                            aTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACggagga                          <  2:1267879/90‑1 (MQ=255)
                                                     ggCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAaaataaat                 <  2:49897/90‑1 (MQ=255)
                                                                    aaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCt   >  1:581382/1‑90 (MQ=255)
                                                                     aGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTa  >  1:208518/1‑90 (MQ=255)
                                                                      |                                                                                        
GTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTA  >  NZ_CP009273/2380272‑2380429

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