Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,287,724 G→A 39.7% A257A (GCG→GCA shiA → shikimate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,287,7240GA39.7% 85.3 / 109.9 161A257A (GCG→GCAshiAshikimate transporter
Reads supporting (aligned to +/- strand):  ref base G (77/20);  new base A (36/28);  total (113/48)
Fisher's exact test for biased strand distribution p-value = 2.59e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 161 total aligned reads displayed.

TCCCGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGCCTTTGCACTTAATTATTCAA  >  minE/1287669‑1287785
                                                       |                                                             
tCCCGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCg                                                  <  1:3008057/69‑1 (MQ=255)
tCCCGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCg                                                  <  1:734557/69‑1 (MQ=255)
tCCCGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCg                                                  <  1:324475/69‑1 (MQ=255)
       tATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGc                                            >  1:2951371/1‑68 (MQ=255)
       tATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCt                                           >  1:906839/1‑69 (MQ=255)
       tATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCt                                           >  1:126785/1‑69 (MQ=255)
       tATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCt                                           >  1:1429858/1‑69 (MQ=255)
       tATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCt                                           >  1:3408213/1‑69 (MQ=255)
       tATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAACATTATTGCGCTACGACTGTGCGAATTGCt                                           >  1:1618955/1‑69 (MQ=255)
        aTCGAAGCTCTGTTACGACATCCCGGTGCTTTTCTGAAGATTATTGCACTACGAC‑‑TGCGAGTTGCTGAc                                        <  1:1323030/69‑1 (MQ=255)
           gAAGCTCTGTTACGACATCCCGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTTAc                                        <  1:1900398/68‑1 (MQ=255)
           gAAGCTCTGTTACGACATCCCGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGAc                                        <  1:2771660/68‑1 (MQ=255)
           gAAGCTCTGTTACGACATCCCGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGAc                                        <  1:2908378/68‑1 (MQ=255)
           gAAGCTCTGTTACGACATCCCGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGAc                                        <  1:1773623/68‑1 (MQ=255)
           gAAGCTCTGTTACGACATCCCGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGAc                                        <  1:1706394/68‑1 (MQ=255)
           gAAGCTCTGTTACGACATCCCGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGAc                                        <  1:1653137/68‑1 (MQ=255)
           gAAGCTCTGTTACGACATCCCGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGAc                                        <  1:2231117/68‑1 (MQ=255)
           gAAGCTCTGTTACGACATCCCGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGAc                                        <  1:2162604/68‑1 (MQ=255)
           gAAGCTCTGTTACGACATCCCGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGAc                                        <  1:2973802/68‑1 (MQ=255)
           gAAGCTCTGTTACGACATCCCGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGAc                                        <  1:1527896/68‑1 (MQ=255)
           gAAGCTCTGTTACGACATCCCGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGAc                                        <  1:3173151/68‑1 (MQ=255)
           gAAGCTCTGTTACGACATCCCGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGAc                                        <  1:1356115/68‑1 (MQ=255)
           gAAGCTCTGTTACGACATCCCGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGAc                                        <  1:3395599/68‑1 (MQ=255)
           gAAGCTCTGTTACGACATCCCGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGAc                                        <  1:391628/68‑1 (MQ=255)
           gAAGCTCTGTTACGACATCCCGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGAc                                        <  1:952098/68‑1 (MQ=255)
           gAAGCTCTGTTACGACATCCCGCTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGAc                                        <  1:2571224/68‑1 (MQ=255)
            aaGCTCTGTTACGACATCCCGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGAc                                        <  1:2445040/67‑1 (MQ=255)
                  tGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTAc                                 <  1:497562/68‑1 (MQ=255)
                      aCGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGt                            >  1:2953656/1‑69 (MQ=255)
                      aCGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGt                            <  1:1048401/69‑1 (MQ=255)
                      aCGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGt                            >  1:2583190/1‑69 (MQ=255)
                      aCGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGt                            >  1:2373768/1‑69 (MQ=255)
                      aCGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGt                            <  1:2200706/69‑1 (MQ=255)
                      aCGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGt                            >  1:975906/1‑69 (MQ=255)
                      aCGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGt                            <  1:1649520/69‑1 (MQ=255)
                      aCGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGt                            >  1:1093778/1‑69 (MQ=255)
                      aCGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGt                            >  1:1552122/1‑69 (MQ=255)
                      aCGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGt                            >  1:1280083/1‑69 (MQ=255)
                      aCGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGt                            >  1:1734952/1‑69 (MQ=255)
                             cccGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGcc                     <  1:1011211/69‑1 (MQ=255)
                             cccGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGcc                     <  1:564238/69‑1 (MQ=255)
                             cccGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGcc                     <  1:3095712/69‑1 (MQ=255)
                              ccGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGa                                                 >  1:3309053/1‑40 (MQ=255)
                              ccGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGt                            >  1:1589053/1‑61 (MQ=255)
                              ccGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGcc                     <  1:921244/68‑1 (MQ=255)
                              ccGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGCCt                    >  1:333139/1‑69 (MQ=255)
                              ccGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGCCt                    >  1:1582824/1‑69 (MQ=255)
                              ccGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGCCt                    >  1:2764276/1‑69 (MQ=255)
                              ccGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGCCt                    >  1:3223799/1‑69 (MQ=255)
                              ccGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGCCt                    >  1:266058/1‑69 (MQ=255)
                              ccGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGCCt                    >  1:2657403/1‑69 (MQ=255)
                              ccGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGCCt                    >  1:837580/1‑69 (MQ=255)
                              ccGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGCCt                    >  1:1947461/1‑69 (MQ=255)
                              ccGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGCCt                    >  1:205185/1‑69 (MQ=255)
                               cGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGACGATGTACAtt                              >  1:2119371/1‑57 (MQ=255)
                               cGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGACGATGTACAtt                              >  1:130898/1‑57 (MQ=255)
                               cGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGACGATGTACAtt                              >  1:509596/1‑57 (MQ=255)
                               cGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGACGATGTACAtt                              >  1:3091864/1‑57 (MQ=255)
                               cGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGACGATGTACAtt                              >  1:719055/1‑57 (MQ=255)
                               cGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGACGATGTACAtt                              >  1:2912251/1‑57 (MQ=255)
                               cGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGACGATGTACAtt                              >  1:2545051/1‑57 (MQ=255)
                               cGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGACGATGTACAtt                              >  1:2028207/1‑57 (MQ=255)
                               cGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGACGATGTACAtt                              >  1:1909335/1‑57 (MQ=255)
                               cGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGACGATGTACAtt                              >  1:2886202/1‑57 (MQ=255)
                               cGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGACGATGTACAtt                              >  1:2489473/1‑57 (MQ=255)
                               cGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGACGATGTACAtt                              >  1:2304281/1‑57 (MQ=255)
                               cGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGACGATGTACAtt                              >  1:2488035/1‑57 (MQ=255)
                               cGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGACGATGTACAtt                              >  1:2815498/1‑57 (MQ=255)
                               cGGTGCTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGACGATGTACAtt                              >  1:2795705/1‑57 (MQ=255)
                                 gTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATg                                    >  1:666939/1‑50 (MQ=255)
                                  tGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGCCTTTGc                <  1:1773871/69‑1 (MQ=255)
                                   gCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTg                       >  1:574781/1‑61 (MQ=255)
                                   gCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTg                       >  1:1194297/1‑61 (MQ=255)
                                   gCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTg                       >  1:861515/1‑61 (MQ=255)
                                   gCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTg                       >  1:212019/1‑61 (MQ=255)
                                   gCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTg                       >  1:2923293/1‑61 (MQ=255)
                                   gCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTg                       >  1:2951651/1‑61 (MQ=255)
                                   gCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTg                       >  1:2717853/1‑61 (MQ=255)
                                   gCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTg                       >  1:3050114/1‑61 (MQ=255)
                                   gCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTg                       >  1:3079129/1‑61 (MQ=255)
                                   gCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTg                       >  1:31071/1‑61 (MQ=255)
                                    cTTTTCTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGACGATGTACATTGTTACTGCCTTTGCa               >  1:1345579/1‑68 (MQ=255)
                                    cTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTg                       >  1:671357/1‑60 (MQ=255)
                                      tttcTGAAGATTATTGCACTACGACTGTGCGAGTTGCTGACGATGTACATTGTTACTGCCTTTGCACtt            >  1:412176/4‑69 (MQ=255)
                                              gATTATTGCACTACGACTGTGCGAGTTGCTGACGATGTACATTGTTACTGCCTTTGCACTTAattat      <  1:577707/67‑1 (MQ=255)
                                              gATTATTGCACTACGACTGTGCGAGTTGCTGACGATGTACATTGTTACTGCCTTTGCACTTAattat      <  1:1709085/67‑1 (MQ=255)
                                              gATTATTGCACTACGACTGTGCGAGTTGCTGACGATGTACATTGTTACTGCCTTTGCACTTAattat      <  1:737554/67‑1 (MQ=255)
                                              gATTATTGCACTACGACTGTGCGAGTTGCTGACGATGTACATTGTTACTGCCTTTGCACTTAattat      <  1:1519027/67‑1 (MQ=255)
                                              gATTATTGCACTACGACTGTGCGAGTTGCTGACGATGTACATTGTTACTGCCTTTGCACTTAattat      <  1:424617/67‑1 (MQ=255)
                                                 tattGCGCTACGACTGTGCGAATTGCTGACGATGTACATCg                             >  1:2835793/1‑41 (MQ=255)
                                                 tattGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGCCTTTGCACTTAATTATTCaa  >  1:1861026/1‑68 (MQ=255)
                                                 tattGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGCCTTTGCACTTAATTATTCaa  >  1:1490488/1‑68 (MQ=255)
                                                 tattGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGCCTTTGCACTTAATTATTCaa  >  1:3115415/1‑68 (MQ=255)
                                                 tattGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGCCTTTGCACTTAATTATTCaa  >  1:1869973/1‑68 (MQ=255)
                                                 tattGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGCCTTTGCACTTAATTATTCaa  >  1:385875/1‑68 (MQ=255)
                                                 tattGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGCCTTTGCACTTAATTATTCaa  >  1:3002134/1‑68 (MQ=255)
                                                 tattGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGCCTTTGCACTTAATTATTCaa  >  1:2890132/1‑68 (MQ=255)
                                                 tattGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGCCTTTGCACTTAATTATTCaa  >  1:132421/1‑68 (MQ=255)
                                                 tattGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGCCTTTGCACTTAATTATTCaa  >  1:1067760/1‑68 (MQ=255)
                                                 tattGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGCCTTTGCACTTAATTATTCaa  >  1:1804552/1‑68 (MQ=255)
                                                       |                                                             
TCCCGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGCCTTTGCACTTAATTATTCAA  >  minE/1287669‑1287785

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: