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

GCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATC  >  NZ_CP009273/2336696‑2336874
                                                                                         |                                                                                         
gCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCa                                                                                           >  1:159957/1‑90 (MQ=255)
gCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCa                                                                                           >  2:218/1‑90 (MQ=255)
               gTATTGCAAAAACCCAGCGTTTTTCCGTCGGAAATATAAAAGGTAACTTCCCTTAAGCCGAGATAAATAAAGCCCCCCGCCGCAGaggaa                                                                            <  1:256290/90‑1 (MQ=255)
                  ttcccAACACCCAGCTTTTTCCCGGCGGAAATTACAAAGGAAACTTCACTTAAGCCGAGATAAAAAACCCCCCCCGCCGCAGaggaagga                                                                         <  1:728974/87‑1 (MQ=255)
                    gCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGAcc                                                                       >  2:265366/1‑90 (MQ=255)
                     ccAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCa                                                                      <  2:541476/90‑1 (MQ=255)
                      cAATAACCAGCTTTTTTCCGTCGGAAATATAAAAGGTAACTTCACTTAAGCCGAGATAAAAAAAGCCCCCCGCCGCAGAGGAAGGACCAt                                                                     <  1:577211/90‑1 (MQ=255)
                      cAATAACCAGCGTTTTTCCGTCGGAAATATAAAAGGTAACTTCCCTTAAGCCGAGATACAAAAAGCCCCCCGCCGCAGAGGAAGGACCAt                                                                     <  2:335853/90‑1 (MQ=255)
                           aCCAGCTTTTTCCCGCCGGCAATATAAAAGGTAACTTCCCTTAAGCCGAAATAAAAAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCa                                                                <  1:507748/90‑1 (MQ=255)
                              aGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAg                                                             >  2:912697/1‑90 (MQ=255)
                              aGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAg                                                             >  2:675275/1‑90 (MQ=255)
                                      cccgccgGCAATAAAAAAGGTAATTTCACTAACGCCGAAATACATAACGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCtgctgc                                                     <  2:750748/89‑1 (MQ=255)
                                                   taaaaGGTACCTTCACTTAAGCCGAAATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCa                                        <  1:74763/88‑1 (MQ=255)
                                                      aaGGTACCTTCACTTAAGCCGAGATACATAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCg                                     <  2:235075/90‑1 (MQ=255)
                                                       aGGTACCTCCACTAAAGCCGAGATACAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGa                                    <  2:268754/90‑1 (MQ=255)
                                                         gaaaCTCCACTAACGCCGAGATACAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACt                                  <  1:565007/88‑1 (MQ=255)
                                                           accTTCACTTAAGCCGAGATACAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTaa                                <  2:211405/88‑1 (MQ=255)
                                                           aaCTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTaa                                >  1:1133147/1‑90 (MQ=255)
                                                                   ttAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGccct                        >  1:268754/1‑88 (MQ=255)
                                                                                       ccACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTg       >  1:691192/1‑87 (MQ=255)
                                                                                         aCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc                                                        >  2:1121117/1‑36 (MQ=255)
                                                                                         aCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc                                                        <  1:1121117/36‑1 (MQ=255)
                                                                                         aCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGtcatc  >  1:1442427/1‑90 (MQ=255)
                                                                                         aCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGtcatc  >  1:1329730/1‑90 (MQ=255)
                                                                                         |                                                                                         
GCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATC  >  NZ_CP009273/2336696‑2336874

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