Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,380,3620CT40.0% 46.5 / 21.8 40intergenic (+52/‑51)yfbL/yfbMM28 family metallopeptidase/YfbM family protein
Reads supporting (aligned to +/- strand):  ref base C (5/19);  new base T (16/0);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 3.39e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATA  >  NZ_CP009273/2380225‑2380500
                                                                                                                                          |                                                                                                                                           
ccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTcc                                                                                                                                            <  1:655290/139‑1 (MQ=255)
      aTACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CAt                                                                                                                                      <  2:296602/138‑1 (MQ=255)
      aTACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CAt                                                                                                                                      <  1:777924/138‑1 (MQ=255)
      aTACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CAt                                                                                                                                      <  1:878169/138‑1 (MQ=255)
      aTACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CAt                                                                                                                                      <  1:290624/138‑1 (MQ=255)
              gacagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAAAGACAAGAACATTTTTTTTTTTCT‑CATCTGTTTTa                                                                                                                              >  2:518702/1‑138 (MQ=255)
                   aTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCC‑GAAAAAAAGACATGCACATTTTATTTTTCCT‑CATATGTTTTTAGAttc                                                                                                                        >  2:213475/1‑137 (MQ=255)
                   aTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CAt                                                                                                                                      <  2:395881/125‑1 (MQ=255)
                   aTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCAAATTTTATTTTTCCT‑CAt                                                                                                                                      >  1:395881/1‑125 (MQ=255)
                         ttATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACa                                                                                                                  <  1:72753/139‑1 (MQ=255)
                              aGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGaaa                                                                                                             <  2:598949/139‑1 (MQ=255)
                                      gCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACg                                                                                                      <  1:518702/138‑1 (MQ=255)
                                          aGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAACGACACGCACATTTTATTTTTCCT‑CATCTGTTTTTAGAAACGAAACATGAAAATCAACggagg                                                                                                  >  1:47943/1‑138 (MQ=255)
                                                      gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCA‑GAAAAAACGACATGCACAATTTATTTTTCCT‑CATCTGTTTTAAGAAACGAAACATAAAATTCCACGGAGGAAaaaaaaattg                                                                                      >  1:270828/1‑136 (MQ=255)
                                                      gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGCTCCAGGAAAAAGGAAAATGC‑CCTTTTATTTTTCCTCCACTTTTTTTAAGTAGGGAAAAAGGAAA‑TCAACGGAAGAAAAAAAAATggg                                                                                     >  2:530469/1‑139 (MQ=255)
                                                                 ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCC‑GAAAAAACGGCATGCACATTTTATTTTTCCT‑CCTCTTTTTTTAGATGCGGAAAATGAAATTTAACCGGGGAAAAAAAAATAGGGATGATTTGCt                                                                          >  1:34799/1‑139 (MQ=255)
                                                                 ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAAAGACAAGCACATTTTTTTTTTCCT‑CATCTGTTTTTAGAAGCGAAAAATGAAAATTAACGGGGGAAAAATAAATTGGAATGAtttggt                                                                          >  1:533725/1‑135 (MQ=255)
                                                                 ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAAAGAAATGCACATTTTATTTTTCCT‑CATCTGTTTTTAGAAGCGAAAAAAGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGc                                                                           >  2:501294/1‑138 (MQ=255)
                                                                 ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAAAGAAAAGAACATTTTTTTTTTCCT‑CATCTGTTTTTAGATGCGAAAAATGAAAATCAAAGGGGGAAAAAAAAATGGGAATGATTGGc                                                                           >  2:186932/1‑138 (MQ=255)
                                                                 ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCA‑GAAAAAAAGAAATGAAAAATTTATTTTTCCT‑CACCTGTTTTTAGAAGGGAAAAAAAAAAATTAAAAAGGGGAAAATAAATGGGGATGAtttgg                                                                           >  2:953469/1‑135 (MQ=255)
                                                                         caacaGTAAATAATAAAAATTTCTCGAGGCGGGGCCA‑GAAAAAAAGACAAGCACATTTTTTTTTCTCT‑CATCTGTTTTTAAAAGGGGAACAAGAAAATCAAAGGGGGGAAAAAAAATGGGGATGATTGGGTAttttt                                                                    >  2:612492/1‑136 (MQ=255)
                                                                          aacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCC‑GAAAAAACGACATGCACAATTTTTTTTTCCT‑CCTCTGTTTTAAGATGCGAAACATGAAAATCAACGGGGGGAAAATAAAATGGAATGATTGGCTAttttt                                                                    >  1:661717/1‑135 (MQ=255)
                                                                              gTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAAt                                                             <  1:938245/139‑1 (MQ=255)
                                                                              gTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAAt                                                             <  2:199268/139‑1 (MQ=255)
                                                                               tAAATAATAAAAATTTCTCGAGGCGGGGCCA‑GAAAAAACGACATGCACATTTTATTTTTCCT‑CATCTTTTTTTAGATGCGGAAAATGAAAATTAAAGGGGGGAAAAAAAATTGGGATGATTGGGTATTTTTTAGaaaat                                                            >  2:938245/1‑137 (MQ=255)
                                                                                 aataatAAAAATTTCTCGAGGCGGGTGCA‑GAAAAAACGACATATACAATTTTTTTTTCCT‑CATCTGTTTTAAGATGAGAAACAAGGAATTTAAAGGGGGGAAAATAAATAGGGATTATTTGGTATTTTGCAGAAAATg                                                           >  1:796043/1‑138 (MQ=255)
                                                                                      taaAAATTTCTCGAGGCGGGTCCA‑GAAAAAACGACATGCACATTTTTTTTTTCCT‑CATCTGTTTTTAGATGCGAAACAAGAAATTTAACGGAGGGAAAATAAATGGGAAAGATTGGGTATTTTGCAGAAAATGACt                                                        >  1:626291/1‑136 (MQ=255)
                                                                                      taaAAATTTCTCGAGGCGGGTCCA‑GAAAAAACGACATGCACAATTTTTTTTTCCT‑CATCTGTTTTTAGATGCGAAAAAAGAAATTTAACGGGGGAAAAAAAAATGGGAATGATTGGCTATTTTTCAGAAAATGACt                                                        >  1:640669/1‑136 (MQ=255)
                                                                                        aaaaTTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa                                                    <  1:381686/138‑1 (MQ=255)
                                                                                               ctcGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGAtt                                             <  1:213475/138‑1 (MQ=255)
                                                                                                            cA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGa                               <  1:501294/139‑1 (MQ=255)
                                                                                                                    aGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc                            >  1:220434/1‑135 (MQ=255)
                                                                                                                    aGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc                            <  2:220434/135‑1 (MQ=255)
                                                                                                                    aGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACt                        >  2:226319/1‑139 (MQ=255)
                                                                                                                     gCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACt                        >  1:693134/1‑138 (MQ=255)
                                                                                                                           tGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAAc                 <  2:364868/139‑1 (MQ=255)
                                                                                                                                attttattttCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTg            <  1:186932/139‑1 (MQ=255)
                                                                                                                                  tttattttCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGAt          >  1:477872/1‑139 (MQ=255)
                                                                                                                                   ttattttCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGAt          <  1:203164/138‑1 (MQ=255)
                                                                                                                                   ttattttCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGAt          <  2:796043/138‑1 (MQ=255)
                                                                                                                                    tattttCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATg         <  2:817772/138‑1 (MQ=255)
                                                                                                                                          cccT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAAcc                                         <  2:51515/100‑1 (MQ=255)
                                                                                                                                          cccT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAAcc                                         >  1:51515/1‑100 (MQ=255)
                                                                                                                                          cccT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCa                                        <  1:147438/101‑1 (MQ=255)
                                                                                                                                          cccT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCa                                        >  2:147438/1‑101 (MQ=255)
                                                                                                                                          cccT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGataata  >  1:408465/1‑139 (MQ=255)
                                                                                                                                          |                                                                                                                                           
CCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATA  >  NZ_CP009273/2380225‑2380500

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