Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,507,263 C→T 32.2% S109S (AGC→AGT yicI → predicted alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,507,2630CT32.2% 89.5 / 51.2 102S109S (AGC→AGTyicIpredicted alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (49/20);  new base T (12/21);  total (61/41)
Fisher's exact test for biased strand distribution p-value = 1.18e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.47e-01

Only 100 of 103 total aligned reads displayed.

GACGTGAAGGTCACAATCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGAACGT  >  minE/2507201‑2507320
                                                              |                                                         
gACGTGAAGGTCACAATCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGTCAACTTAAGCGCGCgt                                                     >  1:3179520/1‑69 (MQ=255)
gACGTGAAGGTCACAATCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGcgcgcg                                                      >  1:1129084/1‑68 (MQ=255)
gACGTGAAGGTCACAATCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCgt                                                     >  1:489876/1‑69 (MQ=255)
gACGTGAAGGTCACAATCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCgt                                                     >  1:1897690/1‑69 (MQ=255)
gACGTGAAGGTCACAATCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCgt                                                     >  1:3284518/1‑69 (MQ=255)
gACGTGAAGGTCACAATCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCgt                                                     >  1:132678/1‑69 (MQ=255)
gACGTGAAGGTCACAATCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCgt                                                     >  1:2162254/1‑69 (MQ=255)
           cacaATCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAg                                           >  1:814985/1‑68 (MQ=255)
           cacaATCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAg                                           >  1:1770684/1‑68 (MQ=255)
           cacaATCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAg                                           >  1:1624888/1‑68 (MQ=255)
           cacaATCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAg                                           >  1:1914528/1‑68 (MQ=255)
           cacaATCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAg                                           >  1:1994488/1‑68 (MQ=255)
           cacaATCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAg                                           >  1:2024738/1‑68 (MQ=255)
           cacaATCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAg                                           >  1:263062/1‑68 (MQ=255)
           cacaATCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAg                                           >  1:2755350/1‑68 (MQ=255)
           cacaATCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAg                                           >  1:2944247/1‑68 (MQ=255)
           cacaATCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAg                                           >  1:1419595/1‑68 (MQ=255)
           cacaATCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAg                                           >  1:1304547/1‑68 (MQ=255)
            acaATCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAg                                                <  1:1452021/62‑1 (MQ=255)
            acaATCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAg                                                >  1:2947342/1‑62 (MQ=255)
            acaATCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAg                                                >  1:309652/1‑62 (MQ=255)
            acaATCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAg                                                >  1:1157769/1‑62 (MQ=255)
            acaATCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAg                                                >  1:432307/1‑62 (MQ=255)
               aTCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAg                                      >  1:2994554/1‑69 (MQ=255)
                tCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGt                                     >  1:1045988/1‑69 (MQ=255)
                 cGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGtt                                    <  1:2531429/69‑1 (MQ=255)
                    aaaCACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTg                                 >  1:1787335/1‑69 (MQ=255)
                    aaaCACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTg                                 >  1:91102/1‑69 (MQ=255)
                    aaaCACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTg                                 >  1:1648499/1‑69 (MQ=255)
                    aaaCACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTg                                 >  1:3100804/1‑69 (MQ=255)
                    aaaCACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTg                                 >  1:1566955/1‑69 (MQ=255)
                    aaaCACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTg                                 >  1:1384797/1‑69 (MQ=255)
                    aaaCACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTg                                 >  1:2880656/1‑69 (MQ=255)
                    aaaCACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTg                                 >  1:1509960/1‑69 (MQ=255)
                       cacaGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTc                              >  1:937209/1‑69 (MQ=255)
                       cacaGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTc                              >  1:821821/1‑69 (MQ=255)
                       cacaGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTc                              >  1:1040025/1‑69 (MQ=255)
                       cacaGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTc                              >  1:330390/1‑69 (MQ=255)
                       cacaGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTc                              >  1:1593936/1‑69 (MQ=255)
                       cacaGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTc                              >  1:2084951/1‑69 (MQ=255)
                                 tATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGAtttt                    >  1:1618986/1‑69 (MQ=255)
                                       gAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGa                        >  1:2125865/1‑59 (MQ=255)
                                       gAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGa                        >  1:3062337/1‑59 (MQ=255)
                                       gAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGa                        >  1:2595002/1‑59 (MQ=255)
                                       gAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTgcgc              <  1:755344/69‑1 (MQ=255)
                                       gAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTgcgc              <  1:813043/69‑1 (MQ=255)
                                       gAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTgcgc              <  1:1155296/69‑1 (MQ=255)
                                        aGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCa             >  1:2395023/1‑69 (MQ=255)
                                        aGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCa             >  1:885636/1‑69 (MQ=255)
                                           ttAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACg          <  1:583862/69‑1 (MQ=255)
                                           ttAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACg          <  1:444888/69‑1 (MQ=255)
                                           ttAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACg          <  1:1167097/69‑1 (MQ=255)
                                            tAAAAGTGGTAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACg          >  1:1445236/1‑68 (MQ=255)
                                            tAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTc                                   >  1:2063608/1‑43 (MQ=255)
                                            tAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACg          >  1:840686/1‑68 (MQ=255)
                                            tAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACg          >  1:39085/1‑68 (MQ=255)
                                            tAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACg          >  1:1715983/1‑68 (MQ=255)
                                                aGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGaa     <  1:2974045/69‑1 (MQ=255)
                                                aGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGaa     <  1:2133229/69‑1 (MQ=255)
                                                aGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGaa     <  1:1873056/69‑1 (MQ=255)
                                                aGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGaa     <  1:718783/69‑1 (MQ=255)
                                                aGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGaa     <  1:1093271/69‑1 (MQ=255)
                                                aGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGaa     <  1:337659/69‑1 (MQ=255)
                                                aGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGaa     <  1:3368444/69‑1 (MQ=255)
                                                aGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGaa     <  1:1208023/69‑1 (MQ=255)
                                                aGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGaa     <  1:2898133/69‑1 (MQ=255)
                                                aGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGaa     <  1:3154166/69‑1 (MQ=255)
                                                aGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACCGCGaa     <  1:766100/69‑1 (MQ=255)
                                                   ggCAACTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGAACGt  <  1:2392055/69‑1 (MQ=255)
                                                   ggCAACTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGAACGt  <  1:2955594/69‑1 (MQ=255)
                                                   ggCAACTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGAACGt  <  1:2102190/69‑1 (MQ=255)
                                                   ggCAACTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGAACGt  <  1:2317878/69‑1 (MQ=255)
                                                   ggCAACTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGAACGt  <  1:834046/69‑1 (MQ=255)
                                                   ggCAACTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGAACGt  <  1:1875381/69‑1 (MQ=255)
                                                   ggCAACTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGAACGt  <  1:814865/69‑1 (MQ=255)
                                                   ggCAACTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGAACGt  <  1:788581/69‑1 (MQ=255)
                                                   ggCAACTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGAACGt  <  1:2345969/69‑1 (MQ=255)
                                                   ggCAACTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGAACGt  <  1:3084847/69‑1 (MQ=255)
                                                   ggCAACTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGAACGt  <  1:312703/69‑1 (MQ=255)
                                                   ggCAACTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGAACGt  <  1:426623/69‑1 (MQ=255)
                                                   ggCAACTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGAACGt  <  1:1455534/69‑1 (MQ=255)
                                                   ggCAACTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGAACGt  <  1:3198452/69‑1 (MQ=255)
                                                   ggCAACTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGAACGt  <  1:388579/69‑1 (MQ=255)
                                                   ggCAACTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGAACGt  <  1:3383150/69‑1 (MQ=255)
                                                   ggCAACTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGAACGt  <  1:1461732/69‑1 (MQ=255)
                                                   ggCAACTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGAACGt  <  1:3341496/69‑1 (MQ=255)
                                                   ggCAACTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGAACGt  <  1:192343/69‑1 (MQ=255)
                                                    gCAACTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGAACGt  <  1:2854851/68‑1 (MQ=255)
                                                      aaCTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTg                          >  1:2327484/1‑42 (MQ=255)
                                                      aaCTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTg                          >  1:1326446/1‑42 (MQ=255)
                                                      aaCTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTg                          >  1:873565/1‑42 (MQ=255)
                                                      aaCTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTg                          >  1:3273358/1‑42 (MQ=255)
                                                      aaCTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTg                          >  1:816881/1‑42 (MQ=255)
                                                      aaCTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTg                          >  1:2323516/1‑42 (MQ=255)
                                                      aaCTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTg                          >  1:810398/1‑42 (MQ=255)
                                                      aaCTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTg                          >  1:1282742/1‑42 (MQ=255)
                                                      aaCTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTg                          >  1:691744/1‑42 (MQ=255)
                                                      aaCTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTg                          >  1:47921/1‑42 (MQ=255)
                                                      aaCTTAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTg                          >  1:2618305/1‑42 (MQ=255)
                                                         ttAAGTGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGAACGt  <  1:2076685/63‑1 (MQ=255)
                                                              |                                                         
GACGTGAAGGTCACAATCGAAAACACAGAACGTTATGCTGAGTTTAAAAGTGGCAACTTAAGCGCGCGTGTCAGCAAAGGTGAGTTCTGGTCACTGGATTTTCTGCGCAACGGCGAACGT  >  minE/2507201‑2507320

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: