Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,910,050 C→T 5.1% L348L (CTG→TTG)  mgtA → magnesium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,910,0500CT 5.1% 303.2 / 4.2 97L348L (CTG→TTG) mgtAmagnesium transporter
Reads supporting (aligned to +/- strand):  ref base C (42/49);  new base T (2/3);  total (44/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

AGGCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGT  >  minE/2909986‑2910119
                                                                |                                                                     
aGGCCTAATGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCg                                                                  >  1:3599319/1‑70 (MQ=255)
aGGCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGAGCGCGGCGCAGTAAAGCTGTCg                                                                  >  1:4120909/1‑70 (MQ=255)
 ggCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGa                                                                 <  1:3307883/70‑1 (MQ=255)
 ggCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGa                                                                 >  1:4252569/1‑70 (MQ=255)
 ggCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGa                                                                 <  1:2850571/70‑1 (MQ=255)
 ggCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGATGTCGa                                                                 <  1:4219598/70‑1 (MQ=255)
  gCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGaaa                                                               <  1:4123414/71‑1 (MQ=255)
   ccTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGaaa                                                               >  1:4643111/1‑70 (MQ=255)
    ttAACGCCAGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGTGGCGCAGTCAAGTTGTCGaaa                                                               >  1:2381115/2‑69 (MQ=255)
    ttAACGCCAGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGTGGCGCAGTCAAGTTGTCGaaa                                                               >  1:2386504/2‑69 (MQ=255)
    cTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACa                                                             <  1:615493/71‑1 (MQ=255)
     tAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGaa                                                                <  1:1104865/67‑1 (MQ=255)
       aCGCCGGAAATGTTGCCGATGATTGTTACCTCTACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGa                                                           <  1:3157112/70‑1 (MQ=255)
       aCGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGaaa                                                               >  1:2026126/1‑66 (MQ=255)
       aCGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGa                                                           >  1:1544116/1‑70 (MQ=255)
       aCGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGa                                                           >  1:4326121/1‑70 (MQ=255)
       aCGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGa                                                           <  1:1009311/70‑1 (MQ=255)
       aCGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGa                                                           <  1:1039953/70‑1 (MQ=255)
          ccGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGt                                                       <  1:3578275/71‑1 (MQ=255)
          ccGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGt                                                       >  1:2545152/1‑71 (MQ=255)
           cGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTg                                                      <  1:3903336/71‑1 (MQ=255)
           cGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTg                                                      <  1:3738283/71‑1 (MQ=255)
           cGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTg                                                      <  1:1216295/71‑1 (MQ=255)
           cGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTg                                                      <  1:2107132/71‑1 (MQ=255)
           cGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTg                                                      <  1:4428809/71‑1 (MQ=255)
            ggAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGa                                                     >  1:3373878/1‑71 (MQ=255)
             gAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGAt                                                    >  1:2362990/1‑71 (MQ=255)
                 tGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGt                                                 >  1:2929383/1‑70 (MQ=255)
                  gTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTc                                                >  1:239013/1‑70 (MQ=255)
                   ttGCCGATGATTGTTACCTCGACGCTGGCGCGCGTCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCa                                               <  1:671322/70‑1 (MQ=255)
                   ttGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCaa                                              <  1:2858664/71‑1 (MQ=255)
                   ttGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCaa                                              <  1:2685087/71‑1 (MQ=255)
                   ttGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCaa                                              <  1:3651492/71‑1 (MQ=255)
                   ttGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCa                                               <  1:463076/70‑1 (MQ=255)
                   ttGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCa                                               <  1:562136/70‑1 (MQ=255)
                     gCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCaaa                                             >  1:3451520/1‑70 (MQ=255)
                      ccGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACa                                           >  1:2720015/1‑71 (MQ=255)
                      ccGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACa                                           <  1:1880165/71‑1 (MQ=255)
                       cGATGATTGTTACCTCGACTCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGa                                                           <  1:2721881/54‑1 (MQ=255)
                       cGATGATTGTTACCTCGACGCTGGCGCGTGGCGCAGTCAAGTTGTCGAAAc                                                              <  1:3306375/51‑1 (MQ=255)
                       cGATGATTGTTACCTCGACGCTGGCGCGTGGCGCAGTCAAGTTGTCGAAAc                                                              <  1:1901069/51‑1 (MQ=255)
                        gatgatTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCaaa                                             <  1:2127847/67‑1 (MQ=255)
                         atgatTGTTACCTCGACGCTGGCGCGTGGCGCAGTCAAGTTGTCGAAAc                                                              <  1:4590905/49‑1 (MQ=255)
                         atgatTGTTACCTCGACGCTGGCGCGCGGCGCATTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCt                                        >  1:861799/1‑71 (MQ=255)
                         atgatTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGa                                                           <  1:705817/52‑1 (MQ=255)
                         atgatTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCt                                        >  1:3359245/1‑71 (MQ=255)
                         atgatTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCt                                        >  1:143791/1‑71 (MQ=255)
                          tgatTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTg                                       <  1:531816/71‑1 (MQ=255)
                          agatTGTTACCTCGACGCTGGCGCGCGGCCCAGTAAAGc                                                                       >  1:2642516/2‑39 (MQ=39)
                           gatTGTTACCTCGACGCTGTCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTg                                       >  1:3108747/1‑70 (MQ=255)
                           gatTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGa                                                                 <  1:3913454/44‑1 (MQ=255)
                           gatTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTg                                       >  1:160976/1‑70 (MQ=255)
                           gatTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTg                                       >  1:4143644/1‑70 (MQ=255)
                           gatTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTg                                       >  1:1388070/1‑70 (MQ=255)
                            atTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGa                                                           <  1:1044164/49‑1 (MQ=255)
                            atTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTgg                                      <  1:1354273/70‑1 (MQ=255)
                            atTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGa                                     <  1:192457/71‑1 (MQ=255)
                             ttGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGaaa                                                         >  1:4532606/1‑50 (MQ=255)
                             ttGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGaaa                                                         >  1:4111310/1‑50 (MQ=255)
                             ttGTTACCTCGACGCTGGCGCGCGGCCCAGTAAAGCTGTCGAAACAGAAAg                                                        <  1:2578572/51‑1 (MQ=255)
                              tGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACAt                                          >  1:4597679/1‑64 (MQ=255)
                               gTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGaa                                                          >  1:2636248/1‑47 (MQ=255)
                               gTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGaa                                                          >  1:1997277/1‑47 (MQ=255)
                               gTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGaa                                                          >  1:4561789/1‑47 (MQ=255)
                               gTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATg                                   <  1:3727291/70‑1 (MQ=255)
                                 tACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGaaa                                                               >  1:3360145/1‑40 (MQ=255)
                                  aCCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTa                                <  1:4013854/70‑1 (MQ=255)
                                  aCCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTAt                               >  1:3142645/1‑71 (MQ=255)
                                   ccTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTAtt                              <  1:1411317/71‑1 (MQ=255)
                                       gACGCTGGCGCGGGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAg                           >  1:2793213/1‑70 (MQ=255)
                                       gACGCTGGCGCGCGGGGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAg                           >  1:1337076/1‑70 (MQ=255)
                                       gACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAg                           >  1:3731149/1‑70 (MQ=255)
                                       gACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAg                           >  1:3375247/1‑70 (MQ=255)
                                        aCGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGa                          <  1:4223449/70‑1 (MQ=255)
                                           cTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACtt                      <  1:2011477/71‑1 (MQ=255)
                                             ggCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTc                                                >  1:3305153/1‑43 (MQ=255)
                                              tcgcgcgGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTg                    <  1:1232690/69‑1 (MQ=255)
                                              gcgcgcGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACttt                     <  1:1311047/69‑1 (MQ=255)
                                              gcgcgcGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTg                    <  1:3789960/70‑1 (MQ=255)
                                              gcgcgcGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTg                    <  1:943801/70‑1 (MQ=255)
                                              gcgcgcGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTg                    <  1:1623193/70‑1 (MQ=255)
                                                 cgcgGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGgcgc                <  1:4517149/71‑1 (MQ=255)
                                                 cgcgGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGgcgc                <  1:4241744/71‑1 (MQ=255)
                                                 cgcgGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGgcgc                <  1:1763110/71‑1 (MQ=255)
                                                  gcggcgCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTgg                                      >  1:1951002/1‑48 (MQ=255)
                                                      cgcAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATgg           <  1:519611/71‑1 (MQ=255)
                                                        cAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGa          >  1:3114228/1‑70 (MQ=255)
                                                        cAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGa          >  1:455486/1‑70 (MQ=255)
                                                        cAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGa          >  1:3057025/1‑70 (MQ=255)
                                                         aGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTg                                       <  1:3932181/40‑1 (MQ=255)
                                                         aGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGata        <  1:312752/71‑1 (MQ=255)
                                                          gTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGatat       <  1:1675205/71‑1 (MQ=255)
                                                             aaGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCaa              >  1:1613730/1‑61 (MQ=255)
                                                             aaGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCaa              >  1:4550335/1‑61 (MQ=255)
                                                              aGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGgcgc                <  1:158849/58‑1 (MQ=255)
                                                              aGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGgcgc                <  1:4519367/58‑1 (MQ=255)
                                                               gCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATg                                   >  1:2332120/1‑38 (MQ=255)
                                                               gCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCtgt  >  1:4570897/1‑71 (MQ=255)
                                                                cTGTCGAAACAGAAAGTGATCGTCAAACATCTGGAt                                    >  1:1502844/1‑36 (MQ=255)
                                                                |                                                                     
AGGCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGT  >  minE/2909986‑2910119

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: ‑