Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,014,583 G→A 43.2% R54R (CGC→CGT exbB ← membrane spanning protein in TonB‑ExbB‑ExbD complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,014,5830GA43.2% 32.7 / 79.3 104R54R (CGC→CGTexbBmembrane spanning protein in TonB‑ExbB‑ExbD complex
Reads supporting (aligned to +/- strand):  ref base G (20/39);  new base A (13/32);  total (33/71)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

GATATCGTTGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAG  >  minE/2014523‑2014646
                                                            |                                                               
gATATCGTTGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGGTCGCGCTTAAGGCGACGCtt                                                         >  1:481791/1‑69 (MQ=255)
gATATCGTTGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGCtt                                                         >  1:1114565/1‑69 (MQ=255)
gATATCGTTGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGCtt                                                         >  1:2898717/1‑69 (MQ=255)
  tatCGTTGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGCTTCt                                                       <  1:285636/69‑1 (MQ=255)
  tatCGTTGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGCTTCt                                                       <  1:2139359/69‑1 (MQ=255)
                 ttAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAActc                                          <  1:1210924/67‑1 (MQ=255)
                 ttAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAActc                                          <  1:3279886/67‑1 (MQ=255)
                 ttAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAActc                                          <  1:2338506/67‑1 (MQ=255)
                 ttAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAActc                                          <  1:912877/67‑1 (MQ=255)
                 ttAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAActc                                          <  1:3186864/67‑1 (MQ=255)
                 ttAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAActc                                          <  1:476836/67‑1 (MQ=255)
                       aaCGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGctct                                  <  1:2832628/69‑1 (MQ=255)
                       aaCGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGctct                                  <  1:467787/69‑1 (MQ=255)
                       aaCGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGctct                                  <  1:2325246/69‑1 (MQ=255)
                        aCGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCtt                                 >  1:804914/1‑69 (MQ=255)
                                cagccagcAGTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTga                            >  1:3212675/1‑66 (MQ=255)
                                cagccagcAGTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTga                            >  1:410608/1‑66 (MQ=255)
                                cagccagcAGTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTga                            >  1:2820659/1‑66 (MQ=255)
                                cagccagcAGTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTga                            >  1:501398/1‑66 (MQ=255)
                                cagccagcAGTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTga                            >  1:512382/1‑66 (MQ=255)
                                cagccagcAGTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTga                            >  1:2556326/1‑66 (MQ=255)
                                cagccagcAGTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTga                            >  1:567286/1‑66 (MQ=255)
                                cagccagcAGTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTga                            >  1:654356/1‑66 (MQ=255)
                                cagccagcAGTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTga                            >  1:2194435/1‑66 (MQ=255)
                                cagccagcAGTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTga                            >  1:465832/1‑66 (MQ=255)
                                cagccagcAGTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTga                            >  1:2038047/1‑66 (MQ=255)
                                cagccagcAGTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTga                            >  1:1815074/1‑66 (MQ=255)
                                  gccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTgaagaa                        <  1:2083515/68‑1 (MQ=255)
                                  gccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTgaagaa                        <  1:3128304/68‑1 (MQ=255)
                                  gccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTgaagaa                        <  1:11357/68‑1 (MQ=255)
                                  gccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTgaagaa                        <  1:3266783/68‑1 (MQ=255)
                                  gccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTgaagaa                        <  1:347442/68‑1 (MQ=255)
                                  gccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTgaagaa                        <  1:835888/68‑1 (MQ=255)
                                  gccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATCGAAGAACTCTACGCTCTTACTgaagaa                        <  1:1909071/68‑1 (MQ=255)
                                   ccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGaagaag                       <  1:1837128/68‑1 (MQ=255)
                                   ccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGaagaag                       <  1:1115103/68‑1 (MQ=255)
                                   ccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGaagaag                       <  1:122761/68‑1 (MQ=255)
                                   ccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGaagaag                       <  1:3341580/68‑1 (MQ=255)
                                   ccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGaagaag                       <  1:1532221/68‑1 (MQ=255)
                                   ccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGaagaag                       <  1:2811025/68‑1 (MQ=255)
                                   ccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGaagaag                       <  1:2282383/68‑1 (MQ=255)
                                   ccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGaagaag                       <  1:1934645/68‑1 (MQ=255)
                                    cagcagTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTgaag                                               >  1:2340690/1‑43 (MQ=255)
                                    cagcagTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTgaag                                               >  1:387281/1‑43 (MQ=255)
                                    cagcagTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTgaag                                               >  1:199851/1‑43 (MQ=255)
                                    cagcagTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTgaag                                               >  1:2856405/1‑43 (MQ=255)
                                    cagcagTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTgaag                                               >  1:3263785/1‑43 (MQ=255)
                                    cagcagTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTgaag                                               >  1:1019327/1‑43 (MQ=255)
                                    cagcagTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGAt                     <  1:1881443/69‑1 (MQ=255)
                                    cagcagTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGAt                     <  1:1763613/69‑1 (MQ=255)
                                    cagcagTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGAt                     <  1:2326508/69‑1 (MQ=255)
                                         gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:2161221/69‑1 (MQ=255)
                                         gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:734079/69‑1 (MQ=255)
                                         gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:828943/69‑1 (MQ=255)
                                         gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:3227497/69‑1 (MQ=255)
                                         gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:439675/69‑1 (MQ=255)
                                         gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:1089154/69‑1 (MQ=255)
                                         gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:2776336/69‑1 (MQ=255)
                                         gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:1801152/69‑1 (MQ=255)
                                         gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:2621651/69‑1 (MQ=255)
                                         gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:2579543/69‑1 (MQ=255)
                                         gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:520904/69‑1 (MQ=255)
                                         gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:2504122/69‑1 (MQ=255)
                                         gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:2372474/69‑1 (MQ=255)
                                          ttGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:2330944/68‑1 (MQ=255)
                                           tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:3421818/69‑1 (MQ=255)
                                           tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:2102210/69‑1 (MQ=255)
                                           tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:3351657/69‑1 (MQ=255)
                                           tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:3397026/69‑1 (MQ=255)
                                           tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:1373434/69‑1 (MQ=255)
                                           tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:1379231/69‑1 (MQ=255)
                                           tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:2897311/69‑1 (MQ=255)
                                           tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:1463468/69‑1 (MQ=255)
                                           tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:1839333/69‑1 (MQ=255)
                                           tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:619444/69‑1 (MQ=255)
                                           tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:1964518/69‑1 (MQ=255)
                                           tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:2106363/69‑1 (MQ=255)
                                           tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTAATGAAGAAGATTGCCCAg              <  1:3177254/69‑1 (MQ=255)
                                            gctgctCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAg              >  1:663262/1‑68 (MQ=255)
                                            gctgctCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAg              >  1:41493/1‑68 (MQ=255)
                                            gctgctCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAg              >  1:1076748/1‑68 (MQ=255)
                                            gctgctCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAg              >  1:2897169/1‑68 (MQ=255)
                                            gctgctCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:2206796/68‑1 (MQ=255)
                                             ctgctCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:398981/65‑1 (MQ=255)
                                                ctCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGaagaag                       <  1:2485577/55‑1 (MQ=255)
                                                  cgcgCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:2334793/60‑1 (MQ=255)
                                                    cgcTTAAGACGACGCTTCTGATTGAAGAACTCCAGGCTCTTACTGAAGAAGATTGCCCAg              <  1:1486925/60‑1 (MQ=255)
                                                     gcTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAg              <  1:3200928/59‑1 (MQ=255)
                                                     gcTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAg              <  1:2139532/59‑1 (MQ=255)
                                                     gcTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAg              <  1:2507309/59‑1 (MQ=255)
                                                     gcTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAg              <  1:2719296/59‑1 (MQ=255)
                                                     gcTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAg              <  1:2736329/59‑1 (MQ=255)
                                                     gcTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAg              <  1:3033336/59‑1 (MQ=255)
                                                      cTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:2147980/58‑1 (MQ=255)
                                                       ttAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTa                                <  1:2147349/39‑1 (MQ=255)
                                                        tAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAg  >  1:2340184/1‑68 (MQ=255)
                                                        tAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAg  >  1:2925206/1‑68 (MQ=255)
                                                        tAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAg  >  1:831126/1‑68 (MQ=255)
                                                        tAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAg  >  1:3125575/1‑68 (MQ=255)
                                                        tAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAg  >  1:1431247/1‑68 (MQ=255)
                                                            |                                                               
GATATCGTTGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAG  >  minE/2014523‑2014646

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: