Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,021,986 C→T 7.5% P842S (CCG→TCG)  PP_2638 → cellulose synthase operon protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,021,9860CT 7.5% 152.6 / 2.4 53P842S (CCG→TCG) PP_2638cellulose synthase operon protein C
Reads supporting (aligned to +/- strand):  ref base C (25/24);  new base T (2/2);  total (27/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 = 1.25e-01

TACGGCTGGCGCGAGGATGACGGTACGCCCGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAGCCTGTTGGCGTTGGGTTGGCAGGGCAGCCGCT  >  NC_002947/3021842‑3022133
                                                                                                                                                |                                                                                                                                                   
tACGGCTGGCGCGAGGATGACGGTACGCCCGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGaa                                                                                                                                                 >  4:86036/1‑149 (MQ=255)
    gCTGGCGCGAGGATGACGGTACGCCCGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGGCTCGAACGac                                                                                                                                             <  2:255718/149‑2 (MQ=255)
     cTGGCGCGAGGATGACGGTACGCCCGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGAt                                                                                                                                             >  3:344744/1‑148 (MQ=255)
        gcgcGAGGATGACGGTACGCCCGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGGCTCGAACGaccc                                                                                                                                           <  8:167363/147‑4 (MQ=255)
         cgcgAGGATGACGGTACGCCCGGTACTTCGGAACTGAGTACCCATTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAgg                                                                                                                                         <  7:355060/148‑1 (MQ=255)
          gcgAGGATGACGGTACGCCCGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAggg                                                                                                                                        >  3:127818/1‑148 (MQ=255)
          gcgAGGATGACGGTACGCCCGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATCCGGGGTCGTTCGAGCCGAACGATCAgg                                                                                                                                         >  8:355060/1‑147 (MQ=255)
            gAGGATGACGGTACGCCCGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGCTCAGGGAAg                                                                                                                                     >  8:254679/1‑149 (MQ=255)
                aTGACGGTACGCCCGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTac                                                                                                                                 <  3:367185/149‑1 (MQ=255)
                   aCGGTACGCCCGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACAcc                                                                                                                              <  6:332788/149‑1 (MQ=255)
                           cccGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCTAGCCGAACGATCAGGGAAGCTACACCCCCGACtt                                                                                                                      <  5:27373/149‑1 (MQ=255)
                            ccGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGTCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTc                                                                                                                     <  2:120881/149‑1 (MQ=255)
                             cGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCg                                                                                                                    <  4:344744/149‑1 (MQ=255)
                             cGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCg                                                                                                                    <  4:127818/149‑1 (MQ=255)
                                       gAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGc                                                                                                            >  2:79093/1‑147 (MQ=255)
                                       gAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCa                                                                                                          <  4:360823/149‑1 (MQ=255)
                                        aaCTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAg                                                                                                         >  8:241865/1‑149 (MQ=255)
                                             aGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTttgt                                                                                                    >  7:337806/1‑149 (MQ=255)
                                              gTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTttgtt                                                                                                   <  5:251560/149‑1 (MQ=255)
                                               tACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCt                                                                                                              <  4:243031/137‑1 (MQ=255)
                                               tACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCt                                                                                                              >  3:243031/1‑137 (MQ=255)
                                               tACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTttgtt                                                                                                   <  3:123850/148‑1 (MQ=255)
                                                   cAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCa                                                                                              <  3:86036/149‑1 (MQ=255)
                                                        gACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCg                                                                                         >  5:303431/1‑149 (MQ=255)
                                                          ccTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCt                                                                                       <  6:303431/149‑1 (MQ=255)
                                                               ttgCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACgg                                                                                  <  4:177474/149‑1 (MQ=255)
                                                                   aTGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCt                                                                                                                                   >  6:288456/1‑96 (MQ=255)
                                                                   aTGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCt                                                                                                                                   <  5:288456/96‑1 (MQ=255)
                                                                    tGTCGACACACCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCAGACCGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCGAGTTTGTTGGCCAGACCGCGGCCGGGaaaa                                                                             >  6:216622/1‑149 (MQ=255)
                                                                           acCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGcc                                                                      >  6:72432/1‑149 (MQ=255)
                                                                             cccTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCcg                                                                    >  1:201402/1‑149 (MQ=255)
                                                                                  ggCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGt                                                                 <  7:241865/147‑1 (MQ=255)
                                                                                    cAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCAc                                                             >  7:159942/1‑149 (MQ=255)
                                                                                     aggaggGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACg                                                            >  5:228952/1‑149 (MQ=255)
                                                                                        aggGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACgggg                                                         <  8:337806/149‑1 (MQ=255)
                                                                                                         cgAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCaaa                                        <  1:79093/149‑1 (MQ=255)
                                                                                                               gagCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAg                                  >  7:2929/1‑149 (MQ=255)
                                                                                                               gagCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAg                                  >  5:269151/1‑149 (MQ=255)
                                                                                                                    cATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAGCCTGt                             >  8:171337/1‑149 (MQ=255)
                                                                                                                      ttGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAGCCTGTTg                           >  3:51648/1‑149 (MQ=255)
                                                                                                                         gCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACggg                                                          >  5:214081/1‑115 (MQ=255)
                                                                                                                         gCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACggg                                                          <  6:214081/115‑1 (MQ=255)
                                                                                                                            tGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAGCCTGTTGGCgttg                     <  4:51648/149‑1 (MQ=255)
                                                                                                                                 gCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAGCCTGTTGGCgttgggttg                >  3:357796/1‑149 (MQ=255)
                                                                                                                                  cGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAGCCTGTTGGCgttgggttgg               <  1:262619/149‑1 (MQ=255)
                                                                                                                                   ggggTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAGCCTGTTGGCgttgggttgg               >  3:308798/1‑148 (MQ=255)
                                                                                                                                   ggggTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAGCCTGTTGGCGTTGGGTTggc              <  2:321845/149‑1 (MQ=255)
                                                                                                                                   ctgatCGTTCGGCTCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAGCCTGTTGGCgttgggttgg               >  2:71972/5‑148 (MQ=255)
                                                                                                                                    tgatCGTTCGGCTCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAGCCTGTTGGCgttgggttgg               >  1:252604/4‑147 (MQ=255)
                                                                                                                                    gggTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAGCCTGTTGGCGTTGGGTTggca             <  6:269151/149‑1 (MQ=255)
                                                                                                                                       tCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAGCCTGTTGGCGTTGGGTTggcaggg          <  7:171337/149‑1 (MQ=255)
                                                                                                                                         gTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAGCCTGTTGGCGTTGGGTTggcaggg          >  4:286909/1‑147 (MQ=255)
                                                                                                                                               gCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAGCCTGTTGGCGTTGGGTTGGCAGGGCAGCCGCt  >  8:127865/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TACGGCTGGCGCGAGGATGACGGTACGCCCGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAGCCTGTTGGCGTTGGGTTGGCAGGGCAGCCGCT  >  NC_002947/3021842‑3022133

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