Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,963,872 C→G 49.1% V50V (GTG→GTC yggD ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,963,8720CG49.1% ‑4.6 / 131.5 104V50V (GTG→GTCyggDpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (28/25);  new base G (30/21);  total (58/46)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

Only 100 of 126 total aligned reads displayed.

CAGTGCAACATCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGTATCCCCTGGAGCAAAATGG  >  minE/1963809‑1963939
                                                               |                                                                   
cAGTGCAACATCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGAt                                                                >  1:1018968/1‑69 (MQ=255)
cAGTGCAACATCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGAt                                                                >  1:3181262/1‑69 (MQ=255)
cAGTGCAACATCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGAt                                                                >  1:1165025/1‑69 (MQ=255)
      aaCATCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTc                                                          <  1:2671052/69‑1 (MQ=255)
      aaCATCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTc                                                          <  1:1054551/69‑1 (MQ=255)
      aaCATCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTc                                                          <  1:486047/69‑1 (MQ=255)
      aaCATCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTc                                                          <  1:609038/69‑1 (MQ=255)
      aaCATCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTc                                                          <  1:738752/69‑1 (MQ=255)
      aaCATCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTc                                                          <  1:2577241/69‑1 (MQ=255)
      aaCATCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTc                                                          <  1:1228146/69‑1 (MQ=255)
      aaCATCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTc                                                          <  1:740531/69‑1 (MQ=255)
       aCATCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTc                                                          <  1:2403085/68‑1 (MQ=255)
       aCATCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTc                                                          <  1:1505912/68‑1 (MQ=255)
       aCATCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCa                                                         >  1:2172412/1‑69 (MQ=255)
          tCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATctct                                                            >  1:250090/1‑63 (MQ=255)
          tCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATc                                                       >  1:1432509/1‑68 (MQ=255)
          tCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATc                                                       >  1:315318/1‑68 (MQ=255)
          tCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATc                                                       >  1:852729/1‑68 (MQ=255)
          tCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATc                                                       >  1:1869197/1‑68 (MQ=255)
               atCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGtt                                                    >  1:380746/1‑66 (MQ=255)
               atCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTAt                                                  >  1:2294240/1‑68 (MQ=255)
               atCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATc                                                 >  1:1984226/1‑69 (MQ=255)
               atCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATc                                                 >  1:1925897/1‑69 (MQ=255)
               atCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATc                                                 >  1:1825878/1‑69 (MQ=255)
               atCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATc                                                 >  1:2003473/1‑69 (MQ=255)
               atCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATc                                                 >  1:1925204/1‑69 (MQ=255)
               atCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATc                                                 >  1:3391765/1‑69 (MQ=255)
               atCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATc                                                 >  1:1594091/1‑69 (MQ=255)
               atCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATc                                                 >  1:1529214/1‑69 (MQ=255)
               atCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATc                                                 >  1:2110837/1‑69 (MQ=255)
               atCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATc                                                 >  1:1463708/1‑69 (MQ=255)
               atCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATc                                                 >  1:2544214/1‑69 (MQ=255)
               atCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATc                                                 >  1:894171/1‑69 (MQ=255)
               atCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTACATCGTTATc                                                 >  1:1099605/1‑69 (MQ=255)
               aacGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATc                                                 >  1:628356/3‑69 (MQ=255)
                tCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATc                                                 >  1:2059762/1‑68 (MQ=255)
                tCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATc                                                 >  1:1009262/1‑68 (MQ=255)
                tCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATc                                                 >  1:765565/1‑68 (MQ=255)
                tCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATc                                                 >  1:1167972/1‑68 (MQ=255)
                tCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATc                                                 >  1:1204747/1‑68 (MQ=255)
                tCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATc                                                 >  1:149616/1‑68 (MQ=255)
                    cAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATCGACaa                                            <  1:1336983/69‑1 (MQ=255)
                    cAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATCGACaa                                            <  1:2460153/69‑1 (MQ=255)
                    cAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATCGACaa                                            <  1:1027184/69‑1 (MQ=255)
                    cAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATCGACaa                                            <  1:2318549/69‑1 (MQ=255)
                    cAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATCGACaa                                            <  1:2909345/69‑1 (MQ=255)
                    cAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATCGACaa                                            <  1:429789/69‑1 (MQ=255)
                    cAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATCGACaa                                            <  1:420700/69‑1 (MQ=255)
                      aGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATCGACaa                                            <  1:3357047/67‑1 (MQ=255)
                          ggCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCg                                      <  1:2305015/69‑1 (MQ=255)
                          ggCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCg                                      <  1:3055037/69‑1 (MQ=255)
                           gCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCg                                      <  1:2757521/68‑1 (MQ=255)
                             cGCTTTGGGCGAGTAACGGCTTCACTGCATATTcc                                                                     >  1:3125079/1‑35 (MQ=255)
                             cGCTTTGGGCGAGTAACGGCTTAACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGAc                                    >  1:924933/1‑68 (MQ=255)
                                tttGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATCGACAAACAGCGACgg                                  >  1:285054/1‑67 (MQ=255)
                                tttGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATCGACAAACAGCGACGGc                                 >  1:2414433/1‑68 (MQ=255)
                                tttGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATCGACAAACAGCGACGGc                                 >  1:16365/1‑68 (MQ=255)
                                tttGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATCGACAAACAGCGACGGc                                 >  1:254648/1‑68 (MQ=255)
                                tttGGGCGAGTAACGGCTTCACTGCATATTCGACGATCTCTTCATCGTTATCGACAAACAGCGACGGc                                 >  1:2581882/1‑68 (MQ=255)
                                    ggCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGa                              <  1:657277/67‑1 (MQ=255)
                                        aGTAACGGCTTCCCTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTgg                        <  1:282261/69‑1 (MQ=255)
                                        aGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTgg                        <  1:1642687/69‑1 (MQ=255)
                                        aGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTgg                        <  1:1480773/69‑1 (MQ=255)
                                        aGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTgg                        <  1:176875/69‑1 (MQ=255)
                                        aGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTgg                        <  1:178894/69‑1 (MQ=255)
                                         gTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTgg                        <  1:2756311/68‑1 (MQ=255)
                                           aaCGGCTTCTCTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGt                     <  1:488084/69‑1 (MQ=255)
                                           aaCGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGt                     <  1:135289/69‑1 (MQ=255)
                                           aaCGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGt                     <  1:2966662/69‑1 (MQ=255)
                                           aaCGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGt                     <  1:1137552/69‑1 (MQ=255)
                                           aaCGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGt                     <  1:2110091/69‑1 (MQ=255)
                                           aaCGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGt                     <  1:2582124/69‑1 (MQ=255)
                                           aaCGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGt                     <  1:424020/69‑1 (MQ=255)
                                            aCGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGTa                    >  1:3083822/1‑69 (MQ=255)
                                            aCGGCTTCACCGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGt                     <  1:510283/68‑1 (MQ=255)
                                                  tCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGt                     <  1:1164488/62‑1 (MQ=255)
                                                        cATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGTATCCCCTGGAGc         >  1:2357770/1‑68 (MQ=255)
                                                           atTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGt                     <  1:1534333/53‑1 (MQ=255)
                                                            ttCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGTATCCCCTGGAGCAAAAt    >  1:879314/1‑69 (MQ=255)
                                                            ttCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGTATCCCCTGGAGCAAAAt    >  1:537536/1‑69 (MQ=255)
                                                            ttCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGTATCCCCTGGAGCAAAAt    >  1:1578701/1‑69 (MQ=255)
                                                            ttCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGAGTATCCCCTGGAGCAAAAt    >  1:2939898/1‑69 (MQ=255)
                                                              ccACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGa                              >  1:1007188/1‑41 (MQ=255)
                                                              ccACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGa                              >  1:2470025/1‑41 (MQ=255)
                                                              ccACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGa                              >  1:1514930/1‑41 (MQ=255)
                                                              ccACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGa                              >  1:1398769/1‑41 (MQ=255)
                                                              ccACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGa                              >  1:1138283/1‑41 (MQ=255)
                                                              ccACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGa                              >  1:153898/1‑41 (MQ=255)
                                                              ccACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGa                              >  1:3262163/1‑41 (MQ=255)
                                                              ccACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGa                              >  1:1695429/1‑41 (MQ=255)
                                                              ccACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGa                              >  1:1884985/1‑41 (MQ=255)
                                                              ccACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGa                              >  1:2063359/1‑41 (MQ=255)
                                                              ccACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGa                              >  1:2718731/1‑41 (MQ=255)
                                                              ccACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGTATCCCCTGGAGCAAAATgg  <  1:407122/69‑1 (MQ=255)
                                                              ccACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGTATCCCCTGGAGCAAAATgg  <  1:579164/69‑1 (MQ=255)
                                                              ccACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGTATCCCCTGGAGCAAAATgg  <  1:1668426/69‑1 (MQ=255)
                                                              ccACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGTATCCCCTGGAGCAAAATgg  <  1:1774353/69‑1 (MQ=255)
                                                              ccACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGTATCCCCTGGAGCAAAATgg  <  1:2832704/69‑1 (MQ=255)
                                                              ccACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGTATCCCCTGGAGCAAAATgg  <  1:2744410/69‑1 (MQ=255)
                                                              ccACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGTATCCCCTGGAGCAAAATgg  <  1:2269564/69‑1 (MQ=255)
                                                               |                                                                   
CAGTGCAACATCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGTATCCCCTGGAGCAAAATGG  >  minE/1963809‑1963939

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: