Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,128,783 G→A 6.0% L83L (TTG→TTA)  yhbW → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,128,7830GA 6.0% 304.4 / 5.6 100L83L (TTG→TTA) yhbWhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (37/57);  new base A (4/2);  total (41/59)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

GGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAAT  >  minE/2128717‑2128851
                                                                  |                                                                    
ggTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCt                                                                   <  1:2933546/70‑1 (MQ=255)
ggTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCt                                                                   <  1:5343502/70‑1 (MQ=255)
ggTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCt                                                                   <  1:471531/70‑1 (MQ=255)
ggTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCt                                                                   >  1:4004748/1‑70 (MQ=255)
ggTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGGCTGGCGGCGTGATGTTGCCTa                                                                  <  1:3758127/71‑1 (MQ=255)
 gTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTaa                                                                 >  1:3096037/1‑71 (MQ=255)
    ttGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCa                                                              >  1:4240667/1‑71 (MQ=255)
     tGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCAc                                                             >  1:1044405/1‑71 (MQ=255)
        tCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAAc                                                                <  1:199427/65‑1 (MQ=255)
        tCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAAc                                                                <  1:150958/65‑1 (MQ=255)
        tCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCa                                                          >  1:3250422/1‑71 (MQ=255)
         cGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCAc                                                         >  1:2852645/1‑71 (MQ=255)
         cGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCAc                                                         >  1:1542697/1‑71 (MQ=255)
          ggCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCAc                                                         <  1:3175396/70‑1 (MQ=255)
           gCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCg                                                       >  1:4809070/1‑71 (MQ=255)
           gCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCg                                                       >  1:5033214/1‑71 (MQ=255)
              aTCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTg                                                    >  1:2881571/1‑71 (MQ=255)
              aTCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTg                                                    >  1:2887994/1‑71 (MQ=255)
                 tGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGt                                                  >  1:5604197/1‑70 (MQ=255)
                  ggcggcGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCa                                                >  1:757079/1‑71 (MQ=255)
                  ggcggcGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCa                                                >  1:1855188/1‑71 (MQ=255)
                  ggcggcGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCa                                                >  1:3135834/1‑71 (MQ=255)
                  ggcggcGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCa                                                >  1:2269511/1‑71 (MQ=255)
                   gcggcgAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTTGTc                                                 <  1:2030040/69‑1 (MQ=255)
                   gcggcgAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGt                                                  <  1:1831747/68‑1 (MQ=255)
                   gcggcaAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTACCAAACCACTCACCGTTGGtg                                                 <  1:2111076/69‑2 (MQ=255)
                   gcggcaAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTACCAAACCACTCACCGTTGGtg                                                 <  1:1155315/69‑2 (MQ=255)
                    cggcgAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCAtt                                              <  1:3205555/71‑1 (MQ=255)
                     ggcgAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTg                                             >  1:898222/1‑71 (MQ=255)
                        gAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAg                                          <  1:2432677/71‑1 (MQ=255)
                        gAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAg                                          <  1:3612360/71‑1 (MQ=255)
                        gAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAg                                          <  1:3690774/71‑1 (MQ=255)
                         aaTACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAg                                          >  1:1195383/1‑70 (MQ=255)
                         aaTACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGa                                         <  1:5525017/71‑1 (MQ=255)
                          aTACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGt                                                  <  1:4987983/61‑1 (MQ=255)
                          aTACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGaa                                        >  1:1282772/1‑71 (MQ=255)
                           tACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGt                                                  <  1:4974268/60‑1 (MQ=255)
                               accacGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGtt                                   >  1:200472/1‑71 (MQ=255)
                                   cGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCgg                                >  1:3442286/1‑70 (MQ=255)
                                        cATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAg                                          <  1:4171005/55‑1 (MQ=255)
                                        cATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAg                                          <  1:2755949/55‑1 (MQ=255)
                                        cATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGc                           <  1:5077610/70‑1 (MQ=255)
                                        cATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGc                           <  1:2109652/70‑1 (MQ=255)
                                        cATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGCTCATTGCAg                                          <  1:2160332/55‑1 (MQ=255)
                                         aTCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCtt                         <  1:5173617/71‑1 (MQ=255)
                                         aTCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCtt                         <  1:4520205/71‑1 (MQ=255)
                                         aTCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCtt                         <  1:5563804/71‑1 (MQ=255)
                                         aTCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCtt                         <  1:1416187/71‑1 (MQ=255)
                                         aTCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCt                          >  1:725394/1‑70 (MQ=255)
                                          tCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTTGTCATTGCAg                                          <  1:820875/53‑1 (MQ=255)
                                          tCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAg                                          <  1:4622542/53‑1 (MQ=255)
                                          tCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAg                                          <  1:1805496/53‑1 (MQ=255)
                                          tCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAg                                          <  1:3186282/53‑1 (MQ=255)
                                          tCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAg                                          <  1:2927221/53‑1 (MQ=255)
                                            tGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAAc                                       >  1:2481454/1‑54 (MQ=255)
                                            tGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCtt                         <  1:4631648/68‑1 (MQ=255)
                                            tGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCtt                         <  1:4250587/68‑1 (MQ=255)
                                             ggggTCTGGCGGCGTGATGTTGCCTAACCACTCAcc                                                        <  1:1948425/36‑1 (MQ=255)
                                             ggggTCTGGCGGCGTGATGTTGCCTAACCACTCAcc                                                        <  1:3989837/36‑1 (MQ=255)
                                             ggggTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCt                          <  1:3158381/66‑1 (MQ=255)
                                              gggTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTaa                       <  1:4569475/68‑1 (MQ=255)
                                              gggTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTaa                       <  1:1666822/68‑1 (MQ=255)
                                              gggTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTaa                       <  1:751515/68‑1 (MQ=255)
                                              gggTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTaa                       <  1:3154073/68‑1 (MQ=255)
                                              gggTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATa                     >  1:4770232/1‑70 (MQ=255)
                                                 tCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGtt                                   <  1:2843948/53‑1 (MQ=255)
                                                 tCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGtt                                   <  1:526108/53‑1 (MQ=255)
                                                 tCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGtt                                   <  1:2297250/53‑1 (MQ=255)
                                                 tCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACAct                 <  1:347671/71‑1 (MQ=255)
                                                   tGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAg                                          <  1:5201139/44‑1 (MQ=255)
                                                   tGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAg                                          <  1:1966137/44‑1 (MQ=255)
                                                   tGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCa                              >  1:386594/1‑56 (MQ=255)
                                                   tGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACActc                >  1:4533210/1‑70 (MQ=255)
                                                   tGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACActc                >  1:1214574/1‑70 (MQ=255)
                                                   tGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGCTCGGCACGCtt                         <  1:3654481/61‑1 (MQ=255)
                                                     gcggcgTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAg                                     >  1:3809273/1‑47 (MQ=255)
                                                     gcggcgTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGtt                                   >  1:5268446/1‑49 (MQ=255)
                                                     gcggcgTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGt                                    >  1:5201349/1‑48 (MQ=255)
                                                     gcggcgTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACActct               <  1:1853371/69‑1 (MQ=255)
                                                     gcggcgTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTa              >  1:3924417/1‑70 (MQ=255)
                                                      cggcgTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATc            >  1:2422984/1‑71 (MQ=255)
                                                       ggcgTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATcc           <  1:4115837/71‑1 (MQ=255)
                                                        gcgTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATcc           <  1:293210/70‑1 (MQ=255)
                                                         cgTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCgg         >  1:2086487/1‑71 (MQ=255)
                                                          gTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCggg        >  1:302183/1‑71 (MQ=255)
                                                           tGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCgggg       >  1:5487209/1‑71 (MQ=255)
                                                            gATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGc      >  1:3159227/1‑71 (MQ=255)
                                                            gATGTTACCAAACCACTCACCGTTGGTGATTGCCGAACAGTTCGGCACGCTCAAt                      >  1:1823771/1‑55 (MQ=37)
                                                            gATGTTACCAAACCACTCACCGTTGGTGATTGCCGAACAGTTCGGCACGCTCAAt                      >  1:4435032/1‑55 (MQ=37)
                                                            gATGTTACCAAACCACTCACCGTTGGTGATTGCCGAACAGTTCGGCACGCTCAAt                      >  1:4944276/1‑55 (MQ=37)
                                                            gATGTTACCAAACCACTCACCGTTGGTGATTGCCGAACAGTTCGGCACGCTCAAt                      >  1:4578327/1‑55 (MQ=37)
                                                              tGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCtt                         <  1:1590835/50‑1 (MQ=255)
                                                              tGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCtt                         <  1:3812922/50‑1 (MQ=255)
                                                              tGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCtt                         <  1:2112398/50‑1 (MQ=255)
                                                              tGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCtt                         <  1:5153639/50‑1 (MQ=255)
                                                              tGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCtt                         <  1:2583454/50‑1 (MQ=255)
                                                                ttGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCtt                         <  1:327661/48‑1 (MQ=255)
                                                                ttGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCtt                         <  1:4672880/48‑1 (MQ=255)
                                                                ttGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATa                     >  1:4714484/1‑52 (MQ=255)
                                                                ttGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAAt  <  1:1097033/71‑1 (MQ=255)
                                                                  |                                                                    
GGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAAT  >  minE/2128717‑2128851

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