Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 703,823 C→G 5.4% R274R (CGG→CGC ymcA ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE703,8230CG 5.4% 228.8 / 2.4 74R274R (CGG→CGCymcAconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (47/23);  new base G (2/2);  total (49/25)
Fisher's exact test for biased strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

TCAGGCGCAGGTTGATATGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTG  >  minE/703755‑703889
                                                                    |                                                                  
tCAGGCGCAGGTTGATATGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCg                                                                  >  1:3250821/1‑71 (MQ=255)
tCAGGCGCAGGTTGATATGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCg                                                                  >  1:3616751/1‑71 (MQ=255)
tCAGGCGCAGGTTGATATGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCg                                                                  >  1:2219542/1‑71 (MQ=255)
tCAGGCGCAGGTTGATATGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCg                                                                  >  1:1480581/1‑71 (MQ=255)
tCAGGCGCAGGTTGATATGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCg                                                                  >  1:4865045/1‑71 (MQ=255)
 cAGGCGCAGGTTGATATGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCg                                                                  >  1:2775724/1‑70 (MQ=255)
   ggCGCAGGTTGATATGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTaa                                                               >  1:4913628/1‑71 (MQ=255)
       cAGGTTGATATGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTg                                                           <  1:5627962/71‑1 (MQ=255)
       cAGGTTGATATGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTg                                                           <  1:2158183/71‑1 (MQ=255)
        aGGTTGATATGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTg                                                           >  1:1840702/1‑70 (MQ=255)
          gTTGATATGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCg                                                                  <  1:4982590/61‑1 (MQ=255)
          gTTGATATGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCg                                                                  <  1:1885572/61‑1 (MQ=255)
          gTTGATATGCCAGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCg                                                                  <  1:4358666/61‑1 (MQ=255)
              atatGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCa                                                     <  1:4355681/70‑1 (MQ=255)
              atatGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCa                                                     <  1:877951/70‑1 (MQ=255)
                atGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCaaa                                                   >  1:2177605/1‑70 (MQ=255)
                atGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCaaa                                                   >  1:5579130/1‑70 (MQ=255)
                   ccGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTg                                                           <  1:4099818/59‑1 (MQ=255)
                   ccGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTg                                                           <  1:394051/59‑1 (MQ=255)
                    cGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATc                                               >  1:1110835/1‑70 (MQ=255)
                    cGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATc                                               >  1:999812/1‑70 (MQ=255)
                    cGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCa                                              <  1:1662944/71‑1 (MQ=255)
                     ggTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCaa                                             >  1:109960/1‑71 (MQ=255)
                     ggTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCaa                                             >  1:3381006/1‑71 (MQ=255)
                        tttGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGa                                                          >  1:4840384/1‑55 (MQ=255)
                        tttGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGa                                                          >  1:2501632/1‑55 (MQ=255)
                            gCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACg                                           <  1:5654407/66‑1 (MQ=255)
                             cGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTAc                                     <  1:2031798/71‑1 (MQ=255)
                               tATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGc                                                       <  1:2895887/51‑1 (MQ=255)
                               tATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACgt                                          <  1:2692703/64‑1 (MQ=255)
                               tATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACgt                                          <  1:2771650/64‑1 (MQ=255)
                                 tCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTaa                                                               >  1:5331475/1‑41 (MQ=255)
                                 tCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACcgcg                                 >  1:5338065/1‑71 (MQ=255)
                                   gcgcAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGtt                               <  1:4107375/71‑1 (MQ=255)
                                   gcgcAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGtt                               <  1:2274816/71‑1 (MQ=255)
                                      cagggcagggCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCAt                            <  1:4726849/71‑1 (MQ=255)
                                       agggcagggcGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCAt                            <  1:1171639/70‑1 (MQ=255)
                                        gggcagggcGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATaa                          >  1:1433550/1‑71 (MQ=255)
                                        gggcagggcGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATaa                          >  1:1789026/1‑71 (MQ=255)
                                          gacgggcGCAGATCGTTAAAGTTGGTGCGCAACGTAACGCCAAACATAAAGGTGTTGCCACGTTCATAGCt                        <  1:1149429/68‑1 (MQ=255)
                                           cagggcGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACt                        >  1:5634674/1‑70 (MQ=255)
                                           acgggcGCAGATCGTTAAAGTTGGTGCGCAACGTAACGCCAAACATAAAGgtg                                         >  1:3028354/3‑53 (MQ=255)
                                            agggcGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACt                        >  1:2925019/1‑69 (MQ=255)
                                             gggcGCAGATCGTTAAAGTTGGTGCGCAACGTAACGCCAAACATaaagg                                           <  1:3648392/49‑3 (MQ=255)
                                               gcgcAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTAcc                                    >  1:3580591/1‑54 (MQ=255)
                                               gcgcAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTAcc                                    >  1:2101668/1‑54 (MQ=255)
                                               gcgcAGATCGTTAAAGTTGGTGCGCAACGTAACGCCAAACATAAAGGTGTTGCCACGTTCATAGCTAAGGt                   >  1:5142751/1‑71 (MQ=255)
                                                 gcAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTg                 <  1:4706880/71‑1 (MQ=255)
                                                  cAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACgtg                                         <  1:3158008/46‑1 (MQ=255)
                                                   aGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTgag               <  1:3768694/71‑1 (MQ=255)
                                                    gATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGtt                                       <  1:2110083/46‑1 (MQ=255)
                                                    gATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGtt                  >  1:1829395/1‑67 (MQ=255)
                                                    gATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGtt                  >  1:2571267/1‑67 (MQ=255)
                                                     aTCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGtt                  >  1:185409/1‑66 (MQ=255)
                                                     aTCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACAGCGTTCATAACTCAGGtt                  >  1:5658889/1‑66 (MQ=255)
                                                       cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCaa                                             >  1:1864379/1‑37 (MQ=255)
                                                       cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCaa                                             >  1:628646/1‑37 (MQ=255)
                                                       cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAAc                                            >  1:921970/1‑38 (MQ=255)
                                                       cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAAc                                            >  1:3437851/1‑38 (MQ=255)
                                                       cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAgg                    >  1:2300236/1‑62 (MQ=255)
                                                       cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAgg                    >  1:1087217/1‑62 (MQ=255)
                                                       cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAgg                    >  1:80547/1‑62 (MQ=255)
                                                       cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAgc                    >  1:76829/1‑61 (MQ=255)
                                                       cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCt           >  1:4514070/1‑71 (MQ=255)
                                                        gTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTga                >  1:3238848/1‑65 (MQ=255)
                                                        gTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTga                >  1:4265829/1‑65 (MQ=255)
                                                          tGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGcc        >  1:3591559/1‑71 (MQ=255)
                                                          tGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGcc        >  1:3713938/1‑71 (MQ=255)
                                                           gAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTga                >  1:5049547/1‑62 (MQ=255)
                                                           gAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTga                >  1:644927/1‑62 (MQ=255)
                                                           gAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTga                >  1:2545410/1‑62 (MQ=255)
                                                                tGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCAt                            >  1:2886105/1‑45 (MQ=255)
                                                                tGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCAt                            >  1:3567367/1‑45 (MQ=255)
                                                                tGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTg  >  1:3031926/1‑71 (MQ=255)
                                                                    |                                                                  
TCAGGCGCAGGTTGATATGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTG  >  minE/703755‑703889

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: