Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,336,7780TA56.2% 8.3 / 19.6 16M291L (ATG→TTG) yfaLAIDA‑I family autotransporter adhesin YfaL/EhaC
Reads supporting (aligned to +/- strand):  ref base T (7/0);  new base A (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: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CAACAGCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGG  >  NZ_CP009273/2336691‑2336860
                                                                                       |                                                                                  
cAACAGCTCCGTCATTCTCTGTATTGCCAATGACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATaa                                                                                  >  2:465423/1‑90 (MQ=255)
                  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)
                                                                                       |                                                                                  
CAACAGCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGG  >  NZ_CP009273/2336691‑2336860

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: