Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 318,136 A→G 35.0% A1044A (GCA→GCG kefA → fused mechanosensitive channel proteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE318,1360AG35.0% 49.9 / 43.9 77A1044A (GCA→GCGkefAfused mechanosensitive channel proteins
Reads supporting (aligned to +/- strand):  ref base A (42/8);  new base G (9/18);  total (51/26)
Fisher's exact test for biased strand distribution p-value = 1.24e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGA  >  minE/318072‑318192
                                                                |                                                        
tGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCAttt                                                       <  1:2349079/68‑1 (MQ=255)
tGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCAttt                                                       <  1:1619706/68‑1 (MQ=255)
       ggcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAg                                               >  1:3388287/1‑69 (MQ=255)
       ggcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAg                                               >  1:2446637/1‑69 (MQ=255)
                     cAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGt                                          >  1:546831/1‑60 (MQ=255)
                     cAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCAt                                  >  1:2467170/1‑68 (MQ=255)
                     cAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCAt                                  >  1:3371087/1‑68 (MQ=255)
                     cAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCAt                                  >  1:1109070/1‑68 (MQ=255)
                     cAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCAt                                  >  1:620774/1‑68 (MQ=255)
                     cAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATg                                 >  1:1139158/1‑69 (MQ=255)
                     cAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATg                                 >  1:485447/1‑69 (MQ=255)
                     cAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATg                                 >  1:2864621/1‑69 (MQ=255)
                     cAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTACGGCATTTGGTGCCAGCACGTTGGATCATGa                                >  1:3141819/1‑69 (MQ=255)
                     cAAGGGTGATGCACGAACCAATGCCGGAAGTATTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATg                                 >  1:2093871/1‑69 (MQ=255)
                          gTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGc                              >  1:2524855/1‑67 (MQ=255)
                          gTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCt                             >  1:798673/1‑68 (MQ=255)
                          gTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTg                            >  1:882343/1‑69 (MQ=255)
                          gTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTg                            >  1:1671953/1‑69 (MQ=255)
                          gTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTg                            >  1:1227181/1‑69 (MQ=255)
                          gTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTg                            >  1:458380/1‑69 (MQ=255)
                                cACGAACCAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTg                      <  1:3172362/69‑1 (MQ=255)
                                cACGAACCAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTg                      <  1:2613415/69‑1 (MQ=255)
                                cACGAACCAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTg                      <  1:507087/69‑1 (MQ=255)
                                cACGAACCAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTg                      <  1:1227523/69‑1 (MQ=255)
                                 aCGAACCAATGCGGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTg                      <  1:2145059/68‑1 (MQ=255)
                                  cGAACCAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTa                    >  1:2799488/1‑69 (MQ=255)
                                  cGAACCAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTa                    >  1:2102143/1‑69 (MQ=255)
                                  cGAACCAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTa                    >  1:2075036/1‑69 (MQ=255)
                                  cGAACCAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTa                    >  1:2218161/1‑69 (MQ=255)
                                  cGAACCAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTa                    >  1:192268/1‑69 (MQ=255)
                                  cGAACCAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTa                    >  1:3435936/1‑69 (MQ=255)
                                  cGAACCAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTa                    >  1:2312442/1‑69 (MQ=255)
                                  cGAACCAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTa                    >  1:1660837/1‑69 (MQ=255)
                                  cGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATg                                 <  1:1528429/56‑1 (MQ=255)
                                  cGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTa                    <  1:811486/69‑1 (MQ=255)
                                  cGAACAAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTa                    >  1:881101/1‑69 (MQ=255)
                                       cAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTAt                   <  1:2937981/65‑1 (MQ=255)
                                       cAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTAt                   <  1:2555346/65‑1 (MQ=255)
                                       cAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTAt                   <  1:2283798/65‑1 (MQ=255)
                                       cAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTAt                   <  1:2420032/65‑1 (MQ=255)
                                       cAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTAt                   <  1:1067407/65‑1 (MQ=255)
                                       cAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTAt                   <  1:1387648/65‑1 (MQ=255)
                                       cAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTAt                   <  1:1781147/65‑1 (MQ=255)
                                       cAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTAt                   <  1:1881160/65‑1 (MQ=255)
                                       cAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTAt                   <  1:2104157/65‑1 (MQ=255)
                                       cAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTAt                   <  1:559837/65‑1 (MQ=255)
                                       cAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGAAt                   <  1:2747754/65‑1 (MQ=255)
                                       cAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGAGCTGCGCCTGTAt                   <  1:437138/65‑1 (MQ=39)
                                       cAATGCCGGAAGTCTTCTTTACGGCGTCTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTAt                   <  1:439957/65‑1 (MQ=255)
                                          tGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTg            <  1:1409641/69‑1 (MQ=255)
                                          tGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTg            <  1:3295974/69‑1 (MQ=255)
                                          tGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTg            <  1:2482647/69‑1 (MQ=255)
                                            ccGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGaa          >  1:1530541/1‑69 (MQ=255)
                                            ccGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGaa          >  1:3254205/1‑69 (MQ=255)
                                                   tcttctTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTa                    <  1:2850216/52‑1 (MQ=255)
                                                    cttcttTACGGCATTTGGTGCCAGCACGTTGGATCATg                                 >  1:2129307/1‑38 (MQ=255)
                                                    cttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTAtgt                 >  1:3139342/1‑54 (MQ=255)
                                                    cttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTg   >  1:2361960/1‑68 (MQ=255)
                                                    cttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGa  >  1:619079/1‑69 (MQ=255)
                                                    cttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGa  >  1:939089/1‑69 (MQ=255)
                                                    cttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGa  >  1:1001225/1‑69 (MQ=255)
                                                    cttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGa  >  1:3385935/1‑69 (MQ=255)
                                                    cttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGa  >  1:3334082/1‑69 (MQ=255)
                                                    cttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGa  >  1:2902967/1‑69 (MQ=255)
                                                    cttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGa  >  1:2836744/1‑69 (MQ=255)
                                                    cttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGa  >  1:271869/1‑69 (MQ=255)
                                                    cttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGa  >  1:2676027/1‑69 (MQ=255)
                                                    cttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGa  >  1:2668170/1‑69 (MQ=255)
                                                    cttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGa  >  1:2593428/1‑69 (MQ=255)
                                                    cttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGa  >  1:2309321/1‑69 (MQ=255)
                                                    cttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGa  >  1:214790/1‑69 (MQ=255)
                                                    cttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGa  >  1:1363212/1‑69 (MQ=255)
                                                    cttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGa  >  1:1300586/1‑69 (MQ=255)
                                                    cttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGa  >  1:1132723/1‑69 (MQ=255)
                                                    cttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGa  >  1:110055/1‑69 (MQ=255)
                                                    cttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGa  >  1:1044177/1‑69 (MQ=255)
                                                    cttcttTACGGCATTTGGTGCCAGAACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGa  >  1:1042299/1‑69 (MQ=255)
                                                                |                                                        
TGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGA  >  minE/318072‑318192

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: