Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,336,7780TA63.6% ‑0.4 / 22.1 22M291L (ATG→TTG) yfaLAIDA‑I family autotransporter adhesin YfaL/EhaC
Reads supporting (aligned to +/- strand):  ref base T (8/0);  new base A (0/14);  total (8/14)
Fisher's exact test for biased strand distribution p-value = 3.13e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-05
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTCAACAGCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTC  >  NZ_CP009273/2336689‑2336852
                                                                                          |                                                                          
gTCAACAGCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACAt                                                                            >  2:419929/1‑90 (MQ=255)
  cAACAGCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATaa                                                                          >  2:456318/1‑90 (MQ=255)
  cAACAGCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATaa                                                                          >  1:574489/1‑90 (MQ=255)
       gCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCa                                                                     >  1:230322/1‑90 (MQ=255)
          ccGTCTTTCTCTGTATTCCCAAAAACCGGTTTTTTCCCGCCGGAAATATCAAAGGTACCTCCATTAAAGCCGAGATAAAAAAACcccccc                                                                  <  2:460842/90‑1 (MQ=255)
          ccGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCAccc                                                                  >  1:524801/1‑90 (MQ=255)
            gTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACccgc                                                                >  1:472378/1‑90 (MQ=255)
                      gTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGaggaa                                                      >  2:205571/1‑90 (MQ=255)
                      gTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGaggaa                                                      >  1:194714/1‑90 (MQ=255)
                         ttGCAAATACCCAGCTTTTTCCCGCCGGAAATAAAAAAGGTAACTTCACTAAAGCCGAAATAAAAAAACCCCCCCGCCGCAGaggaagga                                                   <  2:194714/90‑1 (MQ=255)
                            caaaTACCCAGCTTTTTCCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACAAAAAGCCCCCCGCCGCAGAGGAAGGACCa                                                <  1:451434/88‑1 (MQ=255)
                            caaaTACCCAGCTTTTTCCCGCCGGAAATACAAAAGGAACCTTCCCTTAAGCCGAAATACAAAAAGCCCCCCGCCGCAGAGGAAGGACCa                                                <  2:11133/88‑1 (MQ=255)
                             cAATAACCGGCTTTTTCCCGGCGGAAATATAAAAGGAAACTTCCCTAAAGCCGAAATACAAAAACCCCCCCGCCGCAGAGGAAGGACCAt                                               <  1:303928/90‑1 (MQ=255)
                             aaaTAACCAGCTTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACAAAAAGCCCCCCGCCGCAGAGGAAGGACCAt                                               <  1:511090/89‑1 (MQ=255)
                             aaaTAACCAGCGTTTTCCCGTCGGAAAAATAAAAGGTAACTCCACTAAAGCCGAGAAACAAAAAGCCACCCGCCGCAGAGGAAGGACCAt                                               <  2:602635/89‑1 (MQ=255)
                              aaaaaCCCGGGTTTTTTCGTTGGACAATAAAAAAGTAAACTTCATTAAAGCGGAAATACAAAAACCCCCCCGCCGCAGAGGAAGGACCAt                                               <  2:414020/87‑1 (MQ=255)
                              aaTAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCAtc                                              >  2:556299/1‑90 (MQ=255)
                                          ttttCCGGGGGCAATAAAAAAGGTACTTTCACTAAAGCAGAAAAAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc                                  <  2:591635/90‑1 (MQ=255)
                                             tCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAAATACAAAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCtgctgc                               <  1:205571/90‑1 (MQ=255)
                                                            aaaaGGTAACTTCCCTAAAGCCGAGATACAAAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCAt                 <  1:456318/89‑1 (MQ=255)
                                                                 gTAACTTCACTAAAGCCGAAATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACt            <  1:603867/90‑1 (MQ=255)
                                                                           ttAAGCCGAAATAAAAAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTc  <  1:419929/90‑1 (MQ=255)
                                                                           taaaGCCGAAATACAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTc  <  2:230322/88‑1 (MQ=255)
                                                                                          |                                                                          
GTCAACAGCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTC  >  NZ_CP009273/2336689‑2336852

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: