Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,864,086 C→T 32.8% I129I (ATC→ATT hflK → 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,864,0860CT32.8% 65.9 / 40.3 79I129I (ATC→ATThflKmodulator for HflB protease specific for phage lambda cII repressor
Reads supporting (aligned to +/- strand):  ref base C (29/24);  new base T (6/20);  total (35/44)
Fisher's exact test for biased strand distribution p-value = 8.91e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

CGCTTTGGTAAATTCAGCCATCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAGCCGTGCGTGAACTGGCCGCTTCTGGTGTGATGCTG  >  minE/2864027‑2864149
                                                           |                                                               
cGCTTTGGTAAATTCAGCCATCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGt                                                         >  1:870967/1‑68 (MQ=255)
cGCTTTGGTAAATTCAGCCATCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGt                                                         >  1:775389/1‑68 (MQ=255)
cGCTTTGGTAAATTCAGCCATCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGt                                                         >  1:3063586/1‑68 (MQ=255)
cGCTTTGGTAAATTCAGCCATCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGt                                                         >  1:172016/1‑68 (MQ=255)
cGCTTTGGTAAATTCAGCCATCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGt                                                         >  1:537445/1‑68 (MQ=255)
cGCTTTGGTAAATTCAGCCATCTGGTTGAGCCGGGTCTGAACTCGAAACCGACGTTTATCGACGAAGt                                                         >  1:2789333/1‑68 (MQ=255)
cGCTTTGGTAAATTCAGCCATCTGGTTGAGCCGGGTCTGAACCGGAAACCGACGTTTATCGACGAAGt                                                         >  1:2398421/1‑68 (MQ=255)
 gCTTTGGTAAATTCAGCCATCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGTCa                                                       <  1:2369062/69‑1 (MQ=255)
 gCTTTGGTAAATTCAGCCATCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGTCa                                                       <  1:2318568/69‑1 (MQ=255)
 gCTTTGGTAAATTCAGCCATCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGTCa                                                       <  1:1079286/69‑1 (MQ=255)
  cTTTGGTAAATTCAGCCATCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGTCaa                                                      <  1:1090847/69‑1 (MQ=255)
  cTTTGGTAAATTCAGCCATCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGTCaa                                                      <  1:1654080/69‑1 (MQ=255)
  cTTTGGTAAATTCAGCCATCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGTCaa                                                      <  1:161207/69‑1 (MQ=255)
  cTTTGGTAAATTCAGCCATCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGTCaa                                                      <  1:1264161/69‑1 (MQ=255)
  cTTTGGTAAATTCAGCCATCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGTCa                                                       <  1:2372852/68‑1 (MQ=255)
    ttGGTAAATTCAGCCATCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGTCAAAc                                                    <  1:3215330/69‑1 (MQ=255)
    ttGGTAAATTCAGCCATCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGTCAAAc                                                    <  1:96197/69‑1 (MQ=255)
    ttGGTAAATTCAGCCATCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGTCAAAc                                                    <  1:1200428/69‑1 (MQ=255)
    ttGGTAAATTCAGCCATCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGTCAAAc                                                    <  1:2889108/69‑1 (MQ=255)
    ttGGTAAATTCAGCCATCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGTCAAAc                                                    <  1:456471/69‑1 (MQ=255)
    ttGGTAAATTCAGCCATCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGTCAAAc                                                    <  1:3318094/69‑1 (MQ=255)
                gCCATCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGTCAAACCGGTGAACGTg                                         <  1:1523831/68‑1 (MQ=255)
                    tCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTGGAAGc                                    >  1:1937566/1‑69 (MQ=255)
                    tCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTGGAAGc                                    >  1:732034/1‑69 (MQ=255)
                    tCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTGGAAGc                                    >  1:1429835/1‑69 (MQ=255)
                    tCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTGGAAGc                                    >  1:2085234/1‑69 (MQ=255)
                    tCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTGGAAGc                                    >  1:1252369/1‑69 (MQ=255)
                    tCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTGGAAGc                                    >  1:149741/1‑69 (MQ=255)
                      tGGTTGAGCCGGGTCTGAACTGTAAACCGACGTTTATcgacg                                                             <  1:758266/42‑1 (MQ=255)
                      tGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATcgacg                                                             >  1:3387752/1‑42 (MQ=255)
                      tGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATcgacg                                                             >  1:2991003/1‑42 (MQ=255)
                      tGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATcgacg                                                             >  1:498663/1‑42 (MQ=255)
                      tGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATcgacg                                                             >  1:411717/1‑42 (MQ=255)
                      tGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAgt                                    >  1:2806769/1‑66 (MQ=255)
                      tGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAg                                     >  1:2925939/1‑66 (MQ=255)
                      tGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAg                                     >  1:1436449/1‑66 (MQ=255)
                      tGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAGc                                    >  1:999042/1‑67 (MQ=255)
                      tGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAGc                                    >  1:472120/1‑67 (MQ=255)
                      tGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAGc                                    >  1:951405/1‑67 (MQ=255)
                              ccGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAGCCGTGCGTGa                           <  1:1352107/68‑1 (MQ=255)
                               cGGGTCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTg                                         <  1:127879/53‑1 (MQ=255)
                               cGGGTCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTg                                         <  1:757095/53‑1 (MQ=255)
                               cGGGTCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTg                                         <  1:716989/53‑1 (MQ=255)
                               cGGGTCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTg                                         <  1:570753/53‑1 (MQ=255)
                               cGGGTCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTg                                         <  1:1206645/53‑1 (MQ=255)
                               cGGGTCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTg                                         <  1:486433/53‑1 (MQ=255)
                               cGGGTCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTg                                         <  1:409745/53‑1 (MQ=255)
                               cGGGTCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTg                                         <  1:1480732/53‑1 (MQ=255)
                               cGGGTCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTg                                         <  1:18025/53‑1 (MQ=255)
                               cGGGTCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTg                                         <  1:1917788/53‑1 (MQ=255)
                               cGGGTCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTg                                         <  1:325635/53‑1 (MQ=255)
                               cGGGTCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTg                                         <  1:3240982/53‑1 (MQ=255)
                               cGGGTCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTg                                         <  1:3086585/53‑1 (MQ=255)
                               cGGGTCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTg                                         <  1:2595010/53‑1 (MQ=255)
                               cGGGTCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTg                                         <  1:2080840/53‑1 (MQ=255)
                               cGGGTCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTg                                         <  1:2022776/53‑1 (MQ=255)
                               cGGGTCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTg                                         <  1:2940974/53‑1 (MQ=255)
                               cGGGTCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTg                                         <  1:1959672/53‑1 (MQ=255)
                               cGGGTCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTg                                         <  1:2012019/53‑1 (MQ=255)
                                   tCTGAACTGGAAACCGACGTTTATTGACGAAGTCAAACCGGTGAACGTg                                         <  1:3338854/49‑1 (MQ=255)
                                   tCTGAACTGGAAACCGACGTTTATCGACGAAGTCa                                                       <  1:792807/35‑1 (MQ=255)
                                     tGAACTGGAAACCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAGCcgtgcg                              <  1:347603/58‑1 (MQ=255)
                                     tGAACTGGAAACCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAGCcgtgcg                              <  1:2580605/58‑1 (MQ=255)
                                     tGAACTGGAAACCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAGCcgtgcg                              <  1:516715/58‑1 (MQ=255)
                                         cTGGAAACCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAGCCGTGCGTGAACTGGCCGc                  <  1:244968/66‑1 (MQ=255)
                                              aaCCGACGTTTATCGGCGAAGTCAAACCGGTGAACGTGGAAGCCGTGCGTGAACTGGCCGCTTCTGgt           >  1:903911/1‑68 (MQ=255)
                                              aaCCGACGTTTATCGACGAAGTCAAACCGGTGATCGAGGAAGCCGTGCGTGAACTGGCAGCttgt              >  1:869625/1‑63 (MQ=255)
                                              aaCCGACGTTTATCGACGAAGTCAAACCGGTGAAc                                            >  1:1434932/1‑35 (MQ=255)
                                              aaCCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAGCCGTGCGTGAACTGGCCGc                  >  1:951303/1‑61 (MQ=255)
                                              aaCCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAGCCGTGCGTGAACTGGCCGCTTCTGgt           >  1:2872574/1‑68 (MQ=255)
                                              aaCCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAGCCGTGCGTGAACTGGCCGCTTCTGgt           >  1:1118131/1‑68 (MQ=255)
                                              aaCCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAGCCGTGCGTGAACTGGCCGCTTCTGgt           >  1:590556/1‑68 (MQ=255)
                                              aaCCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAGCCGTGCGTGAACTGGCCGCTTCTGgt           >  1:280664/1‑68 (MQ=255)
                                              aaCCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAGCCGTGCGTGAACTGGCCGCTTCTGgt           >  1:2101531/1‑68 (MQ=255)
                                              aaCCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAGCCGTGCGTGAACTGGCCGCTTCTGgt           >  1:2699235/1‑68 (MQ=255)
                                              aaCCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAGCCGTGCGTGAACTGGCCGCTTCTGgt           >  1:2514709/1‑68 (MQ=255)
                                              aaCCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAGCCGTGCGTGAACTGGCCGCTTCTGgt           >  1:771202/1‑68 (MQ=255)
                                              aaCCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAGCCGTGCGTGAACTGGCCGCTTCTGgt           >  1:2019829/1‑68 (MQ=255)
                                               aCCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAGCCGTGCGTGAACTGGCCGCTTCTGgtgt         <  1:3053067/69‑1 (MQ=255)
                                               aCCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAGCCGTGCGTGAACTGGCCGCTTCTGgtgt         <  1:2342323/69‑1 (MQ=255)
                                               aCCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAGCCGTGCGTGAACTGGCCGCTTCTGgtgt         <  1:3075784/69‑1 (MQ=255)
                                                  gACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAGCCGTGCGTGAACTGGCCGCTTCTGGTGTGAt      >  1:1264708/1‑69 (MQ=255)
                                                      tttATCGACGAAGTCAAACCGGTGAACGTGGAAGCCGTGCGTGAACTGGCCGCTTCTGGTGTGATGCTg  >  1:425399/1‑69 (MQ=255)
                                                      tttATCGACGAAGTCAAACCGGTGAACGTGGAAGCCGTGCGTGAACTGGCCGCTTCTGGTGTGATGCTg  >  1:1473681/1‑69 (MQ=255)
                                                           |                                                               
CGCTTTGGTAAATTCAGCCATCTGGTTGAGCCGGGTCTGAACTGGAAACCGACGTTTATCGACGAAGTCAAACCGGTGAACGTGGAAGCCGTGCGTGAACTGGCCGCTTCTGGTGTGATGCTG  >  minE/2864027‑2864149

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: