Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 526,896 G→T 55.9% T186T (ACG→ACT ybhL → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE526,8960GT55.9% 8.0 / 60.7 59T186T (ACG→ACTybhLpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (7/19);  new base T (7/26);  total (14/45)
Fisher's exact test for biased strand distribution p-value = 7.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGTGGGCAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAACTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCA  >  minE/526836‑526954
                                                            |                                                          
gATGTGGGCAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGaaa                                                      >  1:1920924/1‑67 (MQ=255)
gATGTGGGCAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAAc                                                     >  1:2864460/1‑68 (MQ=255)
gATGTGGGCAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAAc                                                     >  1:1744631/1‑68 (MQ=255)
      ggCAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAACTGaaaaa                                              <  1:2817405/69‑1 (MQ=255)
      ggCAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAACTGaaaaa                                              <  1:2689522/69‑1 (MQ=255)
      ggCAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAACTGaaaaa                                              <  1:2924665/69‑1 (MQ=255)
      ggCAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAACTGaaaaa                                              <  1:2115454/69‑1 (MQ=255)
      ggCAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAACTGaaaaa                                              <  1:765075/69‑1 (MQ=255)
      ggCAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAACTGaaaaa                                              <  1:94883/69‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                             <  1:1394709/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                             <  1:2759885/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                             <  1:275258/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                             <  1:1140355/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                             <  1:1082280/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                             <  1:1159392/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                             <  1:2930796/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                             <  1:2562419/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                             <  1:2459988/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                             <  1:2450352/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                             <  1:3006320/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                             <  1:3393211/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                             <  1:1950469/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                             <  1:1946423/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                             <  1:804850/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                             <  1:716178/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                             <  1:1861131/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                             <  1:1767530/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                             <  1:1282643/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                             <  1:1558315/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                             <  1:1308402/68‑1 (MQ=255)
        cAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAat                                             <  1:147781/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:734837/1‑67 (MQ=255)
             aCCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAACTGAAAAATATgg                                         >  1:967635/1‑67 (MQ=255)
             aCCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAACTGAAAAATATgg                                         >  1:1163332/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:3011769/1‑62 (MQ=255)
                                    tGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAATATGGGTGAGCAGATTGATAccc                       >  1:3095085/1‑62 (MQ=255)
                                    tGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAATATGGGTGAGCAGATTGATAccc                       >  1:1496092/1‑62 (MQ=255)
                                    tGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAATATGGGTGAGCAGATTGATAccc                       >  1:1656880/1‑62 (MQ=255)
                                    tGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAATATGGGTGAGCAGATTGATAccc                       >  1:2062889/1‑62 (MQ=255)
                                    tGTCGGATTGACGGCGTATGACACTCAGAAATTGAAAAATATGGGTGAGCAGATTGATAccc                       >  1:359912/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:1893028/69‑1 (MQ=255)
                                          aTTGACGGCGTATGACACGCAGAAACTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGaa          <  1:2765548/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:2599867/69‑1 (MQ=255)
                                                  cGTATGACACGCAGAAACTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCa  <  1:2610239/69‑1 (MQ=255)
                                                  cGTATGACACGCAGAAACTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCa  <  1:2649658/69‑1 (MQ=255)
                                                            |                                                          
GATGTGGGCAGTTACCTACATCGGCGTGATTGTCTTTGTCGGATTGACGGCGTATGACACGCAGAAACTGAAAAATATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCA  >  minE/526836‑526954

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: