Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,225,9540AG62.5% ‑0.1 / 15.4 16S41G (AGT→GGT) yjbEexopolysaccharide production protein YjbE
Reads supporting (aligned to +/- strand):  ref base A (0/6);  new base G (10/0);  total (10/6)
Fisher's exact test for biased strand distribution p-value = 1.25e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATG  >  NZ_CP009273/4225876‑4226029
                                                                              |                                                                           
gcggcgACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGGGGTTCTCcgcg                                                                  >  2:436370/1‑90 (MQ=255)
    cgACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCa                                                              <  1:565861/90‑1 (MQ=255)
           tGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAg                                                       <  1:486410/90‑1 (MQ=255)
                gggCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGGGTTTCTCCGCGGCCCGCGGCCGCACCg                                                  >  2:340499/1‑90 (MQ=255)
                gggCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGGGGTTCCCCCCGACCAGCGGCAGCAccc                                                  >  2:389788/1‑89 (MQ=255)
                      cACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCt                                            <  1:92588/90‑1 (MQ=255)
                               aGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGGGTTCCTCCGCGACCAGCGGCAGCACCGTAAGCTCGGCggggg                                   >  2:484621/1‑90 (MQ=255)
                                ggTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGGGTTTCTCCCCGGCCCGCGTCAGCACCGGAAACTCCGCgggggg                                  >  2:25638/1‑90 (MQ=255)
                                ggTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGGGTTCCTCCCCGACCCGCGGCAGCACCGTAAGCTCGGCgggggg                                  >  2:577806/1‑90 (MQ=255)
                                  tGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGGGGTCCTCCCCGGCCCGCGGCCGCCCCCCAAGCCCCGCgggggggg                                >  1:592596/1‑88 (MQ=255)
                                  tGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGGGGTCCTCCCCCACCCGCGGCAGCACCCTAAGGTCGGCgggggggg                                >  1:437704/1‑88 (MQ=255)
                                   gggCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGGGTTCCTCCGCGCCCAGCGGCAACACCGTAAGCTCGGCggggggggg                               >  2:143463/1‑87 (MQ=255)
                                   gggCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGGGGTCCTCCCCGACCCGCGTCAGCACCGTAAACTCGGCggggggggg                               >  2:293580/1‑87 (MQ=255)
                                                     ggCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCgg             <  1:500636/90‑1 (MQ=255)
                                                                cGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATg  <  1:484621/90‑1 (MQ=255)
                                                                cGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATg  <  1:487543/90‑1 (MQ=255)
                                                                              |                                                                           
GCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATG  >  NZ_CP009273/4225876‑4226029

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