Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,781,765 T→C 34.4% T72T (ACT→ACC yfjD → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,781,7650TC34.4% 61.1 / 49.2 87T72T (ACT→ACCyfjDpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (36/21);  new base C (4/26);  total (40/47)
Fisher's exact test for biased strand distribution p-value = 1.01e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGAAAAAGTCG  >  minE/1781699‑1781832
                                                                  |                                                                   
cGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACttt                                                                   >  1:2964842/1‑69 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGTCTTTTGTCGt                                                             >  1:3410559/1‑68 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTc                                                               >  1:3086236/1‑66 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCg                                                              >  1:2747033/1‑67 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGt                                                             >  1:3237664/1‑68 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGt                                                             >  1:2098523/1‑68 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGt                                                             >  1:2371700/1‑68 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGt                                                             >  1:2374221/1‑68 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGt                                                             >  1:239761/1‑68 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGt                                                             >  1:1560504/1‑68 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGt                                                             >  1:2786024/1‑68 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGt                                                             >  1:2055274/1‑68 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGt                                                             >  1:1976946/1‑68 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGt                                                             >  1:653992/1‑68 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGt                                                             >  1:730417/1‑68 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGt                                                             >  1:773059/1‑68 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                                                            >  1:1883298/1‑69 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                                                            >  1:1921077/1‑69 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                                                            <  1:1895713/69‑1 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                                                            <  1:1051496/69‑1 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                                                            >  1:2134023/1‑69 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                                                            <  1:2319367/69‑1 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                                                            >  1:3359722/1‑69 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                                                            >  1:42187/1‑69 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                                                            <  1:489925/69‑1 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                                                            <  1:493486/69‑1 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                                                            >  1:887564/1‑69 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                                                            >  1:1386945/1‑69 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                                                            <  1:1417593/69‑1 (MQ=255)
       aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                                                            >  1:1347952/1‑69 (MQ=255)
        attttgGGATGCGTTTGTACGGCGATGCGGGCGTGGCAAATGCGACTGGTGTGCTGACTTTTGTCGTAc                                                           >  1:1753910/4‑69 (MQ=255)
                     tttGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGa                                              >  1:185610/1‑69 (MQ=255)
                     tttGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGa                                              >  1:3132694/1‑69 (MQ=255)
                     tttGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGa                                              >  1:1592228/1‑69 (MQ=255)
                     tttGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGa                                              >  1:2217939/1‑69 (MQ=255)
                     tttGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGa                                              >  1:2201517/1‑69 (MQ=255)
                       tGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAgg                                            <  1:382211/69‑1 (MQ=255)
                       tGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAgg                                            <  1:3331871/69‑1 (MQ=255)
                       tGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAgg                                            <  1:1259206/69‑1 (MQ=255)
                       tGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAgg                                            <  1:535957/69‑1 (MQ=255)
                       tGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAgg                                            <  1:2016217/69‑1 (MQ=255)
                       tGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAgg                                            <  1:1740080/69‑1 (MQ=255)
                       tGTACGGCGATGCGGGCGTGGCAATTGAGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAgg                                            <  1:2555262/69‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:1171893/66‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:1816740/66‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:48408/66‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:3373867/66‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:980345/66‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:926919/66‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:3303441/66‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:3272812/66‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:1462760/66‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:2077059/66‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:1002544/66‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:1543873/66‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:2940175/66‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:2912985/66‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:2796797/66‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:2796098/66‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:1548315/66‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:2759978/66‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:1051553/66‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:1569957/66‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:2729237/66‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:705083/66‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:1727498/66‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:708917/66‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:170288/66‑1 (MQ=255)
                          atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                          <  1:2324799/66‑1 (MQ=255)
                           cGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAGGTAt                                         <  1:2908020/68‑1 (MQ=255)
                           cGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAGGTAt                                         <  1:3166493/68‑1 (MQ=255)
                                     ggCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAAc                               <  1:458224/68‑1 (MQ=255)
                                     ggCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAAc                               <  1:1895410/68‑1 (MQ=255)
                                     ggCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAAc                               <  1:2189055/68‑1 (MQ=255)
                                          ggCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAGGTAt                                         <  1:70601/53‑1 (MQ=255)
                                          ggCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTATTGCCGAAAACCATTGc                         >  1:140396/1‑69 (MQ=255)
                                          ggCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTATTGCCGAAAACCATTGc                         >  1:1980711/1‑69 (MQ=255)
                                          ggCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTATTGCCGAAAACCATTGc                         >  1:1728477/1‑69 (MQ=255)
                                          ggCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTATTGCCGAAAACCATTGc                         >  1:3147088/1‑69 (MQ=255)
                                                 gCGACTGGTGTGCTGACTTTTGTCTTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTg                  >  1:1721356/1‑69 (MQ=255)
                                                 gCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTg                  >  1:1620782/1‑69 (MQ=255)
                                                      tGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAAc                               <  1:3401566/51‑1 (MQ=255)
                                                           tgCTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAAc                               <  1:790615/46‑1 (MQ=255)
                                                             cTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTa                >  1:177341/1‑59 (MQ=255)
                                                             cTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGaaaaa      >  1:1465614/1‑69 (MQ=255)
                                                             cTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGaaaaa      >  1:2975779/1‑69 (MQ=255)
                                                             cTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGaaaaa      >  1:3386557/1‑69 (MQ=255)
                                                             cTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGaaaaa      >  1:918640/1‑69 (MQ=255)
                                                                  ttttGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGAAAAAGTCg  <  1:1197884/68‑1 (MQ=255)
                                                                  ttttGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGAAAAAGTCg  <  1:2703448/68‑1 (MQ=255)
                                                                  ttttGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGAAAAAGTCg  <  1:370971/68‑1 (MQ=255)
                                                                  |                                                                   
CGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGAAAAAGTCG  >  minE/1781699‑1781832

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: