Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,849,476 C→T 71.5% L519L (TTG→TTA yjeP ← predicted mechanosensitive channel

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,849,4760CT71.5% 60.8 / 25.8 63L519L (TTG→TTAyjePpredicted mechanosensitive channel
Reads supporting (aligned to +/- strand):  ref base C (10/8);  new base T (17/28);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTC  >  minE/2849408‑2849532
                                                                    |                                                        
ccGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTgcgc                                                          <  1:2065956/69‑1 (MQ=255)
ccGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTgcgc                                                          <  1:1950881/69‑1 (MQ=255)
ccGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTgcgc                                                          <  1:1372000/69‑1 (MQ=255)
ccGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTgcgc                                                          <  1:261578/69‑1 (MQ=255)
     aCCCAGCGTCATCCACAGCACAGGCCGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:717044/69‑1 (MQ=255)
     aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGAGGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:857916/69‑1 (MQ=255)
     aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:310667/69‑1 (MQ=255)
     aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:2173155/69‑1 (MQ=255)
     aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:2332320/69‑1 (MQ=255)
     aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:2422952/69‑1 (MQ=255)
     aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:2517510/69‑1 (MQ=255)
     aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:2580048/69‑1 (MQ=255)
     aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:2783877/69‑1 (MQ=255)
     aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:2941483/69‑1 (MQ=255)
     aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:2545329/69‑1 (MQ=255)
     aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:315634/69‑1 (MQ=255)
     aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:3223272/69‑1 (MQ=255)
     aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:568812/69‑1 (MQ=255)
     aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:8757/69‑1 (MQ=255)
     aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:893727/69‑1 (MQ=255)
     aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:948170/69‑1 (MQ=255)
     aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:963446/69‑1 (MQ=255)
     aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:1085541/69‑1 (MQ=255)
     aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:153473/69‑1 (MQ=255)
     aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:1332738/69‑1 (MQ=255)
     aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:1885324/69‑1 (MQ=255)
     aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:1945041/69‑1 (MQ=255)
     aCCCAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:1280785/69‑1 (MQ=255)
      cccAGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:2534561/68‑1 (MQ=255)
         aGCGTCATCCACAGCACAGGCAGCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:283303/65‑1 (MQ=255)
                 ccACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTg                                         <  1:1105471/69‑1 (MQ=255)
                 ccACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTg                                         <  1:3018518/69‑1 (MQ=255)
                 ccACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTg                                         <  1:2589377/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:2253769/1‑69 (MQ=255)
                  cacaGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGa                                        >  1:318481/1‑69 (MQ=255)
                  cacaGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGa                                        >  1:1644688/1‑69 (MQ=255)
                  cacaGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGa                                        >  1:3392612/1‑69 (MQ=255)
                  cacaGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGa                                        >  1:1855052/1‑69 (MQ=255)
                  cacaGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGa                                        >  1:642764/1‑69 (MQ=255)
                  cacaGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGa                                        <  1:1011114/69‑1 (MQ=255)
                  cacaGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGa                                        >  1:1911952/1‑69 (MQ=255)
                  cacaGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGa                                        >  1:796199/1‑69 (MQ=255)
                  cacaGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGa                                        >  1:2362629/1‑69 (MQ=255)
                              agcGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGTAACGt                                                     <  1:2184816/42‑1 (MQ=39)
                                                        gAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTc  <  1:2297923/69‑1 (MQ=255)
                                                        gAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTc  <  1:745557/69‑1 (MQ=255)
                                                         aaaaGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTc  <  1:1460667/68‑1 (MQ=255)
                                                          aaaGAGTGCGTAACGTCAGCCAGAAGTGATCCTGAGTCa                              >  1:555701/1‑39 (MQ=255)
                                                          aaaGAGTGCGTAACGTCAGCCAGAAGTGATCCTGAGTCAc                             >  1:2620749/1‑40 (MQ=255)
                                                          aaaGAGTGCGTAACGTCAGCCAGAAGTGATCCTGAGTCAc                             >  1:106271/1‑40 (MQ=255)
                                                          aaaGAGTGCGTAACGTCAGCCAGAAGTGATCCTGAGTCAc                             >  1:1077950/1‑40 (MQ=255)
                                                          aaaGAGTGCGTAACGTCAGCCAGAAGTGATCCTGAGTCAc                             >  1:2058762/1‑40 (MQ=255)
                                                          aaaGAGTGCGTAACGTCAGCCAGAAGTGATCCTGAGTCAc                             >  1:883612/1‑40 (MQ=255)
                                                          aaaGAGTGCGTAACGTCAGCCAGAAGTGATCCTGAGTCAc                             >  1:2428714/1‑40 (MQ=255)
                                                          aaaGAGTGCGTAACGTCAGCCAGAAGTGATCCTGAGTCAc                             >  1:745908/1‑40 (MQ=255)
                                                          aaaGAGTGCGTAACGTCAGCCAGAAGTGATCCTGAGTCAc                             >  1:1554256/1‑40 (MQ=255)
                                                          aaaGAGTGCGTAACGTCAGCCAGAAGTGATCCTGAGTCAc                             >  1:2671802/1‑40 (MQ=255)
                                                          aaaGAGTGCGTAACGTCAGCCAGAAGTGATCCTGAGTCAc                             >  1:613322/1‑40 (MQ=255)
                                                          aaaGAGTGCGTAACGTCAGCCAGAAGTGATCCTGAGTCAc                             >  1:518818/1‑40 (MQ=255)
                                                          aaaGAGTGCGTAACGTCAGCCAGAAGTGATCCTGAGTCAc                             >  1:2852725/1‑40 (MQ=255)
                                                          aaaGAGTGCGTAACGTCAGCCAGAAGTGATCCTGAGTCAc                             >  1:3337707/1‑40 (MQ=255)
                                                          aaaGAGTGCGTAACGTCAGCCAGAAGTGATCCTGAGTCAc                             >  1:3242770/1‑40 (MQ=255)
                                                          aaaGAGTGCGTAACGTCAACCAGAAGTGATCCTGAGTCAc                             >  1:1999490/1‑40 (MQ=255)
                                                          aaaGAGAGCGTAACGTCAGCCAGAAGTGATCCTGAGt                                >  1:1837411/1‑37 (MQ=38)
                                                                    |                                                        
CCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTC  >  minE/2849408‑2849532

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: