Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273721,8000GA57.9% ‑2.4 / 15.4 19V231V (GTC→GTTkdpBpotassium‑transporting ATPase subunit KdpB
Reads supporting (aligned to +/- strand):  ref base G (1/7);  new base A (11/0);  total (12/7)
Fisher's exact test for biased strand distribution p-value = 1.59e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCC  >  NZ_CP009273/721713‑721888
                                                                                       |                                                                                        
gacCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAAAc                                                                                        >  2:681237/1‑90 (MQ=255)
     gcagcGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAAACGAAAg                                                                                   >  1:54911/1‑90 (MQ=255)
         cgcCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAg                                                                               <  1:685249/90‑1 (MQ=255)
          gcCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAgg                                                                              >  2:488126/1‑90 (MQ=255)
           ccaccaGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAggg                                                                             <  1:218307/90‑1 (MQ=255)
                     cGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAAACGATAGTCAGGGCAATcagcag                                                                   >  2:415006/1‑90 (MQ=255)
                        tACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAAACGATAGGCAGGGGAATCAGCCGAAt                                                                >  2:536904/1‑90 (MQ=255)
                           gCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCCGTAAAAAAACCAAAGGCAGGGGCATCAACCGGATggg                                                             >  1:416510/1‑89 (MQ=255)
                            cTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTc                                                            <  1:57902/90‑1 (MQ=255)
                            cTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTc                                                            <  1:199525/90‑1 (MQ=255)
                                          gCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAAACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATttcgtt                                              >  2:181558/1‑90 (MQ=255)
                                          gCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAAACGATAGTCAGGGCAATCAGCAGAATGGGGAGGGCAATttcgtt                                              >  1:106271/1‑90 (MQ=255)
                                          gCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAAACGAAAGTCAGGGCAAACAGCAGAAAGGGCAGGGCAATtttgtt                                              >  1:163692/1‑90 (MQ=255)
                                          gCCCCACGCGGAAAACGACCACGGCGGGGCGGTTGCCCGAAAAAAAACGAAAGCAAGGGGAATCCGGCAAATGGGGAGGGGAAttttttt                                              >  1:664002/1‑86 (MQ=255)
                                                                   gTGGCGGTTGCCAGTAAAAAAACGATAGTCAGGGCAAACAGCAGAAAGGGCAGGGCAATTTCGTTCGGGGGTTTTCGGCGCTGTGCGccc                     >  1:304960/1‑89 (MQ=255)
                                                                      gCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTcc                  <  2:641860/90‑1 (MQ=255)
                                                                           tGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATCTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCAt             <  2:139824/90‑1 (MQ=255)
                                                                           tGCCAGTAAAAAAACGATAGTCAGGGCAATCAGCAGAATGGGCAGGGCAATTTCGTTCGGCGGTTTTCGTCGCTGTGCGCCTTCCACCAt             >  1:220107/1‑90 (MQ=255)
                                                                                      aGACGAAAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATcc  <  1:42370/90‑1 (MQ=255)
                                                                                       |                                                                                        
GACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCC  >  NZ_CP009273/721713‑721888

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