Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,960,632 Δ1 bp 8.2% intergenic (+82/+338) 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,6320A. 8.2% 135.2 / 8.2 49intergenic (+82/+338)yjjV/yjjWpredicted DNase/predicted pyruvate formate lyase activating enzyme
Reads supporting (aligned to +/- strand):  ref base A (21/24);  new base . (2/2);  total (23/26)
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

CGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGA  >  minE/2960568‑2960701
                                                                |                                                                      
cgccgcATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGt                                                                     <  1:3491682/68‑1 (MQ=255)
cgccgcATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTtg                                                                  <  1:3762779/71‑1 (MQ=255)
    gcATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑Gcc                                                             >  1:2242114/1‑71 (MQ=255)
    gcATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑Gc                                                              >  1:2418134/1‑70 (MQ=255)
       tCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGttt                                                                   >  1:978290/1‑63 (MQ=255)
            cAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCg                                                       >  1:209117/1‑69 (MQ=255)
             aGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAg                                                    >  1:1277563/1‑71 (MQ=255)
                tGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAGGcc                                                 >  1:1614761/1‑71 (MQ=255)
                  gcgcACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAGGCCgg                                               <  1:1796414/71‑1 (MQ=255)
                  gcgcACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAGGCCg                                                >  1:983656/1‑70 (MQ=255)
                    gcACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTtgt                                                                 >  1:826109/1‑52 (MQ=255)
                    gcACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAGGCCGGa                                              >  1:2349755/1‑70 (MQ=255)
                    gcACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAGGCCGGAt                                             <  1:3068658/71‑1 (MQ=255)
                     cacaATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAg                                                    >  1:2811138/1‑63 (MQ=255)
                     cacaATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAGGCCGGATa                                            <  1:156772/71‑1 (MQ=255)
                     cacaATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAGGCCGGATa                                            <  1:2087501/71‑1 (MQ=255)
                       caATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAGGCCGGATAAg                                          <  1:2911652/71‑1 (MQ=40)
                         aTGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAGGCCGGATAAGGc                                        >  1:1761796/1‑71 (MQ=39)
                          tGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAGGCCGGATAAGGCg                                       >  1:209598/1‑71 (MQ=39)
                            ccTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCTTAGGCCGGATAAGGCGt                                      <  1:1418000/70‑1 (MQ=255)
                            ccTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAGGCCGGATAAGGCGtt                                     <  1:2587964/71‑1 (MQ=39)
                            ccTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAGGCCGGATAAGGCGt                                      <  1:2483711/70‑1 (MQ=40)
                             cTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAGGCCGGATAAGGCGt                                      <  1:2155124/69‑1 (MQ=40)
                             cTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAGGCCGGATAAGGCGTTc                                    <  1:2588348/71‑1 (MQ=39)
                             cTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAGGCCGGATAAGGCGTTc                                    <  1:1889371/71‑1 (MQ=39)
                                aTGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAGGCCGGATAAGGCGTTCAc                                  <  1:3466625/70‑1 (MQ=255)
                                 tGCGACGCTTAACGCGTCTTATCATGCCTACAGGttt                                                                   <  1:3780405/37‑1 (MQ=255)
                                 tGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAGGcc                                                 >  1:1492276/1‑54 (MQ=255)
                                 tGCGACGCTT‑ACGCGTCTTATCATGCCTAC‑GG‑‑‑GA‑ACCTGACCGTAGGCCGGATAAGGCGTACACGCCGCAt                           <  1:2475241/71‑1 (MQ=2)
                                  gCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAGGCCGGATa                                            <  1:3550161/58‑1 (MQ=255)
                                     aCGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCa                            <  1:3041883/71‑1 (MQ=255)
                                       gCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGa                                                           >  1:1099004/1‑38 (MQ=255)
                                        cTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCg                                                            <  1:3032344/36‑1 (MQ=255)
                                        cTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCg                                                            <  1:1153046/36‑1 (MQ=255)
                                           aaCGCGTCTTATCATGCCTACAGGTTTGT‑GCCGaa                                                          >  1:3250719/1‑35 (MQ=37)
                                           aaCGCGTCTTATCATGCCTACAGGTTTGT‑GCCGaa                                                          >  1:1529913/1‑35 (MQ=37)
                                            aCGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                       >  1:1755987/1‑69 (MQ=37)
                                             cgcgTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAg                    <  1:3904227/71‑1 (MQ=34)
                                              gcgTCTTATCAGGCCTAC‑GGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGt                   >  1:1541134/1‑71 (MQ=2)
                                              gcgTCTTATCAGGCCTAC‑GGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGt                   >  1:448843/1‑71 (MQ=2)
                                               cgTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGtt                  <  1:462366/71‑1 (MQ=34)
                                               cgTCTTATCAGGCCTAC‑GGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGtc                  <  1:1326788/71‑2 (MQ=2)
                                                   ttATCATGCCTACAGGTTTGT‑GCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGgcg              >  1:2420217/1‑71 (MQ=34)
                                                     aTCATGCCTACAGGTTTGT‑GCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGca            >  1:3238637/1‑71 (MQ=34)
                                                        aTGCCTACAGGTTTGT‑GCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGcaca          <  1:2558370/70‑1 (MQ=35)
                                                        aTGCCTACAGGTTTGT‑GCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGcaca          <  1:508048/70‑1 (MQ=35)
                                                             tACAGGTTTGT‑GCCGAACCGTAGGCCGGATAAGGCGTTCAcgccg                              <  1:3639393/45‑1 (MQ=38)
                                                              aCAGGTTTGT‑GCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCt    >  1:2066447/1‑70 (MQ=35)
                                                               cAGGTTTGT‑GCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGa  >  1:3768776/1‑71 (MQ=34)
                                                                |                                                                      
CGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGA  >  minE/2960568‑2960701

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: ‑