Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,131,628 T→A 37.8% V440V (GTA→GTT deaD ← ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,131,6280TA37.8% 60.0 / 77.9 91V440V (GTA→GTTdeaDATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base T (7/49);  new base A (16/18);  total (23/68)
Fisher's exact test for biased strand distribution p-value = 4.26e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

CGCGCGGACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTC  >  minE/2131564‑2131696
                                                                |                                                                    
cgcgcgGACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGTACg                                                                   >  1:233083/1‑68 (MQ=255)
     ggACCACGGTCATCACGGTCACGGAATTCACTTTTCGGACGCATCGGCGCATCTGGCGGAACAATCaga                                                             <  1:681444/69‑1 (MQ=255)
     ggACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGAACAATCaga                                                             <  1:985030/69‑1 (MQ=255)
     ggACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGAACAATCaga                                                             <  1:1003696/69‑1 (MQ=255)
     ggACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGAACAATCaga                                                             <  1:945528/69‑1 (MQ=255)
     ggACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGAACAATCaga                                                             <  1:936805/69‑1 (MQ=255)
     ggACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGAACAATCaga                                                             <  1:810480/69‑1 (MQ=255)
     ggACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGAACAATCaga                                                             <  1:3370109/69‑1 (MQ=255)
     ggACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGAACAATCaga                                                             <  1:308984/69‑1 (MQ=255)
     ggACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGAACAATCaga                                                             <  1:2698451/69‑1 (MQ=255)
     ggACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGAACAATCaga                                                             <  1:2577112/69‑1 (MQ=255)
     ggACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGAACAATCaga                                                             <  1:2562535/69‑1 (MQ=255)
     ggACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGAACAATCaga                                                             <  1:2282725/69‑1 (MQ=255)
     ggACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGAACAATCaga                                                             <  1:2278245/69‑1 (MQ=255)
     ggACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGAACAATCaga                                                             <  1:1949270/69‑1 (MQ=255)
     ggACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGAACAATCaga                                                             <  1:1650883/69‑1 (MQ=255)
     ggACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGAACAATCaga                                                             <  1:135234/69‑1 (MQ=255)
     ggACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGAACAATCaga                                                             <  1:1008025/69‑1 (MQ=255)
      gACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCagag                                                            <  1:2266189/69‑1 (MQ=255)
      gACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCagag                                                            <  1:1262940/69‑1 (MQ=255)
      gACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCagag                                                            <  1:1202477/69‑1 (MQ=255)
      gACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCagag                                                            <  1:2673086/69‑1 (MQ=255)
      gACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCagag                                                            <  1:1964443/69‑1 (MQ=255)
      gACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCagag                                                            <  1:1410176/69‑1 (MQ=255)
      gACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCagag                                                            <  1:2885480/69‑1 (MQ=255)
      gACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCagag                                                            <  1:49667/69‑1 (MQ=255)
      gACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCagag                                                            <  1:2917837/69‑1 (MQ=255)
      gACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCAGCGGCGCATCTGGCGGTACGATCagag                                                            <  1:2382241/69‑1 (MQ=255)
       aCCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCagag                                                            <  1:1560800/68‑1 (MQ=255)
        ccACGGTCCTCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTa                                                          <  1:1590888/69‑1 (MQ=255)
        ccACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTa                                                          <  1:372090/69‑1 (MQ=255)
        ccACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTa                                                          <  1:674513/69‑1 (MQ=255)
        ccACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTa                                                          <  1:1699899/69‑1 (MQ=255)
        ccACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTa                                                          <  1:762488/69‑1 (MQ=255)
        ccACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTa                                                          <  1:785268/69‑1 (MQ=255)
        ccACGGTCATCACGGGCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTa                                                          <  1:1069250/69‑1 (MQ=255)
                           ggAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGAACAATCaga                                                             <  1:333731/47‑1 (MQ=38)
                            gAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCAtt                                        >  1:55222/1‑67 (MQ=255)
                            gAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCAtt                                        >  1:743631/1‑67 (MQ=255)
                            gAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCAtt                                        >  1:1007971/1‑67 (MQ=255)
                            gAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCAtt                                        >  1:1914788/1‑67 (MQ=255)
                                tCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTcagc                                  <  1:2495039/69‑1 (MQ=255)
                                tCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTcagc                                  <  1:1780333/69‑1 (MQ=255)
                                tCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTcagc                                  <  1:1143271/69‑1 (MQ=255)
                                tCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTcagc                                  <  1:1777987/69‑1 (MQ=255)
                                tCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTcagc                                  <  1:1787026/69‑1 (MQ=255)
                                tCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTcagc                                  <  1:872856/69‑1 (MQ=255)
                                tCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTcagc                                  <  1:549695/69‑1 (MQ=255)
                                tCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTcagc                                  <  1:2178042/69‑1 (MQ=255)
                                tCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTcagc                                  <  1:938233/69‑1 (MQ=255)
                                tCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTcagc                                  <  1:1360185/69‑1 (MQ=255)
                                tCACGTTTCGGACGCATCGGCGCATCTGGCGGCACGATCAGAGTACGTTCACCCTGTGCCATTTTcagc                                  <  1:3113697/69‑1 (MQ=255)
                                 cACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCag                                                              <  1:2557811/40‑1 (MQ=255)
                                 cACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCag                                                              <  1:1066656/40‑1 (MQ=255)
                                 cACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCag                                                              <  1:2661269/40‑1 (MQ=255)
                                 cACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCag                                                              <  1:1671903/40‑1 (MQ=255)
                                 cACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCag                                                              <  1:654655/40‑1 (MQ=255)
                                 cACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCag                                                              <  1:2830355/40‑1 (MQ=255)
                                 cACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCag                                                              <  1:2209977/40‑1 (MQ=255)
                                 cACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCag                                                              <  1:321292/40‑1 (MQ=255)
                                 cACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTcagc                                  <  1:921142/68‑1 (MQ=255)
                                   cGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCag                                                              <  1:29367/38‑1 (MQ=255)
                                   cGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTcagcag                                <  1:1520946/68‑1 (MQ=255)
                                   cGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTcagcag                                <  1:1433731/68‑1 (MQ=255)
                                    gTTTCGGACGCATCGGCGCATCTGGCGGAACAATCAGAGTACGTTCACCCTGTGCCATTTTCTGTAGTg                              >  1:2693012/1‑69 (MQ=255)
                                    gTTTCGGACGCATCGGCGCATCTGGCGGAACAATCAGAGTACGTTCACCCTGTGCCATTTTCAGTAGTg                              >  1:852048/1‑69 (MQ=255)
                                    gTTTCGGACGCATCGGCGCATCTGGCGGAACAATCAGAGTACGTTCACCCTGTGCCATTTTCAGTAGTg                              >  1:8987/1‑69 (MQ=255)
                                    gTTTCGGACGCATCGGCGCATCTGGCGGAACAATCAGAGTACGTTCACCCTGTGCCATTTTCAGTAGTg                              >  1:1091566/1‑69 (MQ=255)
                                    gTTTCGGACGCATCGGCGCATCTGGCGGAACAATCAGAGTACGTTCACCCTGTGCCATTTTCAGTAGTg                              >  1:101896/1‑69 (MQ=255)
                                    gTTTCGGACGCATCGGCGCATCTGGCGGAACAATCAGAGTACGTTCACCCTGTGCCATTTTCAGTAGTg                              >  1:1104855/1‑69 (MQ=255)
                                    gTTTCGGACGCATCGGCGCATCTGGCGGAACAATCAGAGTACGTTCACCCTGTGCCATTTTCAGTAGTg                              >  1:3339190/1‑69 (MQ=255)
                                    gTTTCGGACGCATCGGCGCATCTGGCGGAACAATCAGAGTACGTTCACCCTGTGCCATTTTCAGTAGTg                              >  1:2257388/1‑69 (MQ=255)
                                    gTTTCGGACGCATCGGCGCATCTGGCGGAACAATCAGAGTACGTTCACCCTGTGCCATTTTCAGTAGTg                              >  1:2304490/1‑69 (MQ=255)
                                    gTTTCGGACGCATCGGCGCATCTGGCGGAACAATCAGAGTACGTTCACCCTGTGCCATTTTCAGTAGTg                              >  1:2322896/1‑69 (MQ=255)
                                    gTTTCGGACGCATCGGCGCATCTGGCGGAACAATCAGAGTACGTTCACCCTGTGCCATTTTCAGTAGTg                              >  1:304912/1‑69 (MQ=255)
                                    gTTTCGGACGCATCGGCGCATCTGGCGGAACAATCAGAGTACGTTCACCCTGTGCCATTTTCAGTAGTg                              >  1:3003031/1‑69 (MQ=255)
                                    gTTTCGGACGCATCGGCGCATCTGGCGGAACAATCAGAGTACGTTCACCCTGTGCCATTTTCAGTAGTg                              >  1:2347318/1‑69 (MQ=255)
                                    gTTTCGGACGCATCGGCGCATCTGGCGGAACAATCAGAGTACGTTCACCCTGTGCCATTTTCAGTAGTg                              >  1:2009051/1‑69 (MQ=255)
                                    gTTTCGGACGCATCGGCGCATCTGGCGGAACAATCAGAGTACGTTCACCCTGTGCCATTTTCAGTAGTg                              >  1:2449025/1‑69 (MQ=255)
                                    gTTTCGGACGCATCGGCGCATCTGGCGGAACAATCAGAGTACGTTCACCCTGTGCCATTTTCAGTAGTg                              >  1:2087393/1‑69 (MQ=255)
                                                  ggCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCa                     <  1:1650478/64‑1 (MQ=255)
                                                  ggCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCa                     <  1:2096351/64‑1 (MQ=255)
                                                  ggCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCa                     <  1:779512/64‑1 (MQ=255)
                                                  ggCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCa                     <  1:1867179/64‑1 (MQ=255)
                                                  ggCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCa                     <  1:751095/64‑1 (MQ=255)
                                                  ggCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCa                     <  1:2663161/64‑1 (MQ=255)
                                                  ggCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCa                     <  1:2835371/64‑1 (MQ=255)
                                                  ggCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCa                     <  1:1797287/64‑1 (MQ=255)
                                                  ggCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCa                     <  1:1808993/64‑1 (MQ=255)
                                                  ggCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGtt                >  1:1178464/1‑69 (MQ=255)
                                                  ggCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGtt                >  1:188815/1‑69 (MQ=255)
                                                                tACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGctc  <  1:1651324/69‑1 (MQ=255)
                                                                tACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGctc  <  1:2187570/69‑1 (MQ=255)
                                                                tACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGctc  <  1:2257221/69‑1 (MQ=255)
                                                                |                                                                    
CGCGCGGACCACGGTCATCACGGTCACGGAATTCACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTC  >  minE/2131564‑2131696

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: