Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,679,3990GT55.6% 10.7 / 24.7 18N1045K (AAC→AAAbcsCcellulose synthase complex outer membrane protein BcsC
Reads supporting (aligned to +/- strand):  ref base G (5/3);  new base T (0/10);  total (5/13)
Fisher's exact test for biased strand distribution p-value = 6.54e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-03
Rejected as consensus: Frequency below/above cutoff threshold.

GATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGC  >  NZ_CP009273/3679320‑3679487
                                                                               |                                                                                        
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)
                                                                      ccACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCgag          >  1:80678/1‑90 (MQ=255)
                                                                              aGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGc  <  2:80678/90‑1 (MQ=255)
                                                                               |                                                                                        
GATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGC  >  NZ_CP009273/3679320‑3679487

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: