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

GATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAA  >  NZ_CP009273/2380275‑2380451
                                                                                       |                                                                                         
gATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTccc                                                                                         <  1:5433/90‑1 (MQ=255)
  tGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGGAAAAAAGACAAGCACAATTTTTTTTTCCTc                                                                                       >  1:828928/1‑90 (MQ=255)
   ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCAGAAAAAAAGAAAAGAACAATTTTTTTTTCCTct                                                                                      >  2:2201619/1‑89 (MQ=255)
   ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACCTTTTTTTTTTCCTcc                                                                                      >  1:191490/1‑89 (MQ=255)
    gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCAt                                                                                     <  1:24175/90‑1 (MQ=255)
    gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCAt                                                                                     <  2:1125757/90‑1 (MQ=255)
     tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                                                    <  1:868511/90‑1 (MQ=255)
     tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                                                    <  1:76790/90‑1 (MQ=255)
     tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                                                    <  2:1569546/90‑1 (MQ=255)
     tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                                                    <  1:1461846/90‑1 (MQ=255)
     tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                                                    <  2:1033323/90‑1 (MQ=255)
     tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                                                    <  1:2736629/90‑1 (MQ=255)
     tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                                                    <  2:2179289/90‑1 (MQ=255)
     tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                                                    <  1:1990512/90‑1 (MQ=255)
       ttATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCAAAAAAAAGACAAGCAAAATTTTTTTTTTCTCctctt                                                                                  >  1:3376421/1‑89 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAGAAAAGCACAATTTATTTTTTCTCATCTGTTTTTAg                                                                           >  2:254837/1‑90 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCCGAAAAAAAGACAAGAAAAATTTTTTTTTTCCCATCTGtttttaa                                                                           >  2:2138774/1‑87 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCAGAAAAAAAGAAAAGAAAAATTTTTTTTTCCTCATATGtttttaa                                                                           >  2:2624656/1‑87 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCAAAAAAAACGAGACGCGCCTTTTTTTTTTCCTCATCTGtttttta                                                                           >  1:591431/1‑87 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACATGCACATTTTTTTTTTCCCCATCtttttttaa                                                                           >  1:1363496/1‑87 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACAAGCACATTTTATTTTTTCCCATCtttttttta                                                                           >  1:629886/1‑87 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAATCACATTTTATTTTTTCTCCTCtttttttaa                                                                           >  1:166378/1‑87 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGGACATTTTTTTTTTCCTCCTCtttttttaa                                                                           >  1:2234009/1‑87 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCACAATTTTTTTTTCCTCATCTGtttttaa                                                                           >  1:2320218/1‑87 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCACAAAAAACCAAAAGAAAAATTTTTTTTTTCTCATCTGTTTTTAg                                                                           >  2:3205971/1‑90 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGAAAAAAAGAAAAGAACAATTTTTTTTTCCCCCTCtttttttaa                                                                           >  2:1976913/1‑87 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAACGACAAGGACAATTTATTTTTTCTCATTTGTTTTTAg                                                                           >  1:2312657/1‑90 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACATGCACATTTTTTTTTTCCTCATCTGtttttaa                                                                           >  1:2092111/1‑87 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCCCATTTTTTTTTTCCTCCTCtttttttaa                                                                           >  1:1930794/1‑87 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAGGACAATTTTTTTTTTCTCCTCtttttttaa                                                                           >  2:2893049/1‑87 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAACAAATTTTTTTTTTTCTCATCTGtttttat                                                                           >  2:1461846/1‑87 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCAAAAAAAAAAAAAGCAAAATTTTTTTTTTCTCAACTGtttttaa                                                                           >  2:737162/1‑87 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCCTCTGtttttaa                                                                           >  1:1651929/1‑87 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCATCtttttttaa                                                                           >  1:1212746/1‑87 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCACAATTTTTTTTTCCTCATCTGtttttaa                                                                           >  1:586229/1‑87 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACATTTTTTTTTTCCTCTTCTTTTTTTAg                                                                           >  1:632664/1‑90 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACAATTTTTTTTTTCTCTTCtttttttaa                                                                           >  2:2499095/1‑87 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAAACAAGCACAATTTTTTTTTTCTCCTCTGTTTTTAg                                                                           >  2:2736629/1‑90 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAAAAAAAAAGACATGCACCTTTTTTTTTTCCTCCtttttttttaa                                                                           >  1:877309/1‑87 (MQ=255)
               ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAGGACAAGCACAATTTTTTTTTCCTTATCTGtttttaaa                                                                          >  1:672990/1‑86 (MQ=255)
               ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGa                                                                          <  1:1926334/90‑1 (MQ=255)
                    ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGa                                                                     <  2:1579780/90‑1 (MQ=255)
                    ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGa                                                                     <  2:3375584/90‑1 (MQ=255)
                    ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAAAAAAAAGGAAAAGAACATTTTTTTTTTCCTCCTCTTTTTTTAAATGGGa                                                                     >  1:1362347/1‑90 (MQ=255)
                      acaacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGaaa                                                                   <  1:889899/90‑1 (MQ=255)
                       caacaGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAGCACATTTTTTTTTTCCTCCTCTTTTTTTAAATGGGaaaa                                                                  >  1:3384284/1‑89 (MQ=255)
                       caacaGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCACATTTTTTTTTTCCTCCTCTTTTTTTAGATGCGaaaa                                                                  >  1:54168/1‑89 (MQ=255)
                       caacaGAAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAAAGAAAAAGACATTTTTTTTTTCCTCCTCTGTTTTTATATGTGAAAc                                                                  >  2:889899/1‑90 (MQ=255)
                         acaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACAt                                                                <  1:3299695/90‑1 (MQ=255)
                           aGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGa                                                              <  1:2368296/90‑1 (MQ=255)
                           aGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAACGACAAGCACATTTTTTTTTTCCCCCTCTTTTTTTATATGGGGAAAATGa                                                              >  1:1440424/1‑90 (MQ=255)
                              aaaTAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGAAAAGGACATTTTTTTTTTTCTCCTCTGTTTTTATATGGGAAAAAAGaaaa                                                           >  2:784888/1‑89 (MQ=255)
                                 taataaAAATTTCTCGAGGCGGGTCCAAAAAAAACGACATGCACATTTTTTTTTTCCTCCTCTTTTTTTAAATGTGAAAAAAGaaaatta                                                        >  1:2647003/1‑86 (MQ=255)
                                 taataaAAATTTCTCGAGGCGGGGGCCGAAAAAAAGACATGCACAATTTTTTTTTCCTCATCTTTTTTTTTATGGGGAAAATGaaaatta                                                        >  1:1618128/1‑86 (MQ=255)
                                   ataaAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAAc                                                      <  2:3376421/90‑1 (MQ=255)
                                   ataaAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAAc                                                      <  2:191490/90‑1 (MQ=255)
                                   ataaAAATTTCTCGAGGCGGGTCCAGAAAAAAAGAAAAGCACAATTTATTTTTTCTCTTTTGTTTTTAGATGCGAAACATGAAAATCAAc                                                      >  2:940192/1‑90 (MQ=255)
                                     aaaaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACgg                                                    <  2:450656/90‑1 (MQ=255)
                                     aaaaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACgg                                                    <  1:2239405/90‑1 (MQ=255)
                                     aaaaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACgg                                                    <  2:1578909/90‑1 (MQ=255)
                                      aaaaTTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACATTTTTTTTTTCCTCATCTGTTTTTAAATGCGGAAAATGAAAATCaaacgg                                                   >  1:980668/1‑86 (MQ=255)
                                      aaaaTTTCTCGAGGCGGGGCCAGAAAAAACGAAAAGCACATTTTTTTTTTCCTCATCTGTTTTTAAATGCGGAAAATGAAAATTAAcagg                                                   >  1:1095418/1‑87 (MQ=255)
                                         aTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACggagga                                                <  1:1108774/90‑1 (MQ=255)
                                           ttCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGaaa                                              <  2:1882617/90‑1 (MQ=255)
                                           ttCTCGAGGCGGGTCCAGAAAAAAAGACATGCACATTTTTTTTTTCCTCCTCTTTTTTTAGAAGCGAAAAATGAAAATCAACGGGGGaaa                                              >  1:2456598/1‑90 (MQ=255)
                                           ttCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACAATTTTTTTTTTCTCCTCTGTTTTTAAATGCGAAAAAAAAAAATTAAAGGAGGgaa                                              >  2:1851396/1‑90 (MQ=255)
                                            tctcGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGaaaa                                             <  2:1936251/90‑1 (MQ=255)
                                             ctcGAGGCGGATCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGaaaaa                                            <  1:986427/90‑1 (MQ=255)
                                                    cGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATgg                                     <  2:2969731/90‑1 (MQ=255)
                                                       gTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGaa                                  <  1:2840214/90‑1 (MQ=255)
                                                       gTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGaa                                  <  1:1379060/90‑1 (MQ=255)
                                                       gTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGaa                                  <  1:2035320/90‑1 (MQ=255)
                                                        tCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAAt                                 <  2:1840400/90‑1 (MQ=255)
                                                               aaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGc                          <  1:817138/90‑1 (MQ=255)
                                                               aaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGc                          <  2:54168/90‑1 (MQ=255)
                                                                       aTGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGc                  <  1:1976913/90‑1 (MQ=255)
                                                                            caTTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGaaa             >  2:2753684/1‑90 (MQ=255)
                                                                              ttttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAAtt           >  2:1963959/1‑90 (MQ=255)
                                                                                 tattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGAc        >  2:2918189/1‑90 (MQ=255)
                                                                                  attttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACt       <  1:912143/90‑1 (MQ=255)
                                                                                  attttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACt       <  1:2430605/90‑1 (MQ=255)
                                                                                      tCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGTAAAATAAATGGGGATGATTGGCGCTTTTGAGGAAATTGACTCAGa   >  2:2353158/1‑90 (MQ=255)
                                                                                      tCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa   <  2:669957/90‑1 (MQ=255)
                                                                                      tCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa   >  2:679087/1‑90 (MQ=255)
                                                                                      tCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa   >  2:235608/1‑90 (MQ=255)
                                                                                      tCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa   >  1:3191131/1‑90 (MQ=255)
                                                                                      tCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa   >  1:461058/1‑90 (MQ=255)
                                                                                       cccTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGaa  <  2:2956771/90‑1 (MQ=255)
                                                                                       |                                                                                         
GATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAA  >  NZ_CP009273/2380275‑2380451

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