Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,673,5610AC61.9% 14.3 / 18.6 21Y293S (TAC→TCC) ydgHDUF1471 family protein YdgH
Reads supporting (aligned to +/- strand):  ref base A (6/2);  new base C (0/13);  total (6/15)
Fisher's exact test for biased strand distribution p-value = 5.16e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-04
Rejected as consensus: Frequency below/above cutoff threshold.

GCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAA  >  NZ_CP009273/1673489‑1673648
                                                                        |                                                                                       
gCACCTACGGCAACATGACCGAATTCTCCTATCAGGTGGCGAAACGTGCCGCGAAGAAGGGTGATAAGTATTCCCACACCCCCCGCCAGt                                                                        <  1:564729/90‑1 (MQ=255)
         gCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTa                                                                                         >  2:576677/1‑64 (MQ=255)
         gCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTa                                                                                         <  1:576677/64‑1 (MQ=255)
           aaCATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTTTTCCCACATCCCCCGCCAGTGGCAGGAACGt                                                             <  2:686893/90‑1 (MQ=255)
            aCATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACggg                                                            >  1:67890/1‑88 (MQ=255)
            aCATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTg                                                            >  1:462117/1‑90 (MQ=255)
              aTGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGt                                                          >  2:755770/1‑90 (MQ=255)
               tGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTTTTCCCACATCCCCCGCCAGTGGCAGGAACGTGGta                                                         <  2:394019/90‑1 (MQ=255)
                                  ggTTGCGAACCGTGCCGCAAAAAAAGGTGCTAAGTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc                                      <  2:462117/90‑1 (MQ=255)
                                  ggTTGCAAACCGGGCCGGTAAAAAAGGTAGAAAGTTTTCCCTCCCCCCCCGCCAGTGGCAGGAACGTGGTAATAGCCTGACCGTCAgcgc                                      <  1:84173/90‑1 (MQ=255)
                                  ggTTGCAAAACGTGCCGCAAAAAAAGGTGCTAAGTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc                                      <  2:67890/90‑1 (MQ=255)
                                     tGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGGGGTAATAACCTGACCGTCAGCGCCGa                                   >  1:255961/1‑90 (MQ=255)
                                       caaaaCGTGCGGAGAAAAAGGGTGCAAATTTTTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATc                                 <  1:221352/88‑1 (MQ=255)
                                                ccGAGAAAAAAGGTGCTAATTTTTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTataaat                        <  1:517250/90‑1 (MQ=255)
                                                     aagaaAGGTGCTAATTATTACCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTa                   <  2:495515/90‑1 (MQ=255)
                                                            gTGCTAATTTTTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGtt            <  1:542459/90‑1 (MQ=255)
                                                            gTGCTAAGTTTTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGtt            <  1:120554/90‑1 (MQ=255)
                                                            gTGCTAAGTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGtt            <  2:612899/90‑1 (MQ=255)
                                                              gatAAGTTTTCCCACCCCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTAc          <  2:617403/88‑1 (MQ=255)
                                                               acaggTATTCCCCCACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACt         <  2:689092/86‑1 (MQ=255)
                                                                 aaGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc                                      <  2:441660/59‑1 (MQ=255)
                                                                 aaGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc                                      >  1:441660/1‑59 (MQ=255)
                                                                      ttACCACATCACCCGCCAGTGGCAGGAACAAGGGAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACAATTCAGGGACTTCtatta  >  2:402721/1‑88 (MQ=255)
                                                                        |                                                                                       
GCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAA  >  NZ_CP009273/1673489‑1673648

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