Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,725,055 G→A 8.1% G637G (GGC→GGT)  opgH ← glucans biosynthesis glucosyltransferase H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,725,0550GA 8.1% 137.2 / 2.4 49G637G (GGC→GGT) opgHglucans biosynthesis glucosyltransferase H
Reads supporting (aligned to +/- strand):  ref base G (17/28);  new base A (2/2);  total (19/30)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

CACCACAGGAAGCTCGGGTTCAACCACGCCACCAGTGCGGCCCAGGCAATGCCCAGCAGCGTCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACATGGAGAACAGCATCTCCATCAGCATCGACAAGGTGACCTTGACGCGTCCACCGTACTCGACGGCGCCCTTGGCCCAGATCA  >  NC_002947/5724912‑5725197
                                                                                                                                               |                                                                                                                                              
caccacAGGAAGCTCGGGTTCAACCACGCCACCAGTGCGGCCCAGGCAATGCCCAGCAGCGTCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAg                                                                                                                                           <  6:151699/149‑1 (MQ=255)
  ccacAGGAAGCGCGGGTTCAACCACGCCACCAGTGCGGCCCAGGCAATGCCCAGCAGCGTCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGa                                                                                                                                         <  8:263369/149‑1 (MQ=255)
   cacaGGAAGCTCGGGTTCAACCACGCCACCAGTGCGGCCCAGGCAATGCCCAGCAGCGTCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGaa                                                                                                                                        <  3:271250/149‑1 (MQ=255)
   cacaGGAAGCTCGGGTTCAACCACGCCACCAGTGCGGCCCAGGCAATGCCCAGCAGCGTCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGaa                                                                                                                                        >  6:41610/1‑149 (MQ=255)
      aGGAAGCGCGGGTTCAACCACGCCACCAGTGCGGCCCAGGCAATGCCCAGCAGCGGCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGc                                                                                                                                     <  7:243955/149‑1 (MQ=255)
              cGGGTTCAACCACGCCACCAGTGCGGCCCAGGCAATGCCCAGCAGCGTCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCa                                                                                                                             >  1:21245/1‑149 (MQ=255)
                     aaCCACGCCACCAGTGCGGCCCAGGCAATGCCCAGCAGCGTCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCg                                                                                                                      >  1:158328/1‑149 (MQ=255)
                       ccACGCCACCAGTGCGGCCCAGGCAATGCCCAGCAGCGTCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCggg                                                                                                                    >  7:149445/1‑149 (MQ=255)
                        cACGCCACCAGTGCGGCCCAGGCAATGCCCAGCAGCGTCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGgt                                                                                                                   >  7:288732/1‑149 (MQ=255)
                            acaccaTTGCGGCCCAGGCAATGCCCATCAGCGTCTGCTGTCAATGCCGGCGCACCGATTCGCCCCACGTCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCTAAGACCAGCAGATGAATGCGGACAGAACGAAGCGGGTGTgg                                                                                                               <  2:207/148‑1 (MQ=255)
                             caccaGTGCGGCCCAGGCAATGCCCAGCAGCGTCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGa                                                                                                              <  6:317954/149‑1 (MQ=255)
                                caGTGCGGCCCAGGCAATGCCCAGCAGCGTCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGTCCAGCACGAAGCGGGTGTGGAAGa                                                                                                           <  4:30434/149‑1 (MQ=255)
                                  gTGCGGCCCAGGCAATGCCCAGCAGCGTCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGAt                                                                                                          >  6:325948/1‑148 (MQ=255)
                                       gCCCAGGCAATGCCCAGCAGCGTCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCg                                                                                                    >  8:109566/1‑149 (MQ=255)
                                        cccAGGCAATGCCCAGCAGCGTCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCg                                                                                                    >  1:241078/1‑148 (MQ=255)
                                         ccAGGCAATGCCCAGCAGCGTCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCCCGAAGCGGGTGTGGAAGATCATGCGTa                                                                                                  >  5:162326/1‑149 (MQ=255)
                                          cAGGCAATGCCCAGCAGCGTCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTAc                                                                                                 <  5:41610/149‑1 (MQ=255)
                                            ggCAATGCCCAGCAGCGTCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCg                                                                                               <  4:91655/149‑1 (MQ=255)
                                               aaTGCCCAGCAGCGTCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGaa                                                                                                             >  3:115220/1‑132 (MQ=255)
                                               aaTGCCCAGCAGCGTCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGaa                                                                                                             <  4:115220/132‑1 (MQ=255)
                                                aTGCCCAGCAGCGTCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGc                                                                                           <  1:131733/149‑1 (MQ=255)
                                                   cccAGCAGCGTCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCa                                                                                         <  5:325948/148‑1 (MQ=255)
                                                      agcagcGTCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCaa                                                                                     <  4:185772/149‑1 (MQ=255)
                                                         agcGTCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACAt                                                                                  <  4:166834/149‑1 (MQ=255)
                                                         agcGTCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACAt                                                                                  <  8:53401/149‑1 (MQ=255)
                                                                   ggCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACATGGAGAACAGc                                                                        <  2:21245/149‑1 (MQ=255)
                                                                   ggCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACATGGAGAACAGc                                                                        <  3:50739/149‑1 (MQ=255)
                                                                        tGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACATGGAGAACAGCAtctc                                                                   >  1:121229/1‑149 (MQ=255)
                                                                        tGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACATGGAGAACAGCAtctc                                                                   <  7:109566/149‑1 (MQ=255)
                                                                                 aCCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCCAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGCTCATGCGTACCGGTGCCAGCAACATGGAGAACAGCATCTCCATCAGCAt                                                          >  1:302304/1‑149 (MQ=255)
                                                                                  ccGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACATGGAGAACAGCATCTCCATCAGCATc                                                         <  8:137191/149‑1 (MQ=255)
                                                                                  ccGCTTCGCCCCACGGCGTGGAGTCGGCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACATGGAGAACAGCATCTCCATCAGCATc                                                         <  6:162326/149‑1 (MQ=255)
                                                                                          ccccACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACATGGAGAACAGCATCTCCATCAGCATCGACAAGGt                                                 <  2:64348/149‑1 (MQ=255)
                                                                                          ccccACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACATGGAGAACAGCATCTCCATCAGCATCGACAAGGt                                                 <  2:241078/149‑1 (MQ=255)
                                                                                             cACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGGCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACATGGAGAACAGCATCTCCATCAGCATCGACAAGGTGAc                                              <  8:288732/149‑1 (MQ=255)
                                                                                                ggCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACATGGAGAACAGCCTCTCCATCAGCATCGACAAGGTGACCtt                                           >  2:107717/1‑149 (MQ=255)
                                                                                                 gCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACATGGAGAACAGCATCTCCATCAGCATCGACAAGGTGACCTTg                                          <  2:121229/149‑1 (MQ=255)
                                                                                                     ggAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACATGGAGAACAGCATCTCCATCAGCATCGACAAGGTGACCTTGAcgc                                      <  2:158328/149‑1 (MQ=255)
                                                                                                      gAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAACCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACATGGAGAACAGCATCTCCATCAGCATCGACAAGGTGACCTTg                                          >  5:191885/1‑144 (MQ=255)
                                                                                                      gAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAACCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACATGGAGAACAGCATCTCCATCAGCATCGACAAGGTGACCTTg                                          <  6:196326/144‑1 (MQ=255)
                                                                                                      gAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAACCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACATGGAGAACAGCATCTCCATCAGCATCGACAAGGTGACCTTg                                          <  6:191885/144‑1 (MQ=255)
                                                                                                      gAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAACCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACATGGAGAACAGCATCTCCATCAGCATCGACAAGGTGACCTTg                                          >  5:196326/1‑144 (MQ=255)
                                                                                                        gtcgtcgtcACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACATGGAGAACAGCATCTCCATCAGCATCGACAAGGTGACCTTGACGCGTc                                   <  8:149445/149‑1 (MQ=255)
                                                                                                                   gCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACATGGAGAACAGCATCTCCATCAGCATCGACAAGGTGACCTTGACGCGTCCACCGTACTc                         >  3:160995/1‑148 (MQ=255)
                                                                                                                              ttCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACATGGAGAACAGCATCTCCATCAGCATCGACAAGGTGACCTTGACGCGTCCACCGTACTCGACGGCGCCCtt             >  2:51850/1‑149 (MQ=255)
                                                                                                                                ccAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACATGGAGAACAGCATCTCCATCAGCATCGACAAGGTGACCTTGACGCGTCCACCGTACTCGACGGCGCCCTTgg           >  3:228321/1‑149 (MQ=255)
                                                                                                                                 cAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACATGGAGAACAGCATCTCCATCAGCATCGACAAGGTGACCTTGACGCGTCCACCGTACTCGACGGCGCCCTTgg           <  1:51850/148‑1 (MQ=255)
                                                                                                                                     tCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCg                                                                                               <  4:227449/60‑1 (MQ=255)
                                                                                                                                     tCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCg                                                                                               >  3:227449/1‑60 (MQ=255)
                                                                                                                                         agcccagccGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACATGGAGAACATCATCTCCATCAGCATCGACAAGGTGACCTTGACGCGTCCACCGTACTCGACGGCGCCCTTGGCCCAGATCa  <  8:300075/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                              
CACCACAGGAAGCTCGGGTTCAACCACGCCACCAGTGCGGCCCAGGCAATGCCCAGCAGCGTCTGCGGGCCATGCCGGCGCACCGCTTCGCCCCACGGCGTGGAGTCGTCGTCACGCTGCGGCGAGTTCCAGGTCGCAGCCCAGCCGAGGAATGCGGCCAGCACGAAGCGGGTGTGGAAGATCATGCGTACCGGTGCCAGCAACATGGAGAACAGCATCTCCATCAGCATCGACAAGGTGACCTTGACGCGTCCACCGTACTCGACGGCGCCCTTGGCCCAGATCA  >  NC_002947/5724912‑5725197

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