Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,819,2880CT60.0% ‑4.5 / 19.8 25intergenic (‑94/+109)ves/spyenvironmental stress‑induced protein Ves/ATP‑independent periplasmic protein‑refolding chaperone Spy
Reads supporting (aligned to +/- strand):  ref base C (7/3);  new base T (0/15);  total (7/18)
Fisher's exact test for biased strand distribution p-value = 2.50e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGC  >  NZ_CP009273/1819201‑1819371
                                                                                       |                                                                                    
tCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCtt                                                                                    <  2:536778/90‑1 (MQ=255)
tCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCtt                                                                                    <  2:268230/90‑1 (MQ=255)
    ggCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCtttttt                                                                                >  2:56847/1‑90 (MQ=255)
         gTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACaa                                                                           >  2:245358/1‑90 (MQ=255)
         gTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGGGATTTACGTTAGTGGTGATCAGGCTTTTTTTACaa                                                                           >  1:647151/1‑90 (MQ=255)
                             tCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaaga                                                      >  1:171346/1‑90 (MQ=255)
                                                 aTCGGCACCTGGGGTGTTTTACTTTAGTGGTGATCGGTTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCAc                                  <  1:212392/90‑1 (MQ=255)
                                                  ccGGCAACTGGGGGATTTTAAGTTAGGGGTGATCGGGTTTTTTTTAAACCTGTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCAcc                                 <  1:687086/89‑1 (MQ=255)
                                                        aCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTa                           >  2:106117/1‑90 (MQ=255)
                                                              gttaTTTACGTTAGTGTTGACCAGTTTTTTTTTACAACTTTTTGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa                     <  2:119699/87‑1 (MQ=255)
                                                              gtgtTTTACTTTAGTGGTGATCAGGTTTTTTTTACAACTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa                     <  2:638712/90‑1 (MQ=255)
                                                              gtgtTTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa                     <  2:313130/90‑1 (MQ=255)
                                                              gtattttACGTTGGTGGTGATCGGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa                     <  2:601411/86‑1 (MQ=255)
                                                              gTGATTTCCGTAAGGGGTGACCAGCTTTTTTTTACAACTTTTGGAAAAAAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCaa                     <  1:391707/90‑1 (MQ=255)
                                                              gTGATTTACTTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa                     <  2:37933/90‑1 (MQ=255)
                                                              gTGATTTACTTTAGGGGTGATCGGTTTTTTTTTACAATTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa                     <  1:66680/90‑1 (MQ=255)
                                                              gTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa                     <  2:130156/90‑1 (MQ=255)
                                                               tGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAg                    <  2:272135/90‑1 (MQ=255)
                                                                   ttATGTTGGTGGTGATCGGTTTTTTTTTAAATTTGTTGGAAAAAAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg                <  1:606297/90‑1 (MQ=255)
                                                                   ttACTTTGGGGGTGACAGGTTTTTTTTAAAAATTTTTGGAAAAAAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg                <  1:516571/90‑1 (MQ=255)
                                                                   ttACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGGAg                >  2:604359/1‑90 (MQ=255)
                                                                      ctttAGTGGGGATCAGGTTTTTTTTACAACTGTTTGACAAGATAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTgg             <  2:300580/88‑1 (MQ=255)
                                                                          aGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACaa         >  2:545141/1‑90 (MQ=255)
                                                                                 agcAGGTTTTTTTTAAAACTTTTTGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc  <  2:301987/88‑1 (MQ=255)
                                                                                 aTCAGGTTTTTTTAACACTTTTTTGAAAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc  <  2:457805/90‑1 (MQ=255)
                                                                                       |                                                                                    
TCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGC  >  NZ_CP009273/1819201‑1819371

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