Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,920,411 G→A 27.3% V163V (GTC→GTT) ‡ ygeD ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,920,4110GA27.3% 92.9 / 34.4 84V163V (GTC→GTT) ‡ygeDpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (37/24);  new base A (4/19);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 4.99e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

TGGGAATGTAGATATTGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTC  >  minE/1920344‑1920479
                                                                   |                                                                    
tGGGAATGTAGATATTGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGa                                                                     <  1:2423409/69‑1 (MQ=255)
tGGGAATGTAGATATTGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGa                                                                     <  1:2843630/69‑1 (MQ=255)
tGGGAATGTAGATATTGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGa                                                                     <  1:3410844/69‑1 (MQ=255)
      tGTAGATATTGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGcc                                                               >  1:519682/1‑69 (MQ=255)
        tAGATATTGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGCCGAGGACATGCCa                                                              >  1:1798736/1‑68 (MQ=255)
        tAGATATTGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCa                                                              >  1:2458110/1‑68 (MQ=255)
        tAGATATTGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCa                                                              >  1:2863233/1‑68 (MQ=255)
        tAGATATTGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCa                                                              >  1:1321901/1‑68 (MQ=255)
        tAGATATTGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCa                                                              >  1:1592896/1‑68 (MQ=255)
        tAGATATTGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCa                                                              >  1:885825/1‑68 (MQ=255)
            tatTGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTca                                                          <  1:965174/68‑1 (MQ=255)
            tatTGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTca                                                          <  1:422623/68‑1 (MQ=255)
            tatTGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTca                                                          <  1:3271583/68‑1 (MQ=255)
            tatTGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTca                                                          <  1:846421/68‑1 (MQ=255)
            tatTGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTca                                                          <  1:1623960/68‑1 (MQ=255)
              ttGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTcagcc                                                       >  1:1752968/1‑69 (MQ=255)
              ttGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTcagcc                                                       >  1:776595/1‑69 (MQ=255)
              ttGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTcagcc                                                       >  1:3362523/1‑69 (MQ=255)
              ttGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTcagcc                                                       >  1:2824164/1‑69 (MQ=255)
                ggCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTcagccag                                                     >  1:669932/1‑69 (MQ=255)
                ggCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTcagccag                                                     >  1:30587/1‑69 (MQ=255)
                ggCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTcagccag                                                     >  1:3195191/1‑69 (MQ=255)
                ggCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTcagccag                                                     >  1:2281676/1‑69 (MQ=255)
                ggCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTcagccag                                                     >  1:540462/1‑69 (MQ=255)
                ggCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTcagccag                                                     >  1:2128714/1‑69 (MQ=255)
                ggCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTcagccag                                                     >  1:419981/1‑69 (MQ=255)
                ggCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTcagccag                                                     >  1:954076/1‑69 (MQ=255)
                   aaCGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGcac                                                  <  1:482037/69‑1 (MQ=255)
                         cGCACCACCGTAGGCCAGTTCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCAcaccac                                             >  1:288377/1‑68 (MQ=255)
                         cGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGcac                                                  <  1:1820520/63‑1 (MQ=255)
                         cGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACaccacc                                            >  1:2805444/1‑69 (MQ=255)
                         cGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACaccacc                                            >  1:175435/1‑69 (MQ=255)
                               accGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATg                                                                 <  1:1977203/42‑1 (MQ=255)
                               accGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATg                                                                 <  1:599652/42‑1 (MQ=255)
                               accGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATg                                                                 <  1:248103/42‑1 (MQ=255)
                               accGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATg                                                                 <  1:3389164/42‑1 (MQ=255)
                               accGTAGGCCAGTGCGCATGCGGCCAGGGCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTa                                       >  1:2559912/1‑68 (MQ=255)
                               accGTAGGCCAGTGCGCATGCGGCCAGGGCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTa                                       >  1:3017972/1‑68 (MQ=255)
                               accGTAGGCCAGTGCGCATGCGGCCAGGGCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTa                                       >  1:1055469/1‑68 (MQ=255)
                               accGTAGGCCAGTGCGCATGCGGCCAGGGCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTa                                       >  1:1182093/1‑68 (MQ=255)
                                 cGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGa                                                                     <  1:969802/36‑1 (MQ=255)
                                           tGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCaa                          >  1:2356292/1‑69 (MQ=255)
                                           tGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCaa                          >  1:1180425/1‑69 (MQ=255)
                                           tGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCaa                          >  1:1286423/1‑69 (MQ=255)
                                           tGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCaa                          >  1:75424/1‑69 (MQ=255)
                                           tGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCaa                          >  1:1850487/1‑69 (MQ=255)
                                           tGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCaa                          >  1:1900454/1‑69 (MQ=255)
                                           tGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCaa                          >  1:1937302/1‑69 (MQ=255)
                                           tGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCaa                          >  1:2330254/1‑69 (MQ=255)
                                           tGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCaa                          >  1:499807/1‑69 (MQ=255)
                                             cgcATGCGGCCAGGGCGCCGAGGAAATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCaaaa                        >  1:1479506/1‑69 (MQ=255)
                                             cgcATGCGGCCAGGGCGACGAGGAAATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCaaaa                        >  1:1483272/1‑69 (MQ=255)
                                                   cGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGa                                   <  1:251683/52‑1 (MQ=255)
                                                   cGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATcgcc                   >  1:1815907/1‑68 (MQ=255)
                                                   cGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATcgcc                   >  1:724788/1‑68 (MQ=255)
                                                   cGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATcgcc                   >  1:1457340/1‑68 (MQ=255)
                                                   cGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATcgcc                   >  1:2510034/1‑68 (MQ=255)
                                                     gCCAGGGCTACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCt                <  1:1088963/69‑1 (MQ=255)
                                                     gCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATcgccgc                 <  1:1612537/68‑1 (MQ=255)
                                                     gCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATcgccgc                 <  1:2169492/68‑1 (MQ=255)
                                                     gCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCt                <  1:3007670/69‑1 (MQ=255)
                                                     gCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCt                <  1:1161667/69‑1 (MQ=255)
                                                     gCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCt                <  1:500842/69‑1 (MQ=255)
                                                         gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:618782/63‑1 (MQ=255)
                                                         gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:3038177/63‑1 (MQ=255)
                                                         gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:1200869/63‑1 (MQ=255)
                                                         gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:1608864/63‑1 (MQ=255)
                                                         gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:805889/63‑1 (MQ=255)
                                                         gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:28394/63‑1 (MQ=255)
                                                         gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:77658/63‑1 (MQ=255)
                                                         gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:775523/63‑1 (MQ=255)
                                                         gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:1612658/63‑1 (MQ=255)
                                                         gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:728535/63‑1 (MQ=255)
                                                         gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:1966193/63‑1 (MQ=255)
                                                         gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:3294457/63‑1 (MQ=255)
                                                         gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:3219201/63‑1 (MQ=255)
                                                         gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:2283442/63‑1 (MQ=255)
                                                         gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:506493/63‑1 (MQ=255)
                                                         gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:2821113/63‑1 (MQ=255)
                                                         gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:2467480/63‑1 (MQ=255)
                                                         gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:2752026/63‑1 (MQ=255)
                                                         gggCAACCAGAATATGCCAGGCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:2716555/63‑1 (MQ=255)
                                                              acgaGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGaa       <  1:2147071/69‑1 (MQ=255)
                                                              acgaGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGaa       <  1:2396357/69‑1 (MQ=255)
                                                                   gACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTc  <  1:1617154/69‑1 (MQ=255)
                                                                   gACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTc  <  1:2544185/69‑1 (MQ=255)
                                                                   gACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTc  <  1:2641181/69‑1 (MQ=255)
                                                                   gACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTc  <  1:2690746/69‑1 (MQ=255)
                                                                   gACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTc  <  1:3398338/69‑1 (MQ=255)
                                                                   gACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTc  <  1:1602514/69‑1 (MQ=255)
                                                                   gACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTc  <  1:1595725/69‑1 (MQ=255)
                                                                   |                                                                    
TGGGAATGTAGATATTGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTC  >  minE/1920344‑1920479

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: