Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,235,841 G→C 10.2% G318R (GGG→CGG)  uvrB → uvrABC excinuclease protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,235,8410GC10.2% 137.3 / 7.7 49G318R (GGG→CGG) uvrBuvrABC excinuclease protein B
Reads supporting (aligned to +/- strand):  ref base G (21/23);  new base C (2/3);  total (23/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 = 7.43e-01

GCGTCTGGAGTACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGAGTACGGCTTCC  >  NC_002947/2235693‑2235985
                                                                                                                                                    |                                                                                                                                                
gCGTCTGGAGTACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCccggccg                                                                                                                                                  >  8:39685/1‑149 (MQ=255)
gCGTCTGGAGTACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCccggccg                                                                                                                                                  <  7:101739/149‑1 (MQ=255)
   tCTGGAGTACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCgggg                                                                                                                                               <  7:130739/149‑1 (MQ=255)
      ggAGTACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGgcgc                                                                                                                                            >  7:146189/1‑149 (MQ=255)
       gAGTACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGcgcc                                                                                                                                           >  5:241932/1‑149 (MQ=255)
       gAGTACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGcgcc                                                                                                                                           <  6:241932/149‑1 (MQ=255)
        aGTACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCCGGCCGGGcgccc                                                                                                                                          <  6:281069/149‑2 (MQ=255)
             cTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCCGGCGCCCCGGCCGGGcgcccggaca                                                                                                                                     <  2:306352/149‑7 (MQ=255)
               gCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCa                                                                                                                                   >  3:215056/1‑149 (MQ=255)
                cacaAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCAc                                                                                                                                  >  8:112151/1‑149 (MQ=255)
                 acaAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACACGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCAcc                                                                                                                                 >  3:258696/1‑149 (MQ=255)
                  caAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCAccc                                                                                                                                >  2:28418/1‑149 (MQ=255)
                   aaGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCt                                                                                                                               <  3:183634/149‑1 (MQ=255)
                     ggCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGt                                                                                                                             <  3:215963/149‑1 (MQ=255)
                      gCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTa                                                                                                                            <  5:83301/149‑1 (MQ=255)
                       gcaacaaGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTAc                                                                                                                           >  2:263984/2‑149 (MQ=255)
                        caacaaGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTAc                                                                                                                           >  6:265275/1‑148 (MQ=255)
                            aaGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCCGGCGCCCCGGCCGGGcgcccg                                                                                                                                         >  3:141428/1‑128 (MQ=255)
                            aaGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCCGGCGCCCCGGCCGGGcgcccg                                                                                                                                         <  4:141428/130‑3 (MQ=255)
                             aGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTAc                                                                                                                     <  1:305716/149‑1 (MQ=255)
                                  gTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGACGCCGCCGCCCACCCTGTACGATTACCTgc                                                                                                                 <  1:28418/148‑1 (MQ=255)
                                           cAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATgc                                                                                                       >  4:307044/1‑149 (MQ=255)
                                             gCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATgcgc                                                                                                     <  4:215056/149‑1 (MQ=255)
                                                    gAGCAACGCACCCGCTTCGTTCTGGCGATTATCCTCGAGGTGGGATACTGCAACGGCATCTAAAACTACTCGGGTTACCTGTCCGGGCGCCGGGCGGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGt                                                                                              <  4:307422/149‑1 (MQ=255)
                                                    gAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGTCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGctgctg                                                                                                <  4:258696/147‑1 (MQ=255)
                                                    gAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGctgctg                                                                                                <  8:233686/147‑1 (MQ=255)
                                                    gAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGt                                                                                              >  7:261895/1‑149 (MQ=255)
                                                     aGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTg                                                                                             <  1:203754/149‑1 (MQ=255)
                                                      gCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGa                                                                                            <  7:112151/149‑1 (MQ=255)
                                                                     cGATCTGGAGATGATCCTGGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCTCCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACTAATCacac                                                                              <  2:35252/148‑1 (MQ=255)
                                                                        tCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGt                                                                            >  8:2766/1‑147 (MQ=255)
                                                                                tgagCCTCGAGCTGGGCTACTGCACCGGAATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTcc                                                                   <  8:146189/148‑1 (MQ=255)
                                                                                 gatCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCg                                                                  >  3:88710/1‑148 (MQ=255)
                                                                                              ggCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGcc                                                      <  4:88710/147‑1 (MQ=255)
                                                                                               gCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATg                                                   >  6:245428/1‑149 (MQ=255)
                                                                                                      cAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGGCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATgcg                                                                                                      <  1:50761/91‑1 (MQ=255)
                                                                                                      cAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAgg                                            <  7:2766/149‑1 (MQ=255)
                                                                                                      cAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAg                                             <  5:245428/148‑1 (MQ=255)
                                                                                                      cAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCCCCGCCCCCCACCCTCTACGATTACCTGCCGCCCGATgcg                                                                                                      >  2:50761/1‑91 (MQ=255)
                                                                                                             aTCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCg                                     >  8:81126/1‑149 (MQ=255)
                                                                                                               cGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGt                                    >  5:152057/1‑148 (MQ=255)
                                                                                                                  aaaCTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCAc                                >  6:229074/1‑149 (MQ=255)
                                                                                                                  aaaCTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCAc                                >  5:65011/1‑149 (MQ=255)
                                                                                                                        ctCGCGTTACCTGTCCGGGCGCCGGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAg                           <  8:261895/148‑1 (MQ=255)
                                                                                                                        ctCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATcac                                                                                >  2:225024/1‑95 (MQ=255)
                                                                                                                        ctCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATcac                                                                                <  1:225024/95‑1 (MQ=255)
                                                                                                                                aCCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTg                  >  7:42199/1‑149 (MQ=255)
                                                                                                                                aCCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTg                  >  8:16424/1‑149 (MQ=255)
                                                                                                                                 ccTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTg                  <  7:39685/148‑1 (MQ=255)
                                                                                                                                        gggCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGAGTa          <  7:16424/149‑1 (MQ=255)
                                                                                                                                          gcgcCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGGGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGAGTACg        <  1:263984/149‑1 (MQ=255)
                                                                                                                                          gcgcCCCGGCCGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGAGTACg        >  6:193143/1‑149 (MQ=255)
                                                                                                                                              ccggccggGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGAGTACGGCtt    <  3:307044/149‑1 (MQ=255)
                                                                                                                                                ggccggGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGAGTACGGCTTcc  >  6:295194/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GCGTCTGGAGTACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGAGTACGGCTTCC  >  NC_002947/2235693‑2235985

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