Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,819,3180GA57.1% 1.9 / 42.3 35intergenic (‑124/+79)ves/spyenvironmental stress‑induced protein Ves/ATP‑independent periplasmic protein‑refolding chaperone Spy
Reads supporting (aligned to +/- strand):  ref base G (13/2);  new base A (0/20);  total (13/22)
Fisher's exact test for biased strand distribution p-value = 7.11e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.56e-06
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAG  >  NZ_CP009273/1819238‑1819395
                                                                                   |                                                                             
cggTGGTTGGTTATCGGCAACTGGGGTGTTTTACGTTAGTGGTGATCAGGTTTTTTTTACAATTGTTTGAAAAGAAAAAAAAAAAccgcc                                                                      <  1:224875/88‑1 (MQ=255)
cggTGGATGTTTATCGGCAACTGGGGTGATTTACGTTAGTGGTGATCGGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAAccgcc                                                                      <  2:138334/88‑1 (MQ=255)
    ggATGTTTACCGGCAACTGGGGTAATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATc                                                                  <  1:206781/90‑1 (MQ=255)
    ggATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATc                                                                  >  2:110255/1‑90 (MQ=255)
     gATGGTTATCGGCAACTCGGGTGTTTTACGTTAGGGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATcc                                                                 <  2:10839/90‑1 (MQ=255)
         gTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTc                                                             >  2:1457905/1‑90 (MQ=255)
          ttATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTcc                                                            >  2:907313/1‑90 (MQ=255)
            aTCGGCACCGGGGGTGATTTACGTTGGGGGTGATCGGGTTTTTTTTACAATTTTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCAc                                                          <  1:256944/90‑1 (MQ=255)
             tCGGCAACTGGGGTTTTTTACTTTGGTGGTGATCGGGTTTTTTTTACACTTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCAcc                                                         <  2:1392731/90‑1 (MQ=255)
                  agcTCGGTTTATTTACGTTAGTGGTGATCGGGTTTTTTTTACACCTTTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCAtt                                                    <  2:335539/88‑1 (MQ=255)
                  aaCTGGGGTAATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCAtt                                                    <  2:1103786/90‑1 (MQ=255)
                  aaCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCAtt                                                    <  2:204894/90‑1 (MQ=255)
                   aCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTa                                                   >  2:1012247/1‑90 (MQ=255)
                   aCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTa                                                   >  1:169832/1‑90 (MQ=255)
                    cgggggTGTTTTACGTTAGTGGTAATCAGTTTTTTTTAAAAATTTTTTGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTAc                                                  <  1:907313/87‑1 (MQ=255)
                       gggTGATTTCCTTTGGGGGTGACCGGTTTTTTTTAAAACCTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGc                                               <  1:504265/90‑1 (MQ=255)
                         gttttttCCGTTAGTGGTGATCAGGTTTTTTTTACAACTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa                                             <  2:1488369/86‑1 (MQ=255)
                         gttaTTTACGTTAGGGGTGATCGGGTTTTTTTTACAACTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa                                             <  1:1186488/87‑1 (MQ=255)
                         gtgtTTTAGGTTGGTGGTGATCGGCTTTTTTTTAAACTTTTTTGAGAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCaa                                             <  1:1248863/90‑1 (MQ=255)
                         gtgtTTTACGTTAGTGGTGATCAGGTTTTTTTTACAATTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa                                             <  2:991450/90‑1 (MQ=255)
                         gTGATTTACGTTGGGGGTGATCAGTTTTTTTTTACAATTTTTGGAGAAAAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCaa                                             <  1:1012247/90‑1 (MQ=255)
                         aggATTTACGTTAGTGGTGATCGGTTTTTTTTAACAATTTTTTGAAAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCaa                                             <  1:1262781/88‑1 (MQ=255)
                            ttttACGTTAGTGGTGATCGGTTTTTTTTTACACTTTTTGAAAAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGt                                          <  1:1457905/89‑1 (MQ=255)
                            aTTTACTTTAGGGGTCACCAGGTTTTTTTTCCATTTTTTTAAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGt                                          <  1:1262222/90‑1 (MQ=255)
                              ttACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGtct                                        >  1:335539/1‑88 (MQ=255)
                              cagcGTTAGTGGTGATCGGTTTTTTTTTACACCTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg                                        <  1:383823/87‑1 (MQ=255)
                                     aGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAcct                                 >  2:1262781/1‑88 (MQ=255)
                                     aGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACaa                                 >  2:891260/1‑90 (MQ=255)
                                      gtggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAg                                >  1:1204525/1‑90 (MQ=255)
                                           gATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGAcctg                           >  2:383823/1‑88 (MQ=255)
                                                ggCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcgggc                      >  1:1596959/1‑88 (MQ=255)
                                                ggCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTc                      >  1:646040/1‑90 (MQ=255)
                                                     tttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAg                 >  2:1305503/1‑90 (MQ=255)
                                                                 tttttGAAAAAAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCt       <  2:795560/88‑1 (MQ=255)
                                                                      gagaagaaaaAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAg  <  1:1525641/90‑1 (MQ=255)
                                                                             agaaaagaaaACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc                          >  1:107447/1‑60 (MQ=255)
                                                                             aaaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc                          <  2:107447/58‑1 (MQ=255)
                                                                                   |                                                                             
CCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAG  >  NZ_CP009273/1819238‑1819395

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