Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273997,4470TG62.5% 1.3 / 16.1 24V62G (GTA→GGA) elfGfimbrial protein
Reads supporting (aligned to +/- strand):  ref base T (4/5);  new base G (15/0);  total (19/5)
Fisher's exact test for biased strand distribution p-value = 2.96e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CATACTACTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACAATGTCAGAAGTGGCCG  >  NZ_CP009273/997368‑997534
                                                                               |                                                                                       
cATACTTCTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGAc                                                                               <  1:16668/90‑1 (MQ=255)
      actCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGGTAAAAACAAAatt                                                                         >  1:650145/1‑90 (MQ=255)
      actCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAAACAAAttt                                                                         >  2:456345/1‑90 (MQ=255)
      actCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAAAAAAAatt                                                                         >  1:185763/1‑90 (MQ=255)
       ctCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAAtttt                                                                        >  2:52791/1‑90 (MQ=255)
        tCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAATTTTa                                                                       >  1:84196/1‑90 (MQ=255)
        tCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAAATTTa                                                                       >  1:76057/1‑90 (MQ=255)
           ggggAATCCGTATATTTGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAAACAAATTTTATCa                                                                    >  1:441389/1‑90 (MQ=255)
           ggggAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTGGTTAAAGAAAAAATTTATCa                                                                    >  1:339890/1‑90 (MQ=255)
           ggggAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGGTAAAGAAAAATTTTAtta                                                                    >  1:500809/1‑88 (MQ=255)
             ggAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGt                                                                  <  1:156364/90‑1 (MQ=255)
                aTCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAAACAAAATTTTTCAGTGGa                                                               >  1:548943/1‑90 (MQ=255)
                aTCCGTATACTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAATTTTATCAGTGGa                                                               >  2:330486/1‑90 (MQ=255)
                          tGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGAAAAATTTTATCAGTGGAAcgaatcgaa                                                     >  1:55603/1‑90 (MQ=255)
                            gCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATg                                                   >  2:338490/1‑90 (MQ=255)
                              gTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAAGAATCGAATGat                                                 >  2:716639/1‑90 (MQ=255)
                                cAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGGTAAAAAAAAAATTTTTTAATGGGACGAAACGAATGatta                                               >  1:595120/1‑90 (MQ=255)
                                            tAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGAAAAATTTTTTCAAGGGAAAGAAACGAAAGATTATTATGTTTCCtg                                   >  1:710312/1‑90 (MQ=255)
                                                  ccTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTg                             <  1:600468/90‑1 (MQ=255)
                                                           agaagCAAATACGGCAGGGGGAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTTCTGTGATTGCGATAAAGa                    >  1:241761/1‑90 (MQ=255)
                                                             aagCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACa                  <  2:709596/90‑1 (MQ=255)
                                                                        gCAGGGGGAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACAATGTCAGAAGt       >  1:19916/1‑90 (MQ=255)
                                                                         cAGGGGGAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTTCTGTGATTGCGATAAAGACAATGTCAGAAGTg      >  1:284166/1‑90 (MQ=255)
                                                                             ggTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACAATGTCAGAAGTGGCCg  <  1:52791/90‑1 (MQ=255)
                                                                               |                                                                                       
CATACTACTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACAATGTCAGAAGTGGCCG  >  NZ_CP009273/997368‑997534

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