Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,021,985 G→C 10.9% E841D (GAG→GAC) ‡ 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,9850GC10.9% 150.6 / 6.7 53E841D (GAG→GAC) ‡PP_2638cellulose synthase operon protein C
Reads supporting (aligned to +/- strand):  ref base G (21/26);  new base C (3/3);  total (24/29)
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.00e+00

ATGACTACGGCTGGCGCGAGGATGACGGTACGCCCGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAGCCTGTTGGCGTTGGGTTGG  >  NC_002947/3021837‑3022120
                                                                                                                                                      |                                                                                                                                       
aTGACTACGGCTGGCGCGAGGATGACGGTACGCCCGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAg                                                                                                                                         >  4:288656/1‑149 (MQ=255)
 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)
                                                                                                                                                      |                                                                                                                                       
ATGACTACGGCTGGCGCGAGGATGACGGTACGCCCGGTACTTCGGAACTGAGTACCCAGTCGACCTTGTTGCATGTCGACACCCCTTGGCAGGAGGGCACCGCGTTTGCGCGAGTAGAGCGCATTGGCATGGATGCGGGGTCGTTCGAGCCGAACGATCAGGGAAGCTACACCCCCGACTTCGGCAGCTGCCAGTTTGTTGGCCAGACCGCTGACGGCAAAACACTGCCCGCGTGCACGGGGGGATCGCAGACTGCAAACGGCAGCCTGTTGGCGTTGGGTTGG  >  NC_002947/3021837‑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: