Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,910,037 C→T 5.6% R343R (CGC→CGT mgtA → magnesium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,910,0370CT 5.6% 271.4 / 4.3 89R343R (CGC→CGTmgtAmagnesium transporter
Reads supporting (aligned to +/- strand):  ref base C (38/45);  new base T (2/3);  total (41/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 103 total aligned reads displayed.

TGCGCTTTCGGTCGCGGTAGGCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGC  >  minE/2909968‑2910105
                                                                     |                                                                    
tGCGCTTTCGGTCGCGGTAGGCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGcgcgcg                                                                     >  1:128622/1‑71 (MQ=255)
tGCGCTTTCGGTCGCGGTAGGCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGcgcgcg                                                                     >  1:555666/1‑71 (MQ=255)
 gcgcTTTCGGTCGCGGTAGGCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCgcgg                                                                    >  1:649140/1‑71 (MQ=255)
   gcTTTCGGTCGCGGTAGGCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCgcggc                                                                   >  1:2194141/1‑70 (MQ=255)
       tCGGTCGCGGTAGGCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAg                                                               >  1:332666/1‑70 (MQ=255)
       tCGGTCGCGGTAGGCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAg                                                               >  1:1333299/1‑70 (MQ=255)
       tCGGTCGCGGTAGGCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAg                                                               >  1:1882902/1‑70 (MQ=255)
         ggTCGCGGTAGGCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTaa                                                            <  1:3684156/71‑1 (MQ=255)
         ggTCGCGGTAGGCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTa                                                             <  1:3500744/70‑1 (MQ=255)
         ggTCGCGGTAGGCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTa                                                             <  1:4661585/70‑1 (MQ=255)
           tCGCGGTAGGCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAg                                                          <  1:235364/71‑1 (MQ=255)
           tCGCGGTAGGCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAg                                                          <  1:240826/71‑1 (MQ=255)
               ggTAGGCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCgcgg                                                                    <  1:303256/57‑1 (MQ=255)
                  aGGCCTAATGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCg                                                    >  1:3599319/1‑70 (MQ=255)
                  aGGCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGAGCGCGGCGCAGTAAAGCTGTCg                                                    >  1:4120909/1‑70 (MQ=255)
                   ggCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCa                                                                >  1:2641695/1‑57 (MQ=255)
                   ggCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGa                                                   <  1:2850571/70‑1 (MQ=255)
                   ggCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGa                                                   >  1:4252569/1‑70 (MQ=255)
                   ggCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGa                                                   <  1:3307883/70‑1 (MQ=255)
                   ggCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGATGTCGa                                                   <  1:4219598/70‑1 (MQ=255)
                    gCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAg                                                          <  1:644920/62‑1 (MQ=255)
                    gCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAg                                                          <  1:4358515/62‑1 (MQ=255)
                    gCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAg                                                          <  1:2057299/62‑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:1039953/70‑1 (MQ=255)
                         aCGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGa                                             >  1:1544116/1‑70 (MQ=255)
                         aCGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGa                                             >  1:4326121/1‑70 (MQ=255)
                         aCGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGa                                             <  1:1009311/70‑1 (MQ=255)
                            ccGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGt                                         >  1:2545152/1‑71 (MQ=255)
                            ccGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGt                                         <  1:3578275/71‑1 (MQ=255)
                             cGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTg                                        <  1:2107132/71‑1 (MQ=255)
                             cGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTg                                        <  1:4428809/71‑1 (MQ=255)
                             cGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTg                                        <  1:3738283/71‑1 (MQ=255)
                             cGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTg                                        <  1:1216295/71‑1 (MQ=255)
                             cGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTg                                        <  1:3903336/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)
                                     ttGCCGATGATTGTTACCTCGACGCTGGCGCGCGTCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCa                                 <  1:671322/70‑1 (MQ=255)
                                     ttGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTa                                                             <  1:2889043/42‑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:4143644/1‑70 (MQ=255)
                                             gatTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTg                         >  1:1388070/1‑70 (MQ=255)
                                             gatTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTg                         >  1:160976/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:4111310/1‑50 (MQ=255)
                                               ttGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGaaa                                           >  1:4532606/1‑50 (MQ=255)
                                               ttGTTACCTCGACGCTGGCGCGCGGCCCAGTAAAGCTGTCGAAACAGAAAg                                          <  1:2578572/51‑1 (MQ=255)
                                                tGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACAt                            >  1:4597679/1‑64 (MQ=255)
                                                 gTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGaa                                            >  1:1997277/1‑47 (MQ=255)
                                                 gTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGaa                                            >  1:4561789/1‑47 (MQ=255)
                                                 gTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGaa                                            >  1:2636248/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: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:1623193/70‑1 (MQ=255)
                                                                gcgcgcGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTg      <  1:3789960/70‑1 (MQ=255)
                                                                   cgcgGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGgcgc  <  1:1763110/71‑1 (MQ=255)
                                                                   cgcgGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGgcgc  <  1:4241744/71‑1 (MQ=255)
                                                                   cgcgGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGgcgc  <  1:4517149/71‑1 (MQ=255)
                                                                    gcggcgCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTgg                        >  1:1951002/1‑48 (MQ=255)
                                                                     |                                                                    
TGCGCTTTCGGTCGCGGTAGGCCTAACGCCGGAAATGTTGCCGATGATTGTTACCTCGACGCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGC  >  minE/2909968‑2910105

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: