Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,960,710 Δ1 bp 19.6% intergenic (+160/+260) 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,7100T.19.6% 134.8 / 36.5 56intergenic (+160/+260)yjjV/yjjWpredicted DNase/predicted pyruvate formate lyase activating enzyme
Reads supporting (aligned to +/- strand):  ref base T (23/22);  new base . (2/9);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 8.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCC  >  minE/2960644‑2960779
                                                                  |                                                                     
aaCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCt                                                                       >  1:5218969/1‑67 (MQ=33)
aaCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCt                                                                       >  1:3922553/1‑67 (MQ=33)
aaCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCt                                                                       >  1:3547113/1‑67 (MQ=33)
      aGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCt                                                             >  1:1568240/1‑71 (MQ=39)
          cGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTAt                                                          >  1:4005319/1‑70 (MQ=40)
          cGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGGCTTa                                                           >  1:3257840/1‑69 (MQ=255)
           ggATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATc                                                         >  1:4461945/1‑70 (MQ=40)
           ggATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCa                                                        >  1:2901661/1‑71 (MQ=39)
                   cGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTAc                                                 >  1:3276628/1‑70 (MQ=40)
                       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)
                            cgcATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTAc                                                 <  1:605918/61‑1 (MQ=40)
                            cgcATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTAc                                                 <  1:3357733/61‑1 (MQ=40)
                            cgcATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTAc                                                 <  1:5250354/61‑1 (MQ=40)
                              cATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACaa                                               <  1:2874229/61‑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)
                                      aGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGt                                             <  1:68784/55‑1 (MQ=255)
                                        ttGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCAt                                                       <  1:4466510/43‑1 (MQ=40)
                                        ttGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGc                           >  1:3776532/1‑71 (MQ=39)
                                         tGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCa                                                        >  1:5605636/1‑41 (MQ=38)
                                         tGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCa                                                        >  1:5612377/1‑41 (MQ=38)
                                         tGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTa                              >  1:1587003/1‑67 (MQ=255)
                                         tGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGcc                          >  1:5101441/1‑71 (MQ=39)
                                          ggCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCg                         >  1:1665957/1‑71 (MQ=39)
                                           gcgcACAATGCCTGATGCGACGCTTGTCGCGTCTTAt                                                          <  1:2056000/37‑1 (MQ=40)
                                           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:2678661/71‑1 (MQ=255)
                                                                 cTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATcc  <  1:1211336/71‑1 (MQ=255)
                                                                  ttGTCGCGTCTTATCATGCCTACAAGTCTGTGCCg                                     <  1:4764514/35‑1 (MQ=255)
                                                                  |                                                                     
AACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCC  >  minE/2960644‑2960779

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: