Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,673,5610AC56.2% 7.6 / 12.9 16Y293S (TAC→TCC) ydgHDUF1471 family protein YdgH
Reads supporting (aligned to +/- strand):  ref base A (4/3);  new base C (0/9);  total (4/12)
Fisher's exact test for biased strand distribution p-value = 1.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.80e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAA  >  NZ_CP009273/1673493‑1673648
                                                                    |                                                                                       
cTACGGAAACATGCCCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAAAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCa                                                                    <  1:475943/90‑1 (MQ=255)
       aaCATGACCGAATTCTCCTATCAGGTTGCAAAACGGGCCGGGAAAAAGGGTGCTAATTTTTCCCACACCCCCCGCCAGTGGCAGGAACGt                                                             <  2:470997/90‑1 (MQ=255)
       aaCATGACCGAAGTCCCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTCCCACATCACCCGCCAGTGGCAGGAACGt                                                             <  1:235836/90‑1 (MQ=255)
       aaCAGGCCCGAAGTCTCCTATCGGGTTGCGAACCGTCCCGCGAAAAAAGGTGCTAAGTATTCCCACATCCCCCGCCAGTGGCAGGAACGt                                                             <  2:259454/90‑1 (MQ=255)
                         agacAGGTTGCGAAACGGGCCGCGAAGAAAGGTGCTAAGTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc                                           <  2:104325/87‑1 (MQ=255)
                              ggTTGCAAAACGTGCCGAAAAAAAAGGTGCTAAGTTTTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc                                      <  1:278803/90‑1 (MQ=255)
                                 tGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGa                                   >  2:153696/1‑90 (MQ=255)
                                       aCGTGCCGGACAGAAAGGTGCAAAGTTTTCCCACTTCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGta                             <  1:405107/90‑1 (MQ=255)
                                             cgcaAAAAAGGGTGCTAATTTTTCCCCCACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTataaata                       <  2:148753/90‑1 (MQ=255)
                                                  agaaAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTAt                  <  1:153696/90‑1 (MQ=255)
                                                   gaaAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc                                           <  1:247444/64‑1 (MQ=255)
                                                   gaaAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc                                           >  2:247444/1‑64 (MQ=255)
                                                      aGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAg              >  1:148753/1‑90 (MQ=255)
                                                        gTGCTAAGTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGtt            <  2:292333/90‑1 (MQ=255)
                                                          gCTAAGTTTTCCCCCCCCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTAc          <  2:430328/90‑1 (MQ=255)
                                                                  ttACCACATCACCCGCCAGTGGAAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATaa  >  1:293006/1‑90 (MQ=255)
                                                                    |                                                                                       
CTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAA  >  NZ_CP009273/1673493‑1673648

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