Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,781,768 T→C 34.4% F73F (TTT→TTC yfjD → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,781,7680TC34.4% 61.1 / 49.2 87F73F (TTT→TTCyfjDpredicted 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

ATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGAAAAAGTCGCT  >  minE/1781706‑1781834
                                                              |                                                                  
aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGTCTTTTGTCGt                                                               >  1:3410559/1‑68 (MQ=255)
aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTc                                                                 >  1:3086236/1‑66 (MQ=255)
aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCg                                                                >  1:2747033/1‑67 (MQ=255)
aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGt                                                               >  1:2055274/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:3237664/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:1417593/69‑1 (MQ=255)
aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                                                              >  1:1386945/1‑69 (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:1259206/69‑1 (MQ=255)
                tGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAgg                                              <  1:3331871/69‑1 (MQ=255)
                tGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAgg                                              <  1:382211/69‑1 (MQ=255)
                tGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAgg                                              <  1:2016217/69‑1 (MQ=255)
                tGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAgg                                              <  1:535957/69‑1 (MQ=255)
                tGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAgg                                              <  1:1740080/69‑1 (MQ=255)
                tGTACGGCGATGCGGGCGTGGCAATTGAGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAgg                                              <  1:2555262/69‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:1002544/66‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:2940175/66‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:980345/66‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:1462760/66‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:3272812/66‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:3303441/66‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:926919/66‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:3373867/66‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:1051553/66‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:48408/66‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:1171893/66‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:1816740/66‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:708917/66‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:705083/66‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:2912985/66‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:2077059/66‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:1727498/66‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:170288/66‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:2324799/66‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:1569957/66‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:2729237/66‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:1548315/66‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:2759978/66‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:2796098/66‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:2796797/66‑1 (MQ=255)
                   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                            <  1:1543873/66‑1 (MQ=255)
                    cGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAGGTAt                                           <  1:3166493/68‑1 (MQ=255)
                    cGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAGGTAt                                           <  1:2908020/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:3147088/1‑69 (MQ=255)
                                   ggCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTATTGCCGAAAACCATTGc                           >  1:1980711/1‑69 (MQ=255)
                                   ggCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTATTGCCGAAAACCATTGc                           >  1:140396/1‑69 (MQ=255)
                                   ggCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTATTGCCGAAAACCATTGc                           >  1:1728477/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:918640/1‑69 (MQ=255)
                                                      cTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGaaaaa        >  1:2975779/1‑69 (MQ=255)
                                                      cTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGaaaaa        >  1:3386557/1‑69 (MQ=255)
                                                      cTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGaaaaa        >  1:1465614/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)
                                                            tttGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGAAAAAGTCGCt  >  1:3085011/1‑69 (MQ=255)
                                                              |                                                                  
ATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGAAAAAGTCGCT  >  minE/1781706‑1781834

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: