Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,089,943 T→C 7.7% intergenic (+98/‑7) kdgD → / → gudP 5‑dehydro‑4‑deoxyglucarate dehydratase/glucarate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,089,9430TC 7.7% 143.9 / 2.7 52intergenic (+98/‑7)kdgD/gudP5‑dehydro‑4‑deoxyglucarate dehydratase/glucarate transporter
Reads supporting (aligned to +/- strand):  ref base T (24/24);  new base C (2/2);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

AATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGC‑GCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGT  >  NC_002947/4089806‑4090089
                                                                                                                                         |                                                                                                                                                   
aaTGCTCGCAGCGCTTATGGACAAGATGGGCCCGCAGTAAACTCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTTGGCCTAACCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTTGGAGTATTGCCAGCATGc                                                                                                                                           <  3:195676/148‑1 (MQ=255)
  tGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAg                                                                                                                                        <  3:146979/149‑1 (MQ=255)
            gcTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGaa                                                                                                                              <  6:4756/149‑1 (MQ=255)
                   gACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGAc                                                                                                                                    >  3:110636/1‑136 (MQ=255)
                   gACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGAc                                                                                                                                    <  4:110636/136‑1 (MQ=255)
                   gACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCAt                                                                                                                       >  2:214946/1‑149 (MQ=255)
                      aaGATGGGCCCACAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAAATAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCAtgtg                                                                                                                    >  6:174729/1‑149 (MQ=255)
                              cccGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCt                                                                                                            >  4:289940/1‑149 (MQ=255)
                               ccGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCa                                                                                                                                                 >  4:48304/1‑111 (MQ=255)
                               ccGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCa                                                                                                                                                 <  3:48304/111‑1 (MQ=255)
                               ccGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCt                                                                                                            >  2:322944/1‑148 (MQ=255)
                                  cAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATc                                                                                                        <  5:55348/149‑1 (MQ=255)
                                     tAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTg                                                                                                     <  5:97845/149‑1 (MQ=255)
                                         cGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCa                                                                                                 >  3:161001/1‑149 (MQ=255)
                                          gCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCAt                                                                                                <  1:214946/149‑1 (MQ=255)
                                          gCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCAt                                                                                                >  7:86768/1‑149 (MQ=255)
                                          gCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCAt                                                                                                <  3:289940/149‑1 (MQ=255)
                                            tGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGc                                                                                              >  5:325063/1‑149 (MQ=255)
                                              cccGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCt                                                                                             <  6:325063/148‑1 (MQ=255)
                                                    gCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCAGAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTg                                                                                      <  1:69260/149‑1 (MQ=255)
                                                      cTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTggg                                                                                    <  5:236620/149‑2 (MQ=255)
                                                            gACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGc                                                                                                                         >  8:71191/1‑106 (MQ=255)
                                                            gACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGc                                                                                                                         <  7:71191/106‑1 (MQ=255)
                                                               gccgcGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGtt                                                                                                   <  3:77596/125‑1 (MQ=255)
                                                               gccgcGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGtt                                                                                                   >  4:77596/1‑125 (MQ=255)
                                                                  gcgcCCGGAGTGGTCATACCCCACTCCGGATAAAATTCGTTAGCCCGAACAAAAATAACTATTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACTCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCaa                                                                        <  3:258506/149‑1 (MQ=255)
                                                                  gcgcCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGAACAAAAATAACTAGTGGGAGTATTGACAGCATGCAAGCGACTAAGAAAACTCATGTGCGCTACATGATCATGTTCATGCTGTTTCTGGTGACAACGATaaa                                                                        >  4:258506/1‑149 (MQ=255)
                                                                                    ccccACTACGGCTACAATTAGTGAGCCCGCACAAAAATAACTAGAGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGATCATGCTGTTTCTGGTGACCCCGATCAACTATGCCGACCTC‑GCGa                                                      >  6:291836/1‑148 (MQ=255)
                                                                                     cccACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGC‑GCGAc                                                     >  3:122441/1‑148 (MQ=255)
                                                                                     cccACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATCGCCAGCATGCAAGCGACTAAGAAAACGc                                                                                                                         <  7:159730/81‑1 (MQ=255)
                                                                                     cccACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATCGCCAGCATGCAAGCGACTAAGAAAACGc                                                                                                                         >  8:159730/1‑81 (MQ=255)
                                                                                      ccACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGC‑GCGACCa                                                   >  8:70122/1‑149 (MQ=255)
                                                                                       cACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGC‑GCGACCAt                                                  >  3:110761/1‑149 (MQ=255)
                                                                                        aCTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGAGCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGCCCGC‑GCGACCAtt                                                 >  2:42707/1‑149 (MQ=255)
                                                                                         cTCCGGCTACAATTCGTGCGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGC‑GCGACCATTg                                                <  8:86768/149‑1 (MQ=255)
                                                                                                tACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGC‑GCGACCATTGCTATCg                                          >  8:315258/1‑148 (MQ=255)
                                                                                                  cAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGC‑GCGACCATTGCTATCGCa                                        <  1:322944/148‑1 (MQ=255)
                                                                                                   aaTTCGTGAGCCCTCACAAACATAACTAGTGGGAGTATTGCCAGCTTGCAAGCGACTAAGAAAACGCATGTGCGCTACCTTATCCTGTTCATGCTGTTTCTGGTGACCACGATGAACTATGCCGACCGC‑GCGACCATTGCTATCGCAgg                                      <  5:291836/149‑1 (MQ=255)
                                                                                                        gTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGC‑GCGACCATTGCTATCGCAGGCTCCa                                 <  4:122441/149‑1 (MQ=255)
                                                                                                            gCCCGCACAAAAATAACTAGTGGGAGTATCGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGC‑GCGACCATTGCTATCGCAGGCTCCAGCCt                             >  3:22525/1‑149 (MQ=255)
                                                                                                              ccGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATACTGTTTCTGGTGACCACGATCAACTATGCCGACCGC‑GCGACCATTGCTATc                                           <  7:289272/133‑1 (MQ=255)
                                                                                                              ccGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATACTGTTTCTGGTGACCACGATCAACTATGCCGACCGC‑GCGACCATTGCTATc                                           >  8:289272/1‑133 (MQ=255)
                                                                                                                gCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGC‑GCGACCATTGCTATCGCAGGCTCCAGCCTGCAg                         <  7:239018/149‑1 (MQ=255)
                                                                                                                 cacaAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCAGCGACCATTGCTATCGCAGGCTGGAGCCTGcga                         <  1:42707/149‑3 (MQ=255)
                                                                                                                   caAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGC‑GCGACCATTGCTATCGCAGGCTCCAGCCTGCAGaaa                      <  4:110761/149‑1 (MQ=255)
                                                                                                                    aaaaaTAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGC‑GCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAg                     >  3:9379/1‑149 (MQ=255)
                                                                                                                       aaTAACTAGTGGGAGTATCGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGC‑GCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGAcc                  <  4:22525/149‑1 (MQ=255)
                                                                                                                             tAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGC‑GCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCa            >  4:213881/1‑149 (MQ=255)
                                                                                                                                  ggAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGC‑GCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGAc       >  5:14728/1‑149 (MQ=255)
                                                                                                                                  ggAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGC‑GCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGAc       >  5:19744/1‑149 (MQ=255)
                                                                                                                                    aGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGC‑GCGACCATTGCTATc                                           <  5:153660/111‑1 (MQ=255)
                                                                                                                                    aGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGC‑GCGACCATTGCTATc                                           >  6:153660/1‑111 (MQ=255)
                                                                                                                                        ttGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGC‑GCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGt  >  7:107596/1‑148 (MQ=255)
                                                                                                                                         |                                                                                                                                                   
AATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGC‑GCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGT  >  NC_002947/4089806‑4090089

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: ‑