Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,155,591 T→A 9.7% 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 9.7% 166.6 / 8.2 65intergenic (+102/+84)dacB/obgED‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base T (38/19);  new base A (4/2);  total (43/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGT  >  minE/2155522‑2155661
                                                                     |                                                                       
acAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                                                        >  1:1908706/2‑71 (MQ=255)
 cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                                                        <  1:2314111/70‑1 (MQ=255)
 cAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCg                                                                       >  1:320434/1‑71 (MQ=255)
 cAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCg                                                                       >  1:556950/1‑71 (MQ=255)
   gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                     >  1:986703/1‑71 (MQ=255)
     ttgttGACAACGTACGCATTGTTCATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAACAGTGCa                                                                   <  1:3369653/68‑1 (MQ=1)
        ttGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                                >  1:653638/1‑71 (MQ=255)
        ttGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCaaa                                                                 <  1:2719167/70‑1 (MQ=255)
         tGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAAtt                                                               >  1:178360/1‑71 (MQ=255)
           aCAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCaa                                                            >  1:1098698/1‑71 (MQ=11)
           aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                               <  1:1707583/69‑1 (MQ=39)
             acaGTGCGCTTTGTTTATGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCaat                                                           <  1:3620551/69‑2 (MQ=11)
             aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                            >  1:2252239/1‑70 (MQ=25)
             aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                            >  1:2247584/1‑70 (MQ=25)
             aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                            >  1:1649628/1‑70 (MQ=35)
             aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                            >  1:1840022/1‑70 (MQ=35)
                gTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACata                                                        <  1:1502857/71‑1 (MQ=39)
                 tGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACatat                                                       >  1:3163074/1‑71 (MQ=39)
                   cgcTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTg                                                     >  1:592977/1‑71 (MQ=39)
                          tttATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtc                                               <  1:112347/70‑1 (MQ=255)
                            tATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTctct                                           >  1:160016/1‑71 (MQ=255)
                                ccGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                                >  1:1427629/1‑47 (MQ=38)
                                ccGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTa                                       <  1:757168/71‑1 (MQ=255)
                                 cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                              <  1:3889732/48‑1 (MQ=255)
                                  ggATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAg                                      <  1:2820287/70‑1 (MQ=255)
                                   gATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg                                     >  1:378758/1‑71 (MQ=11)
                                     tGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCt                                  >  1:2547970/1‑71 (MQ=255)
                                      gcggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCt                                  >  1:1473641/1‑70 (MQ=255)
                                        ggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                               >  1:924232/1‑40 (MQ=40)
                                        ggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGa                                >  1:1435700/1‑70 (MQ=255)
                                        ggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGAt                               <  1:2525497/71‑1 (MQ=255)
                                        ggcgTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGa                                >  1:3172492/1‑71 (MQ=14)
                                          cgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTg                                                     <  1:204757/48‑1 (MQ=255)
                                          cgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATaa                             <  1:1881545/71‑1 (MQ=255)
                                           gTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                              >  1:1065880/1‑38 (MQ=40)
                                              aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGc                                    >  1:3349055/1‑60 (MQ=255)
                                              aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCg                          >  1:2966096/1‑70 (MQ=255)
                                              aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCg                          >  1:1649692/1‑70 (MQ=255)
                                              aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCg                          >  1:1818822/1‑70 (MQ=255)
                                               aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCCACATATTGCAATTCTCTTGTAGGCCTGATAAGCGtc                        >  1:900512/1‑70 (MQ=255)
                                               aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCg                          >  1:505943/1‑69 (MQ=255)
                                                  ccTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgcg                     >  1:1178601/1‑71 (MQ=255)
                                                   cTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCaa                                                  >  1:634349/1‑42 (MQ=40)
                                                   cTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCgca                        <  1:1723047/68‑3 (MQ=11)
                                                     tATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCAt                  >  1:3104888/1‑71 (MQ=255)
                                                      aTCCGGCCTACAAAATCGTGCAAATTCAACATATTGCTATTCTCTTGTAGGCCTGATAAGCGTAGCGCAt                  >  1:828670/1‑70 (MQ=255)
                                                         cGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAgg              <  1:669121/71‑1 (MQ=255)
                                                           gCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCt            >  1:1363920/1‑71 (MQ=255)
                                                            ccTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGc             >  1:627270/1‑69 (MQ=255)
                                                            ccTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCt            >  1:1067304/1‑70 (MQ=255)
                                                            ccTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCt            >  1:1387815/1‑70 (MQ=255)
                                                             cTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTg                                         <  1:2097864/40‑1 (MQ=255)
                                                             cTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGcaa           >  1:1006971/1‑69 (MQ=2)
                                                                cAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGt                                        >  1:2639210/1‑38 (MQ=255)
                                                                cAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGAt                               >  1:15816/1‑47 (MQ=255)
                                                                cAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAttt       >  1:380661/1‑71 (MQ=255)
                                                                cAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGcaat          <  1:2393310/69‑4 (MQ=0)
                                                                 aaaaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATa                              <  1:3633402/47‑1 (MQ=255)
                                                                 aaaaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCg                          <  1:3844684/51‑1 (MQ=255)
                                                                 aaaaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCa                <  1:3614853/61‑1 (MQ=255)
                                                                 aaaaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCa                <  1:1851692/61‑1 (MQ=255)
                                                                 aaaaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCa                <  1:962781/61‑1 (MQ=255)
                                                                   aaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGc    >  1:227844/1‑71 (MQ=255)
                                                                   aaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGc    >  1:2928902/1‑71 (MQ=255)
                                                                   aaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGc    >  1:983841/1‑71 (MQ=255)
                                                                     tCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGt  >  1:725990/1‑71 (MQ=255)
                                                                     tCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGt  <  1:3388900/71‑1 (MQ=255)
                                                                     tCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGt  >  1:2420039/1‑71 (MQ=255)
                                                                     |                                                                       
TCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGT  >  minE/2155522‑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: