Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 109,697 C→T 57.7% H647H (CAC→CAT aceE → pyruvate dehydrogenase, decarboxylase component E1, thiamin‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE109,6970CT57.7% 21.0 / 59.1 78H647H (CAC→CATaceEpyruvate dehydrogenase, decarboxylase component E1, thiamin‑binding
Reads supporting (aligned to +/- strand):  ref base C (16/17);  new base T (21/24);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCAT  >  minE/109630‑109765
                                                                   |                                                                    
tCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCcaca                                                                     >  1:2141555/1‑69 (MQ=255)
tCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCcaca                                                                     >  1:440283/1‑69 (MQ=255)
               ggcGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGGCACAGCCATATTCAGTCGCTGACTa                                                      <  1:2333988/67‑1 (MQ=255)
               gTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTa                                                      <  1:1782281/69‑1 (MQ=255)
               gTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTa                                                      <  1:2683361/69‑1 (MQ=255)
               gTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTa                                                      <  1:2632199/69‑1 (MQ=255)
               gTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTa                                                      <  1:2623350/69‑1 (MQ=255)
               gTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTa                                                      <  1:2591336/69‑1 (MQ=255)
               gTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTa                                                      <  1:272160/69‑1 (MQ=255)
               gTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTa                                                      <  1:2754872/69‑1 (MQ=255)
               gTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTa                                                      <  1:3042182/69‑1 (MQ=255)
               gTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTa                                                      <  1:1927027/69‑1 (MQ=255)
               gTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTa                                                      <  1:191834/69‑1 (MQ=255)
               gTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTa                                                      <  1:1736459/69‑1 (MQ=255)
               gTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTa                                                      <  1:1548228/69‑1 (MQ=255)
               gTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTa                                                      <  1:1495719/69‑1 (MQ=255)
               gTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTa                                                      <  1:1465699/69‑1 (MQ=255)
               gTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTa                                                      <  1:1155897/69‑1 (MQ=255)
               gTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTa                                                      <  1:800740/69‑1 (MQ=255)
               gTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTa                                                      <  1:1159366/69‑1 (MQ=255)
               gTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTa                                                      <  1:745936/69‑1 (MQ=255)
               gTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTa                                                      <  1:660666/69‑1 (MQ=255)
               gTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGACATATTCAGTCGCTGACTa                                                      <  1:3204994/69‑1 (MQ=255)
                tCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTa                                                      <  1:2104530/68‑1 (MQ=255)
                tCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTa                                                      <  1:1893507/68‑1 (MQ=255)
                          cTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGt                                           >  1:404299/1‑69 (MQ=255)
                            gAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTATCCCGAACTGTa                                          >  1:281871/1‑68 (MQ=255)
                            gAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTATCCCGAACTGTAt                                         >  1:594458/1‑69 (MQ=255)
                            gAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTATCCCGAACTGTAt                                         >  1:1723985/1‑69 (MQ=255)
                                   gAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACt                                                       <  1:1903161/48‑1 (MQ=255)
                                    aaGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACt                                                       <  1:2082999/47‑1 (MQ=255)
                                         cTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGAc                                                        <  1:1247469/41‑1 (MQ=255)
                                         cTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGAc                                                        <  1:827065/41‑1 (MQ=255)
                                         cTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGAc                                                        <  1:821407/41‑1 (MQ=255)
                                         cTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGAc                                                        <  1:3035090/41‑1 (MQ=255)
                                         cTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGAc                                                        <  1:2285185/41‑1 (MQ=255)
                                         cTGCAGCACGAAGATGGTCACAGCCACATTCAGGCGCTGAc                                                        <  1:3430695/41‑1 (MQ=255)
                                          tGCAGCACGAAGATGTTCACAGCCATATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCgg                           >  1:3167310/1‑69 (MQ=255)
                                          tGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCgg                           >  1:995699/1‑69 (MQ=255)
                                          tGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCgg                           >  1:374611/1‑69 (MQ=255)
                                          tGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCgg                           >  1:2805782/1‑69 (MQ=255)
                                          tGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCgg                           >  1:3423400/1‑69 (MQ=255)
                                          tGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCgg                           >  1:496544/1‑69 (MQ=255)
                                          tGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCgg                           >  1:3322214/1‑69 (MQ=255)
                                          tGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCgg                           >  1:2551500/1‑69 (MQ=255)
                                          tGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCgg                           >  1:3226532/1‑69 (MQ=255)
                                          tGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCgg                           >  1:3318216/1‑69 (MQ=255)
                                          tGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCgg                           >  1:560533/1‑69 (MQ=255)
                                          tGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCgg                           >  1:102906/1‑69 (MQ=255)
                                          tGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCgg                           >  1:1387919/1‑69 (MQ=255)
                                          tGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCgg                           >  1:1614922/1‑69 (MQ=255)
                                          tGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCgg                           >  1:1987368/1‑69 (MQ=255)
                                          tGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCgg                           >  1:1230070/1‑69 (MQ=255)
                                          tGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCgg                           >  1:1156079/1‑69 (MQ=255)
                                          tGCAGCACGAAGATGGTCACAGCCATATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCg                            >  1:1381188/1‑68 (MQ=255)
                                              gcaCGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgct                         <  1:2532788/67‑1 (MQ=255)
                                              gcaCGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgct                         <  1:2928049/67‑1 (MQ=255)
                                              gcaCGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgct                         <  1:578417/67‑1 (MQ=255)
                                              gcaCGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgct                         <  1:2119209/67‑1 (MQ=255)
                                              gcaCGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgct                         <  1:1152543/67‑1 (MQ=255)
                                              gcaCGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgct                         <  1:3249697/67‑1 (MQ=255)
                                              gcaCGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgct                         <  1:1727382/67‑1 (MQ=255)
                                              gcaCGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgct                         <  1:1235490/67‑1 (MQ=255)
                                              gcaCGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgct                         <  1:3376792/67‑1 (MQ=255)
                                                 cGAAGATGGTCACAGCCATATTCAGTCGCTGACTa                                                      <  1:19741/35‑1 (MQ=255)
                                                         gTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGCCCCGGCTTACGCTTACGAAg            >  1:2776212/1‑69 (MQ=255)
                                                         gTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAg            >  1:784751/1‑69 (MQ=255)
                                                         gTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAg            >  1:113721/1‑69 (MQ=255)
                                                         gTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAg            >  1:1983368/1‑69 (MQ=255)
                                                         gTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAg            >  1:2248785/1‑69 (MQ=255)
                                                         gTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAg            >  1:2272016/1‑69 (MQ=255)
                                                               gCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCt       >  1:634930/1‑68 (MQ=255)
                                                               gCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCt       >  1:2739205/1‑68 (MQ=255)
                                                               gCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCt       >  1:433437/1‑68 (MQ=255)
                                                               gCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCt       >  1:1523803/1‑68 (MQ=255)
                                                               gCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCt       >  1:1582751/1‑68 (MQ=255)
                                                               gCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCt       >  1:332486/1‑68 (MQ=255)
                                                               gCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCt       >  1:1755/1‑68 (MQ=255)
                                                               gCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCt       >  1:2852582/1‑68 (MQ=255)
                                                               gCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCt       >  1:2393484/1‑68 (MQ=255)
                                                                 cacaTTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGtc    >  1:561868/1‑69 (MQ=255)
                                                                 cacaTTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGtc    >  1:1232214/1‑69 (MQ=255)
                                                                 cacaTTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGtc    >  1:2950396/1‑69 (MQ=255)
                                                                 cacaTTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGtc    >  1:2541471/1‑69 (MQ=255)
                                                                 cacaTTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGtc    >  1:1130720/1‑69 (MQ=255)
                                                                   caTTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGtcat  <  1:1449815/69‑1 (MQ=255)
                                                                   caTTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGtcat  <  1:2753155/69‑1 (MQ=255)
                                                                   |                                                                    
TCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCAT  >  minE/109630‑109765

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: