Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 526,903 C→T 40.7% L189L (CTG→TTG)  ybhL → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE526,9030CT40.7% 35.4 / 55.4 81L189L (CTG→TTG) ybhLpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (29/19);  new base T (7/26);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 6.09e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-03

GATGTGGGCAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAACTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTGG  >  minE/526836‑526970
                                                                   |                                                                   
gATGTGGGCAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAAc                                                                     >  1:1744631/1‑68 (MQ=255)
gATGTGGGCAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAAc                                                                     >  1:2864460/1‑68 (MQ=255)
      ggCAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAACTGaaaaa                                                              <  1:765075/69‑1 (MQ=255)
      ggCAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAACTGaaaaa                                                              <  1:2817405/69‑1 (MQ=255)
      ggCAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAACTGaaaaa                                                              <  1:2115454/69‑1 (MQ=255)
      ggCAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAACTGaaaaa                                                              <  1:2924665/69‑1 (MQ=255)
      ggCAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAACTGaaaaa                                                              <  1:94883/69‑1 (MQ=255)
      ggCAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAACTGaaaaa                                                              <  1:2689522/69‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:1950469/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:1767530/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:1861131/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:275258/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:716178/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:1946423/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:3006320/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:2562419/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:2459988/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:3393211/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:2450352/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:2759885/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:1558315/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:804850/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:147781/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:1394709/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:1308402/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:1282643/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:1159392/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:2930796/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:1140355/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:1082280/68‑1 (MQ=255)
         aGTTACCTACATCTGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:1333427/67‑1 (MQ=255)
         aGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:2322509/67‑1 (MQ=255)
             aCCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAACTGAAAAatat                                                           >  1:2856838/1‑65 (MQ=255)
             aCCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAACTGAAAAATATgg                                                         >  1:967635/1‑67 (MQ=255)
             aCCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAACTGAAAAATATgg                                                         >  1:1163332/1‑67 (MQ=255)
             aCCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAACTGAAAAATATgg                                                         >  1:734837/1‑67 (MQ=255)
                              tGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:3021773/46‑1 (MQ=255)
                               gTCTTTGTCGGATTGACGGCGTATGACACGCAGAAACTGAAAAATATGGGTGAGCAGATTGATACCcgc                                     <  1:1142792/69‑1 (MQ=255)
                               gTCTTTGTCGGATTGACGGCGTATGACACGCAGAAACTGAAAAATATGGGTGAGCAGATTGATACCcgc                                     <  1:847750/69‑1 (MQ=255)
                                tCTTTGTCGGATTGACGGCGTATGACACGCAGAAACTGAAAAATATGGGTGAGCAGATTGATACCcgc                                     <  1:3072082/68‑1 (MQ=255)
                                    tGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAATATGGGTGAGCAGATTGATAccc                                       >  1:2447382/1‑62 (MQ=255)
                                    tGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAATATGGGTGAGCAGATTGATAccc                                       >  1:359912/1‑62 (MQ=255)
                                    tGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAATATGGGTGAGCAGATTGATAccc                                       >  1:2062889/1‑62 (MQ=255)
                                    tGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAATATGGGTGAGCAGATTGATAccc                                       >  1:3095085/1‑62 (MQ=255)
                                    tGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAATATGGGTGAGCAGATTGATAccc                                       >  1:1656880/1‑62 (MQ=255)
                                    tGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAATATGGGTGAGCAGATTGATAccc                                       >  1:1496092/1‑62 (MQ=255)
                                    tGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAATATGGGTGAGCAGATTGATAccc                                       >  1:3011769/1‑62 (MQ=255)
                                    tGTCGGATGGACGGCGTATGACACTCAGAAATTGAAAAat                                                             <  1:3062494/40‑1 (MQ=255)
                                          aTTGACGGCGTATGACACGCAGAAACTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGaa                          <  1:3096986/69‑1 (MQ=255)
                                          aTTGACGGCGTATGACACGCAGAAACTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGaa                          <  1:1261964/69‑1 (MQ=255)
                                          aTTGACGGCGTATGACACGCAGAAACTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGaa                          <  1:2765548/69‑1 (MQ=255)
                                          aTTGACGGCGTATGACACGCAGAAACTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGaa                          <  1:1893028/69‑1 (MQ=255)
                                          aTTGACGGCGTATGACACGCAGAAACTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGa                           <  1:2615663/68‑1 (MQ=255)
                                           ttGACGGCGTATGACACGCAGAAACTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGaa                          <  1:321709/68‑1 (MQ=255)
                                                  cGTATGACACGCAGAAACTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCa                  <  1:3372808/69‑1 (MQ=255)
                                                  cGTATGACACGCAGAAACTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCa                  <  1:2610239/69‑1 (MQ=255)
                                                  cGTATGACACGCAGAAACTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCa                  <  1:2649658/69‑1 (MQ=255)
                                                  cGTATGACACGCAGAAACTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCa                  <  1:2599867/69‑1 (MQ=255)
                                                                 aaCTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCtt    >  1:2193616/1‑68 (MQ=255)
                                                                 aaCTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTg   >  1:2211115/1‑69 (MQ=255)
                                                                 aaCTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTg   >  1:1063602/1‑69 (MQ=255)
                                                                 aaCTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTg   >  1:1212377/1‑69 (MQ=255)
                                                                 aaCTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTg   >  1:1352326/1‑69 (MQ=255)
                                                                 aaCTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTg   >  1:1910581/1‑69 (MQ=255)
                                                                 aaCTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTg   >  1:1913728/1‑69 (MQ=255)
                                                                 aaCTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTg   >  1:597463/1‑69 (MQ=255)
                                                                 aaCTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTg   >  1:458727/1‑69 (MQ=255)
                                                                 aaCTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTg   >  1:2240731/1‑69 (MQ=255)
                                                                 aaCTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTg   >  1:2790937/1‑69 (MQ=255)
                                                                 aaCTGAAAAATAAGGGTGAGCAGATTGATACCcgcg                                    >  1:936057/1‑36 (MQ=255)
                                                                  aCTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTgg  >  1:3241690/1‑69 (MQ=255)
                                                                  aCTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTgg  >  1:3208758/1‑69 (MQ=255)
                                                                  aCTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTgg  >  1:3431636/1‑69 (MQ=255)
                                                                  aCTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTgg  >  1:2173964/1‑69 (MQ=255)
                                                                  aCTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTgg  >  1:3207406/1‑69 (MQ=255)
                                                                  aCTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTgg  >  1:3110920/1‑69 (MQ=255)
                                                                  aCTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTgg  >  1:2412149/1‑69 (MQ=255)
                                                                  aCTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTgg  >  1:1024332/1‑69 (MQ=255)
                                                                  aCTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTgg  >  1:1485769/1‑69 (MQ=255)
                                                                  aCTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTgg  >  1:2678751/1‑69 (MQ=255)
                                                                  aCTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTgg  >  1:2931309/1‑69 (MQ=255)
                                                                  aCTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTg   >  1:2938290/1‑68 (MQ=255)
                                                                  aCTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTg   >  1:1514647/1‑68 (MQ=255)
                                                                   |                                                                   
GATGTGGGCAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAACTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTGG  >  minE/526836‑526970

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: