Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,392,090 G→T 28.8% L655M (CTG→ATG)  yjeP ← predicted mechanosensitive channel

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,392,0900GT28.8% 134.6 / 33.2 59L655M (CTG→ATG) yjePpredicted mechanosensitive channel
Reads supporting (aligned to +/- strand):  ref base G (23/19);  new base T (8/9);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 7.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

ATCATCATGTTCCACAGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCA  >  W3110S.gb/4392026‑4392156
                                                                |                                                                  
atcatcATGTTCCACAGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAAt                                                              >  1:1554045/1‑71 (MQ=255)
atcatcATGTTCCACAGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAAt                                                              >  1:3215082/1‑71 (MQ=255)
 tcatcaTGTTCCACAGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATc                                                             <  1:659270/71‑1 (MQ=255)
   atcatGTTCCACAGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATccc                                                           <  1:4561216/71‑1 (MQ=255)
      atGTTCCACAGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGc                                                        >  1:4819017/1‑71 (MQ=255)
          tCCACAGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCttt                                                    >  1:581360/1‑71 (MQ=255)
           ccACAGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACATCGGAATCCCGGCCttt                                                    <  1:1171679/70‑1 (MQ=255)
           ccACAGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACATCGGAATCCCGGCCttt                                                    <  1:2328590/70‑1 (MQ=255)
            cacaGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACATCGGAATCCCGGCCttttt                                                  <  1:1349482/71‑1 (MQ=255)
            cacaGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACATCGGAATCCCGGCCttttt                                                  <  1:1821421/71‑1 (MQ=255)
                  atatGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCt                                            <  1:2505386/71‑1 (MQ=255)
                   tatGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACATCGGAATCCCGGCCttttt                                                  >  1:730067/1‑64 (MQ=255)
                   tatGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCgg                                                                 <  1:5143324/49‑1 (MQ=255)
                   tatGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCgg                                                                 <  1:2308284/49‑1 (MQ=255)
                     tggttggtAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAg                                         >  1:2111545/1‑71 (MQ=255)
                       gttggtAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAgg                                        >  1:4950680/1‑70 (MQ=255)
                       gttggtAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAgg                                        <  1:2341953/70‑1 (MQ=255)
                        ttggtAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACATCGGAATCCCGGCCTTTTTCAGGCTGAGGGt                                      <  1:1717963/71‑1 (MQ=255)
                        ttggtAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCGTTTTCAGGCTGAGGGt                                      <  1:1392313/71‑1 (MQ=255)
                                gttgtCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGaccac                                >  1:2577017/1‑69 (MQ=255)
                                gttgtCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGaccac                                >  1:2583288/1‑69 (MQ=255)
                                gttgtCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGaccac                                >  1:4810437/1‑69 (MQ=255)
                                gttgtCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGac                                   >  1:1871693/1‑66 (MQ=255)
                                 ttgtCGCCGCTGCCCTCTTTGTTGAGATACATCGGAATCCCGGCCttttt                                                  >  1:4828265/1‑50 (MQ=255)
                                 ttgtCGCCGCTGCCCTCTTTGTTGAGATACATCGGAATCCCGGCCttttt                                                  >  1:3017079/1‑50 (MQ=255)
                                 ttgtCGCCGCTGCCCTCTTTGTTGAGATACATCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGac                                   >  1:260629/1‑65 (MQ=255)
                                 ttgtCGCCGCTGCCCTCTTTGTTGAGATACATCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGc                             >  1:1062857/1‑71 (MQ=255)
                                 ttgtCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGc                             >  1:2877447/1‑71 (MQ=255)
                                 ttgtCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGc                             >  1:425298/1‑71 (MQ=255)
                                   gtCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGc                             >  1:659188/1‑69 (MQ=255)
                                           tGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCttt                                                    <  1:3040668/38‑1 (MQ=255)
                                             cccTCTTTGTTGAGATACATCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCa                 <  1:530133/71‑1 (MQ=255)
                                             cccTCTTTGTTGAGATACATCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCa                 <  1:4258947/71‑1 (MQ=255)
                                             cccTCTTTGTTGAGATACATCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCa                 <  1:4587056/71‑1 (MQ=255)
                                             cccTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCt                                            <  1:1255481/44‑1 (MQ=255)
                                             cccTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCt                                            <  1:1931378/44‑1 (MQ=255)
                                             cccTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGc                  <  1:4951309/70‑1 (MQ=255)
                                               ctctTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGacc                                  <  1:1854818/52‑1 (MQ=255)
                                                tctTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGacc                                  <  1:3778721/51‑1 (MQ=255)
                                                 ctTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAAt              >  1:4625491/1‑70 (MQ=255)
                                                  tttGTTGAGATACATCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAAt              >  1:1486051/1‑69 (MQ=255)
                                                  tttGTTGAGATACATCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAAt              <  1:4541205/69‑1 (MQ=255)
                                                  tttGTTGAGATACAGCGGAATCCCGGCCTTTTTCa                                                >  1:752775/1‑35 (MQ=255)
                                                  tttGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAAt              >  1:3815693/1‑69 (MQ=255)
                                                  tttGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAAt              >  1:2860377/1‑69 (MQ=255)
                                                    tgttgAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGa            <  1:4327752/69‑1 (MQ=255)
                                                           aTACAGGGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATaaag    >  1:658736/1‑70 (MQ=255)
                                                           aTACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCa                      >  1:3225959/1‑52 (MQ=255)
                                                           aTACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATaaag    >  1:4580410/1‑70 (MQ=255)
                                                           aTACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATaaag    >  1:833171/1‑70 (MQ=255)
                                                           aTACAGCGGAATCCCGGCCTATTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATaaag    >  1:2855682/1‑70 (MQ=255)
                                                            tACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATa       <  1:3428480/66‑1 (MQ=255)
                                                            tACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATa       <  1:2399798/66‑1 (MQ=255)
                                                            tACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATa       <  1:2207924/66‑1 (MQ=255)
                                                             aCATCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATaaagca  >  1:4148696/1‑70 (MQ=255)
                                                             aCATCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATaaagca  >  1:3157848/1‑70 (MQ=255)
                                                             aCAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATaaagca  >  1:4006039/1‑70 (MQ=255)
                                                               aGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGAt        <  1:527245/62‑1 (MQ=255)
                                                               aGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGAt        <  1:3444084/62‑1 (MQ=255)
                                                                |                                                                  
ATCATCATGTTCCACAGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCA  >  W3110S.gb/4392026‑4392156

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: