Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,155,591 T→A 5.4% intergenic (+102/+84) dacB → / ← obgE D‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,155,5910TA 5.4% 272.8 / 6.1 96intergenic (+102/+84)dacB/obgED‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base T (51/39);  new base A (3/2);  total (54/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

AGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGT  >  minE/2155524‑2155661
                                                                   |                                                                       
aGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGt                                                                      <  1:1843593/71‑1 (MQ=255)
aGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGt                                                                      >  1:3904086/1‑71 (MQ=255)
  attgttGACAAGGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                     >  1:3074179/1‑70 (MQ=11)
      ttGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                                <  1:3978830/71‑1 (MQ=255)
        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                               >  1:1963225/1‑70 (MQ=40)
         aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                              <  1:2485088/70‑1 (MQ=255)
         aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                              <  1:300113/70‑1 (MQ=39)
         aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                              <  1:411988/70‑1 (MQ=39)
          cAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCaa                                                            >  1:242103/1‑70 (MQ=11)
          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                            <  1:625065/71‑1 (MQ=25)
          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                            >  1:2413976/1‑71 (MQ=35)
           acaGTGCGCTTTGTTTATGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCaa                                                            >  1:3598848/3‑70 (MQ=11)
            aaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACa                                                          >  1:85039/1‑71 (MQ=38)
             aGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACa                                                          >  1:1464419/1‑70 (MQ=39)
                  gcTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGc                                                    <  1:1571595/71‑1 (MQ=39)
                   cTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGc                                                    >  1:187516/1‑70 (MQ=40)
                     ttGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATAtt                                                      >  1:3317833/1‑66 (MQ=40)
                      tGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAAtt                                                >  1:294861/1‑71 (MQ=255)
                       gTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtc                                               >  1:580429/1‑71 (MQ=255)
                       gTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtc                                               >  1:2543323/1‑71 (MQ=255)
                        tttATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCaata                                                          <  1:3599264/59‑3 (MQ=255)
                          tATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACatat                                                       >  1:1554861/1‑60 (MQ=255)
                          tATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtctc                                            <  1:2417866/70‑1 (MQ=255)
                           aTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTctct                                           <  1:109204/70‑1 (MQ=255)
                             gCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTg                                         >  1:2310589/1‑70 (MQ=255)
                              ccGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTg                                         >  1:3414943/1‑69 (MQ=255)
                               cGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtc                                               <  1:1244988/63‑1 (MQ=255)
                                ggATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAgg                                     >  1:1199472/1‑71 (MQ=255)
                                 gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACata                                                        >  1:1385463/1‑52 (MQ=255)
                                 gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTg                                         <  1:2005899/66‑1 (MQ=255)
                                  aTGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCtt                                          <  1:2520956/64‑1 (MQ=255)
                                  aTGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCtt                                          <  1:2042141/64‑1 (MQ=255)
                                  aTGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCtt                                          <  1:2525949/64‑1 (MQ=255)
                                  aTGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTg                                         <  1:3420709/65‑1 (MQ=255)
                                   tGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                              <  1:2928581/44‑1 (MQ=38)
                                   tGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCt                                  >  1:286787/1‑71 (MQ=255)
                                     cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCaaa                                                                 >  1:3315927/1‑39 (MQ=255)
                                     cggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGa                                >  1:53033/1‑71 (MQ=255)
                                     cggcgTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt                                  <  1:1238079/70‑1 (MQ=11)
                                      ggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACa                                                          <  1:858955/45‑1 (MQ=255)
                                      ggcgTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGa                                >  1:2890438/1‑71 (MQ=14)
                                       gcgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGAt                               >  1:924625/1‑70 (MQ=255)
                                        cgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGAt                               >  1:4094102/1‑69 (MQ=255)
                                          tAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACa                                                          <  1:3162022/41‑1 (MQ=255)
                                          tAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACa                                                          <  1:3514097/41‑1 (MQ=255)
                                          tAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACa                                                          <  1:1949693/41‑1 (MQ=255)
                                          tAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGc                           <  1:2910537/71‑1 (MQ=255)
                                            aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACata                                                        >  1:1117306/1‑41 (MQ=255)
                                            aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTctct                                           >  1:3547021/1‑53 (MQ=255)
                                            aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTa                                       >  1:221876/1‑57 (MQ=255)
                                            aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAgga                                    <  1:973272/60‑2 (MQ=255)
                                            aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCg                          >  1:380681/1‑70 (MQ=255)
                                            aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCg                          >  1:1533958/1‑70 (MQ=255)
                                            aaCGCCTTATCCGGCCTACAAAATCCTGCAAATTCAACATATTGCa                                                   >  1:3715942/1‑46 (MQ=255)
                                            aaCGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCg                          >  1:3957879/1‑71 (MQ=11)
                                              cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACata                                                        >  1:383228/1‑39 (MQ=40)
                                              cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGc                                                    <  1:1573077/43‑1 (MQ=38)
                                              cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTg                                         >  1:3306831/1‑53 (MQ=255)
                                                 cTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTa                                       >  1:3476454/1‑52 (MQ=255)
                                                  ttATCCGGCCTACATAATCGTGCAAATTCAACATATTGCAATtctc                                            <  1:1403654/46‑1 (MQ=255)
                                                   tATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCAt                  >  1:1262877/1‑71 (MQ=255)
                                                    aTCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCtt                                          >  1:1560208/1‑46 (MQ=255)
                                                    aTCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCtt                                          >  1:110641/1‑46 (MQ=255)
                                                    aTCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCtt                                          >  1:2851613/1‑46 (MQ=255)
                                                    aTCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGAt                               >  1:2462756/1‑57 (MQ=255)
                                                    aTCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCAt                  >  1:12827/1‑70 (MQ=255)
                                                    aTCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCAt                  >  1:192915/1‑70 (MQ=255)
                                                    aTCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCAt                  >  1:2442764/1‑70 (MQ=255)
                                                     tCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtctc                                            <  1:765947/43‑1 (MQ=255)
                                                     tCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCa                <  1:1051956/71‑1 (MQ=255)
                                                      ccGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtctc                                            <  1:1874577/42‑1 (MQ=255)
                                                       cGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAg               >  1:3777135/1‑70 (MQ=255)
                                                       cGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAg               >  1:3772497/1‑70 (MQ=255)
                                                        ggCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCtt                                          >  1:2745337/1‑42 (MQ=255)
                                                        ggCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGc             >  1:302556/1‑71 (MQ=255)
                                                            tACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgcg                     <  1:519756/59‑1 (MQ=255)
                                                            tACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgcg                     <  1:1384782/59‑1 (MQ=255)
                                                            tACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgcg                     <  1:1273044/59‑1 (MQ=255)
                                                            tACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCa                   <  1:1218725/61‑1 (MQ=255)
                                                            tACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCa                   <  1:1223280/61‑1 (MQ=255)
                                                            tACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAgc                      <  1:2559198/59‑1 (MQ=2)
                                                              cAAAATCGTGCAAATTCAACATATTGCAATTCTCTTg                                         <  1:2486175/37‑1 (MQ=255)
                                                              cAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGt                         <  1:1136094/53‑1 (MQ=255)
                                                              cAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgcgc                    <  1:632625/58‑1 (MQ=255)
                                                              cAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCAt                  >  1:1116796/1‑60 (MQ=255)
                                                              cAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCAt                  >  1:3278021/1‑60 (MQ=255)
                                                              cAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCa                <  1:4085790/62‑1 (MQ=255)
                                                              cAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCa                <  1:2175808/62‑1 (MQ=255)
                                                              cAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCa                <  1:3014952/62‑1 (MQ=255)
                                                              cAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAttt       >  1:1587914/1‑71 (MQ=255)
                                                              cAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAttt       >  1:429785/1‑71 (MQ=255)
                                                              cAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAttt       >  1:91834/1‑71 (MQ=255)
                                                              cAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAttt       >  1:1026083/1‑71 (MQ=255)
                                                                aaaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGc                           >  1:3694941/1‑49 (MQ=255)
                                                                aaaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTgg     >  1:1141643/1‑71 (MQ=255)
                                                                 aaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAttt       <  1:3397377/68‑1 (MQ=255)
                                                                   tCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgc                      >  1:3944330/1‑51 (MQ=255)
                                                                   tCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGt  <  1:1613918/71‑1 (MQ=255)
                                                                   tCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGt  <  1:2526554/71‑1 (MQ=255)
                                                                   tCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGt  <  1:3271417/71‑1 (MQ=255)
                                                                   |                                                                       
AGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGT  >  minE/2155524‑2155661

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: