Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,090,101 C→G 12.8% S36S (TCC→TCG) ‡ gudP → glucarate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,090,1010CG12.8% 103.6 / 8.2 39S36S (TCC→TCG) ‡gudPglucarate transporter
Reads supporting (aligned to +/- strand):  ref base C (21/13);  new base G (3/2);  total (24/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

GCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCTTCAGCATCTTCAC  >  NC_002947/4089962‑4090250
                                                                                                                                           |                                                                                                                                                     
gCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGa                                                                                                                                             >  2:300050/1‑150 (MQ=255)
  acaAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGaaa                                                                                                                                           <  2:192152/150‑1 (MQ=255)
               gTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTg                                                                                                                              >  2:209967/1‑150 (MQ=255)
                aggcaGTATTGCCCGCATGCAAGCGACTAAGAACACGCGTGTGCGCTACCTGATCCGGTTCATGGTGTTTCTGGTGACCACGATCACCTATGCCGAGCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGa                                                                                                                             <  2:30601/146‑1 (MQ=255)
                  ggAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACg                                                                                                                           >  2:381486/1‑150 (MQ=255)
                  ggAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACg                                                                                                                           <  1:300050/150‑1 (MQ=255)
                               cATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTTTTCATTCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTGGAGCCTGcgata                                                                                                                                            <  2:839821/120‑5 (MQ=255)
                               cATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTGGAGCCTGcgata                                                                                                                                            >  1:839821/1‑116 (MQ=255)
                                 tGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGtt                                                                                                            >  1:258034/1‑150 (MQ=255)
                                                    gCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCgg                                                                                         >  1:62961/1‑150 (MQ=255)
                                                        gtgCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGCACTGGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATgg                                                                                     <  1:650139/150‑1 (MQ=255)
                                                         tgCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATggg                                                                                    <  1:209967/150‑1 (MQ=255)
                                                            gcgaCCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCTACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGc                                                                                                                          <  1:149379/106‑1 (MQ=255)
                                                            gcTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGc                                                                                                                          >  2:149379/1‑109 (MQ=255)
                                                                ccTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGt                                                                             >  2:748196/1‑150 (MQ=255)
                                                                  tGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGtgg                                                                           >  2:832922/1‑150 (MQ=255)
                                                                     tCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGc                                                                                   >  1:26109/1‑139 (MQ=255)
                                                                     tCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGc                                                                                   <  2:26109/139‑1 (MQ=255)
                                                                       cTGTTCATGCTGTTTCTGGTGACCACGATCAACTAAGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGGAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGtggctggc                                                                      >  1:274246/1‑150 (MQ=255)
                                                                       cTGTTCATGCTGTTTCTGGCGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGtggctggc                                                                      <  2:563169/150‑1 (MQ=255)
                                                                             aTGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGAtt                                                                >  2:341286/1‑150 (MQ=255)
                                                                                 tGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCtg                                                            >  1:62915/1‑150 (MQ=255)
                                                                                 tGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCtg                                                            >  2:1454/1‑150 (MQ=255)
                                                                                              cACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGc                                                                                   >  2:133528/1‑114 (MQ=255)
                                                                                              cACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGc                                                                                   <  1:133528/114‑1 (MQ=255)
                                                                                                  aTCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGAt                                           <  1:381486/150‑1 (MQ=255)
                                                                                                     aaCTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCt                                       >  2:555706/1‑150 (MQ=255)
                                                                                                             cGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTc                                >  1:212375/1‑150 (MQ=255)
                                                                                                               aCCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGtggctgg                                                    <  1:28576/128‑1 (MQ=255)
                                                                                                               aCCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGtggctgg                                                    >  2:28576/1‑128 (MQ=255)
                                                                                                               aCCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCa                              >  2:379945/1‑150 (MQ=255)
                                                                                                                  gcgcgACCATTGCTATCGCACGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCCTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGa                           >  2:626125/1‑150 (MQ=255)
                                                                                                                      gACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGt                       <  1:341286/150‑1 (MQ=255)
                                                                                                                           ttGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACg                  >  1:170330/1‑150 (MQ=255)
                                                                                                                              tttctgCAGGCTGGAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCt               >  2:433449/6‑150 (MQ=255)
                                                                                                                               tATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCtt              >  2:112032/1‑150 (MQ=255)
                                                                                                                                  tgCAGGCTGGAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGtt                                 <  2:610654/127‑1 (MQ=255)
                                                                                                                                  tgCAGGCTGGAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGtt                                 >  1:610654/2‑128 (MQ=255)
                                                                                                                                        cTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCTTCAGCATCTTc    <  1:555706/150‑1 (MQ=255)
                                                                                                                                           cAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCaa                             <  1:717377/123‑1 (MQ=255)
                                                                                                                                           cAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCaa                             >  2:717377/1‑123 (MQ=255)
                                                                                                                                           cAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCTTCAGCATCTTCAc  <  2:62961/150‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                     
GCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCTTCAGCATCTTCAC  >  NC_002947/4089962‑4090250

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