Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,920,413 C→T 26.5% V163I (GTC→ATC) ‡ ygeD ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,920,4130CT26.5% 100.5 / 31.5 88V163I (GTC→ATC) ‡ygeDpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (35/28);  new base T (4/19);  total (41/47)
Fisher's exact test for biased strand distribution p-value = 2.81e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

TGTAGATATTGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTC  >  minE/1920350‑1920479
                                                               |                                                                  
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:422623/68‑1 (MQ=255)
      tatTGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTca                                                          <  1:846421/68‑1 (MQ=255)
      tatTGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTca                                                          <  1:3271583/68‑1 (MQ=255)
      tatTGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTca                                                          <  1:1623960/68‑1 (MQ=255)
      tatTGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTca                                                          <  1:965174/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:540462/1‑69 (MQ=255)
          ggCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTcagccag                                                     >  1:2281676/1‑69 (MQ=255)
          ggCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTcagccag                                                     >  1:2128714/1‑69 (MQ=255)
          ggCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTcagccag                                                     >  1:419981/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: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:3389164/42‑1 (MQ=255)
                         accGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATg                                                                 <  1:248103/42‑1 (MQ=255)
                         accGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATg                                                                 <  1:599652/42‑1 (MQ=255)
                         accGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATg                                                                 <  1:1977203/42‑1 (MQ=255)
                         accGTAGGCCAGTGCGCATGCGGCCAGGGCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTa                                       >  1:1182093/1‑68 (MQ=255)
                         accGTAGGCCAGTGCGCATGCGGCCAGGGCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTa                                       >  1:2559912/1‑68 (MQ=255)
                         accGTAGGCCAGTGCGCATGCGGCCAGGGCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTa                                       >  1:1055469/1‑68 (MQ=255)
                         accGTAGGCCAGTGCGCATGCGGCCAGGGCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTa                                       >  1:3017972/1‑68 (MQ=255)
                                     tGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCaa                          >  1:499807/1‑69 (MQ=255)
                                     tGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCaa                          >  1:2356292/1‑69 (MQ=255)
                                     tGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCaa                          >  1:2330254/1‑69 (MQ=255)
                                     tGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCaa                          >  1:1937302/1‑69 (MQ=255)
                                     tGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCaa                          >  1:1900454/1‑69 (MQ=255)
                                     tGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCaa                          >  1:1850487/1‑69 (MQ=255)
                                     tGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCaa                          >  1:1286423/1‑69 (MQ=255)
                                     tGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCaa                          >  1:1180425/1‑69 (MQ=255)
                                     tGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCaa                          >  1:75424/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:2510034/1‑68 (MQ=255)
                                             cGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATcgcc                   >  1:1815907/1‑68 (MQ=255)
                                             cGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATcgcc                   >  1:1457340/1‑68 (MQ=255)
                                             cGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATcgcc                   >  1:724788/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:28394/63‑1 (MQ=255)
                                                   gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:77658/63‑1 (MQ=255)
                                                   gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:775523/63‑1 (MQ=255)
                                                   gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:728535/63‑1 (MQ=255)
                                                   gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:1200869/63‑1 (MQ=255)
                                                   gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:618782/63‑1 (MQ=255)
                                                   gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:1608864/63‑1 (MQ=255)
                                                   gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:1612658/63‑1 (MQ=255)
                                                   gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:506493/63‑1 (MQ=255)
                                                   gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:1966193/63‑1 (MQ=255)
                                                   gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:2283442/63‑1 (MQ=255)
                                                   gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:2467480/63‑1 (MQ=255)
                                                   gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:805889/63‑1 (MQ=255)
                                                   gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:3294457/63‑1 (MQ=255)
                                                   gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:3219201/63‑1 (MQ=255)
                                                   gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:2752026/63‑1 (MQ=255)
                                                   gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:2821113/63‑1 (MQ=255)
                                                   gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                  <  1:3038177/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:2641181/69‑1 (MQ=255)
                                                             gACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTc  <  1:3398338/69‑1 (MQ=255)
                                                             gACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTc  <  1:2690746/69‑1 (MQ=255)
                                                             gACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTc  <  1:2544185/69‑1 (MQ=255)
                                                             gACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTc  <  1:1617154/69‑1 (MQ=255)
                                                             gACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTc  <  1:1602514/69‑1 (MQ=255)
                                                             gACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTc  <  1:1595725/69‑1 (MQ=255)
                                                              aCATGCCAGTCAGCCAGCACACCACCGGCTACGAAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTc  <  1:97343/68‑1 (MQ=255)
                                                               |                                                                  
TGTAGATATTGGCAACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTC  >  minE/1920350‑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: