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

AGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGCCCCGGTAACCACTGATACCCCGCCTCCAGCGAGG  >  NZ_CP009273/2334161‑2334317
                                                                                 |                                                                           
aGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGGAAAAACATCCCGCTGCACGCCCCGATaaaac                                                                     >  2:387391/1‑88 (MQ=255)
aGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGGAAAAAAATACTGCTGCCCGCCCCGATaaaac                                                                     >  2:477537/1‑88 (MQ=255)
            atatCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCAcctgcgcctgc                                                         <  2:360446/90‑1 (MQ=255)
                  tcaCCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCa                                                   <  1:647875/90‑1 (MQ=255)
                    aCCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAAACATCCTGCTGCACGCCCCGATAAATCACCTGCGCCCGCGGGTCAAt                                                 >  2:199672/1‑90 (MQ=255)
                       cTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAAACATCCTGCTGCACGCCCCGATAAAACACCTGCGCCTGCGGTTCAAacac                                              >  2:647875/1‑90 (MQ=255)
                       cTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAAACATCCTGCTGCACGCCCCGATAAAACACCCGCGCCCGCGGGTCAATcac                                              >  2:346596/1‑90 (MQ=255)
                       cTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAAACATCCCGCTGCCCGCCCCGATAAAAAACCCGCGCCCGCGGGTTAATcac                                              >  2:150896/1‑90 (MQ=255)
                         gCGATTGTGACACGCGCGCACGGTTAGCGGCGGGAAAAAAATCCTCCTGCCCGCCCCGATAAATTACCTGCGCCCTCGGGGCAATcacca                                            >  1:461765/1‑90 (MQ=255)
                          cGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATcaccac                                           <  2:585806/90‑1 (MQ=255)
                                   cacGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGccc                                  <  1:199672/90‑1 (MQ=255)
                                      gcgcgcACGGTTAGCGGCGGTAAAAACATCCTGCTGCACGCCCCGATAAAACACCCGCGCCCGCGGGTCAATCACCACACCAcccccccg                               >  2:334989/1‑88 (MQ=255)
                                                      gcggTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGCCCCGGTAACCACTGATAcccc               <  2:589100/90‑1 (MQ=255)
                                                         ggaaaaACATCCTGCTGCACCCCCCGATAAAACACCCGCGCCCGCGGGTCAAACACCCCACCCCCCCCCCGCAACCACTGATACcccccc            >  2:458265/3‑90 (MQ=255)
                                                         gTAAAAACATCCTTCTGCACGCCCCGATAAAACACCCGCGCCCGCGGGTCAATCACCACACCCCCCCCCCGGAAACACCGATACcccccc            >  2:322235/1‑90 (MQ=255)
                                                           aaaaTCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGCCCCGGTAACCACTGATACCCCGCCTc          <  1:458265/90‑1 (MQ=255)
                                                                 atcCTGCTGCACGCCCTGATAAATCCCCTGCGCCTGCGGTTCAATCACCACACCCCGCCCCGGTAACCACTGATACCCCGCCTCCAGCGa    <  2:461765/90‑1 (MQ=255)
                                                                   ccTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGCCCCGGTAACCACTGATACCCCGCCTCCAGCGAgg  >  2:255901/1‑90 (MQ=255)
                                                                                 |                                                                           
AGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGCCCCGGTAACCACTGATACCCCGCCTCCAGCGAGG  >  NZ_CP009273/2334161‑2334317

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: