Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,819,3130GA55.9% 9.8 / 41.8 34intergenic (‑119/+84)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/19);  total (13/21)
Fisher's exact test for biased strand distribution p-value = 1.13e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.40e-06
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAA‑C‑TGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAG  >  NZ_CP009273/1819238‑1819395
                                                                              |                                                                                  
cggTGGTTGGTTATCGGCAACTGGGGTGTTTTACGTTAGTGGTGATCAGGTTTTTTTTACAA‑T‑TGTTTG‑AAAAGAAAAAAAAAAAccgcc                                                                      <  1:224875/88‑1 (MQ=255)
cggTGGATGTTTATCGGCAACTGGGGTGATTTACGTTAGTGGTGATCGGGTTTTTTTTACAA‑C‑TGTTTG‑AGAAGAAAAAAAAAAAccgcc                                                                      <  2:138334/88‑1 (MQ=255)
    ggATGTTTACCGGCAACTGGGGTAATTTACGTTAGTGGTGATCAGGCTTTTTTTACAA‑C‑TGTTTG‑AGAAGAGAAAAGAAAACCGCCGATc                                                                  <  1:206781/90‑1 (MQ=255)
    ggATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAA‑C‑TGTTTG‑AGAAGAGAAAAGAAAACCGCCGATc                                                                  >  2:110255/1‑90 (MQ=255)
     gATGGTTATCGGCAACTCGGGTGTTTTACGTTAGGGGTGATCAGGTTTTTTTTACAA‑C‑TGTTTG‑AGAAGAAAAAAAAAAACCGCCGATcc                                                                 <  2:10839/90‑1 (MQ=255)
         gTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAA‑C‑TGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTc                                                             >  2:1457905/1‑90 (MQ=255)
          ttATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAA‑C‑TGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTcc                                                            >  2:907313/1‑90 (MQ=255)
            aTCGGCACCGGGGGTGATTTACGTTGGGGGTGATCGGGTTTTTTTTACAA‑T‑TTTTGG‑AGAAGAAAAAAAAAAACCGCCGATCCTGTCCAc                                                          <  1:256944/90‑1 (MQ=255)
             tCGGCAACTGGGGTTTTTTACTTTGGTGGTGATCGGGTTTTTTTTACAC‑T‑TGTTTG‑AAAAGAAAAAAAAAAACCGCCGATCCTGTCCAcc                                                         <  2:1392731/90‑1 (MQ=255)
                  agcTCGGTTTATTTACGTTAGTGGTGATCGGGTTTTTTTTACAC‑C‑TTTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCAtt                                                    <  2:335539/88‑1 (MQ=255)
                  aaCTGGGGTAATTTACGTTAGTGGTGATCAGGTTTTTTTTACAA‑C‑TTTTTG‑AAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCAtt                                                    <  2:1103786/90‑1 (MQ=255)
                  aaCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAA‑C‑TGTTTG‑AGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCAtt                                                    <  2:204894/90‑1 (MQ=255)
                   aCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAA‑C‑TGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTa                                                   >  2:1012247/1‑90 (MQ=255)
                   aCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAA‑C‑TGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTa                                                   >  1:169832/1‑90 (MQ=255)
                    cgggggTGTTTTACGTTAGTGGTAATCAGTTTTTTTTAAAAA‑T‑TTTTTG‑AGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTAc                                                  <  1:907313/87‑1 (MQ=255)
                       gggTGATTTCCTTTGGGGGTGA‑CCGGTTTTTTTTAAAACC‑TTTTTG‑AAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGc                                               <  1:504265/90‑1 (MQ=255)
                         gttttttCCGTTAGTGGTGATCAGGTTTTTTTTACAA‑C‑TTTTTG‑AAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa                                             <  2:1488369/86‑1 (MQ=255)
                         gttaTTTACGTTAGGGGTGATCGGGTTTTTTTTACAA‑C‑TTTTTG‑AAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa                                             <  1:1186488/87‑1 (MQ=255)
                         gtgtTTTAGGTTGGTGGTGATC‑GGCTTTTTTTTAAA‑CTTTTTTG‑AGAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCaa                                             <  1:1248863/90‑1 (MQ=255)
                         gtgtTTTACGTTAGTGGTGATCAGGTTTTTTTTACAA‑T‑TGTTTG‑AGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa                                             <  2:991450/90‑1 (MQ=255)
                         gTGATTTACGTTGGGGGTGATCAGTTTTTTTTTACAA‑T‑TTTTGG‑AGAAAAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCaa                                             <  1:1012247/90‑1 (MQ=255)
                         aggATTTACGTTAGTGGTGATCGGTTTTTTTTAACAA‑T‑TTTTTG‑AAAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCaa                                             <  1:1262781/88‑1 (MQ=255)
                            ttttACGTTAGTGGTGATCGGTTTTTTTTTAC‑A‑C‑TTTTTGAAAAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGt                                          <  1:1457905/89‑1 (MQ=255)
                            aTTTACTTTAGGGGTCACCAGGTTTTTTTTCCAT‑T‑TTTTTA‑AAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGt                                          <  1:1262222/90‑1 (MQ=255)
                              ttACGTTAGTGGTGATCAGGCTTTTTTTACAA‑C‑TGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGtct                                        >  1:335539/1‑88 (MQ=255)
                              cagcGTTAGTGGTGATCGGTTTTTTTTTACAC‑C‑TTTTTG‑AAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg                                        <  1:383823/87‑1 (MQ=255)
                                     aGTGGTGATCAGGCTTTTTTTACAA‑C‑TGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAcct                                 >  2:1262781/1‑88 (MQ=255)
                                     aGTGGTGATCAGGCTTTTTTTACAA‑C‑TGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACaa                                 >  2:891260/1‑90 (MQ=255)
                                      gtggtgATCAGGCTTTTTTTACAA‑C‑TGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAg                                >  1:1204525/1‑90 (MQ=255)
                                           gATCAGGCTTTTTTTACAA‑C‑TGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGAcctg                           >  2:383823/1‑88 (MQ=255)
                                                ggCTTTTTTTACAA‑C‑TGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcgggc                      >  1:1596959/1‑88 (MQ=255)
                                                ggCTTTTTTTACAA‑C‑TGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTc                      >  1:646040/1‑90 (MQ=255)
                                                     tttttACAA‑C‑TGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAg                 >  2:1305503/1‑90 (MQ=255)
                                                                 tttttG‑AAAAAAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCt       <  2:795560/88‑1 (MQ=255)
                                                                      g‑agaagaaaaAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAg  <  1:1525641/90‑1 (MQ=255)
                                                                             agaaaagaaaACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc                          >  1:107447/1‑60 (MQ=255)
                                                                              |                                                                                  
CCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAA‑C‑TGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAG  >  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 < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑