Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,145,423 A→G 58.9% R428R (CGT→CGC yhbX ← predicted hydrolase, inner membrane

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,145,4230AG58.9% 26.6 / 46.2 68R428R (CGT→CGCyhbXpredicted hydrolase, inner membrane
Reads supporting (aligned to +/- strand):  ref base A (12/16);  new base G (18/22);  total (30/38)
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

GGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTG  >  minE/2145355‑2145491
                                                                    |                                                                    
gggACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCa                                                                      >  1:2338051/1‑69 (MQ=255)
gggACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCa                                                                      >  1:1445993/1‑69 (MQ=255)
    cATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGt                                                                   >  1:1334688/1‑68 (MQ=255)
        atatGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAg                                                              >  1:1520177/1‑69 (MQ=255)
             ccTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGTTCACGTTCCAGGCCGt                                                         >  1:1939936/1‑69 (MQ=255)
             ccTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCg                                                          >  1:2517230/1‑68 (MQ=255)
                          cTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGtt                                                                  >  1:2978084/1‑47 (MQ=255)
                          cTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATa                                            >  1:1535790/1‑69 (MQ=255)
                          cTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATa                                            >  1:3373193/1‑69 (MQ=255)
                          cTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATa                                            >  1:496754/1‑69 (MQ=255)
                               tGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATg                                        <  1:3178175/68‑1 (MQ=255)
                               tGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATg                                        <  1:811022/68‑1 (MQ=255)
                               tGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATg                                        <  1:1311972/68‑1 (MQ=255)
                               tGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATg                                        <  1:2824079/68‑1 (MQ=255)
                               tGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATg                                        <  1:175037/68‑1 (MQ=255)
                               tGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATg                                        <  1:2029968/68‑1 (MQ=255)
                               tGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATg                                        <  1:243379/68‑1 (MQ=255)
                               tGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATg                                        <  1:2954537/68‑1 (MQ=255)
                                        gAAAATAGACGTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTg                                                        >  1:3187981/1‑43 (MQ=255)
                                        gAAAATAGACGTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGt                                                      >  1:3092807/1‑45 (MQ=255)
                                        gAAAATAGACGTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGt                                                      >  1:2930753/1‑45 (MQ=255)
                                        gAAAATAGACGTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGt                                                      >  1:206927/1‑45 (MQ=255)
                                        gAAAATAGACGTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGt                                                      >  1:2275118/1‑45 (MQ=255)
                                        gAAAATAGACGTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGt                                                      >  1:2029014/1‑45 (MQ=255)
                                        gAAAATAGACGTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGt                                                      >  1:2373793/1‑45 (MQ=255)
                                        gAAAATAGACGTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGt                                                      >  1:1928378/1‑45 (MQ=255)
                                        gAAAATAGACGTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGt                                                      >  1:2441623/1‑45 (MQ=255)
                                        gAAAATAGACGTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGt                                                      >  1:1882856/1‑45 (MQ=255)
                                        gAAAATAGACGTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGt                                                      >  1:1752293/1‑45 (MQ=255)
                                        gAAAATAGACGTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGt                                                      >  1:160833/1‑45 (MQ=255)
                                        gAAAATAGACGTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGt                                                      >  1:1515577/1‑45 (MQ=255)
                                        gAAAATAGACGTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGt                                                      >  1:2626975/1‑45 (MQ=255)
                                        gAAAATAGACGTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGt                                                      >  1:953134/1‑45 (MQ=255)
                                        gAAAATAGACGTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGt                                                      >  1:96534/1‑45 (MQ=255)
                                        gAAAATAGACGTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGt                                                      >  1:1098259/1‑45 (MQ=255)
                                        gAAAATAGACGTTCTTCTTCGTAGGGTCGCGTTACAGGCCGTg                                                        >  1:2139477/1‑43 (MQ=255)
                                             tAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGa                         >  1:1567792/1‑69 (MQ=255)
                                             tAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGa                         >  1:981637/1‑69 (MQ=255)
                                                  gTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGTCGGCa                                                 <  1:2425720/40‑1 (MQ=255)
                                                  gTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGTCGGCa                                                 <  1:2782737/40‑1 (MQ=255)
                                                  gTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGTCGGCa                                                 <  1:380538/40‑1 (MQ=255)
                                                  gTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGTCGGCa                                                 <  1:432975/40‑1 (MQ=255)
                                                  gTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGTCGGCa                                                 <  1:802504/40‑1 (MQ=255)
                                                  gTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGTCGGCa                                                 <  1:807696/40‑1 (MQ=255)
                                                  gTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGTCGGCa                                                 <  1:863514/40‑1 (MQ=255)
                                                  gTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGTCGGCa                                                 <  1:978187/40‑1 (MQ=255)
                                                  gTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGTCGGCa                                                 <  1:1966949/40‑1 (MQ=255)
                                                  gTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGTCGGCa                                                 <  1:1080788/40‑1 (MQ=255)
                                                  gTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGTCGGCa                                                 <  1:1084859/40‑1 (MQ=255)
                                                  gTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGTCGGCa                                                 <  1:1250179/40‑1 (MQ=255)
                                                  gTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGTCGGCa                                                 <  1:1468966/40‑1 (MQ=255)
                                                  gTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGTCGGCa                                                 <  1:1631719/40‑1 (MQ=255)
                                                  gTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGTCGGCa                                                 <  1:1645407/40‑1 (MQ=255)
                                                  gTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGTCGGCa                                                 <  1:1655453/40‑1 (MQ=255)
                                                  gTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGTCGGCa                                                 <  1:1918762/40‑1 (MQ=255)
                                                  gTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGTCGGCa                                                 <  1:2061929/40‑1 (MQ=255)
                                                  gTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGTCGGCa                                                 <  1:22104/40‑1 (MQ=255)
                                                  gTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGTCGGCa                                                 <  1:2232075/40‑1 (MQ=255)
                                                  gTTCTTCTTCGTAGGGTCGCGTTCCAGGCCGTGGTCGGCa                                                 <  1:2217809/40‑1 (MQ=255)
                                                  gTTCTTCTTCGTAGGGTCGCGATCCAGGCCGTGGTCGGCa                                                 <  1:41026/40‑1 (MQ=38)
                                                     cttcttcGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCtttt                   <  1:2806559/67‑1 (MQ=255)
                                                     cttcttcGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCtttt                   <  1:2814283/67‑1 (MQ=255)
                                                       tcttcGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTtaat                <  1:2912348/68‑1 (MQ=255)
                                                       tcttcGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTtaat                <  1:192631/68‑1 (MQ=255)
                                                       tcttcGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCCTTtaat                <  1:179079/68‑1 (MQ=255)
                                                       tcttcGTAGGGTCACGTTCCAGGCCGGGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTtaat                <  1:3349559/68‑1 (MQ=255)
                                                        cttcGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCtttt                   <  1:2818164/64‑1 (MQ=255)
                                                        cttcGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTtaat                <  1:3414376/67‑1 (MQ=255)
                                                                  tCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAAc     >  1:291103/1‑68 (MQ=255)
                                                                  tCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAAc     >  1:1224422/1‑68 (MQ=255)
                                                                    aCGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTg  <  1:3424088/69‑1 (MQ=255)
                                                                    aCGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTg  <  1:2403227/69‑1 (MQ=255)
                                                                    aCGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTg  <  1:2737431/69‑1 (MQ=255)
                                                                    aCGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTg  <  1:554980/69‑1 (MQ=255)
                                                                    aCGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTg  <  1:652235/69‑1 (MQ=255)
                                                                    aCGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTg  <  1:105394/69‑1 (MQ=255)
                                                                    aCGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTg  <  1:2475674/69‑1 (MQ=255)
                                                                    aCGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTg  <  1:306468/69‑1 (MQ=255)
                                                                    aCGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTg  <  1:3108163/69‑1 (MQ=255)
                                                                    |                                                                    
GGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTG  >  minE/2145355‑2145491

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