Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,849,3410TC57.1% 11.8 / 20.7 21G109G (GGA→GGGyidHYidH family protein
Reads supporting (aligned to +/- strand):  ref base T (8/1);  new base C (0/12);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 4.42e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.76e-04
Rejected as consensus: Frequency below/above cutoff threshold.

ACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATC  >  NZ_CP009273/3849261‑3849428
                                                                                |                                                                                       
aCGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACg                                                                                >  2:238972/1‑90 (MQ=255)
aCGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACg                                                                                >  1:209338/1‑90 (MQ=255)
           gtgCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCAc                                                                     >  2:1797/1‑90 (MQ=255)
             gcgcTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGa                                                                   >  1:279479/1‑90 (MQ=255)
                  cgggCTGCAGCCCGGGGTGGGAAACGCGTCTGGTTTGGCGGCTACCCGGCATACAACCCCAGCCCCAAAACGATCACCGCCACGACCATc                                                              <  2:75989/88‑1 (MQ=255)
                  ccGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGACCATc                                                              <  1:80764/90‑1 (MQ=255)
                  ccGGCTGCAGTCCGGGGTCGGCAATGCGTTTGGTTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATc                                                              <  2:127302/90‑1 (MQ=255)
                  ccGGCTGAAGTCGGGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACCCCAGCCCAAAACCGATCACCGCCACGACCATc                                                              <  1:238972/90‑1 (MQ=255)
                  ccGGCTGAAGCCCGGGGTAGGCAATGCTTCTGTTTTGGCGGCTACCCGGAAAAAAACACCAGCCCAATAACGACCACCGCCACGACCATc                                                              <  2:188838/90‑1 (MQ=255)
                       ggCAGTCGGGGGTGGGAAATGCTTTTGGTTTGGCGGCTACCCGGAAAACAACACCAGCCCAAAAACGATCACCGCCACGACCATCAatat                                                         <  2:176838/89‑1 (MQ=255)
                                 ggTCGGAAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATg                                               <  1:280490/90‑1 (MQ=255)
                                   ccGGCAATGCTTCTGTCTTGGCGGCTACCCGGAAAACAACCCCAGCCCAAAAACGATCACCGCCACGACCATCAATATTAAGCTAATGAt                                             <  2:159787/89‑1 (MQ=255)
                                     ggCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATta                                           >  1:159787/1‑90 (MQ=255)
                                     ggCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATta                                           >  2:197025/1‑90 (MQ=255)
                                     ggCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATta                                           >  2:112323/1‑90 (MQ=255)
                                          tGCGTCTGGCTTTGCGGCTACCCGGCATACAACACCAGCCCCAAAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAg                                      <  2:87076/90‑1 (MQ=255)
                                               cgggCTTTGCGGCTACCCGGCATAAAACACCAGCCCCAAAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGtt                                 <  1:157930/88‑1 (MQ=255)
                                               cTGGTTTTGCGGCTACCCGGCAAAAAACCCCAGCCCAAAAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGtt                                 <  2:279479/90‑1 (MQ=255)
                                               cTGGCTTTGCGGCTACCCGGCAAACAACACCAGCCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGtt                                 <  1:283415/90‑1 (MQ=255)
                                                               ccGGCAAACACCCCCAGCCCCAAAACGATCCCCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTcc                 <  1:268128/90‑1 (MQ=255)
                                                                              aGTCCCATAACGATCATCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATc  >  1:184941/1‑90 (MQ=255)
                                                                                |                                                                                       
ACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATC  >  NZ_CP009273/3849261‑3849428

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