Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 457,686 C→T 23.2% L306L (CTG→TTG)  ybgK → predicted enzyme subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE457,6860CT23.2% 127.5 / 32.0 69L306L (CTG→TTG) ybgKpredicted enzyme subunit
Reads supporting (aligned to +/- strand):  ref base C (32/21);  new base T (13/3);  total (45/24)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAGGCTGCGCCAGCGACGCAGAGCTATTAACGC  >  minE/457617‑457756
                                                                     |                                                                      
tgtTCACTGGAAGAGGACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTg                                                                      >  1:2129220/1‑71 (MQ=255)
  ttCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTg                                                                      >  1:703415/1‑70 (MQ=255)
        ggAAGAGGCACTGAAAGCACGGCAAGATCAACAACGTTATTTCGAACAATTAGCGTGGCGGCTGCACAAt                                                                >  1:793842/1‑70 (MQ=255)
        ggAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAAt                                                                >  1:2958328/1‑70 (MQ=255)
         gagcttcagtatgatcccaaagtttatgatgaatatcgtaagacgaagcgggcggGGCGGTTGAACAATGa                                                              >  1:915889‑M1/56‑71 (MQ=255)
         gAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCt                                                                       <  1:4350144/62‑1 (MQ=255)
         gAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTAACAATTAGCGTGGCGGCTGCACAATGaa                                                             <  1:1114765/71‑1 (MQ=255)
          agcttcagtatgatcccaaagtttatgatgaatatcgtaagacgaagcgggcggGGCGGTTGAACAATGa                                                              <  1:2236549‑M1/16‑1 (MQ=255)
           gcttcagtatgatcccaaagtttatgatgaatatcgtaagacgaagcgggcggGGCGGTTGAACAATGaaa                                                            >  1:3717041‑M1/54‑71 (MQ=255)
           gcttcagtatgatcccaaagtttatgatgaatatcgtaagacgaagcgggcggGGCGGTTGAACAATGaaa                                                            >  1:2977298‑M1/54‑71 (MQ=255)
           gcttcagtatgatcccaaagtttatgatgaatatcgtaagacgaagcgggcggGGCGGTTGAACAATGaaa                                                            >  1:3710344‑M1/54‑71 (MQ=255)
           gcttcagtatgatcccaaagtttatgatgaatatcgtaagacgaagcgggcggGGCGGTTGAACAATGa                                                              >  1:240379‑M1/54‑69 (MQ=255)
            gagGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGaaaa                                                           <  1:3948132/71‑1 (MQ=255)
            cttcagtatgatcccaaagtttatgatgaatatcgtaagacgaagcgggcggGGCGGTTGAACAATGaaaa                                                           >  1:297011‑M1/53‑71 (MQ=255)
              tcagtatgatcccaaagtttatgatgaatatcgtaagacgaagcgggcggGGCGGTTGAACAATGAAAAt                                                          <  1:5264732‑M1/20‑1 (MQ=255)
               gCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTg                                                        <  1:3112133/71‑1 (MQ=255)
               gCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTg                                                        <  1:1744701/71‑1 (MQ=255)
                 gtatgatcccaaagtttatgatgaatatcgtaagacgaagcgggcggGGCGGTTGAACAATGAAAatc                                                         >  1:2557310‑M1/48‑67 (MQ=255)
                  tatgatcccaaagtttatgatgaatatcgtaagacgaagcgggcggGGCGGTTGaaca                                                                  >  1:2224193‑M1/47‑58 (MQ=255)
                  cTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACaa                                                                 <  1:5033258/59‑1 (MQ=255)
                  cTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGaa                                                             <  1:205400/63‑1 (MQ=255)
                  cTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGaa                                                             <  1:159483/63‑1 (MQ=255)
                   tAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAAt                                                                >  1:5669189/1‑59 (MQ=255)
                   tAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGa                                                       <  1:2727195/68‑1 (MQ=255)
                    aaaaGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGaa                                                             <  1:4995752/61‑1 (MQ=255)
                     aaaGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGAc                                                      >  1:1015479/1‑67 (MQ=255)
                     aaaGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCt                                                    <  1:3095273/69‑1 (MQ=255)
                     aaaGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTg                                                   >  1:856055/1‑70 (MQ=255)
                       aGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAAc                                                <  1:4796028/71‑1 (MQ=255)
                        gCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGaaaa                                                           >  1:2377816/1‑59 (MQ=255)
                           cGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGcc                                             >  1:1941208/1‑70 (MQ=255)
                           cGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGcc                                             >  1:3988515/1‑70 (MQ=255)
                           cGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCg                                            <  1:471770/71‑1 (MQ=255)
                            aagtttatgatgaatatcgtaagacgaagcgggcggGGCGGTTGAACAATGAAAAt                                                          >  1:1064855‑M1/37‑56 (MQ=255)
                              cAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGAc                                                      >  1:4074634/1‑58 (MQ=255)
                               aaGATCAGCAACGTTATTTTGAAGAATTAGCGTGGCGGCTGCACAATGAAAATTGACCt                                                    <  1:1424377/59‑1 (MQ=255)
                               aaGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGAcc                                                     >  1:2918116/1‑58 (MQ=255)
                               aaGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGAcc                                                     >  1:2400825/1‑58 (MQ=255)
                               aaGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGAcc                                                     >  1:170246/1‑58 (MQ=255)
                                aGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTg                                       >  1:4701182/1‑71 (MQ=255)
                                 gATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGa                                                              <  1:5634265/47‑1 (MQ=255)
                                 gATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGa                                                              <  1:1830136/47‑1 (MQ=255)
                                 gATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGa                                                              <  1:3226732/47‑1 (MQ=255)
                                  aTCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGAc                                                      >  1:1799424/1‑54 (MQ=255)
                                   tCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGc                                    >  1:971858/1‑71 (MQ=255)
                                    cAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGc                                    >  1:739114/1‑70 (MQ=255)
                                        aaCGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAgg                               >  1:2989288/1‑71 (MQ=255)
                                           gTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTgg                                      <  1:3048413/61‑1 (MQ=255)
                                           atcgtaagacgaagcgggcggGGCGGTTGAACAATGAAAAt                                                          >  1:2239276‑M1/22‑41 (MQ=255)
                                            ttATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAGGCTgc                           <  1:1072977/71‑1 (MQ=255)
                                            tcgtaagacgaagcgggcggGGCGGTTGAACAATGa                                                              >  1:5139085‑M1/21‑36 (MQ=255)
                                            tcgtaagacgaagcgggcggGGCGGTTGAACAATGa                                                              >  1:3376426‑M1/21‑36 (MQ=255)
                                            tcgtaagacgaagcgggcggGGCGGTTGAACAATGa                                                              >  1:2228061‑M1/21‑36 (MQ=255)
                                              aTTTTGAACAATTAGCGTGGCGGCTGCACAATGaaa                                                            >  1:3144296/1‑36 (MQ=255)
                                              aTTTTGAACAATTAGCGTGGCGGCTGCACAATGaaa                                                            >  1:4838057/1‑36 (MQ=255)
                                              aTTTTGAACAATTAGCGTGGCGGCTGCACAATGaaa                                                            >  1:2802601/1‑36 (MQ=255)
                                              aTTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAGGCTgcgc                         >  1:2433623/1‑71 (MQ=255)
                                              aTTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAGGCTgcgc                         >  1:698919/1‑71 (MQ=255)
                                               ttttGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGa                                                       >  1:4298451/1‑40 (MQ=255)
                                               ttttGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGAc                                                      >  1:5107746/1‑41 (MQ=255)
                                               ttttGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGAc                                                      >  1:5306963/1‑41 (MQ=255)
                                               ttttGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATc                                         >  1:2187942/1‑54 (MQ=255)
                                                aagacgaagcgggcggGGCGGTTGAACAATGAAAAt                                                          <  1:4967640‑M1/20‑1 (MQ=255)
                                                  tGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAgg                               >  1:1472115/1‑61 (MQ=255)
                                                  tGAACAATTAGCGTGGCGGCTCCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAgg                               >  1:3549911/1‑61 (MQ=255)
                                                       aaTTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAGGCt                             <  1:1749936/58‑1 (MQ=255)
                                                       aaTTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAGGCt                             <  1:793860/58‑1 (MQ=255)
                                                          tAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAGGCTGCGCCAGCGACGCa               >  1:3027193/1‑69 (MQ=255)
                                                           aGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAGGCTGCGCCAGCGACGCagag            >  1:3856260/1‑71 (MQ=255)
                                                                   ggcTGCACAATGAAAATTGACCTGAACGCCGATCTGGGGGAAGGCTGCGCCAGCGACGCAGAGCTATTAAc    <  1:1462941/71‑1 (MQ=255)
                                                                     cTGCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAGGCTGCGCCAGCGACGCAGAGCTATTAACGc  <  1:107209/71‑1 (MQ=255)
                                                                     |                                                                      
TGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAGGCTGCGCCAGCGACGCAGAGCTATTAACGC  >  minE/457617‑457756

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: