Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,598,828 A→G 9.5% intergenic (‑67/+64) dctB ← / ← bglX C4‑dicarboxylate transport sensor protein/beta‑D‑glucoside glucohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,598,8280AG 9.5% 109.7 / 2.4 42intergenic (‑67/+64)dctB/bglXC4‑dicarboxylate transport sensor protein/beta‑D‑glucoside glucohydrolase
Reads supporting (aligned to +/- strand):  ref base A (15/23);  new base G (2/2);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCCGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGC  >  NC_002947/1598682‑1598969
                                                                                                                                                  |                                                                                                                                              
cAGCCGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAAACCGCGCCTACAGGTGct                                                                                                                                              >  8:152980/1‑149 (MQ=255)
      gcgcTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑Tgcag                                                                                                                                       <  4:272264/149‑1 (MQ=255)
        gcTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGAGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑Tgcagca                                                                                                                                     >  1:170737/1‑149 (MQ=255)
            aGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGGGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCt                                                                                                                                  <  2:170737/148‑1 (MQ=255)
             gTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTg                                                                                                                                 >  2:59597/1‑148 (MQ=255)
                   tACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAAACCGCGCCTACAGGTGCTG‑CAGAGTACCTGTAGgcgc                                                                                                                          <  7:180364/149‑1 (MQ=37)
                   tACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAAACCGCGCCTACAGGTGCTG‑CAGAGTACCTGTAGgcgc                                                                                                                          <  7:152980/149‑1 (MQ=37)
                     cccAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACCGGTACTC‑TGCAGCACCTGTAGGcgcg                                                                                                                         <  1:59597/148‑1 (MQ=255)
                     cccAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGcgcg                                                                                                                         <  2:93855/148‑1 (MQ=255)
                      ccAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGTTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGt                                                                                                                       <  1:74885/149‑1 (MQ=255)
                       cAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGtt                                                                                                                      >  6:191239/1‑149 (MQ=255)
                         gCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGtttt                                                                                                                    <  1:343845/149‑1 (MQ=255)
                           aCGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTAc                                                                                                                  <  5:191239/149‑1 (MQ=255)
                             gATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCAACACGTGCTC‑TGCCGCACCGGTAGGCGCGGTTTTAcgc                                                                                                                >  5:48584/1‑147 (MQ=255)
                                    agGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGaaaa                                                                                                         <  5:189674/149‑1 (MQ=255)
                                       gTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGc                                                                                                      <  2:107310/149‑1 (MQ=255)
                                           tCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCgg                                                                                                   >  8:223789/1‑148 (MQ=255)
                                           tCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCgg                                                                                                   <  1:134729/148‑1 (MQ=255)
                                             aGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCTGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAAc                                                                                                <  8:239293/149‑1 (MQ=255)
                                                    gACGCAGGGAGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggcggc                                                                                         <  2:341059/149‑1 (MQ=255)
                                                            gcgcGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCa                                                                                 <  1:2895/149‑1 (MQ=255)
                                                             cgcgAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTTCCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCa                                                                                 <  1:102104/149‑1 (MQ=255)
                                                                gAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCg                                                                                                    <  2:31899/126‑1 (MQ=255)
                                                                gAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCg                                                                                                    >  1:31899/1‑126 (MQ=255)
                                                                gAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTc                                                                              >  6:205985/1‑148 (MQ=255)
                                                                  aaGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCa                                                                             >  5:209368/1‑147 (MQ=255)
                                                                   aGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAg                                                                          >  7:42623/1‑149 (MQ=255)
                                                                    ggAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGACACAg                                                                          >  4:251927/1‑148 (MQ=255)
                                                                     gAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGc                                                                         >  1:330455/1‑148 (MQ=255)
                                                                       aaaGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑Tg                                                                                                                                          <  2:132517/81‑1 (MQ=255)
                                                                       aaaGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑Tg                                                                                                                                          >  1:132517/1‑81 (MQ=255)
                                                                       aaaGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGaaa                                                                                                          >  3:251257/1‑113 (MQ=255)
                                                                       aaaGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGaaa                                                                                                          <  4:251257/113‑1 (MQ=255)
                                                                              aTCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCa                                                                                 <  4:172758/131‑1 (MQ=255)
                                                                                        gATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTctgctgctg                                                     >  3:349665/1‑149 (MQ=255)
                                                                                                  tggGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTgg                                           <  2:50416/149‑1 (MQ=255)
                                                                                                  tggGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTgg                                           <  8:42623/149‑1 (MQ=255)
                                                                                                                            gggTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGaa                 <  7:223789/149‑1 (MQ=255)
                                                                                                                                   cccGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCt          >  1:205322/1‑149 (MQ=255)
                                                                                                                                         ccTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCg    <  2:205322/149‑1 (MQ=255)
                                                                                                                                          cTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCgg   >  2:81821/1‑149 (MQ=255)
                                                                                                                                             cAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGc  <  7:37009/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
CAGCCGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTC‑TGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGC  >  NC_002947/1598682‑1598969

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑