Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,380,3620CT53.6% ‑0.8 / 31.2 28intergenic (+52/‑51)yfbL/yfbMM28 family metallopeptidase/YfbM family protein
Reads supporting (aligned to +/- strand):  ref base C (6/7);  new base T (15/0);  total (21/7)
Fisher's exact test for biased strand distribution p-value = 1.45e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-06
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAA  >  NZ_CP009273/2380275‑2380451
                                                                                        |                                                                                         
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:1267879/90‑1 (MQ=255)
                                         aTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACggagga                                                <  2:1117443/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)
                                                                       cATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTg                   >  1:262067/1‑90 (MQ=255)
                                                                       cATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTg                   >  2:751846/1‑90 (MQ=255)
                                                                               ttttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAAtt           <  1:442801/90‑1 (MQ=255)
                                                                                   attttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACt       <  1:893369/90‑1 (MQ=255)
                                                                                       tCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa   >  2:1073466/1‑90 (MQ=255)
                                                                                       tCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa   >  2:375757/1‑90 (MQ=255)
                                                                                        cccTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGaa  <  1:438653/90‑1 (MQ=255)
                                                                                        cccTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGaa  <  2:288160/90‑1 (MQ=255)
                                                                                        cccTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGaa  <  1:986019/90‑1 (MQ=255)
                                                                                        |                                                                                         
GATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAA  >  NZ_CP009273/2380275‑2380451

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