Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,259,5430TG52.9% 15.7 / 18.0 17F68V (TTT→GTT) aphAacid phosphatase AphA
Reads supporting (aligned to +/- strand):  ref base T (3/5);  new base G (9/0);  total (12/5)
Fisher's exact test for biased strand distribution p-value = 9.05e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-03
Rejected as consensus: Frequency below/above cutoff threshold.

GGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAG  >  NZ_CP009273/4259457‑4259612
                                                                                      |                                                                     
ggCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTg                                                                    >  1:874356/1‑90 (MQ=255)
          aCAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATATCGATGa                                                          >  2:402867/1‑90 (MQ=255)
             ggCACCCATTCATTGGGTTTCGGTCGCAAAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGacac                                                       <  2:290378/90‑1 (MQ=255)
                aCCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGt                                                    <  1:882369/90‑1 (MQ=255)
                  ccATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGGGGGGGTTGATATCGATGACACGGGAc                                                  >  1:743479/1‑90 (MQ=255)
                     ttcattGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACttt                                               >  2:639803/1‑90 (MQ=255)
                       cattGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATATCGATGACACGGTACttttt                                             >  2:400232/1‑90 (MQ=255)
                           gggTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATATCGATGACACGGTACTTTTTTCCa                                         >  1:224061/1‑90 (MQ=255)
                           gggTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGGGGGGTTTGATATCGATGACACGGGACTTTTTTTCa                                         >  1:793248/1‑90 (MQ=255)
                                             ttGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGcgcg                       <  1:639803/90‑1 (MQ=255)
                                                     aGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCaaaaaa               >  2:123065/1‑90 (MQ=255)
                                                      gCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATATCGATGACACCGGACTTTTTTTCAGTCCGGGCTTTTGGGCCGGCcaaaaaa              >  1:193320/1‑90 (MQ=255)
                                                              gggCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCg      <  2:874356/90‑1 (MQ=255)
                                                                gCGTCCGCCAATGGCGGTGGGGGTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGcc    >  2:194876/1‑90 (MQ=255)
                                                                 cGTCCGCCAATGGCGGTGGGGGTTTATATCGATGACACGGGACTTTTTTTCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACTTCTCGCCa   >  1:917124/1‑90 (MQ=255)
                                                                  gTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAg  <  1:758729/90‑1 (MQ=255)
                                                                  gTCCGCCAATGGCGGTGGGGGTTGATATCGATGACACGGTACTTTTTTTCAGTCCGGGCTTCTTGCGCGGCAAAAAAACCTTCTCGCCAg  >  1:537034/1‑90 (MQ=255)
                                                                                      |                                                                     
GGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAG  >  NZ_CP009273/4259457‑4259612

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