Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,819,3180GA51.7% 10.8 / 41.5 29intergenic (‑124/+79)ves/spyenvironmental stress‑induced protein Ves/ATP‑independent periplasmic protein‑refolding chaperone Spy
Reads supporting (aligned to +/- strand):  ref base G (13/1);  new base A (0/15);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 2.06e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-05
Rejected as consensus: Frequency below/above cutoff threshold.

AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGC  >  NZ_CP009273/1819232‑1819405
                                                                                       |                                                                                       
aGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaagaaa                                                                                      >  2:627587/1‑90 (MQ=255)
 ggTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAgaaaagaaaa                                                                                     >  2:1454300/1‑90 (MQ=255)
     aCCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccgc                                                                                 >  2:196347/1‑90 (MQ=255)
      cggTGGATGTTTATCGCAACCGGGGGTATTTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTTTTTGAAAAAAAAAAAAAAAAccgcc                                                                                <  1:1454299/88‑1 (MQ=255)
      ccGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccgcc                                                                                >  1:145314/1‑90 (MQ=255)
        ggggATGTTTATCGGCAACTCGGGTGTTTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAGAAAAAAAAACCGCCGa                                                                              <  2:696516/88‑1 (MQ=255)
           gATGTTTATCGGAAACTGGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTTTTTGAAAAGAAAAAAAAAAACCGCCGATcc                                                                           <  2:680315/90‑1 (MQ=255)
                ttACCGCCACCTGGGGTGATTTACGTTTGTGGTGATCGGGTTTTTTTTAAAATTTTTTGAGAAGAAAAAAAAAACCCGCCGATCCTGTcc                                                                      <  1:196347/90‑1 (MQ=255)
                        aaCTCGGGTGTTTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCAtt                                                              <  2:1158425/90‑1 (MQ=255)
                         aCTCGGGTGATTTACGATAGTGGGGATCAGGGTTTTTTTACAACTGTTTGAGAAGGGAAAAGAAAACCGCTGATCATGTTCACCGCagga                                                             >  1:533192/1‑87 (MQ=255)
                               gtgtTTTACGTTGGGGGTGATCAGGTTTTTTTTACAACTGTTTGAAAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCaa                                                       <  2:1886741/90‑1 (MQ=255)
                               gTGATTTACTTTAGTGGTGATCAGTTTTTTTTTAAAACTGTTGGAAAAAAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCaa                                                       <  1:1612435/90‑1 (MQ=255)
                               aggATTTACGTTGGTGGTGATCAGTTTTTTTTTACACCTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa                                                       <  1:1330396/88‑1 (MQ=255)
                                    ttACGTTAGTGGTGATCGGTTTTTTTTTACAATTTTTTGAGAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg                                                  <  1:711193/90‑1 (MQ=255)
                                     tccGTTGGTGGTGACCAGTTTTTTTTTAAAACTTTTTAGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg                                                  <  1:877491/88‑1 (MQ=255)
                                           aGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAcct                                           >  2:1330397/1‑88 (MQ=255)
                                            gtggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAg                                          >  1:296741/1‑90 (MQ=255)
                                            gtggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAg                                          >  1:938649/1‑90 (MQ=255)
                                                 gATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcgg                                     >  1:1269778/1‑90 (MQ=255)
                                                 gATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcgg                                     >  1:959418/1‑90 (MQ=255)
                                                  aTCGGGTTTTTTTTAAACCTGTTTGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc                                    <  2:396624/90‑1 (MQ=255)
                                                  aTCAGTTTTTTTTTACAATTGTTGAAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc                                    <  1:1623216/90‑1 (MQ=255)
                                                    cAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcgg                                  >  1:1277687/1‑90 (MQ=255)
                                                          ttttttAAACTTTTTGGAAAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGATCTTaa                            <  1:1474279/90‑1 (MQ=255)
                                                                     tttttGAGAAAAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCt                 <  1:1249067/88‑1 (MQ=255)
                                                                     tGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTgg                                               >  2:1633255/1‑60 (MQ=255)
                                                                     tGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTgg                                               <  1:1633254/60‑1 (MQ=255)
                                                                      gTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTTTTAAGTTTTTTGGTTg                >  1:365757/1‑90 (MQ=255)
                                                                      gTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTg                >  2:1715817/1‑90 (MQ=255)
                                                                      gTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGGTTTTTGGCTg                >  2:202572/1‑90 (MQ=255)
                                                                           gaaaagagaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAg            <  2:61741/87‑1 (MQ=255)
                                                                                aaaaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGattat       <  2:365757/87‑1 (MQ=255)
                                                                                     aagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGGCTTAAGTTTTTTGGGTGGAAGATTATTcaga  >  1:93210/1‑89 (MQ=255)
                                                                                       |                                                                                       
AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGC  >  NZ_CP009273/1819232‑1819405

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