Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,731,528 G→C 10.7% P781R (CCG→CGG)  PP_1526 ← beta‑(1‑3)‑glucosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,731,5280GC10.7% 153.9 / 9.2 56P781R (CCG→CGG) PP_1526beta‑(1‑3)‑glucosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (30/20);  new base C (3/3);  total (33/23)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.49e-01

TCGCCACCCAGAAGCGCATGTCCGAACTCGGCAGCCCCTGCACCAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCACCCCCACCGCTCGGCGGTACAGGAAGATGATCTT  >  NC_002947/1731382‑1731676
                                                                                                                                                  |                                                                                                                                                    
tCGCCACCCAGAAGCGCATGTCCGAACTCGGCAGCCCCTGCACCAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCCGc                                                                                                                                                    >  6:307755/1‑149 (MQ=255)
    cACCCAGAAGCGCATGTCCGAACTCGGCAGCCCCTGCACCAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGc                                                                                                                                                <  4:10580/149‑1 (MQ=255)
      cccAGAAGCGCATGTCCGAACTCGGCAGCCCCTGCACCAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGc                                                                                                                                                <  8:155801/147‑1 (MQ=255)
      cccAGAAGCGCATGTCCGAACTCGGCAGCCCCTGCACCAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGc                                                                                                                                                >  7:155801/1‑147 (MQ=255)
             gcgcATGTCCGAACTCGGCAGCCCCTGCACCAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCCGCACGCCGAAGAACg                                                                                                                                       <  5:307755/149‑1 (MQ=255)
             gcgcATGTCCGAACTCGGCAGCCCCTGCACCAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCCGCACGCCGAAGAACg                                                                                                                                       <  1:254135/149‑1 (MQ=255)
                cATGTCCGAACTCGGCAGCCCCTGCACCAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCa                                                                                                                                    <  7:216280/149‑1 (MQ=255)
                    tCCGAACTCGGCAGCCCCTGCACCAGGTAGATGCCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCAtgct                                                                                                                                >  7:142047/1‑149 (MQ=255)
                    tCCGAACTCGGCAGCCCCTGCACCAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCAtgct                                                                                                                                >  7:228088/1‑149 (MQ=255)
                    tCCGAACTCGGCAGCCCCTGCACCAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCAtgct                                                                                                                                >  6:30827/1‑149 (MQ=255)
                       gAACTCGGCAGCCCCTGCACCAGGTAGATAGCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCAtgctgc                                                                                                                              >  1:30176/1‑148 (MQ=255)
                       gAACTCGGCAGCCCCTGCACCAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCAtgctgc                                                                                                                              <  2:30176/148‑1 (MQ=255)
                                        cACCAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATa                                                                                                              <  4:162465/147‑1 (MQ=255)
                                           cAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCCGCACGCCGAAGAACg                                                                                                                                       >  2:139414/1‑119 (MQ=255)
                                           cAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCCGCACGCCGAAGAACg                                                                                                                                       <  1:139414/119‑1 (MQ=255)
                                            aGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTAc                                                                                                        >  1:221269/1‑149 (MQ=255)
                                                      ccAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGctgctg                                                                                                                            >  3:227863/1‑119 (MQ=255)
                                                      ccAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGctgctg                                                                                                                            <  4:227863/119‑1 (MQ=255)
                                                      ccAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATAc                                                                                                             >  5:262897/1‑134 (MQ=255)
                                                      ccAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATAc                                                                                                             <  6:262897/134‑1 (MQ=255)
                                                      ccAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGc                                                                                               >  4:138326/1‑148 (MQ=255)
                                                         aCGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATg                                                                                           <  5:30827/149‑1 (MQ=255)
                                                           gccgcGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGgt                                                                                         >  6:301665/1‑149 (MQ=255)
                                                             cgcgCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGgtgt                                                                                       <  2:14079/149‑1 (MQ=255)
                                                              gcgcCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATa                                                                                                              <  2:68283/125‑1 (MQ=255)
                                                              gcgcCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATa                                                                                                              >  1:68283/1‑125 (MQ=255)
                                                                  cccACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGa                                                                                   >  5:213031/1‑148 (MQ=255)
                                                                    cacaGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTAcc                                                                                                       <  3:33008/126‑1 (MQ=255)
                                                                    cacaGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTAcc                                                                                                       >  4:33008/1‑126 (MQ=255)
                                                                    cacaGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAg                                                                                >  3:175550/1‑149 (MQ=255)
                                                                     acaGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGt                                                                               <  8:142047/149‑1 (MQ=255)
                                                                           gcaCCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGGGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGAAGGTGTGCGAAAGTGCCAGc                                                                         >  2:96144/1‑149 (MQ=255)
                                                                                     gaACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGAGACAGTGACAGCCCGGCCAgcg                                                               >  4:116928/1‑149 (MQ=255)
                                                                                       aCAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGt                                                                               >  6:166967/1‑131 (MQ=255)
                                                                                       aCAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGt                                                                               <  5:166967/131‑1 (MQ=255)
                                                                                           ctctTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTg                                                         >  5:186779/1‑149 (MQ=255)
                                                                                             ctTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGc                                                        >  5:294506/1‑148 (MQ=255)
                                                                                               tCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCg                                                       >  1:83628/1‑147 (MQ=255)
                                                                                                cACGGGCTTCGGAAATCGCCACCAGCAACGGTTGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTTAAGAAACCATACAGTACCGCCTTGGCGTTGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGa                                                      <  8:197190/147‑1 (MQ=255)
                                                                                                         cGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCtt                                            <  2:83628/148‑1 (MQ=255)
                                                                                                           gAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAg                                         >  2:161050/1‑149 (MQ=255)
                                                                                                                    accaGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGCAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCAcccc                                >  3:181627/1‑149 (MQ=255)
                                                                                                                             ccATGGCTGTCAGCGTTCTTCGGCGTGCAGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCACCCCCACCGCTcg                       >  1:183018/1‑149 (MQ=255)
                                                                                                                             ccATGGCTGTCAGCGTTCTTCGGCGTGCAGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTCGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCACCCCCACCGCTcg                       >  1:188076/1‑149 (MQ=255)
                                                                                                                               aTGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCACCCCCACCGCTcggc                     >  8:168991/1‑149 (MQ=255)
                                                                                                                                  gCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCACCCCCACCGCTcggcgg                   <  2:221269/148‑1 (MQ=255)
                                                                                                                                   cTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCACCCCCACCGCTCGGCGGTa                 <  3:138326/149‑1 (MQ=255)
                                                                                                                                      tCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCACCCCCACCGCTCGGCGGTACa               <  7:168991/148‑1 (MQ=255)
                                                                                                                                        aGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATcc                                              <  5:33488/115‑1 (MQ=255)
                                                                                                                                        aGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATcc                                              >  6:33488/1‑115 (MQ=255)
                                                                                                                                        aGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCACCCCCACCGCTCGGCGGTACAgg             >  4:44388/1‑148 (MQ=255)
                                                                                                                                        aGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCACCCCCACCGCTCGGCGGTACAgg             >  3:208917/1‑148 (MQ=255)
                                                                                                                                         gCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCACCCCCACCGCTCGGCGGTACAgg             >  4:119288/1‑147 (MQ=255)
                                                                                                                                          cGTTCTTCCGCACGCCGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCACCCCCACCGCTCGGCGGTACAgg             >  6:58970/1‑146 (MQ=255)
                                                                                                                                               ttcGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCACCCCCACCGCTCGGCGGTACAGGAAgatgat     >  2:228504/1‑149 (MQ=255)
                                                                                                                                                 cGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCACCCCCACCGCTCGGCGGTACAGGAAGATGATc    <  8:228088/148‑1 (MQ=255)
                                                                                                                                                  ggCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCACCCCCACCGCTCGGCGGTACAGGAAGATGATCtt  <  4:175550/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TCGCCACCCAGAAGCGCATGTCCGAACTCGGCAGCCCCTGCACCAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCACCCCCACCGCTCGGCGGTACAGGAAGATGATCTT  >  NC_002947/1731382‑1731676

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