Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,865,711 C→T 43.4% R250R (CGC→CGT 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,7110CT43.4% 19.7 / 49.8 69R250R (CGC→CGThflCmodulator for HflB protease specific for phage lambda cII repressor
Reads supporting (aligned to +/- strand):  ref base C (18/21);  new base T (4/26);  total (22/47)
Fisher's exact test for biased strand distribution p-value = 4.37e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTGAGCGTCAGGGCCGCATCATGCG  >  minE/2865645‑2865776
                                                                  |                                                                 
cgcCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTgcgcgc                                                                 <  1:1093366/69‑1 (MQ=255)
 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)
                                                                  |                                                                 
CGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTGAGCGTCAGGGCCGCATCATGCG  >  minE/2865645‑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: