Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,380,3620CT54.2% 4.0 / 32.8 24intergenic (+52/‑51)yfbL/yfbMM28 family metallopeptidase/YfbM family protein
Reads supporting (aligned to +/- strand):  ref base C (5/6);  new base T (13/0);  total (18/6)
Fisher's exact test for biased strand distribution p-value = 3.43e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.68e-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: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                <  2:403187/85‑1 (MQ=255)
                                                                              ttttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAg                >  1:403187/1‑85 (MQ=255)
                                                                              ttttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAAtt           <  1:1827367/90‑1 (MQ=255)
                                                                                   ttttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTc      >  2:1084379/1‑90 (MQ=255)
                                                                                   ttttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTc      >  2:1270137/1‑90 (MQ=255)
                                                                                       cccTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGaa  <  2:1034376/90‑1 (MQ=255)
                                                                                       cccTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGaa  <  1:2052680/90‑1 (MQ=255)
                                                                                       |                                                                                           
GATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAA  >  NZ_CP009273/2380275‑2380451

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