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

CCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCAT  >  NZ_CP009273/2336699‑2336873
                                                                                      |                                                                                        
ccGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCAccc                                                                                       >  2:228327/1‑90 (MQ=255)
            gTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCCCCCGCCGCAGaggaa                                                                           <  2:158830/90‑1 (MQ=255)
                tGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGAc                                                                       >  1:263841/1‑90 (MQ=255)
                 gcaaaTAACCAGCGTTTTTCCGTCGGAAATATCAAAGGTAACTTCACTTAAGCCGAGATACAAAAAGCCCCCCGCCGCAGAGGAAGGAcc                                                                      <  2:296881/87‑1 (MQ=255)
                 gCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGAcc                                                                      >  2:289244/1‑90 (MQ=255)
                   cAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCCCCCGCCGCAGAGGAAGGACCAt                                                                    <  2:52810/90‑1 (MQ=255)
                    aaaaaCCAGCGTTTTCCCGCCGCCAATAACAAGGGTACCTCCCCTTAAGCCGAAATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCAtc                                                                   <  1:264168/87‑1 (MQ=255)
                    aaTACCCAGCTTTTTTCCGCCGGCAATATAAAAGGTAACTTCACTAAAGCCGAGATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCAtc                                                                   <  2:483810/90‑1 (MQ=255)
                    aaTAACCACCTTTTTCCCGTCGGAAATAACAAAGGTAACTTCCCTAAAGCCGAAAAACAAAAACCCCCCCGCCGCAGAGGAAGGACCAtc                                                                   <  1:80011/90‑1 (MQ=255)
                        cccAGCTTTTTCCCGCCGGAAATATCAAAGGTAACTTCACTTAAGCCGAAATACAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCa                                                               <  2:361548/89‑1 (MQ=255)
                           aGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAg                                                            >  1:347092/1‑90 (MQ=255)
                            gCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGc                                                           >  1:343238/1‑90 (MQ=255)
                                ttttCCGCCGGCAATACAAAAGGTAACTTCACTAAAGCCGAGATACAAAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc                                                       <  2:552727/90‑1 (MQ=255)
                                tttcCCGTCGGCAATATAAAAGGTAACTTCCCTAAAGCCGAGATAAAAAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc                                                       <  1:101513/90‑1 (MQ=255)
                                cttcccGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc                                                       <  2:245899/86‑1 (MQ=255)
                                   cccGTCGCCAATAAAAAAGGTAACTTCACTTAAGCGGAGATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCtgctgc                                                    <  1:289244/89‑1 (MQ=255)
                                                      gacaCTTCATTAAACCCGAGAAACATAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACt                                 <  1:122749/87‑1 (MQ=255)
                                                                          aTACAAAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCt             <  2:359399/90‑1 (MQ=255)
                                                                                     cACCCGCCGCGGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGtcat  >  2:429310/1‑90 (MQ=255)
                                                                                     cACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGtcat  >  2:99625/1‑90 (MQ=255)
                                                                                      |                                                                                        
CCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCAT  >  NZ_CP009273/2336699‑2336873

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: