Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,843,500 G→C 8.3% G41A (GGC→GCC)  PP_5122 → GMC family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,843,5000GC 8.3% 145.5 / 2.7 48G41A (GGC→GCC) PP_5122GMC family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (24/20);  new base C (2/2);  total (26/22)
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.85e-01

CAAGAGCATCGCCAATGCCTGTCCCCGACCCGTTTCGCCAAGGCCTCGAACGCGGCTGGATCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGAAGACCAGCAGCGATTTCCACCTGCTGGAGAGCGAAGCCTACGCCAGCCTGTACCAGGAAGGGCTGGGGCGCATGAG  >  NC_002947/5843365‑5843647
                                                                                                                                       |                                                                                                                                                   
caaGAGCATCGCCAATGCCTGTCCCCGACCCGTTTCGCCAAGGCCTCGAACGCGGCTGGATCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGt                                                                                                                                        >  6:138725/1‑149 (MQ=255)
caaGAGCATCGCCAATGCCTGTCCCCGACCCGTTTCGCCAAGGCCTCGAACGCGGCTGGATCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGCGCCGCCGGcaccactg                                                                                                                                          >  4:59073/1‑145 (MQ=255)
caaGAGCATCGCCAATGCCTGTCCCCGACCCGTTTCGCCAAGGCCTCGAACGCGGCTGGATCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGCGCCGCCGGcaccactg                                                                                                                                          <  3:59073/147‑3 (MQ=255)
          gCCAATGCCTGTCCCCGACCCGTTTCGCCAAGGCCTCGAACGCGGCTGGATCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATg                                                                                                                               <  5:155864/148‑1 (MQ=255)
            cAATGCCTGTCCCCGACCCGTTTCGCCAAGGCCTCGAACGCGGCTGGATCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTc                                                                                                                            <  1:53978/149‑1 (MQ=255)
                       cccGACCCGTTTCGCCAAGGCCTCGAACGCGGCTGGATCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCgg                                                                                                                 >  8:19707/1‑149 (MQ=255)
                          gACCCGTTTCGCCAAGGCCTCGAACGCGGCTGGATCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGttt                                                                                                              <  7:60524/149‑1 (MQ=255)
                           aCCCGTTTCGCCAAGGCCTCGAACGCGGCTGGATCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGCGCCGCCGGCACCACCAGTGCGCAGATGCTCAGCGCTGCCGGttt                                                                                                              >  8:23548/1‑148 (MQ=255)
                             ccGTTTCGCCAAGGCCTCGAACGCGGCTGGATCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTc                                                                                                             >  2:217588/1‑147 (MQ=255)
                             ccGTTTCGCCAAGGCCTCGAACGCGGCTGGATCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTc                                                                                                             >  3:257973/1‑147 (MQ=255)
                             ccGTTTCGCCAAGGCCTCGAACGCGGCTGGATCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCaa                                                                                                           <  7:126200/149‑1 (MQ=255)
                               gTTTCGCCAAGGCCTCGAACGCGGCTGGATCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAgg                                                                                                         <  5:138725/149‑1 (MQ=255)
                                tttCCCCAAGGCCTCGAACGCGGCTGGATCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGc                                                                                                                    >  8:79295/1‑137 (MQ=255)
                                tttCCCCAAGGCCTCGAACGCGGCTGGATCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGc                                                                                                                    <  7:79295/137‑1 (MQ=255)
                                  tCGCCAAGGCCTCGAACGCGGCTGGATCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTa                                                                                                       <  4:257973/148‑1 (MQ=255)
                                    gCCAAGGCCTCGAACGCGGCTGGATCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTActg                                                                                                    >  1:69309/1‑149 (MQ=255)
                                     ccAAGGCCTCGAACGCGGCTGGATCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTActgc                                                                                                   <  7:216432/149‑1 (MQ=255)
                                      cAAGGCCTCGAACGCGGCTGGATCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTActgct                                                                                                  >  1:232324/1‑149 (MQ=255)
                                       aaGGCCTCGAACGCGGCTGGATCAACCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTActgct                                                                                                  >  4:272943/1‑148 (MQ=255)
                                        aGGCCTCGAACGCGGCTGGATCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGa                                                                                                <  3:83575/149‑1 (MQ=255)
                                                   gcgGCTGGATCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGCGCCGCCGGCACCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGc                                                                                     <  7:23548/149‑1 (MQ=255)
                                                     ggCTGGATCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGcc                                                                                    >  2:104071/1‑148 (MQ=255)
                                                       cTGGATCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCAc                                                                                 <  7:19707/149‑1 (MQ=255)
                                                         ggATCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACt                                                                                >  6:140857/1‑148 (MQ=255)
                                                           aTCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGaa                                                                             <  3:272943/149‑1 (MQ=255)
                                                                ccACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGAAGAc                                                                          >  1:212945/1‑147 (MQ=255)
                                                                ccACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGAAGAc                                                                          >  6:290146/1‑147 (MQ=255)
                                                                    gacgCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGAAGACcagcag                                                                    <  1:104071/149‑1 (MQ=255)
                                                                       gCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGAAGACCAGCAGCGa                                                                 >  6:259033/1‑149 (MQ=255)
                                                                             cgcCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGAAGACCAGCAGCGATTTCCa                                                           <  2:212945/149‑1 (MQ=255)
                                                                              gcCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATTCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGAATACCAGCAGCGATTTCCAc                                                          <  2:232324/149‑1 (MQ=255)
                                                                                       aCGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCTGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGAAGACCAGCAGCGATTTCCACCTGCTGgag                                                 <  5:290146/149‑1 (MQ=255)
                                                                                        cGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGAAGACCAGCAGCGATTTCCACCTGCTGga                                                  >  5:308577/1‑147 (MQ=255)
                                                                                         gACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTACCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGAAGACCAGCAGCGATTTCCACCTGCTGGagag                                               >  3:135539/1‑149 (MQ=255)
                                                                                              gTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGAAGACCAGCAGCGATTTCCACCTGCTGGAGAGCGAAg                                          >  2:257223/1‑149 (MQ=255)
                                                                                               tttCTGGATGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCGCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGAAGACCAGCAGCGATTTCCACCTGCTGGAGAGCGAAGc                                         <  8:158495/149‑1 (MQ=255)
                                                                                               tttCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGAAGACCAGCAGCGATTTCCACCTGCTGGAGAGCGaa                                           <  5:259033/147‑1 (MQ=255)
                                                                                               tttCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGAAGACCAGCAGCGATTTCCACCTGCTGGAGAGCGAAGc                                         >  4:279421/1‑149 (MQ=255)
                                                                                                  cTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCTACTGAAGACCAGCAGCGATTTCCACCTGCTGGAGAGCGAAGCCt                                       <  1:217588/148‑1 (MQ=255)
                                                                                                  cTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGAAGACCAGCAGCGATTTCCACCTGCTGGAGAGCGAAGCCTa                                      >  6:37877/1‑149 (MQ=255)
                                                                                                     gAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGAAGACCAGCAGCGATTTCCACCTGCTGGAGAGCGAAGCCTACGc                                   >  5:144205/1‑149 (MQ=255)
                                                                                                          cGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGAAGACCAGCAGCGATTTCCACCTGCTGGAGAGCGAAGCCTACGCCAGcc                              >  8:49931/1‑149 (MQ=255)
                                                                                                           gACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGAAGACCAGCAGCGATTTCCACCTGCTGGAGAGCGAAGCCTACGCCAGCCt                             >  8:195222/1‑149 (MQ=255)
                                                                                                              aTTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGAAGACCAGCAGCGATTTCCACCTGCTGGAGAGCGAAGCCTACGCCAGCCTGTa                          >  5:231799/1‑149 (MQ=255)
                                                                                                              aTTGCCGTGATCGGCAGTGGTGCCGGCGCCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGAAGACCAGCAGCGATTTCCACCTGCTGGAGAGCGAAGCCTACGCCAGCCTGTa                          >  1:175686/1‑149 (MQ=255)
                                                                                                                  ccGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTACCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGAAGACCAGCAGCGATTTCCACCTGCTGGAGAGCGAAGCCTACGCCAGCCTGTACCa                       <  4:135539/148‑1 (MQ=255)
                                                                                                                           gCAGTGGTGCCGGCGCCCCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTCCAGAACGAAGAAGGCCCACTGAAGACCAGCAGCGATTTCCACATGCTGGAGAGCGAAGCCTACCCCAGCCTGTACCAGGAAGGGCTg             >  5:175069/1‑149 (MQ=255)
                                                                                                                                 gtgCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGAAGACCAGCAGCGATTTCCACCTGCTGGAGAGCGAAGCCTACGCCAGCCTGTACCAGGAAGGGCTGGGgcg        <  6:144205/148‑1 (MQ=255)
                                                                                                                                     cggcggCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGAAGACCAGCAGCGATTTCCACCTGCTGGAGAGCGAAGCCTACGCCAGCCTGTACCAGGAAGGGCTGGGGCGCATGa   <  2:69309/149‑1 (MQ=255)
                                                                                                                                      ggcggcGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGAAGACCAGCAGCGATTTCCACCTGCTGGAGAGCGAAGCCTACGCCAGCCTGTACCAGGAAGGGCTGGGGCGCATGAg  <  6:308577/149‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                                   
CAAGAGCATCGCCAATGCCTGTCCCCGACCCGTTTCGCCAAGGCCTCGAACGCGGCTGGATCACCCACGACGCCTCGCGCCTGGAGCACGACCTGTTTCTGGAGGCCGACATTGCCGTGATCGGCAGTGGTGCCGGCGGCGCCACCAGTGCGCAGATGCTCAGCGCTGCCGGTTTCAAGGTACTGCTGATCGAAGAAGGCCCACTGAAGACCAGCAGCGATTTCCACCTGCTGGAGAGCGAAGCCTACGCCAGCCTGTACCAGGAAGGGCTGGGGCGCATGAG  >  NC_002947/5843365‑5843647

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