Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,950,841 T→C 7.5% L113S (TTA→TCA)  bglJ → DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,950,8410TC 7.5% 274.7 / 9.7 93L113S (TTA→TCA) bglJDNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (43/43);  new base C (2/5);  total (45/48)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

Only 100 of 105 total aligned reads displayed.

CGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAG  >  minE/2950775‑2950910
                                                                  |                                                                     
cGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAAt                                                                   <  1:3602326/71‑1 (MQ=255)
cGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAAt                                                                   <  1:2810455/71‑1 (MQ=255)
 gCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATg                                                                  >  1:2092794/1‑71 (MQ=255)
 gCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATg                                                                  >  1:468127/1‑71 (MQ=255)
 gCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATg                                                                  >  1:4627822/1‑71 (MQ=255)
 gCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATg                                                                  >  1:1670384/1‑71 (MQ=255)
 gCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATg                                                                  >  1:2963357/1‑71 (MQ=255)
 gCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATg                                                                  >  1:248060/1‑71 (MQ=255)
 gCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATg                                                                  >  1:2518709/1‑71 (MQ=255)
           gtATTAAGTAAAGCGTCAACGCTGGAGATTTTTCCTCAGGAACTCTTTTTGTCATTAAATGGTGTACGt                                                          >  1:1712557/1‑69 (MQ=255)
           gtATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGt                                                          >  1:1885044/1‑69 (MQ=255)
           gtATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGt                                                          >  1:2839312/1‑69 (MQ=255)
           gtATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGt                                                          >  1:4046275/1‑69 (MQ=255)
           gtATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGt                                                          >  1:4836164/1‑69 (MQ=255)
            tATTAAGTAAAGCGTGAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAg                                                       >  1:4766073/1‑71 (MQ=255)
            tATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCa                                                        <  1:2700580/70‑1 (MQ=255)
            tATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAg                                                       >  1:27466/1‑71 (MQ=255)
             aTTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCa                                                        <  1:795633/69‑1 (MQ=255)
              ttAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTc                                                         <  1:4967505/67‑1 (MQ=255)
              ttAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTc                                                         <  1:4757856/67‑1 (MQ=255)
                   tAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGgtg                                                               <  1:2946960/56‑1 (MQ=255)
                   tAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGgtg                                                               <  1:1101097/56‑1 (MQ=255)
                   tAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGgtg                                                               <  1:331970/56‑1 (MQ=255)
                   tAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGgtg                                                               <  1:4083300/56‑1 (MQ=255)
                   tAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGgtg                                                               <  1:3447889/56‑1 (MQ=255)
                    aaaGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGcgaccga                                               >  1:4235681/1‑71 (MQ=255)
                     aaGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGcgaccgac                                              >  1:2954197/1‑71 (MQ=255)
                     aaGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGcgaccgac                                              >  1:4133678/1‑71 (MQ=255)
                      aGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATCAAATGgtg                                                               <  1:2705584/53‑1 (MQ=255)
                      aGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATCAAATGgtg                                                               <  1:2475569/53‑1 (MQ=255)
                      aGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATCAAATGgtg                                                               <  1:1198973/53‑1 (MQ=255)
                      aGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATCAAATGgtg                                                               <  1:3756776/53‑1 (MQ=255)
                        cGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATCAAATGGTGTa                                                             >  1:2244963/1‑53 (MQ=255)
                        cGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATCAAATGGTGTa                                                             >  1:1907135/1‑53 (MQ=255)
                            aaCGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAAt                                                                   >  1:1344068/1‑43 (MQ=255)
                            aaCGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTa                                                             >  1:1282864/1‑49 (MQ=255)
                            aaCGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTa                                                             >  1:4025420/1‑49 (MQ=255)
                            aaCGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACt                                         >  1:2284267/1‑69 (MQ=255)
                            aaCGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTg                                        >  1:4093128/1‑70 (MQ=255)
                               gCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACa                                    <  1:1036311/71‑1 (MQ=255)
                               gCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACa                                    <  1:252139/71‑1 (MQ=255)
                               gCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACa                                    <  1:4651960/71‑1 (MQ=255)
                                  ggAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAAt                                  >  1:4153648/1‑70 (MQ=255)
                                   gagaTTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAAt                                  <  1:3771755/69‑1 (MQ=255)
                                   gagaTTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATc                                 <  1:4541874/70‑1 (MQ=255)
                                   gagaTTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCa                                <  1:1621801/71‑1 (MQ=255)
                                     gaTTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAg                               >  1:4729551/1‑70 (MQ=255)
                                     gaTTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGt                              >  1:1680657/1‑71 (MQ=255)
                                     gaTTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGt                              >  1:2918648/1‑71 (MQ=255)
                                        ttttCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGcga                                                   >  1:4885328/1‑47 (MQ=255)
                                        ttttCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACAGACCGACt                                         >  1:1238059/1‑57 (MQ=255)
                                             atcaGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACAt                       >  1:603723/1‑70 (MQ=255)
                                              tcaGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTgg                            <  1:1793866/64‑1 (MQ=255)
                                              tcaGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTgg                            <  1:1312903/64‑1 (MQ=255)
                                                aGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCTACCGACTGAACAATCAGTg                             <  1:1814883/61‑1 (MQ=255)
                                                aGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTg                             <  1:2371370/61‑1 (MQ=255)
                                                aGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTg                             <  1:1553864/61‑1 (MQ=255)
                                                aGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTg                             <  1:1858290/61‑1 (MQ=255)
                                                aGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTg                             <  1:3316307/61‑1 (MQ=255)
                                                aGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTg                             <  1:3729746/61‑1 (MQ=255)
                                                aGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTaa                    <  1:2887636/70‑1 (MQ=255)
                                                aGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAAc                   >  1:1720923/1‑71 (MQ=255)
                                                  gAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGt                           >  1:2948248/1‑61 (MQ=255)
                                                  gAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGt                           >  1:1278455/1‑61 (MQ=255)
                                                  gAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGt                           >  1:2724928/1‑61 (MQ=255)
                                                  gAACTCTTTTTGTCATTAAATGGGGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGt                           >  1:4757202/1‑61 (MQ=255)
                                                    aCTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCaaa               <  1:4124538/71‑1 (MQ=255)
                                                     ctctTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAAcc                  <  1:2336434/67‑1 (MQ=255)
                                                       ctTTTTGTCATTAAATGGTGTACGTCAGGcgaccgaccg                                            <  1:4988910/39‑1 (MQ=255)
                                                        tttttGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCg           >  1:4152439/1‑71 (MQ=255)
                                                          tttGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCa                                >  1:1910858/1‑48 (MQ=255)
                                                          tttGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCgg          <  1:4268024/70‑1 (MQ=255)
                                                          tttGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCgg          <  1:3941708/70‑1 (MQ=255)
                                                          tttGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCgg          <  1:2443588/70‑1 (MQ=255)
                                                              tCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCgg          <  1:4131597/66‑1 (MQ=255)
                                                              tCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCgg          <  1:2037079/66‑1 (MQ=255)
                                                              tCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCgg          <  1:2292568/66‑1 (MQ=255)
                                                              tCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCgg          <  1:75241/66‑1 (MQ=255)
                                                              tCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCgg          <  1:3166605/66‑1 (MQ=255)
                                                              tCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGtt     <  1:1080822/71‑1 (MQ=255)
                                                                aTTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGAc        >  1:2138498/1‑66 (MQ=255)
                                                                aTTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGAc        >  1:2685079/1‑66 (MQ=255)
                                                                aTTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGAc        >  1:1128400/1‑66 (MQ=255)
                                                                aTTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGAc        >  1:1920912/1‑66 (MQ=255)
                                                                aTTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTa    <  1:1765726/70‑1 (MQ=255)
                                                                aTTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTa    <  1:2833077/70‑1 (MQ=255)
                                                                aTTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTa    <  1:1369486/70‑1 (MQ=255)
                                                                aTTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTa    <  1:3012722/70‑1 (MQ=255)
                                                                aTTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTa    <  1:3199086/70‑1 (MQ=255)
                                                                 ttAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATc                                 <  1:1904640/40‑1 (MQ=255)
                                                                 ttAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATc                                 <  1:573567/40‑1 (MQ=255)
                                                                 ttAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATc                                 <  1:4806688/40‑1 (MQ=255)
                                                                 ttAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATc                                 <  1:2915199/40‑1 (MQ=255)
                                                                 ttAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATc                                 <  1:3416280/40‑1 (MQ=255)
                                                                 ttAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATc                                 <  1:2130092/40‑1 (MQ=255)
                                                                 ttAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATc                                 <  1:2460446/40‑1 (MQ=255)
                                                                 ttAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAg  >  1:704155/1‑71 (MQ=255)
                                                                 ttAAATGGTGTACGGCAGGCGACCGACCGACTGAACAATc                                 <  1:1110379/40‑1 (MQ=255)
                                                                  tAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTa                          >  1:3104496/1‑46 (MQ=255)
                                                                  tAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTa                          >  1:2737633/1‑46 (MQ=255)
                                                                  |                                                                     
CGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAG  >  minE/2950775‑2950910

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: