Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,104,160 G→A 17.0% D54N (GAC→AAC) ‡ kbaZ → tagatose 6‑phosphate aldolase 1, kbaZ subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,104,1600GA17.0% 279.4 / 42.5 123D54N (GAC→AAC) ‡kbaZtagatose 6‑phosphate aldolase 1, kbaZ subunit
Reads supporting (aligned to +/- strand):  ref base G (49/53);  new base A (16/5);  total (65/58)
Fisher's exact test for biased strand distribution p-value = 2.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 128 total aligned reads displayed.

GAAGCTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTC  >  minE/2104091‑2104230
                                                                     |                                                                      
gAAGCTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGa                                                                       >  1:1036656/1‑71 (MQ=255)
 aaGCTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGAc                                                                      <  1:2397512/71‑1 (MQ=255)
    cTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAg                                                                   <  1:5440345/71‑1 (MQ=255)
    cTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAg                                                                   <  1:2181598/71‑1 (MQ=255)
       cAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTc                                                                <  1:2613918/71‑1 (MQ=255)
       cAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTc                                                                >  1:1552048/1‑71 (MQ=255)
       cAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTc                                                                >  1:4476016/1‑71 (MQ=255)
       cAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTc                                                                <  1:1891971/71‑1 (MQ=255)
          tCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTcggc                                                             >  1:3645520/1‑71 (MQ=255)
          tCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTcggc                                                             <  1:183417/71‑1 (MQ=255)
            cGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGAc                                                                      >  1:1308532/1‑60 (MQ=255)
            cGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGAc                                                                      >  1:4230862/1‑60 (MQ=255)
            cGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTcggcgg                                                           >  1:3361451/1‑71 (MQ=255)
             gCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGt                                                          <  1:948892/71‑1 (MQ=255)
             gCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGt                                                          <  1:390094/71‑1 (MQ=255)
               tACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTta                                                        >  1:1221334/1‑71 (MQ=255)
                actccAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTta                                                        >  1:4658344/4‑70 (MQ=255)
                 cGCCAGTGCAAACCAAACGCCGCTACTGATTGAAGCAACCTCCAATCAGGTAGAc                                                                      <  1:2169827/55‑1 (MQ=255)
                   ccAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAg                                                                        >  1:5311209/1‑51 (MQ=255)
                     aGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCgg                                                  <  1:4261789/71‑1 (MQ=255)
                       tGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGaa                                                <  1:2162299/71‑1 (MQ=255)
                         caatccaaacGCCGTTACTGATTGAAGCAACCTCCAATCAGTTAGACCAGTTCGGCGGTTATACCGGAAt                                               <  1:2460635/70‑1 (MQ=255)
                         caaaccaaacGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAAt                                               >  1:1291697/1‑70 (MQ=255)
                          aaaccaaacGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGa                                             <  1:961907/71‑1 (MQ=255)
                           aaccaaacGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAg                                                                        <  1:447841/43‑1 (MQ=255)
                           aaccaaacGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAg                                                                        <  1:833901/43‑1 (MQ=255)
                           aaccaaacGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTta                                                        >  1:1162522/1‑59 (MQ=255)
                              caaacGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGAcgcc                                         >  1:4833883/1‑71 (MQ=255)
                                 acGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTtata                                                      >  1:5337673/1‑55 (MQ=255)
                                  cGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGAcgcccgcc                                     >  1:2849938/1‑71 (MQ=255)
                                  cGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGAcgcccgcc                                     >  1:2964781/1‑71 (MQ=255)
                                   gCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCccgccg                                    <  1:5108881/71‑1 (MQ=255)
                                      gTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAtt                                 <  1:4663429/71‑1 (MQ=255)
                                      gTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAtt                                 <  1:5237288/71‑1 (MQ=255)
                                      gTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAtt                                 <  1:4670184/71‑1 (MQ=255)
                                       ttACTGATTGAAGCAACCTCCAATCAGGTAGACCAg                                                                   <  1:872318/36‑1 (MQ=255)
                                       ttACTGATTGAAGCAACCTCCAATCAGGTAGACCAg                                                                   <  1:819041/36‑1 (MQ=255)
                                       ttACTGATTGAAGCAACCTCCAATCAGGTAGACCAg                                                                   <  1:2823204/36‑1 (MQ=255)
                                       ttACTGATTGAAGCAACCTCCAATCAGGTAGACCAg                                                                   <  1:1508459/36‑1 (MQ=255)
                                       ttACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGt                                                          <  1:5040655/45‑1 (MQ=255)
                                       ttACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGt                                                          <  1:2628446/45‑1 (MQ=255)
                                       ttACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTtata                                                      <  1:5157905/49‑1 (MQ=255)
                                       ttACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTta                                                        <  1:3134208/47‑1 (MQ=255)
                                       ttACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAttt                                <  1:4517395/71‑1 (MQ=255)
                                       ttACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAtt                                 <  1:571206/70‑1 (MQ=255)
                                        tgcTGATTGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTta                                                        <  1:5381246/44‑1 (MQ=255)
                                        tgcTGATTGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGActtt                               >  1:1452742/3‑71 (MQ=255)
                                        tACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTtata                                                      <  1:4972172/48‑1 (MQ=255)
                                        tACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCg                                                   <  1:1515090/51‑1 (MQ=255)
                                         aCTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTcggcg                                                            >  1:3768371/1‑41 (MQ=255)
                                         aCTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCAGCCGATTTTc                              >  1:3781488/1‑71 (MQ=255)
                                          cTGATTGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTtata                                                      >  1:2236531/1‑46 (MQ=255)
                                          cTGATTGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTtata                                                      >  1:4198216/1‑46 (MQ=255)
                                          cTGATTGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTtata                                                      >  1:4144947/1‑46 (MQ=255)
                                          cTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTcg                             >  1:4115025/1‑71 (MQ=255)
                                           tGATTGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGActt                                <  1:4642488/67‑1 (MQ=255)
                                            gATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTcgcg                           >  1:1766241/1‑71 (MQ=255)
                                             aTTGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTcgcg                           >  1:3244603/1‑70 (MQ=255)
                                             aTTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAAt                                               <  1:1889668/50‑1 (MQ=255)
                                             aTTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAAt                                               <  1:4570683/50‑1 (MQ=255)
                                              ttGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTcgcgaa                         >  1:2560392/1‑69 (MQ=255)
                                              ttGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTcgcgaa                         >  1:2350933/1‑69 (MQ=255)
                                              ttGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGc                         >  1:5301465/1‑71 (MQ=255)
                                               tGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCgg                                                  >  1:2220899/1‑45 (MQ=255)
                                               tGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTcgcgaat                        >  1:5112299/1‑68 (MQ=255)
                                               tGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTcgcgaat                        >  1:2493561/1‑68 (MQ=255)
                                               tGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTcgcgaat                        >  1:3330137/1‑68 (MQ=255)
                                               tGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTtat                                                       <  1:4668275/40‑1 (MQ=255)
                                               tGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTtat                                                       <  1:4066572/40‑1 (MQ=255)
                                               tGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACTCCCGCCGATTTTCGCgg                          >  1:189999/1‑69 (MQ=255)
                                                gAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTta                                                        >  1:2316795/1‑38 (MQ=255)
                                                  aGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTta                                                        >  1:1977596/1‑36 (MQ=255)
                                                  aGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTta                                                        >  1:1416292/1‑36 (MQ=255)
                                                    cAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGaa                                                >  1:2724506/1‑42 (MQ=255)
                                                     aaCGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGa                                             >  1:491660/1‑44 (MQ=255)
                                                     aaCGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGa                                             >  1:3224518/1‑44 (MQ=255)
                                                     aaCGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTCGCGAAtttgttt                  <  1:5011078/71‑1 (MQ=255)
                                                     aaCGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTCGCGAAtttgttt                  <  1:4287668/71‑1 (MQ=255)
                                                     aaCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgttt                  <  1:3228111/71‑1 (MQ=255)
                                                     aaCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgttt                  <  1:5089189/71‑1 (MQ=255)
                                                      aCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGAcgc                                          >  1:982795/1‑46 (MQ=255)
                                                      aCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAt                                  <  1:4117836/54‑1 (MQ=255)
                                                      aCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAt                                  <  1:959601/54‑1 (MQ=255)
                                                      aCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgt                    >  1:3450725/1‑68 (MQ=255)
                                                      aCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgt                    >  1:4174387/1‑68 (MQ=255)
                                                       ccTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAAt                                               <  1:5441795/40‑1 (MQ=255)
                                                       ccTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAAt                                               <  1:5025624/40‑1 (MQ=255)
                                                        cTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttg                     <  1:5583957/65‑1 (MQ=255)
                                                        cTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTc               <  1:3011683/71‑1 (MQ=255)
                                                         tCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCt                        >  1:5653725/1‑61 (MQ=255)
                                                             aTCAGGTAGACCAGTTCGGCGGTTATACCGGAATGa                                             >  1:3025590/1‑36 (MQ=255)
                                                             aTCAGGTAGACCAGTTCGGCGGTTATACCGGAATGa                                             >  1:4285320/1‑36 (MQ=255)
                                                              tCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttg                     <  1:1641440/59‑1 (MQ=255)
                                                               cAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGc        >  1:3287018/1‑71 (MQ=255)
                                                                aGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTcgc                            >  1:4649050/1‑50 (MQ=255)
                                                                aGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGcc       >  1:4632385/1‑71 (MQ=255)
                                                                aGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGcc       >  1:2577266/1‑71 (MQ=255)
                                                                aGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGcc       >  1:1520900/1‑71 (MQ=255)
                                                                  gTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgttt                  <  1:1211817/58‑1 (MQ=255)
                                                                   tAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGa     >  1:5044543/1‑70 (MQ=255)
                                                                     |                                                                      
GAAGCTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTC  >  minE/2104091‑2104230

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: