Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,594,9200AC61.1% 9.6 / 12.3 18V176G (GTG→GGG) ftsXpermease‑like cell division protein FtsX
Reads supporting (aligned to +/- strand):  ref base A (4/3);  new base C (0/11);  total (4/14)
Fisher's exact test for biased strand distribution p-value = 1.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.02e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCA  >  NZ_CP009273/3594831‑3595001
                                                                                         |                                                                                 
tccGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATca                                                                                   <  1:439944/90‑1 (MQ=255)
                          tGCGTGATCCGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGttt                                                         <  2:139872/90‑1 (MQ=255)
                               gATACGATCCCGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCCCCACCGCCACTGCCGGAAGCGGGTTttcttc                                                    <  2:9885/90‑1 (MQ=255)
                               gATACGACCCCGCAGCGTATTCAGTGATCCCGCCCCCGGAAAATCGAGTTCCGGGACCCCCACCGCCACTGCCGGAAGCGGGTTttcttc                                                    <  1:862660/90‑1 (MQ=255)
                                         cGCAGCTTATTCAGTGATTCCGCCCCCTGAAAATCGATTTTCGGGATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATcc                                          <  2:87218/90‑1 (MQ=255)
                                             gCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAgcg                                      >  1:666550/1‑90 (MQ=255)
                                               gTATTCAGTGATTCCGCCCCCTGGAAATCGATTTTCGGGATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGca                                    <  2:224971/90‑1 (MQ=255)
                                               gTATTCAGGGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGca                                    >  2:620708/1‑90 (MQ=255)
                                                tATTCAGTGTTTCCGCCCCCTGGAAATCGAGTTTCGGGATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGcac                                   <  2:381206/90‑1 (MQ=255)
                                                tATCCAGTGTTCCCGCCCCCGGAAAATCGTGTTTCGGCCCCCCCCCCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGcac                                   <  1:473644/90‑1 (MQ=255)
                                                       tGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGcac                                   >  1:33961/1‑83 (MQ=255)
                                                       tGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGcac                                   <  2:33961/83‑1 (MQ=255)
                                                         aTTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAac                          >  1:224971/1‑90 (MQ=255)
                                                                ccccTGAAAATCGATTTTCGGGACCCCCCCCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAaccagacca                   <  1:282763/90‑1 (MQ=255)
                                                                ccccTGAAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAaccagacca                   <  2:209943/90‑1 (MQ=255)
                                                                          tCGATTTTGGGGATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAAc         <  2:562838/90‑1 (MQ=255)
                                                                              gTTTCGGGATCCCCCCCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCAc     <  2:248151/90‑1 (MQ=255)
                                                                                ttggggACCCCCCCCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCAccc   <  1:510247/87‑1 (MQ=255)
                                                                                 tggggATCCCCCCCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCa  <  1:160413/88‑1 (MQ=255)
                                                                                         |                                                                                 
TCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCA  >  NZ_CP009273/3594831‑3595001

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