Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,967,352 G→A 42.9% G426G (GGC→GGT tktA ← transketolase 1, thiamin‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,967,3520GA42.9% 56.9 / 107.2 143G426G (GGC→GGTtktAtransketolase 1, thiamin‑binding
Reads supporting (aligned to +/- strand):  ref base G (59/22);  new base A (16/45);  total (76/67)
Fisher's exact test for biased strand distribution p-value = 3.64e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 143 total aligned reads displayed.

CCATACGTACGGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGA  >  minE/1967288‑1967416
                                                                |                                                                
ccATACGTACGGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCAc                                                              <  1:1532550/69‑1 (MQ=255)
ccATACGTACGGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCAc                                                              <  1:1501774/69‑1 (MQ=255)
ccATACGTACGGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCAc                                                              <  1:225568/69‑1 (MQ=255)
  aTACGTACGGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCg                                                            <  1:3375249/69‑1 (MQ=255)
     cgtacgGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGc                                                         >  1:1463850/1‑69 (MQ=255)
     cgtacgGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGc                                                         >  1:1879283/1‑69 (MQ=255)
     cgtacgGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGc                                                         >  1:456071/1‑69 (MQ=255)
     cgtacgGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGc                                                         >  1:481113/1‑69 (MQ=255)
     cgtacgGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGc                                                         >  1:2013352/1‑69 (MQ=255)
     cgtacgGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGc                                                         >  1:1485918/1‑69 (MQ=255)
     cgtacgGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGc                                                         >  1:2834233/1‑69 (MQ=255)
     cgtacgGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGc                                                         >  1:1232464/1‑69 (MQ=255)
     cgtacgGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGc                                                         >  1:774761/1‑69 (MQ=255)
     cgtacgGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGc                                                         >  1:2393314/1‑69 (MQ=255)
     cgtacgGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGc                                                         >  1:2632517/1‑69 (MQ=255)
     cgtacgGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGc                                                         >  1:2714963/1‑69 (MQ=255)
     cgtacgGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGc                                                         >  1:1044099/1‑69 (MQ=255)
     cgtacgGCGTTACGTGCGTATTCCACGAACATAAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGc                                                         >  1:1310410/1‑69 (MQ=255)
         cgGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGc                                                         >  1:72575/1‑65 (MQ=255)
         cgGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCa                                                        >  1:2809697/1‑66 (MQ=255)
         cgGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCa                                                        >  1:3351895/1‑66 (MQ=255)
         cgGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCa                                                        >  1:2391771/1‑66 (MQ=255)
         cgGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCa                                                        >  1:2955332/1‑66 (MQ=255)
         cgGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCa                                                        >  1:194863/1‑66 (MQ=255)
         cgGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCa                                                        >  1:3287133/1‑66 (MQ=255)
         cgGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCa                                                        >  1:1009842/1‑66 (MQ=255)
         cgGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCa                                                        >  1:549293/1‑66 (MQ=255)
         cgGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCa                                                        >  1:1432079/1‑66 (MQ=255)
         cgGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCa                                                        >  1:684896/1‑66 (MQ=255)
         cgGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCa                                                        >  1:624917/1‑66 (MQ=255)
          ggCATTACGTGCGTATTCCACGAACATCAGGAAGGTTGAGGTATACGGCAGGAACCCACCGTGCa                                                        >  1:2062800/1‑65 (MQ=25)
          ggCATTACGTGCGTATTCCACGAACATCAGGAAGGTTGAGGTATACGGCAGGAAACCACCGTGc                                                         >  1:2418803/1‑64 (MQ=21)
          ggCATTACGTGCGTATTCCACGAACATCAGGAAGGTTGAGGTATACGGCAGGAAACCACCGTGCa                                                        >  1:1754035/1‑65 (MQ=25)
          ggCATTACGTGCGTATTCCACGAACATCAGGAAGGTTGAGGTATACGGCAGGAAACCACCGTGCa                                                        >  1:825122/1‑65 (MQ=25)
          ggCATTACGTGCGTATTCCACGAACATCAGGAAGGTTGAGGTATACGGCAGGAAACCACCGTGCa                                                        >  1:1160558/1‑65 (MQ=25)
          ggCATTACGTGCGTATTCCACGAACATCAGGAAGGTTGAGGTATACGGCAGGAAACCACCGTGCa                                                        >  1:1382680/1‑65 (MQ=25)
          ggCATTACGTGCGTATTCCACGAACATCAGGAAGGTTGAGGTATACGGCAGGAAACCACCGTGCa                                                        >  1:2231823/1‑65 (MQ=25)
          ggCATTACGTGCGTATTCCACGAACATCAGGAAGGTTGAGGTATACGGCAGGAAACCACCGTGCa                                                        >  1:2142328/1‑65 (MQ=25)
          ggCATTACGTGCGTATTCCACGAACATCAGGAAGGTTGAGGTATACGGCAGGAAACCACCGTGCa                                                        >  1:3024757/1‑65 (MQ=25)
          ggCATTACGTGCGTATTCCACGAACATCAGGAAGGTTGAGGTATACGGCAGGAAACCACCGTGCa                                                        >  1:3038353/1‑65 (MQ=25)
          ggCATTACGTGCGTATTCCACGAACATCAGGAAGGTTGAGGTATACGGCAGGAAACCACCGTGCa                                                        >  1:3117644/1‑65 (MQ=25)
          ggCATTACGTGCGTATTCCACGAACATCAGGAAGGTTGAGGTATACGGCAGGAAACCACCGTGCa                                                        >  1:3199734/1‑65 (MQ=25)
          ggCATTACGTGCGTATTCCACGAACATCAGGAAGGTTGAGGTATACGGCAGGAAACCACCGTGCa                                                        >  1:478380/1‑65 (MQ=25)
          ggCATTACGTGCGTATTCCACGAACATCAGGAAGGTTGAGGTATACGGCAGGAAACCACCGTGCAggg                                                     >  1:2326529/1‑68 (MQ=25)
           gCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGga                                                    <  1:2258935/68‑1 (MQ=255)
           gCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGga                                                    <  1:942403/68‑1 (MQ=255)
            cGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGga                                                    <  1:1012078/67‑1 (MQ=255)
                     cGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTa                                          >  1:2991911/1‑68 (MQ=255)
                     cGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTa                                          >  1:2942826/1‑68 (MQ=255)
                     cGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTa                                          >  1:730589/1‑68 (MQ=255)
                     cGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTa                                          >  1:2703656/1‑68 (MQ=255)
                        aTTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCaa                                      <  1:509513/69‑1 (MQ=255)
                         ttCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCaa                                      <  1:600615/68‑1 (MQ=255)
                         ttCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAAt                                     >  1:770545/1‑69 (MQ=255)
                         ttCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAAt                                     >  1:2183218/1‑69 (MQ=255)
                                aaCATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCTGGGAGATACCGTTAGCaa                                      <  1:3328905/61‑1 (MQ=255)
                                        gAAGGTTGAGGTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGaa                      <  1:406813/69‑1 (MQ=255)
                                        gAAGGTTGAGGTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGaa                      <  1:944055/69‑1 (MQ=255)
                                        gAAGGTTGAGGTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGaa                      <  1:871819/69‑1 (MQ=255)
                                        gAAGGTTGAGGTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGaa                      <  1:617727/69‑1 (MQ=255)
                                        gAAGGTTGAGGTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGaa                      <  1:1671513/69‑1 (MQ=255)
                                        gAAGGTTGAGGTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGaa                      <  1:1889498/69‑1 (MQ=255)
                                        gAAGGTTGAGGTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGaa                      <  1:3280343/69‑1 (MQ=255)
                                        gAAGGTTGAGGTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGaa                      <  1:2941636/69‑1 (MQ=255)
                                        gAAGGTTGAGGTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGaa                      <  1:2998896/69‑1 (MQ=255)
                                        gAAGGTTGAGGTATACGGCAGAAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGaa                      <  1:1166760/69‑1 (MQ=255)
                                                 ggTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAac             >  1:505274/1‑69 (MQ=255)
                                                 ggTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAac             >  1:2326937/1‑69 (MQ=255)
                                                 ggTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGAATTcacgta              <  1:745490/68‑3 (MQ=21)
                                                 ggTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGAATTcacgta              <  1:541657/68‑3 (MQ=21)
                                                 ggTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGAATTcacgta              <  1:582793/68‑3 (MQ=21)
                                                 ggTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGAATTcacgta              <  1:2900614/68‑3 (MQ=21)
                                                 ggTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGAATTcacgta              <  1:149563/68‑3 (MQ=21)
                                                 ggTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGAATTcacgta              <  1:1363670/68‑3 (MQ=21)
                                                 ggTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGAATTcacgta              <  1:1303948/68‑3 (MQ=21)
                                                 ggTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGAATTcacgta              <  1:1989937/68‑3 (MQ=21)
                                                 ggTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGAATTcacgta              <  1:112568/68‑3 (MQ=21)
                                                 ggTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGAATTcacgta              <  1:829233/68‑3 (MQ=21)
                                                 ggTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGAATTcacgta              <  1:871977/68‑3 (MQ=21)
                                                 ggTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGAATTcacgta              <  1:1060310/68‑3 (MQ=21)
                                                 ggTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGAATTcacgta              <  1:2346268/68‑3 (MQ=21)
                                                 ggTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGAATTcacgta              <  1:1569656/68‑3 (MQ=21)
                                                 ggTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGAATTcacgta              <  1:2072864/68‑3 (MQ=21)
                                                 ggTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGAATTcacgta              <  1:3017136/68‑3 (MQ=21)
                                                 ggTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGAATTcacgta              <  1:3303047/68‑3 (MQ=21)
                                                 ggTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGAATTcacgta              <  1:1768548/68‑3 (MQ=21)
                                                 ggTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGAATTcacgta              <  1:384504/68‑3 (MQ=21)
                                                 ggTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGAATTcacgta              <  1:391417/68‑3 (MQ=21)
                                                 ggTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGAATTcacgta              <  1:1684888/68‑3 (MQ=21)
                                                 ggTATACGGCAGGAAACCACCGTGCAGGGAGATACCGTTGGCAATCGCAGTCATACCGAATTcacgta              <  1:470258/68‑3 (MQ=21)
                                                          cAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTg    <  1:172187/69‑1 (MQ=255)
                                                          cAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTg    <  1:2166205/69‑1 (MQ=255)
                                                          cAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTg    <  1:1967337/69‑1 (MQ=255)
                                                           aGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTg    <  1:268659/68‑1 (MQ=255)
                                                            ggAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGa  >  1:2282483/1‑69 (MQ=255)
                                                            ggAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGa  >  1:201847/1‑69 (MQ=255)
                                                              aaGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGa  >  1:1548391/1‑67 (MQ=255)
                                                              aaGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGa  >  1:3429054/1‑67 (MQ=255)
                                                              aaGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGa  >  1:2862896/1‑67 (MQ=255)
                                                              aaGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGa  >  1:1562746/1‑67 (MQ=255)
                                                                |                                                                
CCATACGTACGGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGA  >  minE/1967288‑1967416

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