Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,267,2830AT55.6% ‑6.4 / 17.2 18V8E (GTA→GAA) chaAsodium‑potassium/proton antiporter ChaA
Reads supporting (aligned to +/- strand):  ref base A (1/7);  new base T (10/0);  total (11/7)
Fisher's exact test for biased strand distribution p-value = 2.51e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.34e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAG  >  NZ_CP009273/1267202‑1267372
                                                                                 |                                                                                           
gtCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCGCccct                                                                                     >  2:353731/1‑90 (MQ=255)
gtCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCCCctct                                                                                     >  2:110492/1‑90 (MQ=255)
gtCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCCCccct                                                                                     >  1:1115836/1‑90 (MQ=255)
  cTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCg                                                                                          >  2:532924/1‑83 (MQ=255)
  cTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCg                                                                                          <  1:532924/83‑1 (MQ=255)
  cTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCGCCCCttt                                                                                   >  1:391051/1‑89 (MQ=255)
  cTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCCCCTCttt                                                                                   >  1:992672/1‑89 (MQ=255)
     gctgctTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTATGGCGGGTTTTTACCGCCTCTTGAGc                                                                                <  1:238846/90‑1 (MQ=255)
         ctTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCAttt                                                                           <  1:594347/90‑1 (MQ=255)
               cacaGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACAtt                                                                     <  1:1216009/90‑1 (MQ=255)
                acaGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTa                                                                    <  1:785428/90‑1 (MQ=255)
                         gCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCCCCCCTTTTGCATTTTGCATTATGGTTTTcc                                                          >  1:94311/1‑90 (MQ=255)
                              aCCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCGCCCCTTGAGCATTTTACATTATTGTTATTCCCttt                                                      >  1:404913/1‑90 (MQ=255)
                                 aGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCCCCCCTTTTGCATTTTTAATTTTTGGTTTTCCCtttcag                                                   >  1:676013/1‑87 (MQ=255)
                                 aGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCCCCCCTTTTGCCTTTTACATTATTGTTTTTCCTTTGcct                                                  >  2:1099238/1‑88 (MQ=255)
                                   cgcCAGAACCGGGGAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGAt                                                <  1:575791/90‑1 (MQ=255)
                                   cgcCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCCCCCCTTTAGGATTTTACATTTTTGGTTTCCCTTTTCAGAt                                                >  1:194947/1‑90 (MQ=255)
                                                 aaaTAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGaaaa                                  <  1:338632/90‑1 (MQ=255)
                                                                                 aCCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAg  <  2:533820/90‑1 (MQ=255)
                                                                                 |                                                                                           
GTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAG  >  NZ_CP009273/1267202‑1267372

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