Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,816,542 G→C 7.0% G67A (GGC→GCC)  soxR → DNA‑binding transcriptional dual regulator, Fe‑S center for redox‑sensing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,816,5420GC 7.0% 410.2 / 15.6 129G67A (GGC→GCC) soxRDNA‑binding transcriptional dual regulator, Fe‑S center for redox‑sensing
Reads supporting (aligned to +/- strand):  ref base G (57/63);  new base C (7/2);  total (64/65)
Fisher's exact test for biased strand distribution p-value = 9.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 132 total aligned reads displayed.

GGCAATCAGCGGCGATATAAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAGGGCATACGTTAAGTGCGAAAG  >  minE/2816472‑2816610
                                                                      |                                                                    
ggCAATCAGCGGCGATATAAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTgg                                                                      <  1:2250602/71‑1 (MQ=255)
  cAATCAGCGGCGATATAAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCa                                                                    >  1:1855299/1‑71 (MQ=255)
    aTCAGCGGCGATATAAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCAtt                                                                  <  1:3404956/71‑1 (MQ=255)
    aTCAGCGGCGATATAAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCAtt                                                                  <  1:2294321/71‑1 (MQ=255)
    aTCAGCGGCGATATAAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCAtt                                                                  <  1:95671/71‑1 (MQ=255)
     tCAGCGGCGATATAAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTc                                                                 >  1:515911/1‑71 (MQ=255)
     tCAGCGGCGATATAAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTc                                                                 >  1:4184713/1‑71 (MQ=255)
        gcggcgATATAAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGc                                                              <  1:3182895/71‑1 (MQ=255)
         cggcgATATAAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCt                                                             <  1:3962990/71‑1 (MQ=255)
         cggcgATATAAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCt                                                             <  1:1739344/71‑1 (MQ=255)
           gcgATATAAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTgg                                                           >  1:993060/1‑71 (MQ=255)
           gcgATATAAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTgg                                                           >  1:3038447/1‑71 (MQ=255)
            cgATATAAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGGAATGCCAATACGCt                                                             <  1:4290721/68‑1 (MQ=37)
            cgATATAAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGGAATGCCAATACGCt                                                             <  1:2426283/68‑1 (MQ=37)
              atatAAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGGAATGCCAATACGCTga                                                           >  1:3644334/1‑67 (MQ=37)
              atatAAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGGAATGCCAATACGCTga                                                           >  1:3387579/1‑67 (MQ=37)
              atatAAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGGAATGCCAATACGCTga                                                           >  1:1348754/1‑67 (MQ=37)
              atatAAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGGAATGCCAATACGCTga                                                           >  1:1452552/1‑67 (MQ=37)
              atatAAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGGAATGCCAATACGCTga                                                           >  1:2180643/1‑67 (MQ=37)
              atatAAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGGAATGCCAATACGCTga                                                           >  1:770757/1‑67 (MQ=37)
               tataAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGAc                                                       <  1:30250/71‑1 (MQ=255)
               tataAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGAc                                                       <  1:1445250/71‑1 (MQ=255)
                ataAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGAc                                                       >  1:2924932/1‑70 (MQ=255)
                ataAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGAc                                                       >  1:4199593/1‑70 (MQ=255)
                ataAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGAc                                                       >  1:2902641/1‑70 (MQ=255)
                 taAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGAc                                                       <  1:188270/69‑1 (MQ=255)
                 taAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGAc                                                       <  1:2869539/69‑1 (MQ=255)
                 taAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGAc                                                       <  1:3927895/69‑1 (MQ=255)
                     cGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTgg                                                 <  1:1481678/71‑1 (MQ=255)
                     cGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTgg                                                 >  1:3490697/1‑71 (MQ=255)
                       tGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGAc                                                       <  1:3361642/63‑1 (MQ=255)
                       tGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTg                                               <  1:3111910/71‑1 (MQ=255)
                       tGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTg                                               <  1:901353/71‑1 (MQ=255)
                       tGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTg                                               <  1:4675041/71‑1 (MQ=255)
                       tGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTg                                               <  1:4437178/71‑1 (MQ=255)
                         aTGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCa                                                     >  1:3543787/1‑63 (MQ=255)
                         aTGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCa                                                     >  1:4597438/1‑63 (MQ=255)
                         aTGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCa                                                     >  1:462443/1‑63 (MQ=255)
                            tgtTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTgg                                                 >  1:2597219/1‑64 (MQ=255)
                            tgtTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGc                                           <  1:3903629/70‑1 (MQ=255)
                            tgtTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGc                                           <  1:1096635/70‑1 (MQ=255)
                            tgtTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCg                                          <  1:4529665/71‑1 (MQ=255)
                                gCGATATGTTGCAATTATCAAAGTTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGttt                                       <  1:1674947/70‑1 (MQ=255)
                                gCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGttt                                       <  1:3252420/70‑1 (MQ=255)
                                gCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGttt                                       <  1:3247130/70‑1 (MQ=255)
                                gCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGttt                                       <  1:1138087/70‑1 (MQ=255)
                                gCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGttt                                       <  1:3718858/70‑1 (MQ=255)
                                gCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTg                                      >  1:2904828/1‑71 (MQ=255)
                                gCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTg                                      >  1:2437854/1‑71 (MQ=255)
                                  gATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCt                                                             <  1:36939/46‑1 (MQ=255)
                                  gATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCt                                                             <  1:1673822/46‑1 (MQ=255)
                                  gATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCt                                                             <  1:2866789/46‑1 (MQ=255)
                                  gATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCt                                                             <  1:4285634/46‑1 (MQ=255)
                                    tatGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCgt                                  <  1:1902572/71‑1 (MQ=255)
                                     atGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGc                                                              >  1:2611073/1‑42 (MQ=255)
                                     atGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGc                                                              >  1:2476349/1‑42 (MQ=255)
                                     atGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGc                                                              >  1:4180570/1‑42 (MQ=255)
                                       gTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGtt                               >  1:922133/1‑71 (MQ=255)
                                        ttGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGtt                               <  1:2802024/70‑1 (MQ=255)
                                         tGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGaa                                             >  1:571953/1‑55 (MQ=255)
                                         tGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGaa                                             >  1:2567705/1‑55 (MQ=255)
                                         tGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGttt                                       <  1:652564/61‑1 (MQ=255)
                                           cAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCt                                                             <  1:3544557/37‑1 (MQ=255)
                                            aaTTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTgg                                     <  1:3089782/60‑1 (MQ=255)
                                             aTTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGtt                                        <  1:2952674/56‑1 (MQ=255)
                                              ttATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCg                                          <  1:2054037/53‑1 (MQ=255)
                                              ttATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGccc                           >  1:500093/1‑68 (MQ=255)
                                              ttATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCg                          >  1:4112973/1‑69 (MQ=255)
                                                aTCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGt                                                >  1:333237/1‑45 (MQ=255)
                                                aTCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGt                                                >  1:2702532/1‑45 (MQ=255)
                                                aTCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGt                                                >  1:1854805/1‑45 (MQ=255)
                                                 tCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAgg                      <  1:3315565/70‑1 (MQ=255)
                                                 tCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAgg                      <  1:1572869/70‑1 (MQ=255)
                                                 tCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAgg                      <  1:4009570/70‑1 (MQ=255)
                                                   aaaaTTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAGGGCa                   >  1:2778356/1‑71 (MQ=255)
                                                   aaaaTTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAGGGCa                   >  1:1431120/1‑71 (MQ=255)
                                                   aaaaTTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAGGGCa                   >  1:222650/1‑71 (MQ=255)
                                                   aaaaTTGCTCAGCGTATTGGCATTCCGCTGGCGACCAATGGTGAAGc                                           >  1:2256273/1‑47 (MQ=255)
                                                       ttGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAGGGCATACg               >  1:3274987/1‑71 (MQ=255)
                                                         gCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAGGGCATACGtt             <  1:3019254/71‑1 (MQ=255)
                                                          cTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAGGGCATACGtt             >  1:776002/1‑70 (MQ=255)
                                                          cTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAGGGCATACGt              <  1:1315135/69‑1 (MQ=255)
                                                                gTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAGGGCATACGTTAAGTGc       >  1:4647424/1‑70 (MQ=255)
                                                                 tATTGGCATTCCGCTGGCGACCATTGGTGAAGCGt                                         <  1:1253436/35‑1 (MQ=255)
                                                                 tATTGGCATTCCGCTGGCGACCATTGGTGAAGCGt                                         <  1:3417716/35‑1 (MQ=255)
                                                                 tATTGGCATTCCGCTGGCGACCATTGGTGAAGCGt                                         <  1:913273/35‑1 (MQ=255)
                                                                 tATTGGCATTCCGCTGGCGACCATTGGTGAAGCGt                                         <  1:195150/35‑1 (MQ=255)
                                                                   ttGGCATTCCGCTGGCGACCATTGGTGAAGCGttt                                       >  1:379828/1‑35 (MQ=255)
                                                                   ttGGCATTCCGCTGGCGACCATTGGTGAAGCGttt                                       >  1:3286082/1‑35 (MQ=255)
                                                                   ttGGCATTCCGCTGGCGACCATTGGTGAAGCGttt                                       >  1:2170362/1‑35 (MQ=255)
                                                                   ttGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAGGGCATACGTTAAGTGCGaaa   <  1:4388410/71‑1 (MQ=255)
                                                                   ttGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAGGGCATACGTTAAGTGCGaaa   <  1:1670615/71‑1 (MQ=255)
                                                                   ttGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAGGGCATACGTTAAGTGCGa     >  1:3171896/1‑69 (MQ=255)
                                                                   ttGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAGGGCATACGTTAAGTGCGa     >  1:9576/1‑69 (MQ=255)
                                                                   ttGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAGGGCATACGTTAAGTGCGa     >  1:2885736/1‑69 (MQ=255)
                                                                   ttGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAGGGCATACGATAAGTGCGa     >  1:386266/1‑69 (MQ=255)
                                                                   ttGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCACGAAGGGCATACGTTAAGTGCGa     >  1:1538908/1‑69 (MQ=255)
                                                                   ttGGCATTCAGCTGGCGACCATTGGTGAAGCGttt                                       >  1:3190322/1‑35 (MQ=255)
                                                                    tGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAGGGCATACGTTAAGTGCGa     >  1:13506/1‑68 (MQ=255)
                                                                     ggCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAGGGCATACGTTAAGTGCGAAag  <  1:1120032/70‑1 (MQ=255)
                                                                      |                                                                    
GGCAATCAGCGGCGATATAAACGTGATGTGTTGCGATATGTTGCAATTATCAAAATTGCTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAGGGCATACGTTAAGTGCGAAAG  >  minE/2816472‑2816610

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: