Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,960,737 C→G 15.5% intergenic (+187/+233) yjjV → / ← yjjW predicted DNase/predicted pyruvate formate lyase activating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,960,7370CG15.5% 171.2 / 19.5 66intergenic (+187/+233)yjjV/yjjWpredicted DNase/predicted pyruvate formate lyase activating enzyme
Reads supporting (aligned to +/- strand):  ref base C (25/29);  new base G (2/8);  total (27/39)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

CACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGC  >  minE/2960667‑2960805
                                                                      |                                                                    
cACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTc                                                                      >  1:2578012/1‑71 (MQ=255)
 aCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTc                                                                      <  1:1839858/70‑1 (MQ=255)
 aCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCt                                                                     <  1:3378591/71‑1 (MQ=40)
 aCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCt                                                                     <  1:3394913/71‑1 (MQ=40)
 aCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCt                                                                     <  1:3399204/71‑1 (MQ=40)
       cATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGcc                                                                >  1:3626803/1‑70 (MQ=40)
       cATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGcc                                                                >  1:4697871/1‑70 (MQ=40)
        aTCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGa                                                              <  1:1469508/71‑1 (MQ=39)
         tCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTc                                                                      <  1:4856407/62‑1 (MQ=40)
          ccGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAAc                                                            >  1:4984778/1‑71 (MQ=39)
           cGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCtg                                                                    >  1:1669761/1‑62 (MQ=40)
           cGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAAcc                                                           >  1:909405/1‑71 (MQ=39)
              cAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCg                                                          <  1:3018763/69‑1 (MQ=40)
                 ttGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGc                                                     >  1:3776532/1‑71 (MQ=39)
                  tGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTa                                                        >  1:1587003/1‑67 (MQ=255)
                  tGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGcc                                                    >  1:5101441/1‑71 (MQ=39)
                   ggCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCg                                                   >  1:1665957/1‑71 (MQ=39)
                    gcgcACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGa                                                              >  1:2958987/1‑59 (MQ=255)
                       atcaATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAg                                             <  1:1476457/69‑1 (MQ=0)
                        tcaATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATa                                               <  1:2266308/67‑1 (MQ=11)
                        tcaATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATa                                               <  1:2289614/67‑1 (MQ=11)
                        tcaATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAgg                                            <  1:2255827/70‑1 (MQ=2)
                        tcaATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAgg                                            <  1:3957540/70‑1 (MQ=2)
                          aaTGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGca                                          >  1:2384491/1‑70 (MQ=2)
                          aaTGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAg                                             >  1:1994611/1‑68 (MQ=11)
                           aTGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGc                                           <  1:3790723/71‑1 (MQ=39)
                           aTGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGAGCCGaa                                                             >  1:4569012/1‑53 (MQ=255)
                           aTGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGcac                                         <  1:4923885/71‑3 (MQ=2)
                            tGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGcact                                        <  1:5339936/71‑4 (MQ=0)
                              ccTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTcg                                      <  1:1454103/71‑2 (MQ=2)
                               cTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTc                                       <  1:2676927/71‑1 (MQ=39)
                               cTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTcgc                                     <  1:3416544/71‑3 (MQ=2)
                                  aTGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAAc                                                            <  1:2326038/47‑1 (MQ=255)
                                  aTGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAc                                     <  1:963089/70‑1 (MQ=255)
                                     cGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGaa                                                             >  1:2202971/1‑43 (MQ=255)
                                     cGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAAc                                                            >  1:724543/1‑44 (MQ=255)
                                         gCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGtt                                        <  1:118522/60‑1 (MQ=255)
                                         gCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATc                             >  1:1882501/1‑71 (MQ=255)
                                          cTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGtt                                        <  1:2036814/59‑1 (MQ=255)
                                          cTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATcc                            <  1:1211336/71‑1 (MQ=255)
                                          cTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATcc                            <  1:2678661/71‑1 (MQ=255)
                                           ttGTCGCGTCTTATCATGCCTACAAGTCTGTGCCg                                                               <  1:4764514/35‑1 (MQ=255)
                                            tGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCa                                      >  1:2900569/1‑59 (MQ=255)
                                            tGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATa                                               <  1:3944611/49‑1 (MQ=11)
                                              tCGCGTCTTATCATGCCTACAAGTCTGTGCCGAAc                                                            <  1:2630404/35‑1 (MQ=255)
                                              tCGCGTCTTATCATGCCTACAAGTCTGTGCCGAAc                                                            <  1:1419202/35‑1 (MQ=255)
                                              tCGCGTCTTATCATGCCTACAAGTCTGTGCCGAAc                                                            <  1:3910395/35‑1 (MQ=255)
                                              tCGCGTCTTATCATGCCTACAAGTCTGTGCCGAAc                                                            <  1:3775098/35‑1 (MQ=255)
                                              tCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCgg                                                  <  1:844988/45‑1 (MQ=255)
                                              tCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCgg                                                  <  1:2710772/45‑1 (MQ=255)
                                              tCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCgg                                                  <  1:352916/45‑1 (MQ=255)
                                              tCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCgg                                                  <  1:157955/45‑1 (MQ=255)
                                                gcgTCTTATCATGCCTACAAGTCTGTGCCGAACCg                                                          >  1:3963214/1‑35 (MQ=40)
                                                gcgTCTTATCATGCCTACAAGTCTGTGCCGAACCg                                                          >  1:1903392/1‑35 (MQ=40)
                                                  gTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCa                                      >  1:5499396/1‑53 (MQ=255)
                                                     ttATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGa                                                 >  1:4321922/1‑39 (MQ=255)
                                                     ttATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGa                                                 >  1:801304/1‑39 (MQ=255)
                                                        tCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCAt              <  1:220514/71‑1 (MQ=255)
                                                          aTGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATa             <  1:941300/70‑1 (MQ=255)
                                                          aTGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATa             <  1:2888511/70‑1 (MQ=255)
                                                          aTGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATa             <  1:1368210/70‑1 (MQ=255)
                                                            gCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATg          >  1:4803741/1‑71 (MQ=255)
                                                             ccTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATg          <  1:3893286/70‑1 (MQ=255)
                                                                 cAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATg          >  1:4228968/1‑66 (MQ=255)
                                                                 cAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTg      <  1:1433105/70‑1 (MQ=255)
                                                                 cAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGa     >  1:899163/1‑71 (MQ=255)
                                                                 cAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGa     >  1:414529/1‑71 (MQ=255)
                                                                    gTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGc  >  1:4845065/1‑71 (MQ=255)
                                                                     tCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgccgc                                <  1:3940811/40‑1 (MQ=37)
                                                                      |                                                                    
CACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGC  >  minE/2960667‑2960805

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: