Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,380,3410GA58.6% ‑4.8 / 30.1 29intergenic (+31/‑72)yfbL/yfbMM28 family metallopeptidase/YfbM family protein
Reads supporting (aligned to +/- strand):  ref base G (5/7);  new base A (17/0);  total (22/7)
Fisher's exact test for biased strand distribution p-value = 5.07e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.75e-05
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACT  >  NZ_CP009273/2380210‑2380478
                                                                                                                                    |                                                                                                                                         
aaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATg                                                                                                                                    <  1:334649/139‑1 (MQ=255)
     tACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAACGACATGcacaa                                                                                                                               >  1:811774/1‑138 (MQ=255)
     tACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCAGAAAAAAAGAAAAGCACat                                                                                                                               >  2:250381/1‑139 (MQ=255)
          cTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCACattttat                                                                                                                          <  2:416823/139‑1 (MQ=255)
             gCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCCGAAAAAAAGAAATGCACattttatttt                                                                                                                       >  1:333255/1‑139 (MQ=255)
                   tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAAAGACATGCACATTTTATTTTTCCTcc                                                                                                                 >  2:66751/1‑138 (MQ=255)
                                          aTCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑TCAGAAAAAACGAAATGCACATTTTATTTTTCCTCCTCTTTTTTAAGATGCGAAACAt                                                                                           >  1:821733/1‑138 (MQ=255)
                                            cAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTATCGAGGCGGGTACAAGAAAAAATGA‑AGACGCACTTTATTTTTCCCCTTCTTTTTTAAGATAGGAAACATGa                                                                                         >  1:419220/1‑138 (MQ=255)
                                                         aGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACggagg                                                                            <  1:134246/138‑1 (MQ=255)
                                                                     gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAACGAAAAGGACATTTTATTTTTTCTCATCTGGTTTTAAATGAGAAAAATGAAATTCAAACGAGGGAAAAAAAATtgg                                                               >  1:250308/1‑139 (MQ=255)
                                                                                ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑TCAGAAAAAAAGACAAGCACATTTTATTTTTCCTTATCTGTTTTTAAATGGGAAAAATGAAAATAAAAGGGGGGAAAAAAAATGGGGATTATTGGc                                                     >  2:547618/1‑138 (MQ=255)
                                                                                ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCCGAAAAAAAGGCAAGCACAATTTTTTTTTCCTCCACTGTTTTTAAAAGCGAAAAAAGAAAATTAAAGGAGGAAAAAAAAATGGGGAAGAtttgg                                                     >  1:511718/1‑135 (MQ=255)
                                                                                ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAACGACATGCACATTTTATTTTTCCTCCTCTGTTTTTAGATGCGAAACATGAAAATTAACGGAGGAAAAATAAAAGGGAATGATTgggt                                                    >  1:67596/1‑137 (MQ=255)
                                                                                ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAAAGACAAGAACATTTTATTTTTCCTCTTCTGTTTTTAGATGCGAAAAAAGAAATTCAACGGGGGAAAAAAAAATGGGAATGATTTGc                                                     >  2:552834/1‑138 (MQ=255)
                                                                                ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑TCAGAAAAAACGACATGAACATTTTTTTTTTCCTCATCTGTTTTTAGATGCGAAACATGAAAATTTACGGAGGAAAAATAAATGGGAATGATTTGCt                                                    >  2:190550/1‑139 (MQ=255)
                                                                                ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCAGAAAAAGCGAAAAGGACATTTTATTTTTCCTCATCTGTTTTAAGATGCGGAACAAGAAAATCAAAAGGGGAAAAATAAATGGGGATGATTGGCt                                                    >  1:621438/1‑139 (MQ=255)
                                                                                  gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCCGAAAAAAAGACCAGCACATTTTTTTTTTCCTCCTCTGTTTTTAGATGCGAAAAATAAAATTTAACGGAGGGAAAAAAAATGGGTATGATTGGCGAt                                                  >  2:315196/1‑139 (MQ=255)
                                                                                           caGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCAGAAAAAACGACATACACATTTTATTTTTCCTCATCTGTTTTAAGATGCGAAAAATTAAATTCAACGGAGGAAAAATAAATGGGGATGATTTGCTAttttttaga                                          >  2:13272/1‑133 (MQ=255)
                                                                                           caGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCAGAAAAAAAGACAAGCACATTTTTTTTTTCCCCCTCTTTTTTTAAATGCGAAACATGAAAATTAACGGGGGAAAAATAAATGGGAATGATTGGCTATTTTTCAGa                                          >  1:37369/1‑138 (MQ=255)
                                                                                             gTAAATAATAAAAATTTCTCGAGGCGGGG‑CCCGGAAAAACAGCAAGGACAATTTATTTTTCCTCCTCTGTTTTTAAAAGCGAAAAAAGAAAATTAACGGAGGGAAAATAAAATGGAAAAATTGGGTCTTTTTCAGcaaa                                       >  2:31365/1‑138 (MQ=255)
                                                                                             gTAAATAATAAAAATTTCTCGAGGCGGGG‑CCAGAAAAAAAGAAATGCACATTTTATTTTTCCTCATCTGTTTTTAGATGCGAAACATGAAATTCAACGGAGGGAAAAAAAATGGGGATGATTGGGTATTTTGCAGaaaa                                       >  2:80641/1‑138 (MQ=255)
                                                                                                  taataaAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACt                                  <  1:31365/139‑1 (MQ=255)
                                                                                                                   ggCGGGT‑CCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCaa                 <  2:511718/139‑1 (MQ=255)
                                                                                                                   ggCGGGT‑CCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCaa                 <  2:621438/139‑1 (MQ=255)
                                                                                                                                   aGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGa                         >  1:528023/1‑116 (MQ=255)
                                                                                                                                   aGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGa                         <  2:528023/116‑1 (MQ=255)
                                                                                                                                   aGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACt  >  2:18474/1‑139 (MQ=255)
                                                                                                                                   aGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACt  >  2:353097/1‑139 (MQ=255)
                                                                                                                                   aGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACt  >  1:430876/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                                         
AACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACT  >  NZ_CP009273/2380210‑2380478

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