Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,673,5690AC63.2% 5.4 / 15.0 19T296P (ACC→CCC) ydgHDUF1471 family protein YdgH
Reads supporting (aligned to +/- strand):  ref base A (6/1);  new base C (0/12);  total (6/13)
Fisher's exact test for biased strand distribution p-value = 2.58e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.51e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCT  >  NZ_CP009273/1673493‑1673655
                                                                            |                                                                                      
cTACGGAAAAATGCCCGAAGTCTCCTATCAGGTTGCGAAACGGGCCGCGAAGAAAGGTGCAAAGTATTACCACACCCCCCGCCAGTGGCa                                                                           <  2:7602/90‑1 (MQ=255)
       accATGACGGACTTCTCCTCTCGGGTGGGGAAGCGGCGCGCGAAAAAAGTTGCTAATTTTTCCCACCCCCCCCGCCAGTGGCAGGAACGt                                                                    <  1:494470/88‑1 (MQ=255)
       aaCATGCCCAAAGTCTCCTATCGGGTTGGAAAACGTCCCGAGAAAAAAGGTGCAAATTTTTCCCACACCCCCCGCCAGTGGCAGGAACGt                                                                    <  1:150352/90‑1 (MQ=255)
           tgcccGAAGTCTCCTATCAGTTTGCGAAACGTGCCGCGAAAAAAGGTGCAAATTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGta                                                                <  2:353922/87‑1 (MQ=255)
            gACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGtaa                                                               >  1:560626/1‑90 (MQ=255)
                         tATCGGGTTCCGAAACGCGCGGAGAAAAAAGGTGCTAAGTTTCCCCCCACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc                                                  <  1:83604/90‑1 (MQ=255)
                         tATCAGGTGGCGAAACGTGCCGAAAAGAAGGGTGCTAAGTTTTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc                                                  <  1:552196/90‑1 (MQ=255)
                              ggTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGGAATAACCTGACCGTCAgcgc                                             >  1:234349/1‑90 (MQ=255)
                              ggTTGCAAAGCGTGCCGGGAAGAAGGGTGCAAAGTTTTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc                                             <  1:414045/90‑1 (MQ=255)
                              ggTTGCAAACCGGGCCGCGAAAAAAGGTGCTAATTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc                                             <  2:212819/90‑1 (MQ=255)
                                     aaaCGTGCCGAGAAGAAGGGTGCTAAGTTTTCCCCCCTCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTg                                      <  1:12374/90‑1 (MQ=255)
                                             cgagaagaaGGGTGCTAAGTTTTCCCCCCCCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTataaata                              <  1:494725/87‑1 (MQ=255)
                                             cgagaaaaaAGTGGCGAAGTATTCCCCCATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTataaata                              <  1:172726/87‑1 (MQ=255)
                                                   aaaaGGTGATAAGTATTCCCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTAtt                        <  2:620705/89‑1 (MQ=255)
                                                           cTAAGTATTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACt                <  2:341896/90‑1 (MQ=255)
                                                                  ttACCACATCACCCGCCAGTGGCAGGAACGTGGGAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATaa         >  1:341896/1‑90 (MQ=255)
                                                                  ttACCACATCACCCGCCAGTGGCAGGAACGTGGGAATAACCTGACCGGCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATaa         >  2:474953/1‑90 (MQ=255)
                                                                       acaTCACCCGCCAGTGGCAGGAACGGGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAggcgg    >  1:55072/1‑90 (MQ=255)
                                                                         aTCACCCGCCAGTGGCAGGAACGTGGGAATAAACTGACCGTCAGCGCAGATCTGTATAAATAAGTATTCAGTTACTCCTATAATGCGGCt  >  1:209912/1‑90 (MQ=255)
                                                                            |                                                                                      
CTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCT  >  NZ_CP009273/1673493‑1673655

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