Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,864,584 C→T 24.1% A295A (GCC→GCT 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,5840CT24.1% 51.3 / 14.6 29A295A (GCC→GCThflKmodulator for HflB protease specific for phage lambda cII repressor
Reads supporting (aligned to +/- strand):  ref base C (4/18);  new base T (4/3);  total (8/21)
Fisher's exact test for biased strand distribution p-value = 6.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

CCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAG  >  minE/2864514‑2864648
                                                                      |                                                                
ccAACGAAGTTCAGCCGCGTGCGAACGGTCAGTCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGcc                                                                  >  1:1576258/1‑71 (MQ=255)
  aaCGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCa                                                                <  1:2082847/71‑1 (MQ=255)
  aaCGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCa                                                                <  1:1787366/71‑1 (MQ=255)
     gAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGAc                                                             >  1:2118823/1‑71 (MQ=255)
       aGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCTCAGa                                                              >  1:388154/1‑68 (MQ=255)
        gTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCTCAGACCAt                                                          <  1:731095/71‑1 (MQ=255)
          tCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATcc                                                        >  1:1614368/1‑71 (MQ=255)
           cAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATcc                                                        <  1:1003668/70‑1 (MQ=255)
                gcgTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCa                                                           <  1:2079136/62‑1 (MQ=255)
                     cGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCa                                              <  1:1855088/70‑1 (MQ=255)
                          ggTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCTCAGACCATCCTGGAAGCTCAGGGTGa                                        <  1:2171794/71‑1 (MQ=255)
                          ggTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTg                                         >  1:106255/1‑70 (MQ=255)
                              aGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTg                                    <  1:1781382/71‑1 (MQ=255)
                              aGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTAAAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTg                                    <  1:1404117/71‑1 (MQ=255)
                                    aaCGTATCCTCGAAGAGTCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGa                                        <  1:2181329/61‑1 (MQ=255)
                                    aaCGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCa                                                                <  1:1994513/37‑1 (MQ=255)
                                    aaCGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGa                                        <  1:250894/61‑1 (MQ=255)
                                            cTCGAAGAGGCGCGTGCGTACAAGGCTCAGACCATcc                                                        >  1:751190/1‑37 (MQ=255)
                                               gAAGAGGCGCGTGCGTACAAGGCTCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAAc                    >  1:1527877/1‑70 (MQ=255)
                                               gAAGAGGCGCGTGCGTACAAGGCTCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACt                   >  1:1679797/1‑71 (MQ=255)
                                               gAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACt                   >  1:205711/1‑71 (MQ=255)
                                                 agagGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTc                 <  1:815369/71‑1 (MQ=255)
                                                        cgtgcgtACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCgg           <  1:1869634/70‑1 (MQ=255)
                                                             gtACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAatat      <  1:1360471/70‑1 (MQ=255)
                                                             gtACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAAtata     <  1:1100685/71‑1 (MQ=255)
                                                              tACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTa                       <  1:226124/52‑1 (MQ=255)
                                                               aCAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTc                 <  1:2072152/57‑1 (MQ=255)
                                                                cAAGGCTCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTa                       <  1:1333468/50‑1 (MQ=255)
                                                                cAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAg  <  1:703891/71‑1 (MQ=255)
                                                                cAAGGCCCAGAACATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAg  <  1:785830/71‑1 (MQ=255)
                                                                      |                                                                
CCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAG  >  minE/2864514‑2864648

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: