Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,850,177 G→A 54.2% Q286* (CAG→TAG)  yjeP ← predicted mechanosensitive channel

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,850,1770GA54.2% 8.2 / 83.3 59Q286* (CAG→TAG) yjePpredicted mechanosensitive channel
Reads supporting (aligned to +/- strand):  ref base G (12/15);  new base A (10/22);  total (22/37)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

GCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTT  >  minE/2850111‑2850233
                                                                  |                                                        
gCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTg                                                       >  1:3618345/1‑70 (MQ=255)
gCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGc                                                      >  1:1730780/1‑71 (MQ=255)
gCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGc                                                      >  1:478974/1‑71 (MQ=255)
gCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGc                                                      >  1:3338752/1‑71 (MQ=255)
 cGACTGTTCACGCCGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGc                                                      <  1:4180325/70‑1 (MQ=255)
 cGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGc                                                      <  1:2646403/70‑1 (MQ=255)
 cGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGc                                                      <  1:3365216/70‑1 (MQ=255)
 cGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGc                                                      <  1:2030038/70‑1 (MQ=255)
 cGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGc                                                      <  1:2093533/70‑1 (MQ=255)
 cGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGc                                                      <  1:1003372/70‑1 (MQ=255)
 cGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGc                                                      <  1:3865448/70‑1 (MQ=255)
   aCTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAg                                                   >  1:3357277/1‑71 (MQ=255)
   aCTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAg                                                   >  1:2604285/1‑71 (MQ=255)
   aCTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAg                                                   >  1:807660/1‑71 (MQ=255)
         cACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGc                                                 >  1:2193155/1‑67 (MQ=255)
         cACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAAc                                              >  1:642924/1‑70 (MQ=255)
         cACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAAc                                              >  1:2340251/1‑70 (MQ=255)
          aCGCAGCGTATTCAACGCCTGGCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACga                                            <  1:809738/71‑1 (MQ=255)
          aCGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACg                                             <  1:2419821/70‑1 (MQ=255)
              aGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTTACGCTGTTTCGAGGCAACGAGatc                                        <  1:2069313/71‑1 (MQ=255)
              aGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGatc                                        <  1:2041719/71‑1 (MQ=255)
              aGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGatc                                        <  1:1997396/71‑1 (MQ=255)
                cGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGatcca                                      >  1:2403419/1‑71 (MQ=255)
                cGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGatcca                                      <  1:3899079/71‑1 (MQ=255)
                cGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGatcca                                      <  1:1025603/71‑1 (MQ=255)
                    ttCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGatccatc                                    >  1:1885775/1‑69 (MQ=255)
                    ttCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGatccatc                                    >  1:3791983/1‑69 (MQ=255)
                    ttCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGatcca                                      <  1:982659/67‑1 (MQ=255)
                    ttCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGatcca                                      <  1:309432/67‑1 (MQ=255)
                    ttCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGatcca                                      <  1:4107742/67‑1 (MQ=255)
                    ttCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGatcca                                      <  1:4140600/67‑1 (MQ=255)
                    ttCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGatcca                                      <  1:255105/67‑1 (MQ=255)
                    ttCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGATCCATCCg                                  >  1:2166332/1‑71 (MQ=255)
                            cTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGATCCATCCGCTGCGCCt                          <  1:4020319/71‑1 (MQ=255)
                            cTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGATCCATCCGCTGCGCCt                          <  1:1640104/71‑1 (MQ=255)
                            cTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGATCCATCCGCTGCGCCt                          <  1:2309024/71‑1 (MQ=255)
                              gCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGc                                                      <  1:289316/41‑1 (MQ=255)
                                     cTGTAACGTCTGGCTTGCAGCCTGACGCTATTGCGAGGCAACGAGatcc                                       <  1:3923796/49‑1 (MQ=255)
                                     cTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGatcc                                       <  1:2024357/49‑1 (MQ=255)
                                      tGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGatcc                                       <  1:568725/48‑1 (MQ=255)
                                        tAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTgcgc                            <  1:3098108/57‑1 (MQ=255)
                                        tAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCaaa               <  1:381688/70‑1 (MQ=255)
                                        tAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCaaa               <  1:1836535/70‑1 (MQ=255)
                                        tAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAg              >  1:930882/1‑71 (MQ=255)
                                        tAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGat                                         >  1:217206/1‑44 (MQ=255)
                                        tAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTc                  <  1:62665/67‑1 (MQ=255)
                                        tAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTc                  <  1:3557462/67‑1 (MQ=255)
                                        tAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTc                  <  1:282300/67‑1 (MQ=255)
                                        tAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAg              >  1:406261/1‑71 (MQ=255)
                                        tAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGATCCATCAGCTGCGCCTGTTGATTCAAAg              >  1:1943040/1‑71 (MQ=255)
                                         aaCGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgc             >  1:2305558/1‑71 (MQ=255)
                                          aCGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgcc            >  1:1302830/1‑71 (MQ=255)
                                             tCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgccgcc         <  1:551267/71‑1 (MQ=255)
                                             tCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgccgcc         <  1:3047657/71‑1 (MQ=255)
                                             tCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgccgcc         <  1:3731123/71‑1 (MQ=255)
                                              ctggctggCAGCCTGACGCTATTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgccgcc         <  1:2540810/70‑1 (MQ=255)
                                               tggctggcAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTc                  >  1:2327813/1‑60 (MQ=255)
                                                 gctggcAGCCTGACGCTATTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGa       >  1:1572663/1‑69 (MQ=255)
                                                     gcAGCCTGACGCTATTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGtt  <  1:677940/70‑1 (MQ=255)
                                                                  |                                                        
GCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTT  >  minE/2850111‑2850233

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: