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

CGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGA  >  minE/2960568‑2960701
                                                                |                                                                      
cgccgcATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGt                                                                     <  1:3491682/68‑1 (MQ=255)
cgccgcATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTtg                                                                  <  1:3762779/71‑1 (MQ=255)
    gcATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGcc                                                             >  1:2242114/1‑71 (MQ=255)
    gcATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGc                                                              >  1:2418134/1‑70 (MQ=255)
       tCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGttt                                                                   >  1:978290/1‑63 (MQ=255)
            cAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCg                                                       >  1:209117/1‑69 (MQ=255)
             aGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAg                                                    >  1:1277563/1‑71 (MQ=255)
                tGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGcc                                                 >  1:1614761/1‑71 (MQ=255)
                  gcgcACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCgg                                               <  1:1796414/71‑1 (MQ=255)
                  gcgcACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCg                                                >  1:983656/1‑70 (MQ=255)
                    gcACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTtgt                                                                 >  1:826109/1‑52 (MQ=255)
                    gcACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGa                                              >  1:2349755/1‑70 (MQ=255)
                    gcACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGAt                                             <  1:3068658/71‑1 (MQ=255)
                     cacaATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAg                                                    >  1:2811138/1‑63 (MQ=255)
                     cacaATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATa                                            <  1:156772/71‑1 (MQ=255)
                     cacaATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATa                                            <  1:2087501/71‑1 (MQ=255)
                       caATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAg                                          <  1:2911652/71‑1 (MQ=40)
                         aTGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGc                                        >  1:1761796/1‑71 (MQ=39)
                          tGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCg                                       >  1:209598/1‑71 (MQ=39)
                            ccTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCTTAGGCCGGATAAGGCGt                                      <  1:1418000/70‑1 (MQ=255)
                            ccTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGtt                                     <  1:2587964/71‑1 (MQ=39)
                            ccTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGt                                      <  1:2483711/70‑1 (MQ=40)
                             cTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGt                                      <  1:2155124/69‑1 (MQ=40)
                             cTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTc                                    <  1:2588348/71‑1 (MQ=39)
                             cTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTc                                    <  1:1889371/71‑1 (MQ=39)
                                aTGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAc                                  <  1:3466625/70‑1 (MQ=255)
                                 tGCGACGCTTAACGCGTCTTATCATGCCTACAGGttt                                                                   <  1:3780405/37‑1 (MQ=255)
                                 tGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGcc                                                 >  1:1492276/1‑54 (MQ=255)
                                 tGCGACGCTTACGCGTCTTATCATGCCTACGG‑‑‑GAACCTGACCGTAGGCCGGATAAGGCGTACACGCCGCAt                           <  1:2475241/71‑1 (MQ=2)
                                  gCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATa                                            <  1:3550161/58‑1 (MQ=255)
                                     aCGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCa                            <  1:3041883/71‑1 (MQ=255)
                                       gCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGa                                                           >  1:1099004/1‑38 (MQ=255)
                                        cTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCg                                                            <  1:3032344/36‑1 (MQ=255)
                                        cTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCg                                                            <  1:1153046/36‑1 (MQ=255)
                                           aaCGCGTCTTATCATGCCTACAGGTTTGTGCCGaa                                                          >  1:3250719/1‑35 (MQ=37)
                                           aaCGCGTCTTATCATGCCTACAGGTTTGTGCCGaa                                                          >  1:1529913/1‑35 (MQ=37)
                                            aCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                       >  1:1755987/1‑69 (MQ=37)
                                             cgcgTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAg                    <  1:3904227/71‑1 (MQ=34)
                                              gcgTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGt                   >  1:1541134/1‑71 (MQ=2)
                                              gcgTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGt                   >  1:448843/1‑71 (MQ=2)
                                               cgTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGtt                  <  1:462366/71‑1 (MQ=34)
                                               cgTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGtc                  <  1:1326788/71‑2 (MQ=2)
                                                   ttATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGgcg              >  1:2420217/1‑71 (MQ=34)
                                                     aTCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGca            >  1:3238637/1‑71 (MQ=34)
                                                        aTGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGcaca          <  1:2558370/70‑1 (MQ=35)
                                                        aTGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGcaca          <  1:508048/70‑1 (MQ=35)
                                                             tACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgccg                              <  1:3639393/45‑1 (MQ=38)
                                                              aCAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCt    >  1:2066447/1‑70 (MQ=35)
                                                               cAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGa  >  1:3768776/1‑71 (MQ=34)
                                                                |                                                                      
CGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGA  >  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: