Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,656,958 G→T 8.0% P263Q (CCG→CAG)  tktA ← transketolase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,656,9580GT 8.0% 151.2 / 3.6 50P263Q (CCG→CAG) tktAtransketolase A
Reads supporting (aligned to +/- strand):  ref base G (24/22);  new base T (2/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

GCGTCGAAGCGCTGGTTCCACTCGGCTTCAACCTTGGCACCGGCGGCCTTGGCATCCCACTCGGCGTAGATGTCGGCCGGGATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCGGGCGATCGCTCTTGCGGGCGGTCTCGATGGCCATCTTGATTTCTTCGGCATCGTGGCCGT  >  NC_002947/5656811‑5657096
                                                                                                                                                   |                                                                                                                                          
gCGTCGAAGCGCTGGTTCCACTCGGCTTCAACCTTGGCACCGGCGGCCTTGGCATCCCACTCGGCGTAGATGTCGGCCGGGATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCggg                                                                                                                                          <  1:319661/150‑1 (MQ=255)
 cGTCGAAGCGCTGGTTCCACTCGGCTTCAACCTTGGCACCGGCGGCCTTGGCATCCCACTCGGCGTAGATGTCGGCCGGGATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCgggg                                                                                                                                         >  2:138965/1‑150 (MQ=255)
   tCGAAGCGCTGGTTCCACTCGGCTTCAACCTTGGCACCGGCGGCCTTGGCATCCCACTCGGCGTAGATGTCGGCCGGGATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCCGCGGGGCa                                                                                                                                       >  2:487518/1‑150 (MQ=255)
     gAAGCGCTGGTTCCACTCGGCTTCAACCTTGGCACCGGCGGCCTTGGCATCCCACTCGGCGTAGATGTCGGCCGGGATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCAcc                                                                                                                                     >  1:455854/1‑150 (MQ=255)
           cTGGTTCCACTCGGCTTCAACCTTGGCACCGGCGGCCTTGGCATCCCACTCGGCGTAGATGTCGGCCGGGATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCgg                                                                                                                                           >  1:670536/1‑138 (MQ=255)
           cTGGTTCCACTCGGCTTCAACCTTGGCACCGGCGGCCTTGGCATCCCACTCGGCGTAGATGTCGGCCGGGATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCgg                                                                                                                                           <  2:670536/138‑1 (MQ=255)
             ggTTCCACTCGGCTTCAACCTTGGCACCGGCGGCCTTGGCATCCCACTCGGCGTAGATGTCGGCCGGGATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGt                                                                                                                             >  2:993396/1‑150 (MQ=255)
                   aCTCGGCTTCAACCTTGGCACCGGCGGCCTTGGCATCCCACTCGGCGTAGATGTCGGCCGGGATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCt                                                                                                                       >  1:274309/1‑150 (MQ=255)
                   aCTCGGCTTCAACCTTGGCACCGGCGGCCTTGGCATCCCACTCGGCGTAGATGTCGGCCGGGATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCt                                                                                                                       >  1:567228/1‑150 (MQ=255)
                          ttCAACCTTGGCACCGGCGGCCTTGGCATCCCACTCGGCGTAGATGTCGGCCGGGATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCCGCTGGGCaacga                                                                                                                                   >  1:1182084/1‑127 (MQ=255)
                          ttCAACCTTGGCACCGGCGGCCTTGGCATCCCACTCGGCGTAGATGTCGGCCGGGATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCCGCTGGGCAACga                                                                                                                                   <  2:1182084/131‑2 (MQ=255)
                             aaCCTTGGCACCGGCGGCCTTGGCATCCCACTCGGCGTAGATGTCGGCCGGGATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGc                                                                                                               >  2:179447/1‑148 (MQ=255)
                             aaCCTTGGCACCGGCGGCCTTGGCATCCCACTCGGCGTAGATGTCGGCCGGGATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGc                                                                                                               <  1:179447/148‑1 (MQ=255)
                                   ggCACCGGCGGCCTTGGCATCCCACTCGGCGTAGATGTCGGCCGGGATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCgg                                                                                                       <  2:274309/150‑1 (MQ=255)
                                                 tGGCATCCCACTCGGCGTAGATGTCGGCCGGGATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGTGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCGTCCTTGCCCTGCTTGTTCGGCGaaccgaaaccga                                                                                         <  2:1071808/150‑1 (MQ=255)
                                                   gCATCCCACTCGGCGTAGATGTCGGCCGGGATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGAtg                                                                                       >  1:311416/1‑150 (MQ=255)
                                                      tCCCACTCGGCGTAGATGTCGGCCGGGATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGAtggtg                                                                                    <  2:311416/150‑1 (MQ=255)
                                                        caaCTCGGCGTAGATGTCGGCCGGGATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATggtggt                                                                                  >  1:814692/3‑150 (MQ=255)
                                                         cACTCGGCGTAGATGTCGGCCGGGATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTc                                                                                 >  2:1149638/1‑150 (MQ=255)
                                                                      tGTCGGCCGGGATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCa                                                                    <  1:960039/150‑1 (MQ=255)
                                                                           gCCGGGATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTgc                                                                             >  2:342421/1‑136 (MQ=255)
                                                                           gCCGGGATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTgc                                                                             <  1:342421/136‑1 (MQ=255)
                                                                             cGGGATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCgg                                                             <  1:138965/150‑1 (MQ=255)
                                                                              gggATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCggg                                                            <  2:455854/150‑1 (MQ=255)
                                                                                 aTTTCGAACGGGCCGTGGTTCCAGTTCATGTCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCGGGCGa                                                         <  1:487518/150‑1 (MQ=255)
                                                                                   ttCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCGGGCGATc                                                       <  1:1149638/150‑1 (MQ=255)
                                                                                   ttCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCGGGCGATc                                                       >  1:355205/1‑150 (MQ=255)
                                                                                            gCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCGGGCGATCGCTCTTgcg                                              >  1:234256/1‑150 (MQ=255)
                                                                                             ccGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCGGGCGATCGCTCTTgcgg                                             <  2:355205/150‑1 (MQ=255)
                                                                                                  gTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCGGGCGATCGCTCTTgcgg                                             >  1:45908/1‑145 (MQ=255)
                                                                                                  gTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCGGGCGATCGCTCTTgcgg                                             <  2:45908/145‑1 (MQ=255)
                                                                                                   ttCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCGGGCGATCGCTCTTGCGGGCGGtc                                       >  2:516692/1‑150 (MQ=255)
                                                                                                   ttCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCGGGCGATCGCTCTTGCGGGCGGtc                                       >  1:950670/1‑150 (MQ=255)
                                                                                                      cAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTc                                                                                                         <  2:987534/81‑1 (MQ=255)
                                                                                                      cAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTc                                                                                                         >  1:987534/1‑81 (MQ=255)
                                                                                                       aGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCGGGCGATCGCTCTTGCGGGCGGTCTCGa                                   >  2:390162/1‑150 (MQ=255)
                                                                                                       aGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCGGGCGATCGCTCTTGCGGGCGGTCTCGa                                   <  2:950670/150‑1 (MQ=255)
                                                                                                        gTTCAGGGCCTGGGGGGTAAGGGAGGTTTAGTCGTTGCCGAGGGGGGCAGCGTGGAAGTCTGCATTGGCGTGATTGTTGGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGAAGGGCGAAGGGTCTTGAGGGCGGTGTCGAt                                  >  2:6839/1‑150 (MQ=255)
                                                                                                         ttCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCGAAATTTTAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCGGGCGATCGCTCTTGCGGGCGGTCTCGATg                                 <  2:338134/150‑1 (MQ=255)
                                                                                                               gCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCGGGCGATCGCTCTTGCGGGCGGTCTCGATGGCCATc                           <  2:567228/150‑1 (MQ=255)
                                                                                                                    gCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGc                                                                                                               >  1:235684/1‑61 (MQ=255)
                                                                                                                    gCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGc                                                                                                               <  2:235684/61‑1 (MQ=255)
                                                                                                                      gggTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCGGGCGATCGCTCTTGCGGGCGGTCTCGATGGCCATCTTGAttt                    >  2:318887/1‑150 (MQ=255)
                                                                                                                         tCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGAtgg                                                                                      <  1:1202181/81‑1 (MQ=255)
                                                                                                                         tCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGAtgg                                                                                      >  2:1202181/1‑81 (MQ=255)
                                                                                                                         tCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCGGGCGATCGCTCTTGCGGGCGGTCTCGATg                                 <  2:718722/134‑1 (MQ=255)
                                                                                                                         tCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCGGGCGATCGCTCTTGCGGGCGGTCTCGATg                                 >  1:718722/1‑134 (MQ=255)
                                                                                                                             ggCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCGGGc                                                           >  1:975115/1‑104 (MQ=255)
                                                                                                                             ggCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCGGGc                                                           <  2:975115/104‑1 (MQ=255)
                                                                                                                                     ccacggtGCCCCGCTGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCGGGCGATCGCTCTTGCGGGCGGTCTCGATGGCCATCTTGATTTCTTCGGCATCGTGGc     >  1:1057285/7‑150 (MQ=255)
                                                                                                                                     ccacggtGCCCCGCTGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCGGGCGATCGCTCTTGCGGGCGGTCTCGATGGCCATCTTGATTTCTTCGGCATCGTGGc     <  2:1057285/144‑1 (MQ=255)
                                                                                                                                      gtcgtTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCGGGAGATCGCTCTTGCGGGCGGTCTCGATGGCCATCTTGa                       <  2:229071/131‑1 (MQ=255)
                                                                                                                                      gtcgtTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCGGGAGATCGCTCTTGCGGGCGGTCTCGATGGCCATCTTGa                       >  1:229071/1‑131 (MQ=255)
                                                                                                                                        cgtTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCGGGCGATCGCTCTTGCGGGCGGTCTCGATGGCCATCTTGATTTCTTCGGCATCGTGGCCGt  <  2:814692/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
GCGTCGAAGCGCTGGTTCCACTCGGCTTCAACCTTGGCACCGGCGGCCTTGGCATCCCACTCGGCGTAGATGTCGGCCGGGATTTCGAACGGGCCGTGGTTCCAGTTCAGGGCCTGGCGGGTCAGGGCGATTTCGTCGTTGCCCAGCGGGGCACCGTGGCAGTCTTCCTTGCCCTGCTTGTTCGGCGAACCGAAACCGATGGTGGTCTTGCAGCAGATCAGGGTCGGGCGATCGCTCTTGCGGGCGGTCTCGATGGCCATCTTGATTTCTTCGGCATCGTGGCCGT  >  NC_002947/5656811‑5657096

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