Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,865,714 G→A 37.9% A251A (GCG→GCA hflC → modulator for HflB protease specific for phage lambda cII repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,865,7140GA37.9% 45.4 / 49.0 79A251A (GCG→GCAhflCmodulator for HflB protease specific for phage lambda cII repressor
Reads supporting (aligned to +/- strand):  ref base G (28/21);  new base A (4/26);  total (32/47)
Fisher's exact test for biased strand distribution p-value = 1.26e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

GCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTGAGCGTCAGGGCCGCATCATGCG  >  minE/2865646‑2865776
                                                                    |                                                              
gcCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGcgcgcg                                                                <  1:3177176/69‑1 (MQ=255)
gcCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGcgcgcg                                                                <  1:603558/69‑1 (MQ=255)
      cGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGcc                                                          >  1:862094/1‑69 (MQ=255)
      cGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGcc                                                          >  1:2698775/1‑69 (MQ=255)
      cGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGcc                                                          >  1:1580532/1‑69 (MQ=255)
      cGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGcc                                                          >  1:1537524/1‑69 (MQ=255)
              ggTAGCGCGTCGTCACCGTTCACAAGGTCAGGACGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGa                                                  <  1:3363708/69‑1 (MQ=255)
              ggTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGa                                                  <  1:1666100/69‑1 (MQ=255)
              ggTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGa                                                  <  1:2566287/69‑1 (MQ=255)
              ggTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGa                                                  <  1:2643381/69‑1 (MQ=255)
              ggTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGa                                                  <  1:2657188/69‑1 (MQ=255)
              ggTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGa                                                  <  1:2308629/69‑1 (MQ=255)
              ggTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGa                                                  <  1:2078322/69‑1 (MQ=255)
              ggTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGa                                                  <  1:1982601/69‑1 (MQ=255)
              ggTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGa                                                  <  1:3026871/69‑1 (MQ=255)
              ggTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGa                                                  <  1:870478/69‑1 (MQ=255)
              ggTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGa                                                  <  1:811138/69‑1 (MQ=255)
              ggTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGa                                                  <  1:78542/69‑1 (MQ=255)
              ggTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGa                                                  <  1:779052/69‑1 (MQ=255)
                  gcgcGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTg                                              <  1:556095/69‑1 (MQ=255)
                  gcgcGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTg                                              <  1:508996/69‑1 (MQ=255)
                  gcgcGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTg                                              <  1:1207248/69‑1 (MQ=255)
                  gcgcGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTg                                              <  1:2758024/69‑1 (MQ=255)
                            aCCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGCCCAGAActc                                    >  1:1049346/1‑67 (MQ=255)
                            aCCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGc                                    >  1:800513/1‑69 (MQ=255)
                            aCCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGc                                    >  1:1556449/1‑69 (MQ=255)
                            aCCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGc                                    >  1:1096832/1‑69 (MQ=255)
                            aCCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGc                                    >  1:412718/1‑69 (MQ=255)
                            aCCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGc                                    >  1:2739139/1‑69 (MQ=255)
                            aCCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGc                                    >  1:2230606/1‑69 (MQ=255)
                            aCCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGc                                    >  1:1517876/1‑69 (MQ=255)
                            aCCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGc                                    >  1:276630/1‑69 (MQ=255)
                            aCCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGc                                    >  1:2922826/1‑69 (MQ=255)
                            aCCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGc                                    >  1:1823540/1‑69 (MQ=255)
                            aCCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGc                                    >  1:1363403/1‑69 (MQ=255)
                            aCCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAAAGc                                    >  1:1943490/1‑69 (MQ=255)
                            aCCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGAACAGAACGc                                    >  1:695219/1‑69 (MQ=255)
                            aCCGTTCACAAGGTCAGCAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACg                                     >  1:3201538/1‑68 (MQ=255)
                                  cacaAGGTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGcc                                                          <  1:3277817/41‑1 (MQ=38)
                                  cacaAGGTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGcc                                                          <  1:939221/41‑1 (MQ=38)
                                  cacaAGGTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGcc                                                          <  1:262872/41‑1 (MQ=38)
                                  cacaAGGTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGcc                                                          <  1:2417047/41‑1 (MQ=38)
                                  cacaAGGTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGcc                                                          <  1:1121494/41‑1 (MQ=38)
                                  cacaAGGTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGcc                                                          <  1:1593263/41‑1 (MQ=38)
                                  cacaAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAg                                         <  1:2690545/58‑1 (MQ=255)
                                  cacaAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAg                                         <  1:1344114/58‑1 (MQ=255)
                                        gTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTg                        <  1:1051969/69‑1 (MQ=255)
                                        gTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTg                        <  1:910505/69‑1 (MQ=255)
                                        gTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTg                        <  1:1142190/69‑1 (MQ=255)
                                        gTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTg                        <  1:1278549/69‑1 (MQ=255)
                                        gTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTg                        <  1:65940/69‑1 (MQ=255)
                                        gTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTg                        <  1:603687/69‑1 (MQ=255)
                                        gTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTg                        <  1:1581911/69‑1 (MQ=255)
                                        gTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTg                        <  1:43306/69‑1 (MQ=255)
                                        gTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTg                        <  1:3374181/69‑1 (MQ=255)
                                        gTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTg                        <  1:3286699/69‑1 (MQ=255)
                                        gTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTg                        <  1:1766242/69‑1 (MQ=255)
                                        gTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTg                        <  1:1863827/69‑1 (MQ=255)
                                        gTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTg                        <  1:2651902/69‑1 (MQ=255)
                                        gTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTg                        <  1:238898/69‑1 (MQ=255)
                                        gTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTg                        <  1:2023285/69‑1 (MQ=255)
                                        gTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTg                        <  1:1942106/69‑1 (MQ=255)
                                        gTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTg                        <  1:199465/69‑1 (MQ=255)
                                        gTCAGGAAGAAGCGGAAAAACTGCGTGCAACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTg                        <  1:3224500/69‑1 (MQ=255)
                                          cAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTGAg                      >  1:2409814/1‑69 (MQ=255)
                                          cAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTGAg                      >  1:2053463/1‑69 (MQ=255)
                                          cAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTGAg                      >  1:1578873/1‑69 (MQ=255)
                                           aGGAAGAAGCGGAAAAACTGCGTGCAACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTg                        <  1:2274666/66‑1 (MQ=255)
                                                 aagCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAg                                         <  1:1582517/43‑1 (MQ=255)
                                                      gAAAAACTGCGTGCAACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTg                        <  1:2848238/55‑1 (MQ=255)
                                                        aaaaCTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTGAGCGTCAGGGCCGcat        <  1:2489220/69‑1 (MQ=255)
                                                              gcgcgcGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAg                           >  1:1891970/1‑44 (MQ=255)
                                                              gcgcgcGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTGAGCGTCAGGGCCGCATCATGCg  >  1:3036044/1‑69 (MQ=255)
                                                              gcgcgcGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTGAGCGTCAGGGCCGCATCATGCg  >  1:659558/1‑69 (MQ=255)
                                                              gcgcgcGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTGAGCGTCAGGGCCGCATCATGCg  >  1:1524815/1‑69 (MQ=255)
                                                              gcgcgcGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTGAGCGTCAGGGCCGCATCATGCg  >  1:344216/1‑69 (MQ=255)
                                                              gcgcgcGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTGAGCGTCAGGGCCGCATCATGCg  >  1:1333052/1‑69 (MQ=255)
                                                              gcgcgcGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTGAGCGTCAGGGCCGCATCATGCg  >  1:3120740/1‑69 (MQ=255)
                                                              gcgcgcGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTGAGCGTCAGGGCCGCATCATGCg  >  1:1731994/1‑69 (MQ=255)
                                                              gcgcgcGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTGAGCGTCAGGGCCGCATCATGCg  >  1:108763/1‑69 (MQ=255)
                                                              gcgcgcGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTGAGCGTCAGGGCCGCATCATGCg  >  1:3282256/1‑69 (MQ=255)
                                                                    |                                                              
GCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTGAGCGTCAGGGCCGCATCATGCG  >  minE/2865646‑2865776

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: