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

GCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATC  >  NZ_CP009273/2336696‑2336874
                                                                                         |                                                                                         
gCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCa                                                                                           >  1:566015/1‑90 (MQ=255)
   ccGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCAccc                                                                                        >  1:464043/1‑90 (MQ=255)
                  ttGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGaggaagga                                                                         >  2:82500/1‑90 (MQ=255)
                    gCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGAcc                                                                       >  2:213518/1‑90 (MQ=255)
                      cAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTACCTTCATTTAAGCCGAGATACAAAACGCCCCCCGCCGCAGAGGAAGGACCAt                                                                     <  2:34435/90‑1 (MQ=255)
                           aCCAGCTTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCa                                                                <  2:61802/90‑1 (MQ=255)
                              ggctttttcCCGCCGGAAATAAAAAAGGTACCTTCACTAAAGCCGAAATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAg                                                             <  1:370378/86‑1 (MQ=255)
                              aGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAg                                                             >  1:591046/1‑90 (MQ=255)
                                   ttttCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACAAAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc                                                        <  2:101295/90‑1 (MQ=255)
                                   ttttCCGCCGCAAATAACAAAGTTACCTCCCCTAAAGCCGAGATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc                                                        <  2:419997/90‑1 (MQ=255)
                                   tttcCCGTCGGCAATTACAAGGGAACCTTCTCTAAACCCGAGATAAAAAAACCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc                                                        <  2:591046/90‑1 (MQ=255)
                                   tttcCCGCCGGAAATAAAAAAGGTAACTCCACTAAAGCCGAGATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc                                                        <  1:82500/90‑1 (MQ=255)
                                        cGTCGCAAATATCAAAGGAACCCTCGCTTAAGCGGGGTCACTTAAGGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTgcgc                                                   <  1:318330/90‑1 (MQ=255)
                                         gTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTgcgcc                                                  >  2:370378/1‑89 (MQ=255)
                                                    cAAAGGTAACTTCACTTAAGCCGAAATAAAAAACCCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCAt                                       <  2:414216/90‑1 (MQ=255)
                                                     aaaGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGaggaagga                                                                         >  1:460237/1‑55 (MQ=255)
                                                     aaaGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGaggaagga                                                                         <  2:460237/55‑1 (MQ=255)
                                                         gTAACTCCACTAAACCCGAGATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACt                                  <  2:162819/90‑1 (MQ=255)
                                                           aaCTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTaa                                >  1:582822/1‑90 (MQ=255)
                                                           aaCTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTaa                                >  2:458297/1‑90 (MQ=255)
                                                              ttCACTAAAGCCGAAAAAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAacac                             <  1:253066/90‑1 (MQ=255)
                                                                        ccGAAATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTc                   <  1:183578/90‑1 (MQ=255)
                                                                            aaaaCAAAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGc               <  1:503708/87‑1 (MQ=255)
                                                                                         aCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGtcatc  >  2:129438/1‑90 (MQ=255)
                                                                                         |                                                                                         
GCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATC  >  NZ_CP009273/2336696‑2336874

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