Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,461,675 T→C 57.0% N305N (AAT→AAC gidA → glucose‑inhibited cell‑division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,461,6750TC57.0% 47.1 / 104.1 137N305N (AAT→AACgidAglucose‑inhibited cell‑division protein
Reads supporting (aligned to +/- strand):  ref base T (40/19);  new base C (10/68);  total (50/87)
Fisher's exact test for biased strand distribution p-value = 2.64e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

Only 100 of 138 total aligned reads displayed.

CGCTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTC  >  minE/2461616‑2461740
                                                           |                                                                 
cgcTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTAt                                                          >  1:2796926/1‑69 (MQ=255)
cgcTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTAt                                                          >  1:1920076/1‑69 (MQ=255)
  cTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATc                                                         >  1:1128729/1‑68 (MQ=255)
  cTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATc                                                         >  1:870113/1‑68 (MQ=255)
  cTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATc                                                         >  1:2449035/1‑68 (MQ=255)
  cTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATc                                                         >  1:2704592/1‑68 (MQ=255)
  cTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATc                                                         >  1:849215/1‑68 (MQ=255)
  cTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATc                                                         >  1:777639/1‑68 (MQ=255)
    tCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGa                                                      >  1:145319/1‑69 (MQ=255)
    tCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGa                                                      >  1:3409260/1‑69 (MQ=255)
     cGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGaa                                                     >  1:513842/1‑69 (MQ=255)
     cGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGaa                                                     >  1:1935978/1‑69 (MQ=255)
     cGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGaa                                                     >  1:857918/1‑69 (MQ=255)
     cGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGaa                                                     >  1:1062831/1‑69 (MQ=255)
     cGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGaa                                                     >  1:1003799/1‑69 (MQ=255)
     cGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACATCTAACGAAATTTATCCGaa                                                     >  1:877584/1‑69 (MQ=255)
        cGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGAACg                                                   <  1:311100/68‑1 (MQ=255)
        cGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGAACg                                                   <  1:3138459/68‑1 (MQ=255)
        cGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGAACg                                                   <  1:314786/68‑1 (MQ=255)
        cGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGAACg                                                   <  1:2870945/68‑1 (MQ=255)
        cGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGAACg                                                   <  1:2847136/68‑1 (MQ=255)
        cGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGAACg                                                   <  1:1296865/68‑1 (MQ=255)
        cGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGAACg                                                   <  1:1064483/68‑1 (MQ=255)
        cGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGAACg                                                   <  1:554866/68‑1 (MQ=255)
        cGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGAACg                                                   <  1:771337/68‑1 (MQ=255)
         gACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGAACg                                                   <  1:407897/67‑1 (MQ=255)
         gACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGt                                                 <  1:3093624/69‑1 (MQ=255)
          aCAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTa                                                <  1:43411/69‑1 (MQ=255)
                         aGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTAtctc                                            >  1:261406/1‑58 (MQ=255)
                         aGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTAtctc                                            >  1:2369010/1‑58 (MQ=255)
                         aGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTAtctc                                            >  1:3189238/1‑58 (MQ=255)
                         aGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTAtctc                                            >  1:3328173/1‑58 (MQ=255)
                         aGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTAtctc                                            >  1:3419445/1‑58 (MQ=255)
                         aGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTAtctc                                            >  1:2790551/1‑58 (MQ=255)
                         aGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTAtc                                              >  1:3293823/1‑56 (MQ=255)
                          gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                                               <  1:3037914/38‑1 (MQ=38)
                          gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                                               <  1:2508339/38‑1 (MQ=38)
                          gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                                               <  1:2219123/38‑1 (MQ=38)
                          gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                                               <  1:2727078/38‑1 (MQ=38)
                          gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                                               <  1:2887380/38‑1 (MQ=38)
                          gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                                               <  1:2887348/38‑1 (MQ=38)
                          gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                                               <  1:2866280/38‑1 (MQ=38)
                          gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                                               <  1:2308313/38‑1 (MQ=38)
                          gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                                               <  1:652263/38‑1 (MQ=38)
                          gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                                               <  1:1374013/38‑1 (MQ=38)
                          gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                                               <  1:1420104/38‑1 (MQ=38)
                          gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                                               <  1:1420283/38‑1 (MQ=38)
                          gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                                               <  1:1902937/38‑1 (MQ=38)
                          gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                                               <  1:1884817/38‑1 (MQ=38)
                          gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                                               <  1:1706440/38‑1 (MQ=38)
                          gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                                               <  1:168197/38‑1 (MQ=38)
                          gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                                               <  1:960816/38‑1 (MQ=38)
                             cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGATCGGTATCTCCACCAGCCTGCCGt                              >  1:1643278/1‑68 (MQ=255)
                             cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTAtc                                              >  1:2943073/1‑52 (MQ=255)
                             cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCCGCCTGCCGt                              >  1:1218204/1‑68 (MQ=255)
                             cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCt                                   >  1:1489215/1‑63 (MQ=255)
                             cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGt                              >  1:466382/1‑68 (MQ=255)
                             cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGt                              >  1:470518/1‑68 (MQ=255)
                             cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGt                              >  1:1004410/1‑68 (MQ=255)
                             cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGt                              >  1:2753029/1‑68 (MQ=255)
                             cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGt                              >  1:113703/1‑68 (MQ=255)
                             cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGt                              >  1:2398750/1‑68 (MQ=255)
                             cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGt                              >  1:2294326/1‑68 (MQ=255)
                             cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGt                              >  1:146878/1‑68 (MQ=255)
                                       ccGGAAGGGCTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAg                   <  1:1045297/69‑1 (MQ=255)
                                       ccGGAAGGGCTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAg                   <  1:1067252/69‑1 (MQ=255)
                                       ccGGAAGGGCTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAg                   <  1:3383815/69‑1 (MQ=255)
                                       ccGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCa                    <  1:2319861/68‑1 (MQ=255)
                                        cGGAAGGGCTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAg                   <  1:1911886/68‑1 (MQ=255)
                                        cGGAAGGGCTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAg                   <  1:3274176/68‑1 (MQ=255)
                                              ggcTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTTCAGATGCaaa            <  1:1617145/67‑1 (MQ=255)
                                              ggcTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATTCaaa            <  1:2098888/67‑1 (MQ=255)
                                              ggcTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCaaa            <  1:1653421/67‑1 (MQ=255)
                                              ggcTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCaaa            <  1:2756091/67‑1 (MQ=255)
                                              ggcTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCaaa            <  1:454294/67‑1 (MQ=255)
                                              ggcTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCaaa            <  1:459561/67‑1 (MQ=255)
                                              ggcTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCaaa            <  1:1454396/67‑1 (MQ=255)
                                              ggcTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCaaa            <  1:1327405/67‑1 (MQ=255)
                                              ggcTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCaaa            <  1:1256793/67‑1 (MQ=255)
                                              ggcTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCaaa            <  1:8500/67‑1 (MQ=255)
                                              ggcTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCaaa            <  1:1168427/67‑1 (MQ=255)
                                              ggcTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCaaa            <  1:3317589/67‑1 (MQ=255)
                                              ggcTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCaaa            <  1:2124274/67‑1 (MQ=255)
                                              ggcTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCaaa            <  1:3292455/67‑1 (MQ=255)
                                              ggcTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCaaa            <  1:2746966/67‑1 (MQ=255)
                                              ggcTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCaaa            <  1:222870/67‑1 (MQ=255)
                                              ggcTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCaaa            <  1:1839684/67‑1 (MQ=255)
                                              ggcTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCaaa            <  1:2544780/67‑1 (MQ=255)
                                                cTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAg                   >  1:2779356/1‑60 (MQ=255)
                                                 tGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAg                   <  1:1516176/59‑1 (MQ=255)
                                                     ctctAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTc       <  1:70731/67‑1 (MQ=255)
                                                     ctctAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTc       <  1:3026910/67‑1 (MQ=255)
                                                     ctctAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTc       <  1:1338766/67‑1 (MQ=255)
                                                     ctctAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTc       <  1:3306857/67‑1 (MQ=255)
                                                     ctctAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCg     <  1:1296113/69‑1 (MQ=255)
                                                     ctctAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCg     <  1:2979650/69‑1 (MQ=255)
                                                      tctAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCg     <  1:2610912/68‑1 (MQ=255)
                                                          aTGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGctc  <  1:1292337/67‑1 (MQ=255)
                                                          aTGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGctc  <  1:2158749/67‑1 (MQ=255)
                                                          aTGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGctc  <  1:982116/67‑1 (MQ=255)
                                                           |                                                                 
CGCTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTC  >  minE/2461616‑2461740

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: