Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,679,3880TC56.2% 6.4 / 14.4 16E1049G (GAG→GGG) bcsCcellulose synthase complex outer membrane protein BcsC
Reads supporting (aligned to +/- strand):  ref base T (4/3);  new base C (0/9);  total (4/12)
Fisher's exact test for biased strand distribution p-value = 1.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.80e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCC  >  NZ_CP009273/3679320‑3679471
                                                                    |                                                                                   
gATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTgc                                                                >  2:74232/1‑90 (MQ=255)
  ttCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTACCACGACCAGTTTTCCGTgcgc                                                              <  1:171298/90‑1 (MQ=255)
               aaCGCGGCATGTTTTTGGTGCGTGAAGGCGCCCACGAGCCAGACGCCCCCAGCCCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCa                                                 <  1:102467/90‑1 (MQ=255)
                   gggCATGGTTTTGGTGCGTGAATGCGCCCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCCCATCACCgg                                             <  2:356535/89‑1 (MQ=255)
                      cATGTTTTTGGTGCGTGAATGCGCCCCCGAGCCAGACGCCCCCAGCCCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTAt                                          <  1:1986/90‑1 (MQ=255)
                        tGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATgg                                        >  1:87420/1‑90 (MQ=255)
                                    gTGAGTGCCCCCCCGACCCAGCCCCCCCCAGCCCCCCCGCCCTTTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGt                            <  2:87420/90‑1 (MQ=255)
                                    gTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGt                            >  1:192624/1‑90 (MQ=255)
                                               cccGAGCCAGACGCCCCCCGCCCCCACACCTTTTTTTCCGGGCGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGa                 <  2:375481/88‑1 (MQ=255)
                                                ccgcgCCAGCCCCCCCCACCCCCCCCGCCCTTTTTTCCGGGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGAc                <  2:66942/86‑1 (MQ=255)
                                                ccGACCCAGACCCCCCCCGCCCCCACGACCGGTTTTCCGGGCCCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGAc                <  2:358820/89‑1 (MQ=255)
                                                aCGAGCCAGACGCCCCCAGCCCCCACGACCATTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGAc                <  1:70626/90‑1 (MQ=255)
                                                          cccccccACCCCCCACGACCGGTTTCCGGTGCTCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAAcc      <  1:350962/86‑1 (MQ=255)
                                                            cACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGa                                 <  1:64632/61‑1 (MQ=255)
                                                            cACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGa                                 >  2:64632/1‑61 (MQ=255)
                                                              cccACCCCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGccc  <  1:157787/90‑1 (MQ=255)
                                                                    |                                                                                   
GATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCC  >  NZ_CP009273/3679320‑3679471

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: