Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 229,430 C→T 72.5% G147G (GGC→GGT)  proB → gamma‑glutamate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE229,4300CT72.5% 58.7 / 22.1 58G147G (GGC→GGT) proBgamma‑glutamate kinase
Reads supporting (aligned to +/- strand):  ref base C (11/5);  new base T (15/27);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 3.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.04e-01

TCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTTGCTG  >  minE/229363‑229490
                                                                   |                                                            
tCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGc                                                              <  1:2574743/68‑1 (MQ=255)
tCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGc                                                              <  1:2216513/68‑1 (MQ=255)
tCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGc                                                              <  1:456772/68‑1 (MQ=255)
tCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGc                                                              <  1:555387/68‑1 (MQ=255)
tCGATAACAATATCGTTCCCGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGc                                                              <  1:2667654/68‑1 (MQ=255)
 cGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcg                                                             >  1:2163933/1‑68 (MQ=255)
      aCAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgat                                                           <  1:3232953/65‑1 (MQ=255)
      aCAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgat                                                           <  1:2454528/65‑1 (MQ=255)
      aCAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgat                                                           <  1:2504043/65‑1 (MQ=255)
      aCAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgat                                                           <  1:2150652/65‑1 (MQ=255)
      aCAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgat                                                           <  1:2051700/65‑1 (MQ=255)
      aCAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgat                                                           <  1:162178/65‑1 (MQ=255)
      aCAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgat                                                           <  1:1501230/65‑1 (MQ=255)
      aCAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgat                                                           <  1:985456/65‑1 (MQ=255)
      aCAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgat                                                           <  1:759707/65‑1 (MQ=255)
          taacGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgataacgata                                                    <  1:2633523/65‑1 (MQ=255)
           atCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgataactata                                                    <  1:1102225/67‑1 (MQ=255)
           atCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgataacgata                                                    <  1:48655/67‑1 (MQ=255)
           atCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgataacgata                                                    <  1:855063/67‑1 (MQ=255)
           atCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgataacgata                                                    <  1:1076240/67‑1 (MQ=255)
           atCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgataacgata                                                    <  1:2064955/67‑1 (MQ=255)
           atCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgataacgata                                                    <  1:1291564/67‑1 (MQ=255)
           atCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgataacgata                                                    <  1:2049991/67‑1 (MQ=255)
           atCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgataacgata                                                    <  1:2759486/67‑1 (MQ=255)
           atCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgataacgata                                                    <  1:2903636/67‑1 (MQ=255)
           atCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgataacgata                                                    <  1:3121183/67‑1 (MQ=255)
           atCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgataacgata                                                    <  1:3242908/67‑1 (MQ=255)
              gTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGCGATAACGATAACCtt                                               >  1:346607/1‑69 (MQ=255)
                    gTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgataacgata                                                    <  1:2036773/58‑1 (MQ=255)
                      aaTCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGCGATAACGATAACCTTTCTGCGCt                                       >  1:2520103/1‑69 (MQ=255)
                      aaTCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGCGATAACGATAACCTTTCTGCGCt                                       >  1:1061109/1‑69 (MQ=255)
                      aaTCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGCGATAACGATAACCTTTCTGCGCt                                       >  1:1629721/1‑69 (MQ=255)
                        ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGtgataa                                                         >  1:3299447/3‑49 (MQ=255)
                        ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                                     >  1:2786599/3‑69 (MQ=255)
                        ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                                     >  1:2908842/3‑69 (MQ=255)
                        ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                                     >  1:583756/3‑69 (MQ=255)
                        ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                                     >  1:2529399/3‑69 (MQ=255)
                        ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                                     >  1:570630/3‑69 (MQ=255)
                        ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                                     >  1:2566539/3‑69 (MQ=255)
                        ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                                     >  1:1012135/3‑69 (MQ=255)
                        ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                                     >  1:235567/3‑69 (MQ=255)
                        ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                                     >  1:1934235/3‑69 (MQ=255)
                        ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                                     >  1:1761417/3‑69 (MQ=255)
                        ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                                     >  1:1732580/3‑69 (MQ=255)
                        ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                                     >  1:1974659/3‑69 (MQ=255)
                        ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                                     >  1:1276936/3‑69 (MQ=255)
                        ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                                     >  1:2141648/3‑69 (MQ=255)
                        ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTg                                      >  1:940989/3‑68 (MQ=255)
                                   gATGCTGTCGCTTCTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:307558/69‑2 (MQ=255)
                                   gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCTATTCtg                          <  1:36493/69‑2 (MQ=255)
                                   gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:1351477/69‑2 (MQ=255)
                                   gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:1489476/69‑2 (MQ=255)
                                   gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:1116829/69‑2 (MQ=255)
                                   gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:766102/69‑2 (MQ=255)
                                   gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:1076215/69‑2 (MQ=255)
                                   gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:1783840/69‑2 (MQ=255)
                                   gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:3258907/69‑2 (MQ=255)
                                   gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:934070/69‑2 (MQ=255)
                                   gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:1814150/69‑2 (MQ=255)
                                   gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:3192344/69‑2 (MQ=255)
                                   gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:245115/69‑2 (MQ=255)
                                   gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:3160617/69‑2 (MQ=255)
                                   gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:2733478/69‑2 (MQ=255)
                                   gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:2427776/69‑2 (MQ=255)
                                   gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:225412/69‑2 (MQ=255)
                                   gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:2615371/69‑2 (MQ=255)
                                   gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:1834553/69‑2 (MQ=255)
                                   gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:2776216/69‑2 (MQ=255)
                                   gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:1969324/69‑2 (MQ=255)
                                   gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:281555/69‑2 (MQ=255)
                                   gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:2926836/69‑2 (MQ=255)
                                   gATGCTGTCGCTACGGCAGAGATTAAGGTCGGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtt                          >  1:2507037/1‑69 (MQ=255)
                                   gATGCTGTCGCTACGGCAGAGATTAAGGTCGGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtt                          >  1:14850/1‑69 (MQ=255)
                                   gATGCTGTCGCTACGGCAGAGATTAAGGTCGGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtt                          >  1:3355163/1‑69 (MQ=255)
                                    aTGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:1893231/68‑2 (MQ=255)
                                           cGCTACTGCAGAGATTAAGGTCGGTTATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:645252/61‑2 (MQ=255)
                                               aCGGCAGAGATTAAGGTCGGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGa               >  1:22259/1‑68 (MQ=255)
                                                  gCAGAGATTAAGGTCGGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATa             <  1:3335526/67‑1 (MQ=255)
                                                  gCAGAGATTAAGGTCGGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATa             <  1:285695/67‑1 (MQ=255)
                                                  gCAGAGATTAAGGTCGGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATa             <  1:229234/67‑1 (MQ=255)
                                                  gCAGAGATTAAGGTCGGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATa             <  1:2670386/67‑1 (MQ=255)
                                                  gCAGAGATTAAGGTCGGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATa             <  1:688764/67‑1 (MQ=255)
                                                           aaGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:598454/45‑2 (MQ=255)
                                                           aaGGTCGGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTtgctg  >  1:1809991/1‑69 (MQ=255)
                                                           aaGGTCGGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTtgctg  >  1:101109/1‑69 (MQ=255)
                                                           aaGGTCGGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTtgct   >  1:3022902/1‑68 (MQ=255)
                                                                cGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg                          <  1:2862837/40‑2 (MQ=39)
                                                                   |                                                            
TCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTTGCTG  >  minE/229363‑229490

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