Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,673,5610AC60.0% 9.3 / 13.8 15Y293S (TAC→TCC) ydgHDUF1471 family protein YdgH
Reads supporting (aligned to +/- strand):  ref base A (5/1);  new base C (0/9);  total (5/10)
Fisher's exact test for biased strand distribution p-value = 2.00e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCC  >  NZ_CP009273/1673476‑1673643
                                                                                     |                                                                                  
aTCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACcac                                                                                >  1:467900/1‑90 (MQ=255)
            ggCAACTACGGAAACATGACCGAAGTCCCCTATCAGGTTGCGAAACGTGCCGAGAAAAAGGGTGCTAATTATTCCCCCACCCCCCGCCAg                                                                    <  1:613608/90‑1 (MQ=255)
                              aCCGAATTCTCCTATCAGGTGGCAAACCGTGCCGCGAAAAAAGGTGATAAGTATTCCCACACCCCCCGCCAGTGGCAGGAACGTGGtaat                                                  <  1:1008218/90‑1 (MQ=255)
                                 gAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAAc                                               >  2:1008218/1‑90 (MQ=255)
                                               ggTTGCAAACCGTGCCGGAAAAAAGGGTGCTAAGTTTTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc                                 <  2:642975/90‑1 (MQ=255)
                                               ggTTGCAAAACGTCCCGCAAAAAAAGGTGCTAATTTTTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc                                 <  1:287837/90‑1 (MQ=255)
                                                 ttGCGAACCGTGCCGCGAAAAAAGGTGCTAAGTATTACCCCATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAg                               <  2:712575/90‑1 (MQ=255)
                                                  tGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGGACGTGGTAATAACCTGACCGTCAGCGCAGa                              >  1:909641/1‑90 (MQ=255)
                                                  tGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGGAATAACCTGACCGTCAGCGCAGa                              >  1:712575/1‑90 (MQ=255)
                                                    cGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTCCCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATc                            <  2:406584/90‑1 (MQ=255)
                                                    cGAAACGTCCCGCGAAAAAGGGTTAAAAGTTTTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATc                            <  1:1038565/90‑1 (MQ=255)
                                                                    aaagggTGCAAAGTATTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTAtt            <  2:467900/86‑1 (MQ=255)
                                                                         gTGCTAATTATTCCCACACCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGtt       <  2:970081/90‑1 (MQ=255)
                                                                           gatAAGTTTTCCCCCACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTAc     <  1:1010974/88‑1 (MQ=255)
                                                                              aaGTATTACCACATCACCCGCCAGTGGCAGGAACGGGGTAATAAACTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTcc  >  1:880799/1‑90 (MQ=255)
                                                                                     |                                                                                  
ATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCC  >  NZ_CP009273/1673476‑1673643

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: