Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,064,095 C→G 5.3% T112T (ACC→ACG rpoD → RNA polymerase, sigma 70 (sigma D) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,064,0950CG 5.3% 452.3 / 10.1 132T112T (ACC→ACGrpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base C (88/36);  new base G (3/4);  total (91/41)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

Only 100 of 137 total aligned reads displayed.

TGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTG  >  minE/2064026‑2064165
                                                                     |                                                                      
tGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGAccc                                                                       >  1:1193875/1‑71 (MQ=255)
  aaaTCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCcgc                                                                     >  1:656230/1‑71 (MQ=255)
   aaTCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCcgc                                                                     >  1:3110790/1‑70 (MQ=255)
      cGGCCGCACCACCGACCCGGTCCGCATGTACATGCGCGAGATGGGCGCGGTGTCGCTGCTGACGCGCGAAg                                                                 >  1:4729303/1‑71 (MQ=255)
       gggCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAgg                                                                >  1:5016737/1‑71 (MQ=255)
         gcgcACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGc                                                               <  1:6257689/70‑1 (MQ=255)
          cgcACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGa                                                             <  1:6080671/71‑1 (MQ=255)
             aCGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAAt                                                          >  1:4572189/1‑71 (MQ=255)
               gACTGACCCGGTACGCATGTACATGCGTGAAATGGGTACCGTTGAACTGTTGACCCGCGAAGGCGAAATTg                                                        <  1:3514912/71‑1 (MQ=255)
               gACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTg                                                        <  1:1976413/71‑1 (MQ=255)
                 cTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGa                                                       >  1:2579748/1‑70 (MQ=255)
                 cTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGa                                                       >  1:5769790/1‑70 (MQ=255)
                 cTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGa                                                       >  1:1989583/1‑70 (MQ=255)
                 cTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGa                                                       >  1:1950799/1‑70 (MQ=255)
                 cTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGa                                                       >  1:1759616/1‑70 (MQ=255)
                 cTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGa                                                       >  1:5067565/1‑70 (MQ=255)
                 cTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGAc                                                      <  1:6385041/71‑1 (MQ=255)
                     cccGGTCCGCATGTACATGCGCGAGATGGGCGCGGTGTCGCTGCTGACGCGCGAAGGCGagatcgagatcg                                                  <  1:3969888/71‑1 (MQ=255)
                     cccGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCg                                                  <  1:5218307/71‑1 (MQ=255)
                      ccGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGc                                                 >  1:5151413/1‑71 (MQ=255)
                        ggTCCGCATGTATATGCGCGAGATGGGCTCGGTCGAGCTGCTCACGCGCGAAGGCGAGATCGAGATCGcga                                               <  1:4710248/71‑3 (MQ=255)
                                 gTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTAtt                                       >  1:4802033/1‑70 (MQ=255)
                                  tACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGaaa                                                           >  1:112708/1‑49 (MQ=255)
                                  tACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGaaa                                                           >  1:1163336/1‑49 (MQ=255)
                                  tACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGaaa                                                           >  1:6827173/1‑49 (MQ=255)
                                  tACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGaaa                                                           >  1:5745828/1‑49 (MQ=255)
                                  tACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTg                                      <  1:5396024/70‑1 (MQ=255)
                                   aCATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTg                                      >  1:2275490/1‑69 (MQ=255)
                                   aCATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTg                                      >  1:550127/1‑69 (MQ=255)
                                   aCATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTg                                      >  1:1748459/1‑69 (MQ=255)
                                    cATGCGTGAAATGGGCACCGTTGAACTGTTGACCcgcg                                                                    <  1:2722520/38‑1 (MQ=255)
                                    cATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTa                                         >  1:3862035/1‑65 (MQ=255)
                                       gCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACg                                >  1:4387304/1‑71 (MQ=255)
                                       gCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACg                                >  1:3237163/1‑71 (MQ=255)
                                        cGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAg                                                                 >  1:5343440/1‑37 (MQ=255)
                                        cGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAg                                                                 >  1:1877126/1‑37 (MQ=255)
                                        cGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAg                                                                 >  1:2489841/1‑37 (MQ=255)
                                        cGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAg                                                                 >  1:594988/1‑37 (MQ=255)
                                         gTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGa                                                       <  1:2718488/46‑1 (MQ=255)
                                         gTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACggg                              <  1:4918071/71‑1 (MQ=255)
                                         gTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACggg                              <  1:4926132/71‑1 (MQ=255)
                                         gTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACgg                               >  1:5131990/1‑70 (MQ=255)
                                         gTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACgg                               >  1:2610697/1‑70 (MQ=255)
                                           gAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACa                                                     >  1:5451769/1‑46 (MQ=255)
                                           gAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACa                                                     >  1:4439339/1‑46 (MQ=255)
                                            aaaTGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGAt                            >  1:4077446/1‑70 (MQ=255)
                                            aaaTGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGAt                            >  1:2570542/1‑70 (MQ=255)
                                             aaTGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATc                           >  1:503257/1‑70 (MQ=255)
                                             aaTGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATc                           >  1:3640458/1‑70 (MQ=255)
                                             aaTGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATc                           >  1:6543466/1‑70 (MQ=255)
                                               tGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCa                          <  1:6633275/69‑1 (MQ=255)
                                                      cGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAg                                             <  1:680362/43‑1 (MQ=255)
                                                      cGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAg                                             <  1:94000/43‑1 (MQ=255)
                                                      cGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAg                                             <  1:941492/43‑1 (MQ=255)
                                                      cGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAg                                             <  1:2488474/43‑1 (MQ=255)
                                                      cGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAg                                             <  1:2907331/43‑1 (MQ=255)
                                                      cGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTc                 <  1:2521589/71‑1 (MQ=255)
                                                      cGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCATATTGAAGACGGGATCAACCAGGTTc                 <  1:4017084/71‑1 (MQ=255)
                                                       gTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCg                                                  >  1:1630734/1‑37 (MQ=255)
                                                       gTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGaa                                    >  1:4003481/1‑51 (MQ=255)
                                                       gTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGaa                                    >  1:217083/1‑51 (MQ=255)
                                                       gTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGaa                                    >  1:3639039/1‑51 (MQ=255)
                                                       gTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGaa                                    >  1:1570786/1‑51 (MQ=255)
                                                       gTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTc                 >  1:4283105/1‑70 (MQ=255)
                                                       gTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTc                 >  1:3469940/1‑70 (MQ=255)
                                                        ttGAACTGTTGACGCGCGAAGGCGAAATTGACATCGCTAAGCGTa                                         >  1:3018553/1‑45 (MQ=255)
                                                        ttGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTa                                         >  1:5389163/1‑45 (MQ=255)
                                                        ttGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTa                                         >  1:6113002/1‑45 (MQ=255)
                                                        ttGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTa                                         >  1:2264239/1‑45 (MQ=255)
                                                        ttGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTa                                         >  1:2869096/1‑45 (MQ=255)
                                                        ttGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTa                                         >  1:6889997/1‑45 (MQ=255)
                                                          gAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCa                      >  1:6033410/1‑62 (MQ=255)
                                                          gAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATg             >  1:341719/1‑71 (MQ=255)
                                                          gAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATg             >  1:1185318/1‑71 (MQ=255)
                                                            aCTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCg                                           >  1:6697374/1‑39 (MQ=255)
                                                            aCTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCg                                           >  1:3274140/1‑39 (MQ=255)
                                                            aCTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCg                                           >  1:1279028/1‑39 (MQ=255)
                                                            aCTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCg                                           >  1:6852131/1‑39 (MQ=255)
                                                            aCTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCg                                           >  1:6921258/1‑39 (MQ=255)
                                                             cTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCaa                         >  1:4637324/1‑56 (MQ=255)
                                                             cTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCaa                         >  1:1919442/1‑56 (MQ=255)
                                                             cTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCaa                         >  1:5951432/1‑56 (MQ=255)
                                                             cTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCaa                         >  1:5829399/1‑56 (MQ=255)
                                                             cTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCaa                         >  1:3000200/1‑56 (MQ=255)
                                                             cTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCaa                         >  1:3654748/1‑56 (MQ=255)
                                                             cTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCa                          >  1:2725767/1‑55 (MQ=255)
                                                             cTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTc          <  1:4384115/71‑1 (MQ=255)
                                                              tGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTAt                                        <  1:6214249/40‑1 (MQ=255)
                                                              tGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTAt                                        <  1:466031/40‑1 (MQ=255)
                                                              tGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCaa               >  1:6375583/1‑65 (MQ=255)
                                                              tGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCaa               >  1:6383277/1‑65 (MQ=255)
                                                              tGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCaa               >  1:3684364/1‑65 (MQ=255)
                                                              tGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTcc         >  1:6135116/1‑71 (MQ=255)
                                                               gTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCg        >  1:1074117/1‑71 (MQ=255)
                                                                ttGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGt       <  1:1000001/71‑1 (MQ=255)
                                                                  gACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGa                                  >  1:770767/1‑42 (MQ=255)
                                                                  gACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGa                                  >  1:5646005/1‑42 (MQ=255)
                                                                  gACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGtt      <  1:6564360/70‑1 (MQ=255)
                                                                   aCCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATc                           <  1:3916321/48‑1 (MQ=255)
                                                                     ccGCGAAGGCGAAATTACATCGCTAAGCGTATTGaa                                    <  1:5638773/36‑1 (MQ=38)
                                                                     |                                                                      
TGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTG  >  minE/2064026‑2064165

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