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

GCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTG  >  NZ_CP009273/2336696‑2336869
                                                                                          |                                                                                    
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)
                                                                                     aaGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTg  >  2:110314/1‑90 (MQ=255)
                                                                                          aCCCGCCGCAGAGGAAGGACCATCTCCATAGCCtgctgc                                                <  2:169794/39‑1 (MQ=255)
                                                                                          aCCCGCCGCAGAGGAAGGACCATCTCCATAGCCtgctgc                                                >  1:169794/1‑39 (MQ=255)
                                                                                          aCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCAt                                  <  2:148795/53‑1 (MQ=255)
                                                                                          aCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCAt                                  >  1:148795/1‑53 (MQ=255)
                                                                                          |                                                                                    
GCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTG  >  NZ_CP009273/2336696‑2336869

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: