Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,864,584 C→T 14.0% 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,5840CT14.0% 174.0 / 17.1 64A295A (GCC→GCThflKmodulator for HflB protease specific for phage lambda cII repressor
Reads supporting (aligned to +/- strand):  ref base C (24/31);  new base T (4/5);  total (28/36)
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

CCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCG  >  minE/2864514‑2864651
                                                                      |                                                                   
ccAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGcc                                                                     <  1:2760352/71‑1 (MQ=255)
  aaCGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCa                                                                   <  1:2064772/71‑1 (MQ=255)
           cAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCt                                                          >  1:3344935/1‑71 (MQ=255)
           cAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCt                                                          >  1:1285499/1‑71 (MQ=255)
              ccGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTgg                                                        >  1:3370000/1‑70 (MQ=255)
              ccGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTgg                                                        >  1:2184930/1‑70 (MQ=255)
              ccGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTgg                                                        >  1:1441571/1‑70 (MQ=255)
                   tGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTc                                                  <  1:2158384/71‑1 (MQ=255)
                   tGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTc                                                  <  1:370184/71‑1 (MQ=255)
                   tGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTc                                                  <  1:4013714/71‑1 (MQ=255)
                   tGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTc                                                  <  1:1995691/71‑1 (MQ=255)
                     cGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCa                                                 <  1:532261/70‑1 (MQ=255)
                         cGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCTCAGACCATCCTGGAAGCTCAGGGTg                                            <  1:4040995/71‑1 (MQ=255)
                         cGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCTCAGACCATCCTGGAAGCTCAGGGTg                                            <  1:1952255/71‑1 (MQ=255)
                         cGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCTCAGACCATCCTGGAAGCTCAGGGTg                                            <  1:1483321/71‑1 (MQ=255)
                         cGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCTCAGACCATCCTGGAAGCTCAGGGTg                                            <  1:1065807/71‑1 (MQ=255)
                         cGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCTCAGACCATCCTGGAAGCTCAGGGTg                                            <  1:2078050/71‑1 (MQ=255)
                         cGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTg                                            <  1:879235/71‑1 (MQ=255)
                          ggTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGa                                           <  1:1381319/71‑1 (MQ=255)
                            tCAGGCGCAACGTATCCTCGAAGCGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAgg                                               >  1:715833/1‑65 (MQ=255)
                            tCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAgg                                               >  1:1068174/1‑65 (MQ=255)
                            tCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAgg                                               >  1:2490352/1‑65 (MQ=255)
                            tCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAgg                                               >  1:75109/1‑65 (MQ=255)
                            tCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAgg                                               >  1:1554845/1‑65 (MQ=255)
                            tCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAg                                         >  1:3232403/1‑71 (MQ=255)
                                      cGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCa                                                              >  1:2542778/1‑40 (MQ=255)
                                      cGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTgg                                                        >  1:3103890/1‑46 (MQ=255)
                                      cGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAg                                                >  1:2476451/1‑54 (MQ=255)
                                      cGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTg                                            <  1:3111884/58‑1 (MQ=255)
                                      cGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCtt                               <  1:849785/71‑1 (MQ=255)
                                      cGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCtt                               <  1:3840120/71‑1 (MQ=255)
                                      cGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCtt                               <  1:2889028/71‑1 (MQ=255)
                                      cGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCtt                               <  1:4022322/71‑1 (MQ=255)
                                      cGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCtt                               <  1:2042699/71‑1 (MQ=255)
                                        tATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGaa                                          >  1:2538457/1‑58 (MQ=255)
                                        tATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGaa                                          >  1:3996691/1‑58 (MQ=255)
                                        tATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGaa                                          >  1:3535501/1‑58 (MQ=255)
                                        tATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGaa                                          >  1:1036279/1‑58 (MQ=255)
                                        tATCCTCGAAGAGGAGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAggg                                              >  1:3844344/1‑54 (MQ=255)
                                          tCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGa                                           <  1:1453353/55‑1 (MQ=255)
                                          tCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGa                                           <  1:3643733/55‑1 (MQ=255)
                                            cTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCTCGCTTTGCTa                          <  1:1367011/70‑1 (MQ=255)
                                            cTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTa                          <  1:1707157/70‑1 (MQ=255)
                                               gAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAAc                       <  1:2678219/70‑1 (MQ=255)
                                               gAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACt                      <  1:4039567/71‑1 (MQ=255)
                                               gAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACt                      <  1:913597/71‑1 (MQ=255)
                                               gAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACt                      <  1:1057811/71‑1 (MQ=255)
                                                          tgcgtACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGaa                                          <  1:2401120/40‑1 (MQ=255)
                                                          tgcgtACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGaa                                          <  1:1633751/40‑1 (MQ=255)
                                                          tgcgtACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGaa                                          <  1:1641408/40‑1 (MQ=255)
                                                          tgcgtACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGaa                                          <  1:1662008/40‑1 (MQ=255)
                                                          tgcgtACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGaa                                          <  1:555550/40‑1 (MQ=255)
                                                          tgcgtACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGaa                                          <  1:1675678/40‑1 (MQ=255)
                                                          tgcgtACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGaa                                          <  1:2606053/40‑1 (MQ=255)
                                                          tgcgtACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGaa                                          <  1:3875039/40‑1 (MQ=255)
                                                          tgcgtACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGaa                                          <  1:3713176/40‑1 (MQ=255)
                                                          tgcgtACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCtt                               >  1:3682442/1‑51 (MQ=255)
                                                          tgcgtACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGaa            >  1:2942485/1‑70 (MQ=255)
                                                            cgtACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAata          >  1:2418086/1‑70 (MQ=255)
                                                                cAAGGCTCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAg     >  1:4053434/1‑71 (MQ=255)
                                                                cAAGGCTCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAg     >  1:2926478/1‑71 (MQ=255)
                                                                cAAGGCTCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAg     >  1:1190016/1‑71 (MQ=255)
                                                                 aaGGCTCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGc    >  1:1839895/1‑71 (MQ=255)
                                                                   ggCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTg                             >  1:3664645/1‑44 (MQ=255)
                                                                   ggCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCcg  >  1:3517229/1‑71 (MQ=255)
                                                                      |                                                                   
CCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCG  >  minE/2864514‑2864651

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: