Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,296,484 T→G 22.6% R261S (AGA→AGC)  U712_16965 ← 2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,296,4840TG22.6% 60.4 / 8.9 43R261S (AGA→AGC) U712_169652,3‑bisphosphoglycerate‑independent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base T (14/19);  new base G (8/2);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 6.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACATAGCCGTCAACTGTTTCACTGAAGTGAGTCAGGCAGACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAATCATCAACAACATCCAGGGCG  >  CP006881/3296350‑3296622
                                                                                                                                      |                                                                                                                                          
aCATAGCCGTCAACTGTTTCACTGAAGTGAGTCAGGCAGACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTGTGTGTTTGTGTACGTGTTGGTAATATGGATGGCGgggtgcggtcgg                                                                                                                                         >  3:298517/1‑136 (MQ=255)
       cGTCAACTGTTTGACTGAAGTGAGTCTGGCAGACTAAATACAGGTTCTTTTGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGGCTGAAATtata                                                                                                                                 <  5:28760/139‑1 (MQ=255)
         tCAACTGTTTCACTGAAGTGAGTCAGGCAGACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTGGGTGTTTGTGGTCGTGGTGGGAATATGGATGGCGCGGGCAGGGCTGAtattattt                                                                                                                                >  2:366931/1‑136 (MQ=255)
         tCAACTGTTTCACTGAAGTGAGTCAGGCAGACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCGTTGTTTGTGAACGTGTTTGGAATCTGGATGGCGCGGGCCGGGCTGAAATTATAg                                                                                                                                >  2:3679/1‑138 (MQ=255)
         tCAACTGTTTCACTGAAGTGAGTCAGGCAGACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCGTTGGTTGTGTTCGTGTTGGAAATGTGGATGGCGGGGGCCGGGCTGTAATtattg                                                                                                                                >  3:326046/1‑136 (MQ=255)
           aaCTGTTTCACTGAAGTGAGTCAGGCAGACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTGGTTGTTTGTGGACGTGTTGGAAATGTGGGTTGGGCGGTGCGGGCTGGAATTATAGaa                                                                                                                              >  1:74806/1‑138 (MQ=255)
           aaCTGTTTCACTGAAGTGAGTCAGGCAGACGAAATACAGGTTCTTTGGATAATTATCGCCGCGGTCGAAGTCACGGAAGTCGGTTTTTTTTTAAGTGTTGTATATATGGATTGGGCGGGACGGTCTGTAAATATATaa                                                                                                                              >  7:345338/1‑138 (MQ=255)
                tttCACTGAAGTGAGTCAGGCAGACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGCACTCGTTGTTTGGGAGCGGGTTGGTAAAGAGGTTGGGGCGGGGCGGTGTGATATTATTTATAATca                                                                                                                         >  7:386622/1‑138 (MQ=255)
                tttCACTGAAGTGAGTCAGGCAGACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTGGTTGTTTGTTGTCGGGGGTGGAATCTGGGTGGGGGGGGGCGGGCGGAAATTTTTGAAATTcac                                                                                                                        >  2:70741/1‑139 (MQ=255)
                tttCACTGAAGTGAGTCAGGCAGACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCGTTGTTTGTGTACGTGTTGGAAATCTGGATGGCGCGGGCCGGTCTGAAATTATAGAAAATca                                                                                                                         >  1:22269/1‑138 (MQ=255)
                                cAGGCCGAAGTAATACAGGGTCTGTGTAAAATTCTCGGAGCGGTCGAAGTCAAGGAAGGGGTTGTTTGTGACCGTGTTGGAAATCTGGCTGGCGCGGTCCGGTCTGAAATTATAGAAACTCACAGAATCGCCGTCTTgg                                                                                                        <  2:378532/139‑1 (MQ=255)
                                  ggCATACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCGTTGTTTTTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGAt                                                                                                      >  8:107048/1‑139 (MQ=255)
                                  ggCAGACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTGTTTGTTTGTGAACGTGTTGGCAATCTGGAGGGCGCGGTGCGGGGGGTAATTTTAGAAAATGCCAGGATCGCCGTCTTGGag                                                                                                      >  4:99213/1‑138 (MQ=255)
                                          aaaTACAGGTTCTTTGGATAATTCTCGCCGCGGGCGAAGTCCCGGCAGTCTTTGTTGGGGCACGTGTGGGAAATCTGGATGGCGCGGTCCGGGCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCaa                                                                                              <  5:106453/139‑1 (MQ=255)
                                             tACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCTTTGTTTGTGAACGTGTTGGAAATATGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAg                                                                                           >  5:201440/1‑139 (MQ=255)
                                                  gTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCa                                                                                      >  6:193203/1‑139 (MQ=255)
                                                     cTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCTTTGTTTTTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATATAAAATCACAGAATCGCCGTCGTGTCTTTTAGCAACAGGCGCACCg                                                                                   >  7:188140/1‑139 (MQ=255)
                                                     cTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCGTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGGCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTAGCAACAGGCTCCCCg                                                                                   >  5:151036/1‑139 (MQ=255)
                                                     cTTTGGATAATTCTCGCCGCGGACGAAGTCACGGAAGTGTTAGGTTGTTACCGAGGTGGAAATCTGGATGGCGCGGTCCGGTAAGAAATAATAGAAAATCAACGAAACGCCGTCTTCGATATTCGCAACCGGCTCACCg                                                                                   >  8:267528/1‑139 (MQ=255)
                                                      tttGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCGTTGTTTGTGAACGTGTTGGTAATCTGGATGGGGCGGTCCGGTCTGAAATTATTGAAAATCACAGAATCGCCGGGTTGGATTTTCGCACACGGCGCACCGt                                                                                  >  6:170561/1‑139 (MQ=255)
                                                                    gccgcGGTCGAAGTCACGGAAGTCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTtgtg                                                                       <  5:193203/136‑1 (MQ=255)
                                                                          gTCGAAGTCACGGAAGTGGTTGTTTGTGTACGTGGTGGAAATCGGGATGGCGGGGTGCGGGGTGAAATTATTGGAAATGACAGAATCGCCGGGGTGGAGTTTGGCCACAGGCGCAGCGGTTTGTGTTGTGTTGTCAgtg                                                              >  6:366920/1‑137 (MQ=255)
                                                                                         gTCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCa                                               <  1:237501/139‑1 (MQ=255)
                                                                                          tCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCa                                               <  4:84091/138‑1 (MQ=255)
                                                                                                tttgtGAACGTGTTTGAAATCTGGATGGCGCTGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTATTTTGTGATGACAGTTGGAATCACGAATTCATCGTAGa                                        <  7:267528/139‑1 (MQ=255)
                                                                                                tttgtGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGa                                        <  3:59484/139‑1 (MQ=255)
                                                                                                tttgtGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGa                                        <  1:237253/139‑1 (MQ=255)
                                                                                                   gtgAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCGTTTGTGTTGCCAGCTGGAATCACGAATTGAGCGTAGATAc                                     >  8:331685/1‑139 (MQ=255)
                                                                                                    tgAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCGCCCCGTGTTCTTTTGTGATGCCAGATGGAATCCCGAATTCATCGTTGAGAcc                                    >  6:315561/1‑139 (MQ=255)
                                                                                                      accGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCAt                                  <  8:89277/137‑1 (MQ=255)
                                                                                                      aaCGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCAt                                  <  6:201440/139‑1 (MQ=255)
                                                                                                       aCGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCAtt                                 <  5:366920/139‑1 (MQ=255)
                                                                                                          tgtTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTg                               >  5:53106/1‑138 (MQ=255)
                                                                                                          tgtTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTg                               >  5:174961/1‑138 (MQ=255)
                                                                                                                     tGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCTCCGGTTTCTTTTGTGATGACAGATGGAATCACGATTTGATCGTAGATACCATTTGGATATGaataat                   >  8:20620/1‑139 (MQ=255)
                                                                                                                     agacaggCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAatcat                   <  5:253863/134‑1 (MQ=255)
                                                                                                                       gATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAatcat                   <  2:74806/137‑1 (MQ=255)
                                                                                                                       gATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAatcat                   <  8:345338/137‑1 (MQ=255)
                                                                                                                       gATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAatc                     <  2:22269/135‑1 (MQ=255)
                                                                                                                         tGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGTCAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAATCATcaac               <  3:346536/139‑1 (MQ=255)
                                                                                                                         tGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAATCATcaac               <  1:122976/139‑1 (MQ=255)
                                                                                                                          ggCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAATCATcaaca              <  3:99213/139‑1 (MQ=255)
                                                                                                                           gcgcGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAATCATcaaca              <  7:107048/138‑1 (MQ=255)
                                                                                                                           gcgcGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAATCATcaaca              <  4:326046/138‑1 (MQ=255)
                                                                                                                                 tccggtcTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCAc                                                      >  7:147765/1‑92 (MQ=255)
                                                                                                                                 tccggtcTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCAc                                                      <  8:147765/92‑1 (MQ=255)
                                                                                                                                  ccggtcTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAatcat                   <  7:118335/126‑1 (MQ=255)
                                                                                                                                  ccggtcTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAatcat                   >  8:118335/1‑126 (MQ=255)
                                                                                                                                    ggtcTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAATCATCAACAACATCCAggg    <  7:221477/139‑1 (MQ=255)
                                                                                                                                      tcTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCAt                                  >  3:313892/1‑107 (MQ=255)
                                                                                                                                      tcTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCAt                                  <  4:313892/107‑1 (MQ=255)
                                                                                                                                      tcTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAATCATCAACAACATCCAGGgcg  <  1:366931/139‑1 (MQ=255)
                                                                                                                                      tcTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAATCATCAACAACATCCAGGgcg  >  3:173617/1‑139 (MQ=255)
                                                                                                                                      tcTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAATCATCAACAACATCCAGGgcg  <  5:44748/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                          
ACATAGCCGTCAACTGTTTCACTGAAGTGAGTCAGGCAGACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAATCATCAACAACATCCAGGGCG  >  CP006881/3296350‑3296622

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