Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,679,3880TC64.7% ‑5.1 / 11.9 17E1049G (GAG→GGG) bcsCcellulose synthase complex outer membrane protein BcsC
Reads supporting (aligned to +/- strand):  ref base T (4/2);  new base C (0/11);  total (4/13)
Fisher's exact test for biased strand distribution p-value = 6.30e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.55e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCC  >  NZ_CP009273/3679305‑3679470
                                                                                   |                                                                                  
aGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACg                                                                              <  2:355819/90‑1 (MQ=255)
               gATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTgc                                                               >  1:40403/1‑90 (MQ=255)
               gATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTgc                                                               >  2:271714/1‑90 (MQ=255)
                          ggATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCcac                                                    >  2:21629/1‑90 (MQ=255)
                                cgcgGTTTGTTTTGGGCGCGTGAAGGCCCCCCCGGCCCAGCCCCCCCCGGCCCCCCCGACCATTTTCCCGTGCGCCCCCGCCACATCAcc                                              <  1:233999/90‑1 (MQ=255)
                                cgcgGCTTGTTTTGGGTGGGTGAATGCCCCCCCGAGCCAGCCCCCCCCAGCCCCCCCGACCATTTTTCCGTGCCCCCCCGCCACATCAcc                                              <  1:138849/90‑1 (MQ=255)
                                cgcgGCTTGGTTTGGGTGGGTGAAGGCGCCCACGACCCAGCCCCCCCCAGCCCCCCCGCCCATTTTTCCGTGCGCTCCCGCCACATCAcc                                              <  1:21629/90‑1 (MQ=255)
                                          ttttGGGGCGTGAAGGCGCCCGCGACCCAGACGCCCCCGGCCCCCACGCCAAGTTTTCCGTGCGCCCCCGCCACATCACCGGTATGGCaa                                    <  1:184007/90‑1 (MQ=255)
                                              ggTGCGTGAATGCGACCACGAGCCAGCCGCCCCCAGCTCCCACGACCATTTTCCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGa                                <  2:73144/90‑1 (MQ=255)
                                                   gtaaaTGCCCCCACGAGCCAGACCCCCCCCGCCCCCCCGCCCATTTTTCCGTGCCCCCCCGCCACATCACCGGTATGGCAAACGACAGGt                           <  2:385317/87‑1 (MQ=255)
                                                   gTGAATGCGCCCACGAGCCAGCCGCCCCCCCCCCCCACGACCATTTTTCCGTGCCCTCCCGCCACATCACCGGTATGGCAAACGACAGGt                           <  2:313339/90‑1 (MQ=255)
                                                          cgcccACGCCCCAGCCCCCCCCCGCCCCCACGCCCATTTTCCCGTGCGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTg                    <  1:126324/87‑1 (MQ=255)
                                                          cgcccACGAGCCAGACGCCCCCAGCCCCCACGACCATTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTg                    <  2:40403/87‑1 (MQ=255)
                                                               ccgcgCCAGACGCACCCAGCCCCCACGCCCAGTTTTCCGGGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGAc               <  2:195337/86‑1 (MQ=255)
                                                               ccGACCCAGCCCCCCCCAGCCCCCACGACCATTTTTCCGTGCGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGAc               <  2:372184/89‑1 (MQ=255)
                                                               aCGACCCAGCCCCCCCCAGCCCCCCCGCCCGTTTTCCCGGCCCCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGAc               <  1:483367/90‑1 (MQ=255)
                                                                            aCCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGcc  >  1:159322/1‑90 (MQ=255)
                                                                                   |                                                                                  
AGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCC  >  NZ_CP009273/3679305‑3679470

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