Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,279,641 A→G 28.3% A175A (GCA→GCG trmA → tRNA (uracil‑5‑)‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,279,6410AG28.3% 97.7 / 40.9 95A175A (GCA→GCGtrmAtRNA (uracil‑5‑)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (32/36);  new base G (17/10);  total (49/46)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

GGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACGCGGCGATGAATA  >  minE/2279576‑2279705
                                                                 |                                                                
ggTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAggg                                                               <  1:3119116/69‑1 (MQ=255)
ggTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAggg                                                               <  1:479626/69‑1 (MQ=255)
ggTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAggg                                                               <  1:1196361/69‑1 (MQ=255)
  tCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGa                                                              <  1:2767015/68‑1 (MQ=255)
  tCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGa                                                              <  1:1711634/68‑1 (MQ=255)
  tCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGa                                                              <  1:3263605/68‑1 (MQ=255)
  tCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGa                                                              <  1:2025096/68‑1 (MQ=255)
  tCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGa                                                              <  1:1164268/68‑1 (MQ=255)
       cAACGAAAACCAAAATCGAG‑‑GGATCAGGATTACATCGATGAACGTCTGCCGGCCGCAGGGAAAGAgatg                                                      <  1:2707685/69‑1 (MQ=255)
         tcGAAAACCAAAATCGAGCTGGATCAGTATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatg                                                      <  1:2276098/68‑1 (MQ=255)
         aCGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatg                                                      <  1:2349738/69‑1 (MQ=255)
         aCGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatg                                                      <  1:740521/69‑1 (MQ=255)
         aCGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatg                                                      <  1:2843835/69‑1 (MQ=255)
         aCGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatg                                                      <  1:3339138/69‑1 (MQ=255)
         aCGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatg                                                      <  1:2296320/69‑1 (MQ=255)
         aCGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatg                                                      <  1:1587833/69‑1 (MQ=255)
         aCGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatg                                                      <  1:3245187/69‑1 (MQ=255)
           gAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatga                                                     >  1:2507901/1‑68 (MQ=255)
           gAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatga                                                     >  1:801299/1‑68 (MQ=255)
           gAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatga                                                     >  1:2103285/1‑68 (MQ=255)
           gAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatga                                                     >  1:2998198/1‑68 (MQ=255)
           gAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatga                                                     >  1:2274507/1‑68 (MQ=255)
           gAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatga                                                     >  1:3157646/1‑68 (MQ=255)
           gAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatga                                                     >  1:1213381/1‑68 (MQ=255)
           gAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatga                                                     >  1:1445521/1‑68 (MQ=255)
                       cGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatg                                                      <  1:1750704/55‑1 (MQ=255)
                         aGCTGGATCAGGATTACATCGATGAACGGCTGCCGGTCGCAGGGa                                                              <  1:3295365/45‑1 (MQ=255)
                                 cAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGc                              >  1:3122466/1‑69 (MQ=255)
                                 cAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGc                              >  1:1150674/1‑69 (MQ=255)
                                 cAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGc                              >  1:2223441/1‑69 (MQ=255)
                                 cAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGc                              >  1:1166928/1‑69 (MQ=255)
                                 cAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGc                              >  1:717285/1‑69 (MQ=255)
                                 cAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGc                              >  1:392006/1‑69 (MQ=255)
                                 cAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGc                              >  1:284245/1‑69 (MQ=255)
                                 cAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGc                              >  1:1332376/1‑69 (MQ=255)
                                 cAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGc                              >  1:3051259/1‑69 (MQ=255)
                                 cAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGc                              >  1:3072781/1‑69 (MQ=255)
                                 cAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGc                              >  1:1809221/1‑69 (MQ=255)
                                 cAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGc                              >  1:1690140/1‑69 (MQ=255)
                                    gATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGaa                                    <  1:547130/60‑1 (MQ=255)
                                        aCATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTAccc                       >  1:3315001/1‑69 (MQ=255)
                                        aCATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTAccc                       >  1:3046368/1‑69 (MQ=255)
                                        aCATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTAccc                       >  1:826839/1‑69 (MQ=255)
                                        aCATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTAccc                       >  1:110680/1‑69 (MQ=255)
                                        aCATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTAccc                       >  1:2173223/1‑69 (MQ=255)
                                        aCATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTAccc                       >  1:119000/1‑69 (MQ=255)
                                        aCATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTAccc                       >  1:2104411/1‑69 (MQ=255)
                                        aCATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTAccc                       >  1:1810688/1‑69 (MQ=255)
                                        aCATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTAcc                        >  1:1104047/1‑68 (MQ=255)
                                             gATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCg                  <  1:2110621/69‑1 (MQ=255)
                                               cgAACGTCTGCTGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAACCg                  <  1:1343934/66‑1 (MQ=255)
                                               cgAACGTCTGCCGATCGCGGGGAAAGATATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAACCg                  <  1:759102/66‑1 (MQ=255)
                                               cgAACGTCTGCCGATCGCGGGGAAAGATATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAACCg                  <  1:3283211/66‑1 (MQ=255)
                                               cgAACGTCTGCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAACCg                  <  1:2173323/66‑1 (MQ=255)
                                               cgAACGTCTGCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAACCg                  <  1:1209287/66‑1 (MQ=255)
                                               cgAACGTCTGCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAACCg                  <  1:50498/66‑1 (MQ=255)
                                               cgAACGTCTGCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAACCg                  <  1:3410884/66‑1 (MQ=255)
                                               cgAACGTCTGCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAACCg                  <  1:1658936/66‑1 (MQ=255)
                                               cgAACGTCTGCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAACCg                  <  1:2460577/66‑1 (MQ=255)
                                               cgAACGTCTGCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAACCg                  <  1:292582/66‑1 (MQ=255)
                                                     tCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACgcggc          <  1:2066573/69‑1 (MQ=255)
                                                     tCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACgcggc          <  1:2680264/69‑1 (MQ=255)
                                                     tCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACgcggc          <  1:2131855/69‑1 (MQ=255)
                                                     tCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACgcggc          <  1:2725067/69‑1 (MQ=255)
                                                     tCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACgcggc          <  1:1585579/69‑1 (MQ=255)
                                                     tCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACgcggc          <  1:3274886/69‑1 (MQ=255)
                                                     tCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACgcggc          <  1:3160311/69‑1 (MQ=255)
                                                     tCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACgcggc          <  1:1744560/69‑1 (MQ=255)
                                                     tCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACgcggc          <  1:1951421/69‑1 (MQ=255)
                                                      gtGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACgcggc          <  1:1252140/67‑1 (MQ=255)
                                                      cTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAATCAGCTTTACCCAGCCGAACgcggc          <  1:3410202/68‑1 (MQ=255)
                                                      cTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACgcggc          <  1:728304/68‑1 (MQ=255)
                                                      cTGCCGGTCGCAGGGAAAGAGATGATCTACAGTCAGGTAGAAAACAGCTTTACCCAGCCGAACgcggc          <  1:2077560/68‑1 (MQ=255)
                                                        gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCttt                           >  1:2216504/1‑49 (MQ=255)
                                                        gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCttt                           >  1:1678920/1‑49 (MQ=255)
                                                        gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCttt                           >  1:246952/1‑49 (MQ=255)
                                                        gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCttt                           >  1:2663768/1‑49 (MQ=255)
                                                        gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCttt                           >  1:481833/1‑49 (MQ=255)
                                                        gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCttt                           >  1:417863/1‑49 (MQ=255)
                                                        gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCttt                           >  1:35576/1‑49 (MQ=255)
                                                        gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCttt                           >  1:355534/1‑49 (MQ=255)
                                                        gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCttt                           >  1:347400/1‑49 (MQ=255)
                                                        gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCttt                           >  1:3095114/1‑49 (MQ=255)
                                                        gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCttt                           >  1:150881/1‑49 (MQ=255)
                                                        gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCttt                           >  1:3109791/1‑49 (MQ=255)
                                                        gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCttt                           >  1:3294880/1‑49 (MQ=255)
                                                        gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCttt                           >  1:1609505/1‑49 (MQ=255)
                                                        gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCtt                            >  1:1274302/1‑48 (MQ=255)
                                                        gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCtt                            >  1:2660784/1‑48 (MQ=255)
                                                        gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCtt                            >  1:3042887/1‑48 (MQ=255)
                                                         ccGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCTAACGCGGCGATg      >  1:562083/1‑69 (MQ=255)
                                                           ggTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACGCGGCGATGaa    <  1:3277204/69‑1 (MQ=255)
                                                            gTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACGCGGCGATGaa    <  1:1207143/68‑1 (MQ=255)
                                                              cGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACGCGGCGATGAata  >  1:1998805/1‑68 (MQ=255)
                                                              cGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACGCGGCGATGAata  >  1:2108742/1‑68 (MQ=255)
                                                                 |                                                                
GGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACGCGGCGATGAATA  >  minE/2279576‑2279705

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