Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,850,177 G→A 36.8% 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,1770GA36.8% 51.3 / 58.7 57Q286* (CAG→TAG) yjePpredicted mechanosensitive channel
Reads supporting (aligned to +/- strand):  ref base G (21/15);  new base A (11/10);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 7.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

TGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATT  >  minE/2850110‑2850244
                                                                   |                                                                   
tGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTg                                                                  >  1:3621087/1‑71 (MQ=255)
 gCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGc                                                                 >  1:1075304/1‑71 (MQ=255)
 gCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGc                                                                 >  1:490937/1‑71 (MQ=255)
 gCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGc                                                                 >  1:1835419/1‑71 (MQ=255)
 gCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGc                                                                 >  1:259373/1‑71 (MQ=255)
 gCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGc                                                                 >  1:2488568/1‑71 (MQ=255)
       gTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGc                                                            >  1:560946/1‑70 (MQ=255)
       gTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGc                                                            >  1:190428/1‑70 (MQ=255)
       gTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGc                                                            >  1:3868450/1‑70 (MQ=255)
        ttCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGc                                                            <  1:1940961/69‑1 (MQ=255)
        ttCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCaa                                                          >  1:564117/1‑71 (MQ=255)
        ttCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGc                                                            <  1:1335028/69‑1 (MQ=255)
        ttCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGc                                                            <  1:549357/69‑1 (MQ=255)
        ttCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGc                                                            <  1:3584436/69‑1 (MQ=255)
          cACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAAc                                                         >  1:1581767/1‑70 (MQ=255)
           aCGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACg                                                        <  1:993786/70‑1 (MQ=255)
           aCGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACga                                                       <  1:788153/71‑1 (MQ=255)
           aCGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACga                                                       <  1:2577874/71‑1 (MQ=255)
           aCGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACga                                                       <  1:3815338/71‑1 (MQ=255)
           aCGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACga                                                       <  1:1691244/71‑1 (MQ=255)
              cAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGat                                                    <  1:436966/71‑1 (MQ=255)
               aGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGat                                                    <  1:150940/70‑1 (MQ=255)
                gCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGatcc                                                  <  1:744112/71‑1 (MQ=255)
                   tATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGatc                                                   <  1:3831861/67‑1 (MQ=255)
                     ttCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCg                                             >  1:525007/1‑71 (MQ=255)
                       cAACGCGTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCt                                           >  1:530837/1‑71 (MQ=255)
                       cAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGc                                            >  1:2775253/1‑70 (MQ=255)
                       cAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCt                                           >  1:1812018/1‑71 (MQ=255)
                       cAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCt                                           >  1:3595580/1‑71 (MQ=255)
                       cAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCt                                           >  1:2981497/1‑71 (MQ=255)
                       cAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCt                                           >  1:2725512/1‑71 (MQ=255)
                       cAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGATCCATCCGCt                                           >  1:3671733/1‑71 (MQ=255)
                         aCGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTgc                                         <  1:625222/71‑1 (MQ=255)
                         aCGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTgc                                         <  1:2611159/71‑1 (MQ=255)
                         aCGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTgc                                         <  1:2165789/71‑1 (MQ=255)
                                  ggACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGa                                <  1:947085/71‑1 (MQ=255)
                                    aCCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGatc                                                   <  1:3205760/50‑1 (MQ=255)
                                        gTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCa                            >  1:1715366/1‑69 (MQ=255)
                                        gTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCa                            >  1:1983068/1‑69 (MQ=255)
                                           aCGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgcc                       >  1:1001070/1‑71 (MQ=255)
                                           aCGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgcc                       >  1:3857520/1‑71 (MQ=255)
                                           aCGTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgcc                       >  1:2411858/1‑71 (MQ=255)
                                             gTCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGAt                               >  1:2311450/1‑61 (MQ=255)
                                              tCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgccgcc                    <  1:394436/71‑1 (MQ=255)
                                              tCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgccgcc                    <  1:2399883/71‑1 (MQ=255)
                                              tCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgccgcc                    <  1:8452/71‑1 (MQ=255)
                                              tCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgccgcc                    <  1:814750/71‑1 (MQ=255)
                                              tCTGGCTGGCAGCCTGACGCTATTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgccgcc                    <  1:3750766/71‑1 (MQ=255)
                                                tggctggcAGCCTGACGCTATTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgccgc                     >  1:201801/1‑68 (MQ=255)
                                                  gctggcAGCCTGACGCTGTTGCGAGGCAAGGAGATCCATCCGCTGCGCCTGTTGATTCAAAGccgccg                   <  1:417307/68‑1 (MQ=255)
                                                  gctggcAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGccgccg                   <  1:3913805/68‑1 (MQ=255)
                                                   ctggcAGCCTGACGCTGTTGCGAGGCAACGAGatccat                                                >  1:2153457/1‑38 (MQ=255)
                                                      gcAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGa                  >  1:2019684/1‑65 (MQ=255)
                                                      gcAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGa                  >  1:3386023/1‑65 (MQ=255)
                                                      gcAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGa                  >  1:2638596/1‑65 (MQ=255)
                                                                gCTATTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAAtt  >  1:3648183/1‑71 (MQ=255)
                                                                gCTATTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAAtt  >  1:3799361/1‑71 (MQ=255)
                                                                   |                                                                   
TGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATT  >  minE/2850110‑2850244

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