Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,848,531 C→T 45.3% V834V (GTG→GTA yjeP ← predicted mechanosensitive channel

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,848,5310CT45.3% 16.5 / 69.3 86V834V (GTG→GTAyjePpredicted mechanosensitive channel
Reads supporting (aligned to +/- strand):  ref base C (11/36);  new base T (27/12);  total (38/48)
Fisher's exact test for biased strand distribution p-value = 2.77e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACA  >  minE/2848466‑2848595
                                                                 |                                                                
tCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAgcac                                                                >  1:3394587/1‑68 (MQ=255)
tCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAgcac                                                                >  1:3093044/1‑68 (MQ=255)
         gcgggcggGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAg                                                       >  1:2540/1‑68 (MQ=255)
         gcgGGCAGGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGa                                                        >  1:2521609/1‑67 (MQ=255)
         gcgGGCAGGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAg                                                       >  1:2605357/1‑68 (MQ=255)
         gcgGGCAGGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAg                                                       >  1:147849/1‑68 (MQ=255)
         gcgGGCAGGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAg                                                       >  1:1563571/1‑68 (MQ=255)
         gcgGGCAGGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAg                                                       >  1:1580682/1‑68 (MQ=255)
         gcgGGCAGGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAg                                                       >  1:160038/1‑68 (MQ=255)
         gcgGGCAGGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAg                                                       >  1:164836/1‑68 (MQ=255)
         gcgGGCAGGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAg                                                       >  1:2049023/1‑68 (MQ=255)
         gcgGGCAGGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAg                                                       >  1:2085673/1‑68 (MQ=255)
         gcgGGCAGGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAg                                                       >  1:2432070/1‑68 (MQ=255)
         gcgGGCAGGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAg                                                       >  1:1017290/1‑68 (MQ=255)
         gcgGGCAGGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAg                                                       >  1:2815643/1‑68 (MQ=255)
         gcgGGCAGGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAg                                                       >  1:2905661/1‑68 (MQ=255)
         gcgGGCAGGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAg                                                       >  1:2913996/1‑68 (MQ=255)
         gcgGGCAGGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAg                                                       >  1:3086624/1‑68 (MQ=255)
         gcgGGCAGGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAg                                                       >  1:3367930/1‑68 (MQ=255)
         gcgGGCAGGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAg                                                       >  1:899463/1‑68 (MQ=255)
         gcgGGCAGGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAg                                                       >  1:916327/1‑68 (MQ=255)
         gcgGGCAGGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAg                                                       >  1:133252/1‑68 (MQ=255)
         gcgGGCAGGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAg                                                       >  1:1242471/1‑68 (MQ=255)
         gcgGGCAGGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAg                                                       >  1:1043184/1‑68 (MQ=255)
         gcgGGCAGGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAg                                                       >  1:1461870/1‑68 (MQ=255)
            ggCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTa                                                   <  1:1288726/69‑1 (MQ=255)
            ggCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTa                                                   <  1:1049679/69‑1 (MQ=255)
                   tGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAGGGTAATCGGc                                             >  1:3409575/1‑68 (MQ=255)
                   tGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAGGGTAATCGGc                                             >  1:1002674/1‑68 (MQ=255)
                   tGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAGGGTAATCGGc                                             >  1:1186792/1‑68 (MQ=255)
                   tGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAGGGTAATCGGc                                             >  1:360464/1‑68 (MQ=255)
                        aCCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGa                                       <  1:2755/69‑1 (MQ=255)
                        aCCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGa                                       <  1:1798601/69‑1 (MQ=255)
                        aCCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGa                                       <  1:651918/69‑1 (MQ=255)
                        aCCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGa                                       <  1:807179/69‑1 (MQ=255)
                        aCCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGa                                       <  1:3209838/69‑1 (MQ=255)
                        aCCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGa                                       <  1:2509663/69‑1 (MQ=255)
                        aCCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGa                                       <  1:3196451/69‑1 (MQ=255)
                        aCCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGa                                       <  1:2492933/69‑1 (MQ=255)
                          cAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACt                                     <  1:105877/69‑1 (MQ=255)
                          cAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACt                                     <  1:2614914/69‑1 (MQ=255)
                          cAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACt                                     <  1:109800/69‑1 (MQ=255)
                          cAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACt                                     <  1:567886/69‑1 (MQ=255)
                          cAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACt                                     <  1:647648/69‑1 (MQ=255)
                           aGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACt                                     <  1:1411105/68‑1 (MQ=255)
                               gCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAgtgt                                                    <  1:1156966/49‑1 (MQ=255)
                               gCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGt                                                    <  1:3201034/49‑1 (MQ=255)
                               gCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGt                                                    <  1:952259/49‑1 (MQ=255)
                               gCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGt                                                    <  1:3380309/49‑1 (MQ=255)
                               gCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGt                                                    <  1:1602552/49‑1 (MQ=255)
                               gCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGt                                                    <  1:2573449/49‑1 (MQ=255)
                                      gatgatAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTg                         >  1:196470/1‑69 (MQ=255)
                                          atAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTAcc                     >  1:789629/1‑69 (MQ=255)
                                          atAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTAcc                     >  1:3224917/1‑69 (MQ=255)
                                            aaaCACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGt                   <  1:812738/69‑1 (MQ=255)
                                            aaaCACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGt                   <  1:448923/69‑1 (MQ=255)
                                            aaaCACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGt                   <  1:2145824/69‑1 (MQ=255)
                                            aaaCACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGt                   <  1:2250571/69‑1 (MQ=255)
                                            aaaCACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGt                   <  1:225535/69‑1 (MQ=255)
                                            aaaCACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGt                   <  1:2839088/69‑1 (MQ=255)
                                                    gAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGa           <  1:3175436/69‑1 (MQ=255)
                                                       tGGCTATCAGTACCGCACCGAGGGTAATCGGCTCCAGGCTTTc                                  <  1:483355/43‑1 (MQ=25)
                                                       tGGCAATCAGTACCGCACCGAGGGTAATCGGCTCCAGGCTTTc                                  <  1:1939040/43‑1 (MQ=37)
                                                       tGGCAATCAGTACCGCACCGAGGGTAATCGGCTCCAGGCTTTc                                  <  1:2584699/43‑1 (MQ=37)
                                                       tGGCAATCAGTACCGCACCGAGGGTAATCGGCTCCAGGCTTTc                                  <  1:1609684/43‑1 (MQ=37)
                                                       tGGCAATCAGTACCGCACCGAGGGTAATCGGCTCCAGGCTTTc                                  <  1:1347142/43‑1 (MQ=37)
                                                       tGGCAATCAGTACCGCACCGAGGGTAATCGGCTCCAGGCTTTc                                  <  1:62344/43‑1 (MQ=37)
                                                       tGGCAATCAGTACCGCACCGAGGGTAATCGGCTCCAGGCTTTc                                  <  1:50455/43‑1 (MQ=37)
                                                       tGGCAATCAGTACCGCACCGAGGGTAATCGGCTCCAGGCTTTc                                  <  1:391969/43‑1 (MQ=37)
                                                       tGGCAATCAGTACCGCACCGAGGGTAATCGGCTCCAGGCTTTc                                  <  1:3355190/43‑1 (MQ=37)
                                                       tGGCAATCAGTACCGCACCGAGGGTAATCGGCTCCAGGCTTTc                                  <  1:2287889/43‑1 (MQ=37)
                                                       tGGCAATCAGTACCGCACCGAGGGTAATCGGCTCCAGGCTTTc                                  <  1:3086932/43‑1 (MQ=37)
                                                          cAATCAGTACCGCACCGAGGGTAATCGGCTCCAGGCTTTc                                  <  1:1271763/40‑1 (MQ=37)
                                                              cAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCcaca  >  1:1335482/1‑68 (MQ=255)
                                                              cAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCcaca  >  1:1878512/1‑68 (MQ=255)
                                                              cAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCcaca  >  1:1404143/1‑68 (MQ=255)
                                                              cAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCcaca  >  1:1400146/1‑68 (MQ=255)
                                                              cAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCcaca  >  1:3222526/1‑68 (MQ=255)
                                                              cAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCcaca  >  1:1160019/1‑68 (MQ=255)
                                                              cAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCcaca  >  1:1290977/1‑68 (MQ=255)
                                                              cAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCcaca  >  1:2641474/1‑68 (MQ=255)
                                                               aGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTTACa         <  1:1840017/60‑1 (MQ=255)
                                                               aGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACa         <  1:1527349/60‑1 (MQ=255)
                                                               aGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACa         <  1:897231/60‑1 (MQ=255)
                                                               aGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACa         <  1:2338328/60‑1 (MQ=255)
                                                               aGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACa         <  1:117581/60‑1 (MQ=255)
                                                               aGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACa         <  1:2640119/60‑1 (MQ=255)
                                                               aGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACACCCTGTACCGTGGAGGTGACa         <  1:2784814/60‑1 (MQ=255)
                                                                 |                                                                
TCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACA  >  minE/2848466‑2848595

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