Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,336,7850AC56.2% 12.4 / 20.1 16G288G (GGT→GGGyfaLAIDA‑I family autotransporter adhesin YfaL/EhaC
Reads supporting (aligned to +/- strand):  ref base A (7/0);  new base C (0/9);  total (7/9)
Fisher's exact test for biased strand distribution p-value = 8.74e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-04
Rejected as consensus: Frequency below/above cutoff threshold.

CTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTC  >  NZ_CP009273/2336709‑2336871
                                                                            |                                                                                      
ctGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGagg                                                                           >  1:211462/1‑90 (MQ=255)
  gttttGCCAAAAACCAGCTTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATAAAAAAACCCCCCCGCCGCAGaggaa                                                                         <  1:138434/87‑1 (MQ=255)
  gTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGaggaa                                                                         >  1:476844/1‑90 (MQ=255)
       gCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGAcc                                                                    >  1:363829/1‑90 (MQ=255)
       gCCAAAACCCAGTTTTTTCCCGGCGGCAATATCAAAGGAAATTTCCCTTAACCCGAAAAACAAAAACCCCCCCGCCGCAGAGGAAGGAcc                                                                    <  1:465423/90‑1 (MQ=255)
         aaaaaaCCAGCTTTTTCCCGCCGGCAATAGAAAAGGAACCTTCACTAAACCCGAAATAAAAAACCCCCCCCGCCGCAGAGGAAGGACCAt                                                                  <  1:535778/86‑1 (MQ=255)
          aaTAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCAtc                                                                 >  2:535778/1‑90 (MQ=255)
              aCCGGCTTTTTCCCGTCGGAAATATAAAAGGTAACTTCACTTAAGCCGAGATACAAAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCa                                                             <  2:476844/90‑1 (MQ=255)
                      ttttCCGTCGACAATAACAAAGGAAACTTCCCTAAAGCCGAGAAACAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc                                                     <  1:94033/90‑1 (MQ=255)
                      ttttCCGCCGGAAATAACAAAGGTAACTTCCCTAAAGCCGAAATACAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc                                                     <  1:553252/90‑1 (MQ=255)
                      tttcCCGTCGGAAATATAAAAGGTAACTTCATTTAAGCCGAAATACAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc                                                     <  2:363829/90‑1 (MQ=255)
                      tttcCCGCCGGAAATACCAAAGGTACCTCCACTAAAGCCAAAAAACAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc                                                     <  1:519147/90‑1 (MQ=255)
                       tttCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCctgct                                                    >  2:491890/1‑90 (MQ=255)
                         cccgccgGAAATATAAAAGGAAACTCCACTAAACCCGAAATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCtgctgc                                                  <  1:257800/89‑1 (MQ=255)
                                                              agaTACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTgg             >  2:171528/1‑90 (MQ=255)
                                                                         gCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGtc  >  1:26651/1‑90 (MQ=255)
                                                                            aCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTgcgc                                                >  1:534066/1‑41 (MQ=255)
                                                                            |                                                                                      
CTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTC  >  NZ_CP009273/2336709‑2336871

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