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

GGCAGTGTCATACTACTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACA  >  NZ_CP009273/997360‑997518
                                                                                       |                                                                       
ggCAGTGTCATACTACTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAg                                                                       >  1:122650/1‑90 (MQ=255)
         aTACTACTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACa                                                              <  2:104246/90‑1 (MQ=255)
              actCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGGTTAAGAAAAAttt                                                         >  2:579988/1‑90 (MQ=255)
                 cAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAAACAAAATTTAt                                                      >  2:113565/1‑90 (MQ=255)
                  aGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATc                                                     <  2:504691/90‑1 (MQ=255)
                   ggggAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTTAAGACAAAATTTATCa                                                    >  2:127641/1‑90 (MQ=255)
                        aTCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAATTTTTTCAATGGa                                               >  1:550548/1‑90 (MQ=255)
                        aTCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGAAAAATTTTATCAGTGGa                                               >  1:389564/1‑90 (MQ=255)
                        aTCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGGTAAAGAAAAATTTTTTCAGGGGa                                               >  2:141156/1‑90 (MQ=255)
                                   ggCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAcgaatcgaat                                    <  1:205260/90‑1 (MQ=255)
                                                                ggaAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATaa       <  1:155029/90‑1 (MQ=255)
                                                                    gaagCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGAc   <  2:102901/90‑1 (MQ=255)
                                                                     aagCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACa  <  2:39164/90‑1 (MQ=255)
                                                                                       |                                                                       
GGCAGTGTCATACTACTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACA  >  NZ_CP009273/997360‑997518

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