Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,380,3620CT54.8% 7.3 / 37.5 31intergenic (+52/‑51)yfbL/yfbMM28 family metallopeptidase/YfbM family protein
Reads supporting (aligned to +/- strand):  ref base C (2/12);  new base T (17/0);  total (19/12)
Fisher's exact test for biased strand distribution p-value = 6.45e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.78e-06
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACAT‑GAAATTCAACGGAGG‑AAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGG  >  NZ_CP009273/2380225‑2380494
                                                                                                                                         |                                                                                                                                      
ccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTcc                                                                                                                                       <  1:724557/139‑1 (MQ=255)
                   aTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGACAAGAAAAAAAGAGAAGAGAACTTTATTTTTTCTCATCTGTTTTTAGaaac                                                                                                                    >  2:384360/1‑136 (MQ=255)
                       aaTTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGaa                                                                                                                 <  1:384360/138‑1 (MQ=255)
                           aTCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAGGAAAAGCATATTTTTTTTTTTTTCATATGTTTTTAAATGGGaaaaaa                                                                                                             >  1:335052/1‑135 (MQ=255)
                                                tGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCCGAAAAAACGAACTGCACATTTTTTTTTTCCTCCTTTGTTTTTAAATGCGAAAAAT‑GAAATTTAAAGGAGG‑AAaaat                                                                                      >  1:587175/1‑138 (MQ=255)
                                                      gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAAGACATGCAAATTTTTTTTTTTCTCCTCTGTTTTTAAATGCGAAAAAA‑GAAAATTAACGGAGG‑GAAAATAAATggg                                                                               >  2:578410/1‑139 (MQ=255)
                                                      gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGGCAAGCACATTTTTTTTTTCCTCATCTGTTTTTAGATGAGAAAAAT‑GAAAATCAACGGAGG‑GAAAATAAATtgg                                                                               >  2:373154/1‑139 (MQ=255)
                                                      gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACATTTTATTTTTCCTCCTCTGTTTTTAAATGCGAAACAT‑GAAAATCAAAGGAGG‑AAAAATAAATggg                                                                               >  1:284269/1‑139 (MQ=255)
                                                      gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCACAATTTATTTTTCCTCATCTGTTTTTAGATGCGGAAAAA‑GAAAATTAACGGGGG‑GAAAATAAATggg                                                                               >  1:192900/1‑139 (MQ=255)
                                                                 ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAAGAAATGCACATTTTATTTTTCCTCTTCTGTTTTTAAATGCGAAACAT‑GAAATTCAAAGGGGG‑GAAAATAAATGGGGATGATTgggt                                                                    >  2:227743/1‑137 (MQ=255)
                                                                 ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCAGAAAAAACGACATGCACAATTTTTTTTTCCTCATCTGTTTTAAGATGCGGAAAAT‑GAAATTCAAAGGGGG‑GAAAATAAATTGGGATGATTGGc                                                                     >  1:445024/1‑138 (MQ=255)
                                                                 ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGCACATTTTTTTTTTCCTCATCTGTTTTTAGATTCGAAAAAT‑TAAATTTAACGGAGG‑AAAAATAAATGGGGATGATTGGCt                                                                    >  1:29633/1‑139 (MQ=255)
                                                                 ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTACAGAAAAAGAGACATGCACATTTTATTTTTCCTCATCTGTTTTAAGAAGCGAAACAA‑GAAATTCAACGGGGG‑AAAAATAAATGGGAATGGTTggg                                                                     >  2:741702/1‑137 (MQ=255)
                                                                 ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCATCTGTTTTTAGATGCGGAAAAA‑GAAATTCAAAGGAGG‑AAAAATAAATGGGAATGATTTGCt                                                                    >  1:132765/1‑139 (MQ=255)
                                                                       taCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACAT‑GAAATTCAACGGAGG‑AAAAATAAATGGGAATGATTGGCTATTTTg                                                              <  1:227743/139‑1 (MQ=255)
                                                                       acagacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACAT‑GAAATTCAACGGAGG‑AAAAATAAATGGGAATGATTGGCTATTTTg                                                              <  1:88041/135‑1 (MQ=255)
                                                                            caGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACAATTTATTTTTCCCCATCTGTTTTAAAAAGCGAAAAAT‑GAAAATTAAAGGAGG‑AAAAATAAATGGGAATGATTTGCTTTTTTTCaaa                                                          >  1:291808/1‑136 (MQ=255)
                                                                              gTAAATAATAAAAATTTCTCGAGGCGGGTCCCAAAAAAACAAGAATAACAATTTATTTTTTTTCATCTGTTTTTAAAAGCGAAAAAA‑GAAACTCAAAGGAGG‑GAAAATAAAAAGGAATGGTTGGGGATTTTTCAGaaaa                                                       >  2:18383/1‑138 (MQ=255)
                                                                               tAAATAATAAAAATTTCTCGAGGCGGGTGCAGAAAAAACGACAAGCACATTTTTTTTTTCCTCATCTGTTTTTAGATGCGAAAAAT‑GAAAATTAACGGGGG‑GAAAATAAAAGGGGATGATTTGCTATTTTTCAGAAAtt                                                      >  1:717756/1‑139 (MQ=255)
                                                                                 aataatAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGCACATTTTATTTTTCCTCATCTGTTTTTAGATGCGGAAAAT‑GAAATTTAACGGAGG‑AAAAATAAATGGGGATGATTTGCTATTTTTCAGAAATTg                                                     >  1:268545/1‑138 (MQ=255)
                                                                                  ataataAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACAT‑GAAATTCAACGGAGG‑AAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGAc                                                   <  2:717756/139‑1 (MQ=255)
                                                                                      taaAAATTTCTCGAGGCGGGTAGCACAAAAGCGAGACAAAAAATTTTTTTTCCATCATCAGTTTTATGATGCAGAACATGGAAATTC‑ACGGAGGTAAAAAAAAATGGG‑TTGATTTGGTATTTTTTAGaaactatct                                                  >  1:546457/1‑130 (MQ=255)
                                                                                      taaAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCACATTTTTTTTTTCCTCATCTGTTTTAAGATGCGGAACAT‑GAAATTCAACGGAGG‑GAAAAAAAAATGGAATGATTGGGTATTTTTCAGAAATTGaactgt                                               >  2:88041/1‑134 (MQ=255)
                                                                                        aaaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACAT‑GAAATTCAACGGAGG‑AAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGaa                                             <  2:29633/139‑1 (MQ=255)
                                                                                                  gAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACAT‑GAAATTCAACGGAGG‑AAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAAcc                                   <  2:368202/139‑1 (MQ=255)
                                                                                                    ggCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACAT‑GAAATTCAACGGAGG‑AAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCaa                                 <  2:262283/139‑1 (MQ=255)
                                                                                                                         aTGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACAT‑GAAATTCAACGGAGG‑AAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGaaaaa            >  1:723555/1‑139 (MQ=255)
                                                                                                                                  ttattttCCCTCATCTGTTTTAAGATGCGAAACAT‑GAAATTCAACGGAGG‑AAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGAt    <  2:268545/138‑1 (MQ=255)
                                                                                                                                  ttattttCCCTCATCTGTTTTAAGATGCGAAACAT‑GAAATTCAACGGAGG‑AAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGAt    <  1:578410/138‑1 (MQ=255)
                                                                                                                                  ttattttCCCTCATCTGTTTTAAGATGCGAAACAT‑GAAATTCAACGGAGG‑AAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGAt    <  2:445024/138‑1 (MQ=255)
                                                                                                                                  ttattttCCCTCATCTGTTTTAAGATGCGAAACAT‑GAAATTCAACGGAGG‑AAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGAt    <  2:674329/138‑1 (MQ=255)
                                                                                                                                    attttCCCTCATCTGTTTTAAGATGCGAAACAT‑GAAATTCAACGGAGG‑AAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATgg  <  2:192900/138‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                      
CCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACAT‑GAAATTCAACGGAGG‑AAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGG  >  NZ_CP009273/2380225‑2380494

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑