Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,441,545 C→G 40.7% P251P (CCC→CCG ilvY → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,441,5450CG40.7% 43.3 / 78.5 81P251P (CCC→CCGilvYDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (18/30);  new base G (6/27);  total (24/57)
Fisher's exact test for biased strand distribution p-value = 8.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GGTATCGATGGTGGCACTCGGCTGTGGCGTGGCATTGTTGCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAGCGATGAGAAAACGCCG  >  minE/2441479‑2441599
                                                                  |                                                      
ggTATCGATGGTGGCACTCGGCTGTGGCGTGGCATTGTTGCCGGAAGTGGTGCTGGAAAACAGCCCCGa                                                      <  1:447929/69‑1 (MQ=255)
 gTATCGATGGTGGCACTCGGCTGTGGCGTGGCATTGTTGCCGGAAGTGGTGCTGGAAAACAGCCCCGa                                                      <  1:1662394/68‑1 (MQ=255)
   aTCGATGGTGGCACTCGGCTGTGGCGTGGCATTGTTGCCGGAAGTGGTGCTGGAAAACAGCCCCGAAcc                                                   <  1:2724311/69‑1 (MQ=255)
   aTCGATGGTGGCACTCGGCTGTGGCGTGGCATTGTTGCCGGAAGTGGTGCTGGAAAACAGCCCCGAAcc                                                   <  1:708/69‑1 (MQ=255)
          gtggCACTCGGCTGTGGCGTGGCATTGTTGCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGt                                            <  1:1026634/69‑1 (MQ=255)
          gtggCACTCGGCTGTGGCGTGGCATTGTTGCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGt                                            <  1:2729279/69‑1 (MQ=255)
          gtggCACTCGGCTGTGGCGTGGCATTGTTGCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGt                                            <  1:2589971/69‑1 (MQ=255)
          gtggCACTCGGCTGTGGCGTGGCATTGTTGCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGt                                            <  1:2773875/69‑1 (MQ=255)
          gtggCACTCGGCTGTGGCGTGGCATTGTTGCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGt                                            <  1:2843084/69‑1 (MQ=255)
          gtggCACTCGGCTGTGGCGTGGCATTGTTGCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGt                                            <  1:2409656/69‑1 (MQ=255)
          gtggCACTCGGCTGTGGCGTGGCATTGTTGCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGt                                            <  1:2888802/69‑1 (MQ=255)
          gtggCACTCGGCTGTGGCGTGGCATTGTTGCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGt                                            <  1:2336081/69‑1 (MQ=255)
          gtggCACTCGGCTGTGGCGTGGCATTGTTGCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGt                                            <  1:3087962/69‑1 (MQ=255)
          gtggCACTCGGCTGTGGCGTGGCATTGTTGCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGt                                            <  1:1655432/69‑1 (MQ=255)
          gtggCACTCGGCTGTGGCGTGGCATTGTTGCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGt                                            <  1:1568117/69‑1 (MQ=255)
          gtggCACTCGGCTGTGGCGTGGCATTGTTGCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGt                                            <  1:1401346/69‑1 (MQ=255)
          gtggCACTCGGCTGTGGCGTGGCATTGTTGCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGt                                            <  1:1182029/69‑1 (MQ=255)
           tggCACTCGGCTGTGGCGTGGCATTGTTGCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGt                                            <  1:709655/68‑1 (MQ=255)
            ggCACTCGGCTGTGGCGTGGCATTGTTGCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTaa                                          >  1:3192732/1‑69 (MQ=255)
            ggCACTCGGCTGTGGCGTGGCATTGTTGCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTa                                           <  1:2858154/68‑1 (MQ=255)
                cTCGGCTGTGGCGTGGCATTGTTGCCGGAAGTGGTGCTGGAAAACAGCCCCGa                                                      <  1:1989006/53‑1 (MQ=255)
                  cGGCTGTGGCGTGGCATTGTTGCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGt                                    >  1:1701887/1‑69 (MQ=255)
                                  ttgttgCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGTGATGATTTTAGAgcgc                    >  1:1801307/1‑69 (MQ=255)
                                  ttgttgCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGTGATGATTTTAGAgcgc                    >  1:3184378/1‑69 (MQ=255)
                                  ttgttgCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGTGATGATTTTAGAgcgc                    >  1:2737910/1‑69 (MQ=255)
                                  ttgttgCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGTGATGATTTTAGAgcgc                    >  1:80768/1‑69 (MQ=255)
                                  ttgttgCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGTGATGATTTTAGAgcgc                    >  1:2703777/1‑69 (MQ=255)
                                  ttgttgCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGTGATGATTTTAGAgcg                     >  1:2211844/1‑68 (MQ=255)
                                   tgttgCCGGAAGTGGTGCTGTAAAACAGCCCCGAAcc                                                   <  1:145610/37‑1 (MQ=255)
                                    tttgCCGGAAGCGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAg                  <  1:3104395/68‑1 (MQ=255)
                                    gttgCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCa                   >  1:2861258/1‑68 (MQ=255)
                                    gttgCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCa                   >  1:1319843/1‑68 (MQ=255)
                                    gttgCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCa                   >  1:2849882/1‑68 (MQ=255)
                                    gttgCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCa                   >  1:2253960/1‑68 (MQ=255)
                                    gttgCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCa                   >  1:2161613/1‑68 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGCGCGCAg                  <  1:439388/67‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAg                  <  1:79437/67‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAg                  <  1:873578/67‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAg                  <  1:95335/67‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAg                  <  1:971724/67‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAg                  <  1:2981827/67‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAg                  <  1:3136864/67‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAg                  <  1:3175526/67‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAg                  <  1:3230430/67‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAg                  <  1:568680/67‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAg                  <  1:3350064/67‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAg                  <  1:3404830/67‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAg                  <  1:2512092/67‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAg                  <  1:1131230/67‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAg                  <  1:1494706/67‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAg                  <  1:1530992/67‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAg                  <  1:1560304/67‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAg                  <  1:1039122/67‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAg                  <  1:166405/67‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAg                  <  1:2280802/67‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAg                  <  1:2315472/67‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAg                  <  1:1030911/67‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAg                  <  1:2444301/67‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAg                  <  1:2358525/67‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTGACCGCGTGATGATTTTAGAGCGCAGc                 <  1:2503969/68‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAGc                 <  1:1954722/68‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAGc                 <  1:1559192/68‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAGc                 <  1:1604588/68‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAGc                 <  1:1657360/68‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAGc                 <  1:1637823/68‑1 (MQ=255)
                                      tgCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAAACGCGTGATGATTTTAGAGCGCAGc                 <  1:2345216/68‑1 (MQ=255)
                                       gCCGGAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAg                  <  1:1916081/66‑1 (MQ=255)
                                       gCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAGCg                >  1:2878622/1‑68 (MQ=255)
                                       gCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAGCg                >  1:1842918/1‑68 (MQ=255)
                                       gCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAGCg                >  1:941343/1‑68 (MQ=255)
                                       gCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAGCg                >  1:2966672/1‑68 (MQ=255)
                                       gCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAGCg                >  1:1301268/1‑68 (MQ=255)
                                          ggAAGTGGTGCTGGAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGAAg                  <  1:1115680/63‑1 (MQ=255)
                                             aGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAGCGATGAGaa         <  1:567387/69‑1 (MQ=255)
                                             aGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAGCGATGAGaa         <  1:752510/69‑1 (MQ=255)
                                                       gAAAACAGCCCGGACCCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAGtgatga            >  1:599073/1‑56 (MQ=255)
                                                       gAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAGTGATGAAAAAACGCCg  >  1:1004873/1‑66 (MQ=255)
                                                       gAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAGTGATGAAAAAACGCCg  >  1:3308596/1‑66 (MQ=255)
                                                       gAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAGTGATGAAAAAACGCCg  >  1:2238468/1‑66 (MQ=255)
                                                       gAAAACAGCCCGGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAGTGATGAAAAAACGCCg  >  1:2468619/1‑66 (MQ=255)
                                                       gAAAACAGCCCGGAACCGGTGCGTAACCGCGGGATGATTTTAGAGCGCAGTGATGAAAAAACGCCg  >  1:2529895/1‑66 (MQ=255)
                                                                  |                                                      
GGTATCGATGGTGGCACTCGGCTGTGGCGTGGCATTGTTGCCGGAAGTGGTGCTGGAAAACAGCCCCGAACCGGTGCGTAACCGCGTGATGATTTTAGAGCGCAGCGATGAGAAAACGCCG  >  minE/2441479‑2441599

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: