Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,761,357 T→A 7.5% L118Q (CTG→CAG)  gpmI → 2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,761,3570TA 7.5% 165.8 / 4.7 53L118Q (CTG→CAG) gpmI2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base T (27/22);  new base A (2/2);  total (29/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 = 7.79e-01

ATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAGTCGTCCATCGAAC  >  NC_002947/5761209‑5761506
                                                                                                                                                    |                                                                                                                                                     
atGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGActg                                                                                                                                                      >  2:982518/1‑150 (MQ=255)
      cTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTcc                                                                                                                                                >  2:972757/1‑150 (MQ=255)
       tCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCg                                                                                                                                               >  2:239185/1‑150 (MQ=255)
         ggCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAt                                                                                                                                             <  1:123197/150‑1 (MQ=255)
            gcAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtggt                                                                                                                                          >  1:879198/1‑150 (MQ=255)
               ggTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATggtggt                                                                                                                                       <  1:239185/150‑1 (MQ=255)
                     gtcgtGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTacac                                                                                                                                 >  1:752523/1‑150 (MQ=255)
                         tgtATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGcc                                                                                                                             >  1:63793/1‑150 (MQ=255)
                         tgtATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACAGCAGTCCGag                                                                                                                                             <  1:624974/134‑2 (MQ=255)
                         tgtATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACAGCAGTCCGag                                                                                                                                             >  2:624974/1‑133 (MQ=255)
                             tCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtggtg                                                                                                                                         >  2:369758/1‑134 (MQ=255)
                             tCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtggtg                                                                                                                                         <  1:369758/134‑1 (MQ=255)
                                ggACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGa                                                                                                                      <  2:1324420/150‑1 (MQ=255)
                                   cTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtggtg                                                                                                                                         <  1:954492/128‑1 (MQ=255)
                                   cTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtggtg                                                                                                                                         >  2:954492/1‑128 (MQ=255)
                                    ttCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCAc                                                                                                                  <  1:1043882/150‑1 (MQ=255)
                                           gggTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTCCTGTCCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTtggtgg                                                                                                           <  1:87310/150‑1 (MQ=255)
                                           gggTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTtggtgg                                                                                                           >  1:355983/1‑150 (MQ=255)
                                              tGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCa                                                                                                        >  1:698321/1‑150 (MQ=255)
                                                  cAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACAGCAGTCCGa                                                                                                                                              >  2:86277/1‑108 (MQ=255)
                                                  cAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACAGCAGTCCGa                                                                                                                                              <  1:86277/108‑1 (MQ=255)
                                                   aaGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGcc                                                                                                   >  1:1381710/1‑150 (MQ=255)
                                                             gcgACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGgccg                                                                                         >  1:171642/1‑150 (MQ=255)
                                                                   gCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGc                                                                                   >  1:503873/1‑150 (MQ=255)
                                                                     gAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGt                                                                                 >  1:344495/1‑150 (MQ=255)
                                                                      aGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTg                                                                                <  2:344495/150‑1 (MQ=255)
                                                                        ttcttcGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATgg                                                                                                     <  2:1227852/127‑1 (MQ=255)
                                                                        ttcttcGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATgg                                                                                                     >  1:1227852/1‑127 (MQ=255)
                                                                        ttcttcGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGgc                                                                              <  2:171642/150‑1 (MQ=255)
                                                                           ttcGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGa                                                                                                                                              >  1:1115141/1‑83 (MQ=255)
                                                                           ttcGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGa                                                                                                                                              <  2:1115141/83‑1 (MQ=255)
                                                                           ttcGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGa                                                                                                                      <  2:1289189/107‑1 (MQ=255)
                                                                           ttcGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGa                                                                                                                      >  1:1289189/1‑107 (MQ=255)
                                                                                    ccGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATc                                                                  >  2:1339907/1‑150 (MQ=255)
                                                                                       gTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTAc                                                               <  1:697853/150‑1 (MQ=255)
                                                                                       gTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGCACTGCTGTCCGATGGTGTTGTACTCAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTTGCCGCCCATCGTGGCGCTGAAAAGATCTAc                                                               <  2:78073/150‑1 (MQ=255)
                                                                                           tGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGc                                                           >  2:1181201/1‑150 (MQ=255)
                                                                                             aCCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCACTGAAAAGATCTACCTGCAc                                                         >  1:880420/1‑150 (MQ=255)
                                                                                                 gTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCt                                                     >  1:643191/1‑150 (MQ=255)
                                                                                                     tgtgGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCt                                                 >  2:440543/1‑150 (MQ=255)
                                                                                                      gtgGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAg                                                                                                                        <  1:485932/78‑1 (MQ=255)
                                                                                                      gtgGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAg                                                                                                                        >  2:485932/1‑78 (MQ=255)
                                                                                                      gtgGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTc                                                <  1:440543/150‑1 (MQ=255)
                                                                                                       tgGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCAGAAAAGATCTACCTGCACGCCTTCCTCg                                               >  1:1071975/1‑150 (MQ=255)
                                                                                                            aaGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGc                                          >  2:980000/1‑150 (MQ=255)
                                                                                                              ggCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGgccgc                                                                                        <  1:945414/102‑1 (MQ=255)
                                                                                                              ggCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGgccgc                                                                                        >  2:945414/1‑102 (MQ=255)
                                                                                                                  gTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGac                                    >  1:1129160/1‑150 (MQ=255)
                                                                                                                   tCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGaca                                   <  1:1339907/150‑1 (MQ=255)
                                                                                                                    cGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGacac                                  <  1:982518/150‑1 (MQ=255)
                                                                                                                                      gCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCa                <  1:288492/150‑1 (MQ=255)
                                                                                                                                          aTCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAgtcg            >  1:845046/1‑150 (MQ=255)
                                                                                                                                             cTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAGTCGTcc         <  2:698321/150‑1 (MQ=255)
                                                                                                                                                ggACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAGTCGTCCATc      <  2:1234820/150‑1 (MQ=255)
                                                                                                                                                    tgctgTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAGTCGTCCATCGAAc  >  2:1339016/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAGTCGTCCATCGAAC  >  NC_002947/5761209‑5761506

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