Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,021,986 C→T 9.6% 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 9.6% 143.3 / 4.2 52P842S (CCG→TCG) PP_2638cellulose synthase operon protein C
Reads supporting (aligned to +/- strand):  ref base C (21/26);  new base T (3/2);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

TGACTACGGCTGGCGCGAGGATGACGGTACGCCCGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAGCCTGTTGGCGTTGGGTTGG  >  NC_002947/3021838‑3022120
                                                                                                                                                      |                                                                                                                                      
tGACTACGGCTGGCGCGAGGATGACGGTACGCCCGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGc                                                                                                                                        <  2:141045/149‑1 (MQ=255)
          tGGCGCGAGGATGACGGTACGCCCGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCa                                                                                                                              <  3:288656/149‑1 (MQ=255)
              gcgAGGATGACGGTACGCCCGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAgg                                                                                                                            <  8:308759/147‑1 (MQ=255)
              gcgAGGATGACGGTACGCCCGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAgg                                                                                                                            >  7:308759/1‑147 (MQ=255)
               cgAGGATGACGGTACGCCCGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAg                                                                                                                             >  5:341540/1‑145 (MQ=255)
               cgAGGATGACGGTACGCCCGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAg                                                                                                                             <  6:341540/145‑1 (MQ=255)
                              gCCCGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGTATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACt                                                                                                          >  8:39140/1‑149 (MQ=255)
                                ccGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTc                                                                                                        >  8:213656/1‑149 (MQ=255)
                                 cGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCTAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTc                                                                                                        >  1:270440/1‑149 (MQ=255)
                                  gggaCTTCGGAACTGAGTAACCATTCGACCTTGTTGCATGGCGACACCCCTTGGCACGAGGGGACCGCGTTTGGGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGACCCGTACGTTCAGGGAAGCTACACCCCCGACTTCgg                                                                                                      <  2:3923/146‑1 (MQ=255)
                                       ttCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGGCTCGAACGacc                                                                                                                               <  3:84780/119‑3 (MQ=255)
                                       ttCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGGCTCGAACGacc                                                                                                                               >  4:84780/1‑117 (MQ=255)
                                       ttCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCt                                                                                                 >  6:358738/1‑149 (MQ=255)
                                         cGGAACTGAGTACCCAGTCGACCTTGTTGCATTTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGc                                                                                               >  6:49291/1‑149 (MQ=255)
                                         cGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCATGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGATCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCt                                                                                                 <  1:82280/148‑1 (MQ=255)
                                          ggAACTGAGTACCCAGTCGACCTTGTTGCATTTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGc                                                                                               <  5:49291/148‑1 (MQ=255)
                                                   tACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTa                                                                                                                     >  3:114083/1‑117 (MQ=255)
                                                   tACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTa                                                                                                                     <  4:114083/117‑1 (MQ=255)
                                                     cccAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGc                                                                                   <  2:270440/149‑1 (MQ=255)
                                                      ccAGTCGACCTTGTTGTATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGcc                                                                                  <  7:39140/149‑1 (MQ=255)
                                                      ccAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGGCTCGAACGacc                                                                                                                               >  3:255827/1‑102 (MQ=255)
                                                      ccAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGGCTCGAACGacc                                                                                                                               <  4:255827/104‑3 (MQ=255)
                                                              ccTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACtt                                                                                                         >  6:139842/1‑118 (MQ=255)
                                                              ccTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCa                                                                                             >  5:285648/1‑130 (MQ=255)
                                                              ccTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCa                                                                                             <  6:285648/130‑1 (MQ=255)
                                                               cTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTacac                                                                                                                  <  5:222257/108‑1 (MQ=255)
                                                               cTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTacac                                                                                                                  >  6:222257/1‑108 (MQ=255)
                                                                  gttgCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGAc                                                                       <  7:213656/148‑1 (MQ=255)
                                                                  gttgCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGAc                                                                       <  3:282550/148‑1 (MQ=255)
                                                                            gACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAg                                                                                                                             >  2:156684/1‑84 (MQ=255)
                                                                            gACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAg                                                                                                                             <  1:156684/84‑1 (MQ=255)
                                                                                 cccTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGc                                                                    >  7:137941/1‑136 (MQ=255)
                                                                                 cccTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGc                                                                    <  8:137941/136‑1 (MQ=255)
                                                                                 cccTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCcg                                                       >  1:188242/1‑149 (MQ=255)
                                                                                  ccTTGGCAGGAGGGCACCGCGTTTGCGCGAGTTGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACtt                                                                                                         <  4:309371/98‑1 (MQ=255)
                                                                                  ccTTGGCAGGAGGGCACCGCGTTTGCGCGAGTTGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACtt                                                                                                         >  3:309371/1‑98 (MQ=255)
                                                                                                  aCCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGc                                                                                   <  1:332401/105‑1 (MQ=255)
                                                                                                  aCCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGc                                                                                   >  2:332401/1‑105 (MQ=255)
                                                                                                       gTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCa                                    >  6:104432/1‑147 (MQ=255)
                                                                                                       gTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCa                                    <  5:104432/147‑1 (MQ=255)
                                                                                                        tttGCGAGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGAc                                 <  2:188242/149‑1 (MQ=255)
                                                                                                                    gagGGCATTGGGATGGATGGGGTGTCGTTCGAGCCGAACGTTCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACAGTGGCCGCGTGCACTGGGGGATCGCAGACTGCAAACGGCa                      <  3:167958/148‑1 (MQ=255)
                                                                                                                          aTTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTg                                                           >  5:309748/1‑105 (MQ=255)
                                                                                                                          aTTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTg                                                           <  6:309748/105‑1 (MQ=255)
                                                                                                                           ttGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAGCCTGTTg              <  5:358738/149‑1 (MQ=255)
                                                                                                                              gCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAg                                                                                                   <  2:47826/61‑1 (MQ=255)
                                                                                                                                        cctgatCGTTCGGCTCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAGCCTGTTGGCgttgggttgg  >  3:93315/6‑149 (MQ=255)
                                                                                                                                           gTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGa                                  <  7:264949/113‑1 (MQ=255)
                                                                                                                                           gTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGa                                  >  8:264949/1‑113 (MQ=255)
                                                                                                                                           gTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAGCCTGTTGGCgttgggtt    >  8:137977/1‑143 (MQ=255)
                                                                                                                                           gTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAGCCTGTTGGCgttgggtt    <  7:137977/143‑1 (MQ=255)
                                                                                                                                                    aGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAGCCTGtt               >  3:31190/1‑124 (MQ=255)
                                                                                                                                                    aGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAGCCTGtt               <  4:31190/124‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                      
TGACTACGGCTGGCGCGAGGATGACGGTACGCCCGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAGCCTGTTGGCGTTGGGTTGG  >  NC_002947/3021838‑3022120

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