Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,563,3980GC38.5% 30.1 / 11.5 28S779C (TCC→TGC) dosPoxygen sensor, c‑di‑GMP phosphodiesterase, heme‑regulated; cold‑ and stationary phase‑induced bioflim regulator
Reads supporting (aligned to +/- strand):  ref base G (10/6);  new base C (0/10);  total (10/18)
Fisher's exact test for biased strand distribution p-value = 2.55e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGT  >  NC_000913/1563334‑1563462
                                                                |                                                                
tCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGaaa                                                               >  2:372900/1‑68 (MQ=255)
         ctgCGGTTACACGCTGCTCATCCAGCCCGGAATTTCTTCGGCGGGTTGGTGGCGGTCAAAGTATTATc                                                      <  2:377875/66‑1 (MQ=255)
         cAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATc                                                      <  2:543269/68‑1 (MQ=255)
         cAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGCAAAAGAAATATc                                                      <  2:668277/68‑1 (MQ=255)
             ggcaaCACGCTGCTCATCCAGCCTGGAACTTCTTCGGCGGGTAGGGGGCGGCAAAAGAAATATCCCTg                                                  <  2:1177834/65‑1 (MQ=255)
                   acGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAAc                                            >  1:581689/1‑68 (MQ=255)
                   acGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGCAAAAGAAATATCCCTGAATAAc                                            <  1:216350/68‑1 (MQ=255)
                     gctgctCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAAcgc                                          >  2:937518/1‑68 (MQ=255)
                     gctgctCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAAcgc                                          >  2:879372/1‑68 (MQ=255)
                        gctCATCCTGCCCGGAATTTCTTCGGCGGGTAGGGGGCGGCAAAAGAAATATCCCTGAATAACGCGAc                                       <  1:1104175/68‑1 (MQ=255)
                        gctCATCCAGCCTGGAATTTCTTCGGCGTGTAGGGGGCGGCAAAAGAAATATCCCTGAATAACGCGAc                                       <  2:703205/68‑1 (MQ=255)
                        gctCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGAc                                       <  2:932511/68‑1 (MQ=255)
                        gctCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGAc                                       <  2:173441/68‑1 (MQ=255)
                        gctCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAACAAATATCCCTGAATAACGCGAc                                       <  1:1032728/68‑1 (MQ=255)
                        gctCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGCAAAAGAAATATCCCTGAATAACGCGAc                                       <  2:818084/68‑1 (MQ=255)
                        gctCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGCAAAAGAAATATCCCTGAATAACGCGAc                                       <  2:670657/68‑1 (MQ=255)
                        gctCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGCAAAAGAAATATCCCTGAATAACGCGAc                                       <  2:546587/68‑1 (MQ=255)
                          tCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAg                                     >  2:1073721/1‑68 (MQ=255)
                           cATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGCAAAAGAAATATCCCTGAATAACGCGACAGt                                    <  2:570957/68‑1 (MQ=255)
                            aTCCTGCCTGGAATTTCTTCTGCGTGTAGGGGTCGGGAAAAGAAATATCCCTGAATAACGCGACAGTg                                   <  2:581689/68‑1 (MQ=255)
                               cAGCCTGGAATTTCTTCGGCGGGTAGGGTGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGAt                                <  1:1105219/68‑1 (MQ=255)
                               cAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGAt                                >  1:98528/1‑68 (MQ=255)
                                  ccTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCtt                             >  1:16152/1‑68 (MQ=255)
                                   cTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATACCTGAATAACGCGACAGTGGATCTTg                            >  2:525738/1‑68 (MQ=255)
                                             ttcGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCAtct                  >  2:913715/1‑68 (MQ=255)
                                                        gggggCGGCAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGctct       <  1:1038456/68‑1 (MQ=255)
                                                          gggCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCttt     >  1:7592/1‑68 (MQ=255)
                                                             cGGTAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGt  <  1:340668/68‑1 (MQ=255)
                                                                |                                                                
TCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGT  >  NC_000913/1563334‑1563462

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: