Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,334,2420TC62.5% ‑0.6 / 11.6 16Q1136R (CAG→CGG) yfaLAIDA‑I family autotransporter adhesin YfaL/EhaC
Reads supporting (aligned to +/- strand):  ref base T (1/5);  new base C (10/0);  total (11/5)
Fisher's exact test for biased strand distribution p-value = 1.37e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.47e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGCCCCGGTAACCACTGATACCCCGCCTCCAGCGAGGCGATAATCCCCGAC  >  NZ_CP009273/2334160‑2334331
                                                                                  |                                                                                         
cAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAAt                                                                                    <  1:580528/90‑1 (MQ=255)
 aGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAAAAATCCTGCTGCACGCCCCGATaaaac                                                                                   >  1:601858/1‑88 (MQ=255)
           gAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGGAAAAACATCCCGCTGCACCCCCCGATAAAACAcctgcgcct                                                                         >  2:274146/1‑90 (MQ=255)
                   tcaCCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCCGCTGCCCCCCCCGGTAAAACACCCCCGCCCGCggggca                                                                 >  2:288700/1‑86 (MQ=255)
                        cTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGGAAAAACATCCCCCCGCACCCCCCGATAAAACACCCGCGCCTCGGGTTCAATaac                                                            >  2:558008/1‑90 (MQ=255)
                         tGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAAACATCCTGCTGCACCCCCCGATAAAACACCCGCGCCTGCGGGTCAATcacc                                                           >  1:41048/1‑90 (MQ=255)
                               tgtgACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCCCATAAATAACCTGCCCCCGCGGTTCAATcaccacacca                                                     >  2:762863/1‑90 (MQ=255)
                                gtgACACGCGCGCACGGTTAGCGGCGGGAAAAACATCCTCCTGCACCCCCCGATAAAACACCCGCGCCTGCGGGTCAATCaccacaccac                                                    >  2:59768/1‑90 (MQ=255)
                                gtgACACGCGCGCACGGTTAGCGGCGGGAAAAACATCCCGCTGCGCGCCCCGATAAAAAACCCCCGCCCGCGGGTCAATCaccacaccac                                                    >  2:59754/1‑90 (MQ=255)
                                    cacGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGccc                                                <  2:109901/90‑1 (MQ=255)
                                          cgcACGGTTAGCGGCGGTAAAAACATCCTGCTGCACGCCCCGGTAAAACACCCGCGCCCGCGGGTCAATCACCACACCACCCCCCCGTaa                                          >  2:763818/1‑90 (MQ=255)
                                          cgcACGGTTAGCGGCGGGAAAAACATCCCTCGGCACGCCCTTCTAAAAAACCCCCGCCCGCGGGTCCATCACCACACCCCCCCCCGGCaa                                          >  2:761182/1‑90 (MQ=255)
                                            cACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGCCCCGGTAAcc                                        <  1:403158/90‑1 (MQ=255)
                                                        cggTAAAATCATCCTGCTGCCCGCCCCGGTAAATCACCCGCGCCCGCGGGTCAAACACCACACCCCCCCCCGGTAACCCCTAGTACCCCg                            >  2:162319/1‑90 (MQ=255)
                                                                  atcCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGCCCCGGTAACCACTGATACCCCGCCTCCAGCGa                  <  1:7721/90‑1 (MQ=255)
                                                                           cACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGCCCCGGTAACCACTGATACCCCGCCTCCAGCGAGGCGATAAt         <  1:233872/90‑1 (MQ=255)
                                                                                  tGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGCCCCGGTAACCACTGATACCCCGCCTCCAGCGAGGCGATAATCCCcgac  >  2:343514/1‑90 (MQ=255)
                                                                                  tGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGCCCCGGTAACCACTGATACCCCGCCTCCAGCGAGGCGATAATCCCcgac  >  2:69030/1‑90 (MQ=255)
                                                                                  |                                                                                         
CAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGCCCCGGTAACCACTGATACCCCGCCTCCAGCGAGGCGATAATCCCCGAC  >  NZ_CP009273/2334160‑2334331

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