Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,237,336 C→A 6.9% R395R (CGG→CGT edd ← 6‑phosphogluconate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,237,3360CA 6.9% 281.5 / 3.7 84R395R (CGG→CGTedd6‑phosphogluconate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (39/39);  new base A (3/3);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGC  >  minE/1237266‑1237401
                                                                      |                                                                 
caccaTGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCTAGTGATTTTTCCGCCCCTTcc                                                                   >  1:2823788/1‑71 (MQ=255)
 accaTGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTccc                                                                  >  1:466266/1‑71 (MQ=255)
   caTGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGc                                                                >  1:688716/1‑71 (MQ=255)
   caTGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGc                                                                >  1:1984271/1‑71 (MQ=255)
        tgaGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGccagtc                                                           <  1:2514781/71‑1 (MQ=255)
        tgaGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGccagtc                                                           <  1:1414921/71‑1 (MQ=255)
          agagaAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGccagtcc                                                          <  1:999423/70‑1 (MQ=255)
              aaaGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGccagtc                                                           <  1:1033728/65‑1 (MQ=255)
               aaGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTccc                                                                  >  1:2154920/1‑57 (MQ=255)
               aaGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTccc                                                                  >  1:3097116/1‑57 (MQ=255)
                aGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCAc                                                   <  1:2308415/71‑1 (MQ=255)
                aGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCAc                                                   >  1:2948898/1‑71 (MQ=255)
                 ggTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGccagtccagt                                                       <  1:366194/66‑1 (MQ=255)
                  gttgttCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCa                                                 >  1:1579423/1‑71 (MQ=255)
                   ttgttCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCat                                                <  1:2029205/71‑1 (MQ=255)
                     gttCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCatt                                               <  1:1249041/70‑1 (MQ=255)
                     gttCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCACGCCAGTCCAGTTCAc                                                   >  1:3731427/1‑66 (MQ=255)
                     gttCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCACGCCAGTCCAGTTCAc                                                   >  1:1652616/1‑66 (MQ=255)
                     attCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCACGCCAGTCCAGTTCAc                                                   >  1:3362817/2‑66 (MQ=255)
                      ttCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCattat                                             <  1:2193351/71‑1 (MQ=255)
                        cGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCattatt                                            >  1:1314794/1‑70 (MQ=255)
                        cGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCattatt                                            >  1:2461531/1‑70 (MQ=255)
                         gtaggaagCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTc                                           <  1:1441107/68‑1 (MQ=255)
                         gaaggaagCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTcc                                                                   >  1:3340608/1‑46 (MQ=255)
                         gaaggaagCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTcc                                                                   >  1:3344950/1‑46 (MQ=255)
                         gaaggaagCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTcc                                                                   >  1:3431853/1‑46 (MQ=255)
                         gaaggaagCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTcc                                                                   >  1:3457617/1‑46 (MQ=255)
                         gaaggaagCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTcc                                                                   >  1:2681757/1‑46 (MQ=255)
                         gaaggaagCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTcc                                                                   >  1:417922/1‑46 (MQ=255)
                         gaaggaagCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTcc                                                                   >  1:2521589/1‑46 (MQ=255)
                         gaaggaagCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGcc                                                               <  1:3243630/50‑1 (MQ=255)
                         gaaggaagCGATCACATTGCTGTCGAGTGATTTTTCCGCCACTTCCCGCCAGTCCAGTTCACAattat                                             >  1:715132/1‑68 (MQ=255)
                           aggaagCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGc                                        >  1:2199466/1‑71 (MQ=255)
                            ggaagCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGc                                        >  1:2579679/1‑70 (MQ=255)
                              aagCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCa                                      >  1:2182376/1‑70 (MQ=255)
                              aagCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCAt                                     >  1:3044979/1‑71 (MQ=255)
                              aagCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCAt                                     <  1:2635462/71‑1 (MQ=255)
                               agCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATg                                    <  1:1836239/71‑1 (MQ=255)
                                gCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATgg                                   <  1:2336635/71‑1 (MQ=255)
                                gCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATgg                                   >  1:3302944/1‑71 (MQ=255)
                                  gATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGccagtccag                                                        <  1:1200541/48‑1 (MQ=255)
                                  gATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGtt                                 >  1:458273/1‑71 (MQ=255)
                                    tCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCTTGGTTCa                               <  1:3582607/71‑1 (MQ=255)
                                    tCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCa                               <  1:1718725/71‑1 (MQ=255)
                                    tCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCa                               <  1:1089694/71‑1 (MQ=255)
                                     cacaTTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTc                                >  1:792630/1‑69 (MQ=255)
                                      acaTTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCaa                              <  1:3504351/70‑1 (MQ=255)
                                       caTTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAgg                            >  1:2521633/1‑71 (MQ=255)
                                           gCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGta                         >  1:1924082/1‑70 (MQ=255)
                                            cTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCa                                                    >  1:81502/1‑42 (MQ=255)
                                            cTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCa                                                    >  1:495992/1‑42 (MQ=255)
                                            cTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGtat                        <  1:3335312/70‑1 (MQ=255)
                                              gTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCa                               >  1:142090/1‑61 (MQ=255)
                                              gTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAAc                     >  1:1259517/1‑71 (MQ=255)
                                              gTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAAc                     >  1:127340/1‑71 (MQ=255)
                                               tCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCAt                                     <  1:417274/54‑1 (MQ=255)
                                               tCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCAt                                     <  1:2540341/54‑1 (MQ=255)
                                               tCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCAt                                     <  1:639425/54‑1 (MQ=255)
                                               tCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAAc                     >  1:2259768/1‑70 (MQ=255)
                                                 gAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACga                   <  1:2599606/70‑1 (MQ=255)
                                                 gAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTAAGCCATGGTTCAAGGGTATAACga                   <  1:1962653/70‑1 (MQ=255)
                                                  aGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACgaga                 <  1:3424941/71‑1 (MQ=255)
                                                  aGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACgaga                 <  1:330235/71‑1 (MQ=255)
                                                  aGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACgaga                 <  1:1388142/71‑1 (MQ=255)
                                                   gTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCa                                                 <  1:1830897/38‑1 (MQ=255)
                                                   gTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCa                                                 <  1:1029398/38‑1 (MQ=255)
                                                   gTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCa                                                 <  1:2887866/38‑1 (MQ=255)
                                                   gTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGagac                <  1:1082549/71‑1 (MQ=255)
                                                   gTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGagac                <  1:2608550/71‑1 (MQ=255)
                                                   gTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGagac                <  1:998473/71‑1 (MQ=255)
                                                      aTTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGagacag              <  1:2005944/70‑1 (MQ=255)
                                                       tttttCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCa                               <  1:2905954/52‑1 (MQ=255)
                                                        ttttCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGagacag              >  1:3472023/1‑68 (MQ=255)
                                                        ttttCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGagacag              >  1:1668406/1‑68 (MQ=255)
                                                         tttCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATgg                                   >  1:2984062/1‑46 (MQ=255)
                                                         tttCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAg                             >  1:2694221/1‑52 (MQ=255)
                                                          ttCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCa          >  1:1767186/1‑70 (MQ=255)
                                                             cGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAAc      <  1:2090558/71‑1 (MQ=255)
                                                             cGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAAc      <  1:2605200/71‑1 (MQ=255)
                                                               ccccTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTc                                >  1:690126/1‑43 (MQ=255)
                                                               ccccTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTc                                >  1:243553/1‑43 (MQ=255)
                                                               ccccTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTc                                >  1:2389516/1‑43 (MQ=255)
                                                               ccccTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTc                                >  1:2960466/1‑43 (MQ=255)
                                                               ccccTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACg                    >  1:163093/1‑55 (MQ=255)
                                                               ccccTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCaaaa       >  1:2233887/1‑68 (MQ=255)
                                                               ccccTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCt    <  1:3517602/71‑1 (MQ=255)
                                                                cccTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAGCCTg   >  1:1654202/1‑71 (MQ=255)
                                                                cccTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTg   >  1:1143629/1‑71 (MQ=255)
                                                                 ccTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGc  >  1:3080421/1‑71 (MQ=255)
                                                                  cTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGc  >  1:2617856/1‑70 (MQ=255)
                                                                  cTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGc  >  1:233988/1‑70 (MQ=255)
                                                                  cTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGc  >  1:487431/1‑70 (MQ=255)
                                                                  cTTCACGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGagaca               <  1:2506093/57‑1 (MQ=255)
                                                                  cTTCACGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGagaca               <  1:722335/57‑1 (MQ=255)
                                                                  cTTCACGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGagaca               <  1:2538471/57‑1 (MQ=255)
                                                                      ccGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGtat                        <  1:2478103/44‑1 (MQ=255)
                                                                      ccGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGtat                        <  1:784267/44‑1 (MQ=255)
                                                                      ccGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGtat                        <  1:489064/44‑1 (MQ=255)
                                                                      ccGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGtat                        <  1:3058395/44‑1 (MQ=255)
                                                                      |                                                                 
CACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGC  >  minE/1237266‑1237401

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: