Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,926,225 G→A 65.6% G56G (GGG→GGA lysR → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,926,2250GA65.6% 62.2 / 51.5 90G56G (GGG→GGAlysRDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base G (11/20);  new base A (38/21);  total (49/41)
Fisher's exact test for biased strand distribution p-value = 1.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

GTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGG  >  minE/1926160‑1926282
                                                                 |                                                         
gTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCg                                                         >  1:1056198/1‑68 (MQ=255)
gTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCg                                                         >  1:1080474/1‑68 (MQ=255)
gTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCg                                                         >  1:3310460/1‑68 (MQ=255)
gTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCg                                                         >  1:2054584/1‑68 (MQ=255)
gTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCg                                                         >  1:1561928/1‑68 (MQ=255)
  cAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGAt                                                       >  1:391255/1‑68 (MQ=255)
  cAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGAt                                                       >  1:1658232/1‑68 (MQ=255)
  cAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGAt                                                       >  1:885075/1‑68 (MQ=255)
     ccGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTAc                                                    <  1:335243/68‑1 (MQ=255)
     ccGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACa                                                   <  1:77723/69‑1 (MQ=255)
      cgcgAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACAt                                                  <  1:987488/69‑1 (MQ=255)
      cgcgAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACAt                                                  >  1:1216001/1‑69 (MQ=255)
       gcgAACTTGCGCGCTTTGAGAAGGTGATCTGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACAt                                                  <  1:3080090/68‑1 (MQ=255)
       gcgAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACAt                                                  <  1:2274830/68‑1 (MQ=255)
             tggcgctcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGAtt                                                      <  1:2456249/56‑1 (MQ=255)
             tggcgcgcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGAtt                                                      <  1:702964/56‑1 (MQ=255)
             tggcgcgcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGAtt                                                      <  1:513153/56‑1 (MQ=255)
             tggcgcgcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGAtt                                                      <  1:448398/56‑1 (MQ=255)
             tggcgcgcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGAtt                                                      <  1:2167019/56‑1 (MQ=255)
             tggcgcgcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGAtt                                                      <  1:69235/56‑1 (MQ=255)
             tggcgcgcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGAtt                                                      <  1:945830/56‑1 (MQ=255)
             tggcgcgcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGAtt                                                      <  1:700926/56‑1 (MQ=255)
             tggcgcgcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGAtt                                                      <  1:1440642/56‑1 (MQ=255)
             tggcgcgcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGAtt                                                      <  1:2527831/56‑1 (MQ=255)
             tggcgcgcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGAtt                                                      <  1:2532581/56‑1 (MQ=255)
             tggcgcgcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGAtt                                                      <  1:1734373/56‑1 (MQ=255)
             tggcgcgcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGAtt                                                      <  1:1705307/56‑1 (MQ=255)
             tggcgcgcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGAtt                                                      <  1:1142609/56‑1 (MQ=255)
             tggcgcgcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGAtt                                                      <  1:1452194/56‑1 (MQ=255)
             tggcgcgcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGAtt                                                      <  1:280233/56‑1 (MQ=255)
             tggcgcgcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGAtt                                                      <  1:3009250/56‑1 (MQ=255)
             tggcgcgcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGAt                                                       <  1:2717935/55‑1 (MQ=255)
             tggcgcgcTTTGAGAAGGTGATCGGGCTGAAATTGTTGAGCGCATACGTGGACGAtt                                                      <  1:3076157/55‑1 (MQ=255)
              ggcgcgcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGAtt                                                      <  1:642743/56‑1 (MQ=255)
                   gcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAg                                     <  1:1508491/69‑1 (MQ=255)
                   gcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAg                                     <  1:373279/69‑1 (MQ=255)
                   gcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAg                                     <  1:1317028/69‑1 (MQ=255)
                   gcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAg                                     <  1:1497441/69‑1 (MQ=255)
                    cTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAg                                     <  1:3022359/68‑1 (MQ=255)
                    cTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAg                                     <  1:1512357/68‑1 (MQ=255)
                    cTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAg                                     <  1:1636784/68‑1 (MQ=255)
                    cTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAg                                     <  1:2726996/68‑1 (MQ=255)
                      ttGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGAtt                                                      <  1:325286/49‑1 (MQ=39)
                       tGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAg                                     >  1:3264517/1‑65 (MQ=255)
                       tGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACt                                 >  1:2320828/1‑69 (MQ=255)
                       tGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACt                                 >  1:3145605/1‑69 (MQ=255)
                       tGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACt                                 >  1:3050184/1‑69 (MQ=255)
                       tGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACt                                 >  1:3048806/1‑69 (MQ=255)
                       tGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACt                                 >  1:149535/1‑69 (MQ=255)
                       tGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACt                                 >  1:3223679/1‑69 (MQ=255)
                       tGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACt                                 >  1:2775312/1‑69 (MQ=255)
                       tGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACt                                 >  1:3307109/1‑69 (MQ=255)
                       tGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACt                                 >  1:904886/1‑69 (MQ=255)
                       tGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACt                                 >  1:2608789/1‑69 (MQ=255)
                       tGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACt                                 >  1:1293555/1‑69 (MQ=255)
                       tGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACt                                 >  1:849064/1‑69 (MQ=255)
                       tGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACt                                 >  1:25324/1‑69 (MQ=255)
                       tGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACt                                 >  1:187868/1‑69 (MQ=255)
                       tGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACt                                 >  1:2524650/1‑69 (MQ=255)
                       tGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACt                                 >  1:2066875/1‑69 (MQ=255)
                       tGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACt                                 >  1:2498113/1‑69 (MQ=255)
                       tGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACt                                 >  1:617717/1‑69 (MQ=255)
                               tGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTg                          <  1:1062775/68‑1 (MQ=255)
                               tGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTg                          <  1:2746768/68‑1 (MQ=255)
                               tGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTg                          <  1:1782623/68‑1 (MQ=255)
                               tGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTg                          <  1:2305661/68‑1 (MQ=255)
                               tGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTg                          <  1:2346847/68‑1 (MQ=255)
                                    gggCTGGAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTga                     >  1:437693/1‑68 (MQ=255)
                                    gggCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGt                                          >  1:2466603/1‑47 (MQ=39)
                                    gggCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTga                     >  1:2417679/1‑68 (MQ=255)
                                    gggCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTga                     >  1:1121565/1‑68 (MQ=255)
                                    gggCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTga                     >  1:1401018/1‑68 (MQ=255)
                                    gggCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTga                     >  1:1444412/1‑68 (MQ=255)
                                    gggCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTga                     >  1:155363/1‑68 (MQ=255)
                                    gggCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTga                     >  1:1889476/1‑68 (MQ=255)
                                    gggCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTga                     >  1:735632/1‑68 (MQ=255)
                                    gggCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTga                     >  1:1970789/1‑68 (MQ=255)
                                    gggCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTga                     >  1:2017561/1‑68 (MQ=255)
                                    gggCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTga                     >  1:3066329/1‑68 (MQ=255)
                                    gggCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTga                     >  1:250551/1‑68 (MQ=255)
                                    gggCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTga                     >  1:3317154/1‑68 (MQ=255)
                                    gggCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTga                     >  1:2630477/1‑68 (MQ=255)
                                    gggCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTga                     >  1:2706575/1‑68 (MQ=255)
                                    gggCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTga                     >  1:2728534/1‑68 (MQ=255)
                                    gggCTGAAATTGTTTGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTga                     >  1:3180815/1‑68 (MQ=255)
                                    gggCTGAAATTGTATGAGCGCATACGTGGACGATTACATCCTACCGTGCAAGGACTGCGTCTg                          >  1:626387/1‑63 (MQ=255)
                                               gTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACg         >  1:210690/1‑69 (MQ=255)
                                               gTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACg         >  1:2473669/1‑69 (MQ=255)
                                                      cgcgTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTgg  <  1:2583304/69‑1 (MQ=255)
                                                      cgcgTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTgg  <  1:2789735/69‑1 (MQ=255)
                                                                 |                                                         
GTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGG  >  minE/1926160‑1926282

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: