Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,849,434 C→A 61.1% P533P (CCG→CCT yjeP ← predicted mechanosensitive channel

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,849,4340CA61.1% 41.7 / 81.6 90P533P (CCG→CCTyjePpredicted mechanosensitive channel
Reads supporting (aligned to +/- strand):  ref base C (12/23);  new base A (28/27);  total (40/50)
Fisher's exact test for biased strand distribution p-value = 1.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGA  >  minE/2849380‑2849494
                                                      |                                                            
cAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGAcgc                                                <  1:3294596/69‑1 (MQ=255)
           gCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCg                                     >  1:3010554/1‑69 (MQ=255)
           gCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCg                                     >  1:394533/1‑69 (MQ=255)
                  cTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGaa                              <  1:2941652/69‑1 (MQ=255)
                  cTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGaa                              <  1:2878756/69‑1 (MQ=255)
                  cTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGaa                              <  1:2789904/69‑1 (MQ=255)
                  cTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGaa                              <  1:1209160/69‑1 (MQ=255)
                  cTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGaa                              <  1:526849/69‑1 (MQ=255)
                  cTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGaa                              <  1:2005999/69‑1 (MQ=255)
                  cTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGaa                              <  1:758859/69‑1 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGaaaa                            <  1:1010508/69‑1 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGaaaa                            <  1:2830446/69‑1 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGaaaa                            <  1:2155012/69‑1 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGaaaa                            <  1:835744/69‑1 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGAcg                                                 >  1:2375218/1‑48 (MQ=37)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGcgac                                             >  1:3251745/1‑52 (MQ=37)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGcgac                                             >  1:2247250/1‑52 (MQ=37)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGaa                                        >  1:191439/1‑57 (MQ=38)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGaa                              >  1:3206656/1‑67 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGaa                              >  1:2623081/1‑67 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGaa                              >  1:102835/1‑67 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGaa                              >  1:3228621/1‑67 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGaa                              >  1:675827/1‑67 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGaa                              >  1:935083/1‑67 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGaa                              >  1:3304974/1‑67 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGaa                              >  1:3099172/1‑67 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGaa                              >  1:900757/1‑67 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGaa                              >  1:2927585/1‑67 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGaa                              >  1:1021139/1‑67 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGaa                              >  1:577664/1‑67 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGaa                              >  1:1047415/1‑67 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGaa                              >  1:696125/1‑67 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGaa                              >  1:1815414/1‑67 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGaa                              >  1:1631528/1‑67 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGaa                              >  1:1607266/1‑67 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGaa                              >  1:726438/1‑67 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGaa                              >  1:193267/1‑67 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGaa                              >  1:1956539/1‑67 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGaa                              >  1:1974755/1‑67 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGaa                              >  1:2069313/1‑67 (MQ=255)
                    cgcgCAAGCCGTAACCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGa                               >  1:3307516/1‑66 (MQ=255)
                    cgcgCAAGCCGTAAACCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGaa                              >  1:1924270/1‑67 (MQ=255)
                            ccGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTgcgc                    <  1:1950881/69‑1 (MQ=255)
                            ccGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTgcgc                    <  1:261578/69‑1 (MQ=255)
                            ccGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTgcgc                    <  1:2065956/69‑1 (MQ=255)
                            ccGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTgcgc                    <  1:1372000/69‑1 (MQ=255)
                                 aCCCAGCGTCATCCACAGCACAGGCCGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:717044/69‑1 (MQ=255)
                                 aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGAGGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:857916/69‑1 (MQ=255)
                                 aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:963446/69‑1 (MQ=255)
                                 aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:948170/69‑1 (MQ=255)
                                 aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:310667/69‑1 (MQ=255)
                                 aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:315634/69‑1 (MQ=255)
                                 aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:3223272/69‑1 (MQ=255)
                                 aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:893727/69‑1 (MQ=255)
                                 aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:568812/69‑1 (MQ=255)
                                 aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:8757/69‑1 (MQ=255)
                                 aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:2941483/69‑1 (MQ=255)
                                 aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:1885324/69‑1 (MQ=255)
                                 aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:1945041/69‑1 (MQ=255)
                                 aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:2173155/69‑1 (MQ=255)
                                 aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:153473/69‑1 (MQ=255)
                                 aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:2332320/69‑1 (MQ=255)
                                 aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:1332738/69‑1 (MQ=255)
                                 aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:1280785/69‑1 (MQ=255)
                                 aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:2422952/69‑1 (MQ=255)
                                 aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:2517510/69‑1 (MQ=255)
                                 aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:2545329/69‑1 (MQ=255)
                                 aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:2580048/69‑1 (MQ=255)
                                 aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:1085541/69‑1 (MQ=255)
                                 aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:2783877/69‑1 (MQ=255)
                                  cccAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:2534561/68‑1 (MQ=255)
                                     aGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:283303/65‑1 (MQ=255)
                                             ccACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTg   <  1:2589377/69‑1 (MQ=255)
                                             ccACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTg   <  1:1105471/69‑1 (MQ=255)
                                             ccACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTg   <  1:3018518/69‑1 (MQ=255)
                                             ccACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt               <  1:1879777/57‑1 (MQ=255)
                                              cacaGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGGCAGAAGTGa  <  1:2263383/69‑1 (MQ=255)
                                              cacaGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTg   >  1:2635357/1‑68 (MQ=255)
                                              cacaGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGa  >  1:1644688/1‑69 (MQ=255)
                                              cacaGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGa  >  1:1855052/1‑69 (MQ=255)
                                              cacaGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGa  <  1:1011114/69‑1 (MQ=255)
                                              cacaGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGa  >  1:796199/1‑69 (MQ=255)
                                              cacaGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGa  >  1:1911952/1‑69 (MQ=255)
                                              cacaGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGa  >  1:642764/1‑69 (MQ=255)
                                              cacaGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGa  >  1:2253769/1‑69 (MQ=255)
                                              cacaGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGa  >  1:2362629/1‑69 (MQ=255)
                                              cacaGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGa  >  1:3392612/1‑69 (MQ=255)
                                              cacaGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGa  >  1:318481/1‑69 (MQ=255)
                                                caGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAag                           <  1:2930818/42‑1 (MQ=255)
                                                    aCCGGTAACGGTGACGCGACGAGAATCGACCAGaaaa                            <  1:2766777/37‑1 (MQ=255)
                                                      |                                                            
CAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGA  >  minE/2849380‑2849494

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: