Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,279,636 G→A 29.2% V174I (GTC→ATC)  trmA → tRNA (uracil‑5‑)‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,279,6360GA29.2% 94.4 / 39.4 92V174I (GTC→ATC) trmAtRNA (uracil‑5‑)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (31/34);  new base A (17/10);  total (48/44)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

CTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACGCGGCGATGAA  >  minE/2279570‑2279703
                                                                  |                                                                   
cTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGt                                                                    >  1:3164953/1‑68 (MQ=255)
      ggTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAggg                                                             <  1:3119116/69‑1 (MQ=255)
      ggTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAggg                                                             <  1:479626/69‑1 (MQ=255)
      ggTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAggg                                                             <  1:1196361/69‑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:2767015/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:740521/69‑1 (MQ=255)
               aCGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatg                                                    <  1:2843835/69‑1 (MQ=255)
               aCGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatg                                                    <  1:2296320/69‑1 (MQ=255)
               aCGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatg                                                    <  1:3339138/69‑1 (MQ=255)
               aCGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatg                                                    <  1:1587833/69‑1 (MQ=255)
               aCGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatg                                                    <  1:2349738/69‑1 (MQ=255)
               aCGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatg                                                    <  1:3245187/69‑1 (MQ=255)
                 gAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatga                                                   >  1:2274507/1‑68 (MQ=255)
                 gAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatga                                                   >  1:2998198/1‑68 (MQ=255)
                 gAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatga                                                   >  1:2103285/1‑68 (MQ=255)
                 gAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatga                                                   >  1:3157646/1‑68 (MQ=255)
                 gAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatga                                                   >  1:2507901/1‑68 (MQ=255)
                 gAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatga                                                   >  1:801299/1‑68 (MQ=255)
                 gAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatga                                                   >  1:1445521/1‑68 (MQ=255)
                 gAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatga                                                   >  1:1213381/1‑68 (MQ=255)
                             cGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatg                                                    <  1:1750704/55‑1 (MQ=255)
                               aGCTGGATCAGGATTACATCGATGAACGGCTGCCGGTCGCAGGGa                                                            <  1:3295365/45‑1 (MQ=255)
                                       cAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGc                            >  1:1332376/1‑69 (MQ=255)
                                       cAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGc                            >  1:1150674/1‑69 (MQ=255)
                                       cAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGc                            >  1:1166928/1‑69 (MQ=255)
                                       cAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGc                            >  1:717285/1‑69 (MQ=255)
                                       cAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGc                            >  1:2223441/1‑69 (MQ=255)
                                       cAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGc                            >  1:284245/1‑69 (MQ=255)
                                       cAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGc                            >  1:1690140/1‑69 (MQ=255)
                                       cAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGc                            >  1:3051259/1‑69 (MQ=255)
                                       cAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGc                            >  1:392006/1‑69 (MQ=255)
                                       cAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGc                            >  1:3072781/1‑69 (MQ=255)
                                       cAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGc                            >  1:3122466/1‑69 (MQ=255)
                                       cAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGc                            >  1:1809221/1‑69 (MQ=255)
                                          gATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGaa                                  <  1:547130/60‑1 (MQ=255)
                                              aCATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTAccc                     >  1:3046368/1‑69 (MQ=255)
                                              aCATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTAccc                     >  1:3315001/1‑69 (MQ=255)
                                              aCATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTAccc                     >  1:826839/1‑69 (MQ=255)
                                              aCATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTAccc                     >  1:110680/1‑69 (MQ=255)
                                              aCATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTAccc                     >  1:1810688/1‑69 (MQ=255)
                                              aCATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTAccc                     >  1:119000/1‑69 (MQ=255)
                                              aCATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTAccc                     >  1:2104411/1‑69 (MQ=255)
                                              aCATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTAccc                     >  1:2173223/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:1209287/66‑1 (MQ=255)
                                                     cgAACGTCTGCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAACCg                <  1:50498/66‑1 (MQ=255)
                                                     cgAACGTCTGCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAACCg                <  1:2173323/66‑1 (MQ=255)
                                                     cgAACGTCTGCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAACCg                <  1:2460577/66‑1 (MQ=255)
                                                     cgAACGTCTGCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAACCg                <  1:292582/66‑1 (MQ=255)
                                                     cgAACGTCTGCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAACCg                <  1:3410884/66‑1 (MQ=255)
                                                     cgAACGTCTGCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAACCg                <  1:1658936/66‑1 (MQ=255)
                                                           tCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACgcggc        <  1:1585579/69‑1 (MQ=255)
                                                           tCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACgcggc        <  1:1744560/69‑1 (MQ=255)
                                                           tCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACgcggc        <  1:3160311/69‑1 (MQ=255)
                                                           tCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACgcggc        <  1:3274886/69‑1 (MQ=255)
                                                           tCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACgcggc        <  1:2066573/69‑1 (MQ=255)
                                                           tCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACgcggc        <  1:2680264/69‑1 (MQ=255)
                                                           tCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACgcggc        <  1:2725067/69‑1 (MQ=255)
                                                           tCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACgcggc        <  1:2131855/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:355534/1‑49 (MQ=255)
                                                              gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCttt                         >  1:2663768/1‑49 (MQ=255)
                                                              gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCttt                         >  1:2216504/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:347400/1‑49 (MQ=255)
                                                              gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCttt                         >  1:150881/1‑49 (MQ=255)
                                                              gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCttt                         >  1:1609505/1‑49 (MQ=255)
                                                              gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCttt                         >  1:3294880/1‑49 (MQ=255)
                                                              gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCttt                         >  1:3095114/1‑49 (MQ=255)
                                                              gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCttt                         >  1:246952/1‑49 (MQ=255)
                                                              gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCttt                         >  1:3109791/1‑49 (MQ=255)
                                                              gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCttt                         >  1:1678920/1‑49 (MQ=255)
                                                              gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCtt                          >  1:3042887/1‑48 (MQ=255)
                                                              gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCtt                          >  1:1274302/1‑48 (MQ=255)
                                                              gCCGATCGCGGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCtt                          >  1:2660784/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)
                                                                  |                                                                   
CTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACGCGGCGATGAA  >  minE/2279570‑2279703

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: