Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,334,2240TA52.9% 1.9 / 19.2 17D1142V (GAT→GTT) yfaLAIDA‑I family autotransporter adhesin YfaL/EhaC
Reads supporting (aligned to +/- strand):  ref base T (1/7);  new base A (9/0);  total (10/7)
Fisher's exact test for biased strand distribution p-value = 4.11e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGTACGCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGCCCCGGTAACC  >  NZ_CP009273/2334136‑2334293
                                                                                        |                                                                     
ggTACGCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAatc                                                                      >  1:181850/1‑90 (MQ=255)
   aCGCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAgtcatc                                                                   >  1:487389/1‑90 (MQ=255)
   aCGCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAaacatc                                                                   >  1:3053/1‑90 (MQ=255)
      ccATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCctg                                                                <  2:181850/90‑1 (MQ=255)
      ccATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGGAAAAACATcccc                                                                >  1:515189/1‑88 (MQ=255)
        aTTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGGAAAAACATccctct                                                              >  2:99906/1‑86 (MQ=255)
                tgtAAACTCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGGAAAAACATCCTGCTGCACcccc                                                      >  2:643617/1‑90 (MQ=255)
                    aaCCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGAt                                                  <  1:710137/90‑1 (MQ=255)
                         aGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAAAAATCCTGCTGCCCCCCCCGATaaaac                                             >  2:492347/1‑88 (MQ=255)
                                        tcatcaCCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAAACATCCCGCTGCACCCCCCGATAAAACACCCGCGCCTGCggg                              >  2:113399/1‑89 (MQ=255)
                                             aCCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAAACATCCCTCTGCACCCCCTGGTAAAACACCTGCGCCTGCGGGTGAAt                         >  1:343992/1‑90 (MQ=255)
                                                     aTTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATcaccacac                 <  2:343992/90‑1 (MQ=255)
                                                        gtgACACGCGCGCACGGTTAGCGGCGGTAAAAACATCCTGCTGCACGCCCCGATAAAACACCCGCGCCCGCGGGTCAATCaccacaccac              >  2:83950/1‑90 (MQ=255)
                                                            cacGCGCGCACGGTTAGCGACGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGccc          <  2:142563/90‑1 (MQ=255)
                                                               gcgcgcACGGTTAGCGGCGGTAAAAACATCCTTCTGCACGCCCCGATAAAACACCCCCGCCGGCGGGTCAATCACCCCACcccccccccg       >  2:487714/1‑88 (MQ=255)
                                                               gcgcgcACGGTTAGCGGCGGTAAAAACATCCTGCTGCACCCCCCCATAAAACACCTGCGCCCGCGGGTCAATCACCACACCAcccccccg       >  2:489975/1‑88 (MQ=255)
                                                                    cACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGCCCCGGTAAcc  <  1:83950/90‑1 (MQ=255)
                                                                    cACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGCCCCGGTAAcc  <  2:56777/90‑1 (MQ=255)
                                                                                cggTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGtt                             >  1:603272/1‑51 (MQ=255)
                                                                                cggTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGtt                             <  2:603272/51‑1 (MQ=255)
                                                                                        |                                                                     
GGTACGCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGCCCCGGTAACC  >  NZ_CP009273/2334136‑2334293

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