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

TGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAA  >  NZ_CP009273/2380280‑2380451
                                                                                  |                                                                                          
tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                                                     <  2:141720/90‑1 (MQ=255)
tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                                                     <  1:522876/90‑1 (MQ=255)
tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                                                     <  2:114973/90‑1 (MQ=255)
tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                                                     <  2:806919/90‑1 (MQ=255)
         tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAACGACAAGAAAATTTTATTTTTCCTTTTCTG‑TTTTTAg                                                                           >  1:512627/1‑90 (MQ=255)
         tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACAAGAACATTTTATTTTTTCTCCTCTG‑tttttaa                                                                           >  1:1031404/1‑87 (MQ=255)
         tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAGCACATTTTTTTTTTCCTCATCtt‑tttttaa                                                                           >  2:1084658/1‑87 (MQ=255)
         tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGGACATTTTTTTTTTCCTCATCTT‑TTTTTAg                                                                           >  2:700585/1‑90 (MQ=255)
         tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCACAATTTTTTTTTTCTCATCTG‑TTTTTAg                                                                           >  2:835294/1‑90 (MQ=255)
         tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAACAAAAGAAAATTTTATTTTTTCTCAGCTG‑GTTTTAg                                                                           >  1:979837/1‑90 (MQ=255)
         tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCATCtt‑tttttaa                                                                           >  2:432432/1‑87 (MQ=255)
         tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCACAATTTATTTTTCCTCATCTG‑tttttaa                                                                           >  2:446295/1‑87 (MQ=255)
         tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACATGCACATTTTATTTTTCCTCCtttt‑tttttaa                                                                           >  2:483802/1‑87 (MQ=255)
         tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACAATTTTTTTTTTCTCCTCTG‑tttttaa                                                                           >  1:314224/1‑87 (MQ=255)
           tGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCATCTG‑TTTTTAGAt                                                                         >  2:106107/1‑90 (MQ=255)
                  caacaGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAAACAAGATCAATTTTTTTTTTCTCCTCTGTTTTTAAAATGGGaaa                                                                   >  1:633039/1‑90 (MQ=255)
                         aaaTAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACATGCACATTTTTTTTTTCCCCCTCTT‑TTTTTAAAAGCGGAAAAAGaaaa                                                           >  2:240837/1‑89 (MQ=255)
                            taataaAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACATTTTTTTTTTCCTCATCTT‑TTTTTAAATGCGAAAAAAGaaaatta                                                        >  1:69546/1‑86 (MQ=255)
                              ataaAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGCACATTTTTTTTTTCCCCCTCTG‑TTTTTAGATGCGAAAAATGAAAATCAAc                                                      >  2:281878/1‑90 (MQ=255)
                              ataaAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAGCACATTTTTTTTTTCCTCATCTT‑TTTTTTGATGGGAAAAAAGAAAATTAAc                                                      >  1:132677/1‑90 (MQ=255)
                                                      aGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGAt                              <  1:206830/90‑1 (MQ=255)
                                                                  aTGCACATTTTATTTTCCCTCATCTG‑TTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGc                  <  1:893241/90‑1 (MQ=255)
                                                                     cacaTTTTATTTTCCCTCATCTG‑TTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGa               >  1:1072985/1‑90 (MQ=255)
                                                                              ttttCCCTCATCTG‑TTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTc      <  2:447507/90‑1 (MQ=255)
                                                                                 tCCCTCATCTG‑TTTTAAGATGCGAAACATGAAATTCAACGGAGGGAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa   >  1:66968/1‑90 (MQ=255)
                                                                                 tCCCTCATCTG‑TTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa   >  2:180008/1‑90 (MQ=255)
                                                                                 tCCCTCATCTG‑TTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa   >  1:421853/1‑90 (MQ=255)
                                                                                 tCCCTCATCTG‑TTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa   >  2:958884/1‑90 (MQ=255)
                                                                                  cccTCATCTG‑TTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGGTTGGCTATTTTGCAGAAATTGACTCAGaa  <  2:282895/90‑1 (MQ=255)
                                                                                  |                                                                                          
TGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAA  >  NZ_CP009273/2380280‑2380451

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: ‑