Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273997,4470TG56.2% 12.2 / 18.7 16V62G (GTA→GGA) elfGfimbrial protein
Reads supporting (aligned to +/- strand):  ref base T (1/6);  new base G (9/0);  total (10/6)
Fisher's exact test for biased strand distribution p-value = 8.74e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.98e-03
Rejected as consensus: Frequency below/above cutoff threshold.

ACTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACAATGTCAGAAGTG  >  NZ_CP009273/997374‑997530
                                                                         |                                                                                   
actCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTGGTTAAAAACAAAatt                                                                     >  2:583343/1‑90 (MQ=255)
actCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTTAAAAAAAAatt                                                                     >  2:238192/1‑90 (MQ=255)
actCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAAttt                                                                     >  1:15821/1‑90 (MQ=255)
actCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGGTAAAGAAAAAatt                                                                     >  2:634998/1‑90 (MQ=255)
    aGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAACCCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATc                                                                 <  1:634998/90‑1 (MQ=255)
     ggggAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAAACAAAATTTATCa                                                                >  2:91531/1‑90 (MQ=255)
          aTCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTTAAAACAAAATTTATCAGTGGa                                                           >  2:286856/1‑90 (MQ=255)
          aTCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAAAAAAATTTTTTCAGTGGa                                                           >  2:602013/1‑90 (MQ=255)
                   ttGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAAAAAAAATTTATCAGGGGGAagaaacga                                                  >  2:538202/1‑90 (MQ=255)
                     ggCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAcgaatcgaat                                                <  1:91531/90‑1 (MQ=255)
                                                 aggaAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATa                    <  1:538202/90‑1 (MQ=255)
                                                    aagaagCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGAt                     <  1:179668/86‑1 (MQ=255)
                                                    aagaagCAAATACGGCAGGGGGAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGAt                     >  2:179668/1‑86 (MQ=255)
                                                       aagCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACa              <  1:286856/90‑1 (MQ=255)
                                                              tACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACAATGTCAg       <  2:689598/90‑1 (MQ=255)
                                                                   cAGGGGGAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACAATGTCAGAAGTg  >  2:346969/1‑90 (MQ=255)
                                                                         |                                                                                   
ACTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACAATGTCAGAAGTG  >  NZ_CP009273/997374‑997530

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