Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,563,3980GC34.6% 30.1 / 15.5 27S779C (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 (11/6);  new base C (0/9);  total (11/16)
Fisher's exact test for biased strand distribution p-value = 2.41e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTG  >  NC_000913/1563331‑1563460
                                                                   |                                                              
ttGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCggg                                                                >  1:25473/1‑68 (MQ=255)
ttGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGCAGGGGGCggg                                                                <  2:897609/68‑1 (MQ=255)
   tCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGaaa                                                             >  2:773899/1‑68 (MQ=255)
       aTTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGCAAAAGaa                                                         <  2:1460385/68‑1 (MQ=255)
            cAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGCAAAACAAATATc                                                    <  2:835863/68‑1 (MQ=255)
                 gTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGa                                               >  1:1443538/1‑68 (MQ=255)
                     cacGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATaa                                           >  1:1144171/1‑68 (MQ=255)
                       cGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAAcg                                         >  2:465494/1‑68 (MQ=255)
                        gctgctCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAAcgc                                        >  1:74611/1‑68 (MQ=255)
                        gctgctCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGAGGGAAAAGAAATATCCCTGAATAAcga                                        >  1:1137996/1‑67 (MQ=255)
                           gctCATCCATCCTGGAATTTCTTCGGCGGGTAGGGGGCGGTAAAATAAATATCCCTGAATAACGCGAc                                     <  1:576891/68‑1 (MQ=255)
                           gctCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGAc                                     <  1:576893/68‑1 (MQ=255)
                           gctCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGCAAAAGAAATATCCCTGAATAACGCGAc                                     <  1:1425781/68‑1 (MQ=255)
                           gctCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGCAAAACAAATATCCCTGAATAACGCGAc                                     <  1:885140/68‑1 (MQ=255)
                             tCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGCAAAAGAAATATCCCTGAATAACGCGACAg                                   <  2:1436774/68‑1 (MQ=255)
                             tCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGCAAAACAAATATCCCTGAATAACGCGACAg                                   <  2:918593/68‑1 (MQ=255)
                                cccAGCCTGGAATTTCTTCGGCGGTTAGGGGGCGGGAAATGAAATATCCCTGAATAACGCGACAGTgg                                <  2:1311050/67‑1 (MQ=255)
                                 ccAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGa                               <  1:1005851/68‑1 (MQ=255)
                                  cAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACTGTGGAt                              >  1:462620/1‑68 (MQ=255)
                                    gCCTGGAATTTCTTCGGCGGGTAGGGGGCGGCAAAAGAAATATCCCTGAATAACGCGACAGCGGATCt                            <  2:683039/68‑1 (MQ=255)
                                       tgtaaTTCCTTCGGCGGGTAGGGGGCGTCAAAACAAATATCCCTGAATAACGCGACAGTGGATCTTGc                         <  1:1424809/65‑1 (MQ=255)
                                           aTTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAg                     >  1:111810/1‑68 (MQ=255)
                                            cttcttcGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGc                    <  2:1443533/67‑1 (MQ=255)
                                              tcttcGGCGGGTAGGGGGCGGCAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCAt                  <  1:543018/68‑1 (MQ=255)
                                                        gTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGc        <  2:25473/68‑1 (MQ=255)
                                                             gggCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCttt   <  1:36036/68‑1 (MQ=255)
                                                              ggCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTg  >  2:795473/1‑68 (MQ=255)
                                                              ggCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTg  >  2:1267487/1‑68 (MQ=255)
                                                              ggCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTg  >  2:944254/1‑68 (MQ=255)
                                                                   |                                                              
TTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTG  >  NC_000913/1563331‑1563460

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: