Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,319,9290TG51.8% 36.5 / 119.4 85G114G (GGT→GGG) atoBacetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (2/39);  new base G (44/0);  total (46/39)
Fisher's exact test for biased strand distribution p-value = 3.27e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-16
Rejected as consensus: Frequency below/above cutoff threshold.

TAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAAT  >  NZ_CP009273/2319842‑2320018
                                                                                       |                                                                                         
tAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTAt                                                                                         <  2:2460974/90‑1 (MQ=255)
  aaGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATgg                                                                                       <  2:2917526/90‑1 (MQ=255)
     tATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGgaa                                                                                    >  1:1934227/1‑90 (MQ=255)
     tATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGaaa                                                                                    >  1:1408245/1‑90 (MQ=255)
       tgtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGgaaaa                                                                                  >  2:2516287/1‑89 (MQ=255)
       tgtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGgaaaa                                                                                  >  2:321240/1‑89 (MQ=255)
       tgtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGgaaaa                                                                                  >  1:1929126/1‑89 (MQ=255)
       tgtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGgaaaa                                                                                  >  2:1688815/1‑89 (MQ=255)
          ggTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATg                                                                               <  1:1257356/90‑1 (MQ=255)
           gTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGa                                                                              >  1:1340370/1‑90 (MQ=255)
           gTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGa                                                                              >  1:1380946/1‑90 (MQ=255)
            ttCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATGAg                                                                             >  2:675161/1‑90 (MQ=255)
             tCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATGAGt                                                                            >  2:730749/1‑90 (MQ=255)
              cGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGaaa                                                                                    <  2:478222/81‑1 (MQ=255)
              cGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGgaa                                                                                    >  1:478222/1‑81 (MQ=255)
              cGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAggt                                                                           >  1:3299697/1‑88 (MQ=255)
              cGGGTCTTAAAAGTGTGGCGCTTACCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGtt                                                                           >  1:2531306/1‑90 (MQ=255)
               gggTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGttt                                                                          <  2:877969/90‑1 (MQ=255)
               gggTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGttt                                                                          <  1:2516287/90‑1 (MQ=255)
                ggTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATGAGGTTa                                                                         >  2:192864/1‑90 (MQ=255)
                ggTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGTTTa                                                                         >  1:2565112/1‑90 (MQ=255)
                   cTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGcc                                                                      <  2:427001/90‑1 (MQ=255)
                    ttAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGccc                                                                     <  1:1633025/90‑1 (MQ=255)
                    ttAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGccc                                                                     <  1:1613356/90‑1 (MQ=255)
                    ttAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGccc                                                                     <  1:2386990/90‑1 (MQ=255)
                    ttAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGccc                                                                     <  2:2726003/90‑1 (MQ=255)
                     tAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGcccc                                                                    <  1:2708374/90‑1 (MQ=255)
                          gtgtGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCcact                                                               >  1:2347758/1‑90 (MQ=255)
                          gtgtGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTTGCCCCCcact                                                               >  1:1257375/1‑90 (MQ=255)
                          gtgtGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGGTTTGCCCCCcact                                                               >  1:1928500/1‑90 (MQ=255)
                           tgtgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGGTTAGCCCCCcactt                                                              >  1:1586873/1‑90 (MQ=255)
                            gtgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGttt                                                                          <  2:689245/77‑1 (MQ=255)
                            gtgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTTTCCCCCcacttt                                                             >  1:2944391/1‑89 (MQ=255)
                            gtgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGttt                                                                          >  1:689245/1‑77 (MQ=255)
                              ggCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGGTTAGCCCCCcacttact                                                           >  1:2058774/1‑90 (MQ=255)
                               gcgcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGGGGGAAAAATGAGTTTAGGCCCCCACCTtctc                                                          >  2:1028257/1‑90 (MQ=255)
                                cgcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCg                                                         <  1:362297/90‑1 (MQ=255)
                                 gcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATGAGGTTTGCCCCCCACTTACTCGa                                                        >  1:2861171/1‑90 (MQ=255)
                                 gcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGGGGGAAAAATTAGTTTTGCCCCCCCCTTTCTCGa                                                        >  2:831272/1‑90 (MQ=255)
                                  cTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAAAGGGGTTTGCCCCCCCCTTACTCggt                                                       >  2:884656/1‑88 (MQ=255)
                                    tGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATTAGTTTTGCCCCCCACTTAATCGATGc                                                     >  1:508884/1‑90 (MQ=255)
                                           cAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCa                                              <  1:2705026/90‑1 (MQ=255)
                                            aGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCAc                                             <  1:521235/90‑1 (MQ=255)
                                             ggCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGTTGGGAAATATGAGGTTTTGCCCCTACTTTCTCCATGGCAAAACACg                                            >  1:1443952/1‑90 (MQ=255)
                                             ggCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGGTTAGCCCCCCACTTACTCGATGCAAAAACACg                                            >  1:2763729/1‑90 (MQ=255)
                                              gCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGc                                           <  1:831272/90‑1 (MQ=255)
                                              gCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGc                                           <  1:2964590/90‑1 (MQ=255)
                                              gCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGc                                           <  1:2319415/90‑1 (MQ=255)
                                              gCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGc                                           <  1:3169792/90‑1 (MQ=255)
                                               ccATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGct                                          <  2:3299697/90‑1 (MQ=255)
                                                 aTTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGttt                                                                          <  2:3351543/56‑1 (MQ=255)
                                                 aTTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGttt                                                                          >  1:3351543/1‑56 (MQ=255)
                                                  ttCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTTGCCCCCCACTTACTCGATGCAAAAACACGCTCTg                                       >  1:3284359/1‑90 (MQ=255)
                                                  ttCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGTTTAGCCCCCCACTTACTCGATGCAAAAACACGCTCTg                                       >  1:151715/1‑90 (MQ=255)
                                                  ttCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAAAGAGGTTTGCCCCCCCCCTACTCGATGCAAAAAGACGCCCTg                                       >  2:2964590/1‑90 (MQ=255)
                                                      ggcaggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTa                                   <  2:945376/90‑1 (MQ=255)
                                                       gcaggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTAt                                  <  1:675161/90‑1 (MQ=255)
                                                       gcaggTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTAGCCCCCCACTTTCTCCATGCAAAAACACGCCCTGGTTAt                                  >  1:78331/1‑90 (MQ=255)
                                                       gcaggTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAAAGAGTTTTGCCCCCCCCTTTCTCGATGCAAAAGAACGCCCTGGGTAt                                  >  2:2764071/1‑90 (MQ=255)
                                                         aggTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATTAGTTTAGCCCCCCCCTTACTCGATGCAAAAACACGCCCTGGGTAttt                                >  2:2736431/1‑88 (MQ=255)
                                                         aggTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTAAGCCCCCCCTTATTCCAAGCCAAAAGACGCTCTGGGTAtag                                >  2:2028233/1‑88 (MQ=255)
                                                         aggTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAACGAGTTTTGCCCCCCCCTTACTCCATGCAAAAACACGCTCTGGTTATCg                                >  1:1527512/1‑90 (MQ=255)
                                                                  gcAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGgaga                       <  2:3166682/90‑1 (MQ=255)
                                                                    agagCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGcc                                                                      <  2:3304189/41‑1 (MQ=255)
                                                                    agagCATTGTGGCGGGGGGGATGGGAAAAATGAGGTTAGCCCCCCACTTACTCGATGGAAAAAAACGCCCTGGGTATCGTCTTGGAGACg                     >  2:398824/1‑90 (MQ=255)
                                                                    agagCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGcc                                                                      >  1:3304189/1‑41 (MQ=255)
                                                                    agagCATTGTGGCGGGGGGGATGGAAAATATAAGCTTAGCCCCCCACTTCTTCAATGCAAAAACCCGCGCTGGTTATCGTATTGgggagg                     >  1:2195300/1‑88 (MQ=255)
                                                                       gCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGAc                  <  1:2552584/90‑1 (MQ=255)
                                                                       gCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGAc                  <  1:2543759/90‑1 (MQ=255)
                                                                       gCATTGTGGCGGGGGGTATGGAAAATATGAGGTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGAc                  <  1:2446619/90‑1 (MQ=255)
                                                                        cATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACa                 <  1:2781174/90‑1 (MQ=255)
                                                                        cATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACa                 <  2:741625/90‑1 (MQ=255)
                                                                        cATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACa                 <  2:3020565/90‑1 (MQ=255)
                                                                         aTTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAg                <  1:1274507/90‑1 (MQ=255)
                                                                          ttGTGGCGGGGGGGATGGAAAATATTAGTTTAGCCCCCCACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTTGAGACGGACAgg               >  1:655896/1‑90 (MQ=255)
                                                                            gtgGCGGGGGGTATGGAAAATATGAGTTTAG‑CCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGttt            <  2:32071/90‑1 (MQ=255)
                                                                            gtgGCGGGGGGGATGGAAAATATGAGTTTTGCCCCCCACTTACTCGATGCAAAAACACGCTCTGGTTATCGTCTTTGGGACGGACAgggt             >  1:593506/1‑88 (MQ=255)
                                                                               gCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTAt          <  1:2972071/90‑1 (MQ=255)
                                                                               gCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTAt          <  1:2003922/90‑1 (MQ=255)
                                                                               gCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTAt          <  2:410445/90‑1 (MQ=255)
                                                                               gCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTAt          <  2:2306524/90‑1 (MQ=255)
                                                                                 ggggggGATGGAAAAAATGAGTTTAGCCCCCCACTTACTCGATGCAAAAGGACGCTCTGGTTATCGTCTTTGGGACGGACAGGGTTTTGa        >  2:1993615/1‑90 (MQ=255)
                                                                                      gTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACa                 >  1:2095436/1‑76 (MQ=255)
                                                                                      gTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACa                 <  2:2095436/76‑1 (MQ=255)
                                                                                      gTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTaa   >  1:876964/1‑90 (MQ=255)
                                                                                       tATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAAt  <  2:644512/90‑1 (MQ=255)
                                                                                       |                                                                                         
TAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAAT  >  NZ_CP009273/2319842‑2320018

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