Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 760,905 C→T 39.0% V26I (GTC→ATC)  ymfA ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE760,9050CT39.0% 45.0 / 57.3 87V26I (GTC→ATC) ymfApredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (29/24);  new base T (20/14);  total (49/38)
Fisher's exact test for biased strand distribution p-value = 8.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

GGTATCAACAATGACGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTT  >  minE/760840‑760964
                                                                 |                                                           
ggTATCAACAATGACGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGacaac                                                          >  1:1781249/1‑69 (MQ=255)
ggTATCAACAATGACGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGacaac                                                          >  1:2771751/1‑69 (MQ=255)
ggTATCAACAATGACGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGacaac                                                          >  1:54468/1‑69 (MQ=255)
      aacaaTGACGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGa                                                    <  1:3137884/69‑1 (MQ=255)
      aacaaTGACGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGa                                                    <  1:970202/69‑1 (MQ=255)
      aacaaTGACGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGa                                                    <  1:1973984/69‑1 (MQ=255)
        caaTGACGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGaa                                                   <  1:1754103/68‑1 (MQ=255)
        caaTGACGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGaa                                                   <  1:1900879/68‑1 (MQ=255)
        caaTGACGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGaa                                                   <  1:2105497/68‑1 (MQ=255)
             aCGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGATCTGGTTATAAAAGACAACAACAGAAATGAgg                                             <  1:2548116/69‑1 (MQ=255)
             aCGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAgt                                             >  1:425696/1‑68 (MQ=255)
             aCGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAgg                                             >  1:2512062/1‑69 (MQ=255)
             aCGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAgg                                             >  1:2270114/1‑69 (MQ=255)
             aCGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAgg                                             >  1:214914/1‑69 (MQ=255)
             aCGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAgg                                             <  1:373526/69‑1 (MQ=255)
             aCGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAgg                                             <  1:632649/69‑1 (MQ=255)
             aCGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAgg                                             <  1:259726/69‑1 (MQ=255)
             aCGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAgg                                             >  1:2824203/1‑69 (MQ=255)
             aCGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAgg                                             >  1:1012484/1‑69 (MQ=255)
             aCGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAgg                                             <  1:1708014/69‑1 (MQ=255)
             aCGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAgg                                             >  1:853210/1‑69 (MQ=255)
             aCGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAgg                                             <  1:1622480/69‑1 (MQ=255)
             aCGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAgg                                             <  1:1720283/69‑1 (MQ=255)
             aCGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAgg                                             >  1:1493153/1‑69 (MQ=255)
             aCGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAg                                              >  1:2233059/1‑68 (MQ=255)
                      gTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTAc                                     <  1:2827708/68‑1 (MQ=255)
                      gTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTAc                                     <  1:2073648/68‑1 (MQ=255)
                             tCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGa                             <  1:1255668/69‑1 (MQ=255)
                             tCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGa                             <  1:326120/69‑1 (MQ=255)
                             tCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGa                             <  1:2567284/69‑1 (MQ=255)
                             tCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGa                             <  1:907896/69‑1 (MQ=255)
                             tCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGa                             <  1:1484659/69‑1 (MQ=255)
                             tCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGa                             <  1:2343063/69‑1 (MQ=255)
                             tCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGa                             <  1:294047/69‑1 (MQ=255)
                             tCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGa                             <  1:297928/69‑1 (MQ=255)
                             tCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGa                             <  1:3177548/69‑1 (MQ=255)
                             tCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGa                             <  1:755740/69‑1 (MQ=255)
                             tCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGa                             <  1:1652331/69‑1 (MQ=255)
                             tCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGa                             <  1:1916973/69‑1 (MQ=255)
                              cgcgTCTTTGCTATAGGTGAACTGGTTATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGa                             <  1:2818738/68‑1 (MQ=255)
                              cgcgTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGa                                                    <  1:2099993/45‑1 (MQ=255)
                              cgcgTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGAcc                               >  1:2217852/1‑66 (MQ=255)
                              cgcgTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGAcc                               >  1:1800411/1‑66 (MQ=255)
                              cgcgTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGAcc                               >  1:810911/1‑66 (MQ=255)
                              cgcgTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGAcc                               >  1:3046032/1‑66 (MQ=255)
                              cgcgTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGAcc                               >  1:1701764/1‑66 (MQ=255)
                              cgcgTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGAcc                               >  1:1466482/1‑66 (MQ=255)
                              cgcgTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGAcc                               >  1:1698003/1‑66 (MQ=255)
                              cgcgTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGAcc                               >  1:1642101/1‑66 (MQ=255)
                                cgTCTTTGCTATAGGTGAACTGGTTATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGa                             <  1:3173417/66‑1 (MQ=255)
                                    tttGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGa                             <  1:719770/62‑1 (MQ=255)
                                              gTGAACTGGTTATAAAAGATAACAACAGAAATGAGGATCAGTACCAGAc                                >  1:2497700/1‑49 (MQ=255)
                                              gTGAACTGGTTATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAAGAAtata               >  1:284785/1‑66 (MQ=255)
                                              gTGAACTGGTTATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAAGAATATAc              >  1:2478266/1‑67 (MQ=255)
                                              gTGAACTGGTTATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAAGAATATAc              >  1:2813336/1‑67 (MQ=255)
                                              gTGAACTGGTTATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAAGAATATAc              >  1:2011592/1‑67 (MQ=255)
                                              gTGAACTGGTTATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAAGAATATAc              >  1:1209451/1‑67 (MQ=255)
                                              gTGAACTGGTTATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAAGAATATAc              >  1:1830459/1‑67 (MQ=255)
                                              gTGAACTGGTTATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAAGAATATAc              >  1:90706/1‑67 (MQ=255)
                                                 aaCTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTa                                      >  1:3145005/1‑40 (MQ=255)
                                                 aaCTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTAc                                     >  1:963095/1‑41 (MQ=255)
                                                 aaCTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTAc                                     >  1:517193/1‑41 (MQ=255)
                                                 aaCTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTAc                                     >  1:1041162/1‑41 (MQ=255)
                                                 aaCTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTAc                                     >  1:525352/1‑41 (MQ=255)
                                                 aaCTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTAc                                     >  1:3086519/1‑41 (MQ=255)
                                                 aaCTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTAc                                     >  1:2093221/1‑41 (MQ=255)
                                                 aaCTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTAc                                     >  1:2149336/1‑41 (MQ=255)
                                                 aaCTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTAc                                     >  1:241652/1‑41 (MQ=255)
                                                   cTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACgaaagaa       <  1:1210041/69‑1 (MQ=255)
                                                   cTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACgaaagaa       >  1:1290738/1‑69 (MQ=255)
                                                   cTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACgaaagaa       >  1:828274/1‑69 (MQ=255)
                                                   cTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACgaaagaa       <  1:2527273/69‑1 (MQ=255)
                                                      gTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAAc    >  1:708029/1‑69 (MQ=255)
                                                       ttATAAAAGATAACAACAGAAATGAGGATCAGTAc                                     >  1:1921693/1‑35 (MQ=255)
                                                       ttATAAAAGATAACAACAGAAATGAGGATCAGTACCATACCGATACCGAAGAATATACGAAAGAAAAc    >  1:2651397/1‑68 (MQ=255)
                                                       ttATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGAt                            >  1:1883052/1‑44 (MQ=255)
                                                       ttATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAAGAATATACGAAAGAAAAc    >  1:183070/1‑68 (MQ=255)
                                                       ttATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAAGAATATACGAAAGAAAAc    >  1:3026882/1‑68 (MQ=255)
                                                       ttATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAAGAATATACGAAAGAAAAc    >  1:925101/1‑68 (MQ=255)
                                                       ttATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAAGAATATACGAAAGAAAAc    >  1:1496993/1‑68 (MQ=255)
                                                       ttATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAAGAATATACGAAAGAAAAc    >  1:3216115/1‑68 (MQ=255)
                                                       ttATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAAGAATATACGAAAGAAAAc    >  1:3368655/1‑68 (MQ=255)
                                                       ttATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAAGAATATACGAAAGAAAAc    >  1:3279513/1‑68 (MQ=255)
                                                       ttATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAAGAATATACGAAAGAAAAc    >  1:2514642/1‑68 (MQ=255)
                                                       ttATAAAAGATAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAAGAATATACGAAAGAAAAc    >  1:415337/1‑68 (MQ=255)
                                                        tataAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGTAATATACGAAAGAAAACtt  <  1:3177710/69‑1 (MQ=255)
                                                        tataAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACtt  <  1:1820609/69‑1 (MQ=255)
                                                        tataAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACtt  <  1:322086/69‑1 (MQ=255)
                                                        tataAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACtt  <  1:1546582/69‑1 (MQ=255)
                                                        tataAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACtt  <  1:1278866/69‑1 (MQ=255)
                                                                 |                                                           
GGTATCAACAATGACGCCTTCGGTATGAATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTT  >  minE/760840‑760964

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: