Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,380,3570AT55.0% 8.1 / 22.2 20intergenic (+47/‑56)yfbL/yfbMM28 family metallopeptidase/YfbM family protein
Reads supporting (aligned to +/- strand):  ref base A (4/5);  new base T (11/0);  total (15/5)
Fisher's exact test for biased strand distribution p-value = 8.13e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATT  >  NZ_CP009273/2380275‑2380442
                                                                                  |                                                                                       
gATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTccc                                                                                  <  1:923485/90‑1 (MQ=255)
    gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACATGCACAATTTTTTTTTCCTcct                                                                              >  1:405067/1‑88 (MQ=255)
    gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACATTTTTTTTTTCCTcct                                                                              >  1:824881/1‑88 (MQ=255)
     tgtTATACCTTTGTTATCCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                                            <  2:1454876/90‑1 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCCGAAAAAAAGAAAAGCAAAATTTTTTTTTCCTCCTCTTTTTTTAg                                                                  >  2:503491/1‑90 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAGCACATTTTATTTTTCCTCCTCTGtttttaa                                                                  >  1:919985/1‑87 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGGGCTGCAAAATTTTTTTTTTCCCCTCTGtttttaa                                                                  >  1:1455589/1‑87 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCCTCTGtttttaa                                                                  >  1:717551/1‑87 (MQ=255)
                    ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCACATTTTTTTTTTCCTCATCTGTTTTTAAATGCGa                                                            >  1:120157/1‑90 (MQ=255)
                     taCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGaa                                                           <  1:1998839/90‑1 (MQ=255)
                         acaGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAAAAAAGCACATTTTTTTTTTTCTCATCTGTTTTTTAATGCGaaaaaa                                                       >  2:1998839/1‑87 (MQ=255)
                           aGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACAATTTTTTTTTCCTCCTCCTTTTTTAAAAGGGGaaaaaga                                                      >  1:1636643/1‑86 (MQ=255)
                                 taataaAAATTTCTCGAGGCGGGGCCCGAAAAAACGACAAGCACATTTTTTTTTTTCTCCTCTTTTTTTATATGGGGAACAAGaaaataa                                               >  2:550314/1‑86 (MQ=255)
                                 taataaAAATTTCTCGAGGCGGGCACCGAAAAAACGACAAGGAAATTTTTTTTTTTTTCCTCTGTTTTTAAGATGGAAAAAAGaaaataa                                               >  2:1211845/1‑86 (MQ=255)
                                   ataaAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCACATTTTATTTTTCCTCATCTGTTTTTAAATGCGAAAAATGAAAATCAAc                                             >  1:1034376/1‑90 (MQ=255)
                                   ataaAAATTTCTCGAGGCGGGGCCCGAAAAAAGGAAAGGCACATTTTTTTTTTCCTCTTCTGTTTTTAAATGGGAAAAAAGAAATTTaaa                                             >  1:1845404/1‑89 (MQ=255)
                                                 aGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAaaataaa                               <  2:1636643/90‑1 (MQ=255)
                                                  ggCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAaaataaat                              <  2:919985/90‑1 (MQ=255)
                                                          cAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGa                      >  1:231221/1‑90 (MQ=255)
                                                                      cATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTg          >  1:1840177/1‑90 (MQ=255)
                                                                              ttttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAg       >  1:403187/1‑85 (MQ=255)
                                                                              ttttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAg       <  2:403187/85‑1 (MQ=255)
                                                                              ttttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAAtt  <  1:1827367/90‑1 (MQ=255)
                                                                                  |                                                                                       
GATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATT  >  NZ_CP009273/2380275‑2380442

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