Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,920,429 C→A 35.3% V157V (GTG→GTT ygeD ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,920,4290CA35.3% 63.7 / 46.5 65V157V (GTG→GTTygeDpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (25/17);  new base A (4/19);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 1.54e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

AACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTCCATTAAACC  >  minE/1920363‑1920488
                                                                  |                                                           
aaCGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGcac                                                           <  1:482037/69‑1 (MQ=255)
      cGCACCACCGTAGGCCAGTTCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCAcaccac                                                      >  1:288377/1‑68 (MQ=255)
      cGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGcac                                                           <  1:1820520/63‑1 (MQ=255)
      cGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACaccacc                                                     >  1:175435/1‑69 (MQ=255)
      cGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACaccacc                                                     >  1:2805444/1‑69 (MQ=255)
            accGTAGGCCAGTGCGCATGCGGCCAGGGCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTa                                                >  1:1055469/1‑68 (MQ=255)
            accGTAGGCCAGTGCGCATGCGGCCAGGGCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTa                                                >  1:2559912/1‑68 (MQ=255)
            accGTAGGCCAGTGCGCATGCGGCCAGGGCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTa                                                >  1:1182093/1‑68 (MQ=255)
            accGTAGGCCAGTGCGCATGCGGCCAGGGCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTa                                                >  1:3017972/1‑68 (MQ=255)
                        tGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCaa                                   >  1:2330254/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:75424/1‑69 (MQ=255)
                        tGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCaa                                   >  1:1286423/1‑69 (MQ=255)
                        tGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCaa                                   >  1:1180425/1‑69 (MQ=255)
                        tGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCaa                                   >  1:499807/1‑69 (MQ=255)
                        tGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCaa                                   >  1:2356292/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:1457340/1‑68 (MQ=255)
                                cGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATcgcc                            >  1:724788/1‑68 (MQ=255)
                                cGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATcgcc                            >  1:1815907/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:500842/69‑1 (MQ=255)
                                  gCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCt                         <  1:1161667/69‑1 (MQ=255)
                                      gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                           <  1:2752026/63‑1 (MQ=255)
                                      gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                           <  1:2821113/63‑1 (MQ=255)
                                      gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                           <  1:28394/63‑1 (MQ=255)
                                      gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                           <  1:3294457/63‑1 (MQ=255)
                                      gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                           <  1:2467480/63‑1 (MQ=255)
                                      gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                           <  1:3219201/63‑1 (MQ=255)
                                      gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                           <  1:3038177/63‑1 (MQ=255)
                                      gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                           <  1:506493/63‑1 (MQ=255)
                                      gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                           <  1:805889/63‑1 (MQ=255)
                                      gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                           <  1:1200869/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:1608864/63‑1 (MQ=255)
                                      gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                           <  1:1612658/63‑1 (MQ=255)
                                      gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                           <  1:618782/63‑1 (MQ=255)
                                      gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                           <  1:1966193/63‑1 (MQ=255)
                                      gggCAACCAGAATATGCCAGTCAGCCAGAACACCACCAGCTACGGAACCGAGCAAAATcgccg                           <  1:2283442/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:2690746/69‑1 (MQ=255)
                                                gACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTc           <  1:1617154/69‑1 (MQ=255)
                                                gACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTc           <  1:2544185/69‑1 (MQ=255)
                                                gACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTc           <  1:1602514/69‑1 (MQ=255)
                                                gACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTc           <  1:3398338/69‑1 (MQ=255)
                                                gACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTc           <  1:1595725/69‑1 (MQ=255)
                                                 aCATGCCAGTCAGCCAGCACACCACCGGCTACGAAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTc           <  1:97343/68‑1 (MQ=255)
                                                      ccAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCt                         >  1:157762/1‑49 (MQ=255)
                                                      ccAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTa                        >  1:1142935/1‑50 (MQ=255)
                                                      ccAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTAt                       >  1:2925361/1‑51 (MQ=255)
                                                      ccAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTAt                       >  1:840263/1‑51 (MQ=255)
                                                      ccAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTAt                       >  1:1238103/1‑51 (MQ=255)
                                                      ccAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTAt                       >  1:1605094/1‑51 (MQ=255)
                                                      ccAGTCAGCCAGCACACCACCGGCTACGGAACAGAGCAAAATCGCCGCTAt                       >  1:2188868/1‑51 (MQ=255)
                                                           cagccagcACACCACCGGCTACGGAACCGAGCAAAATcgccgc                          <  1:2188270/43‑1 (MQ=255)
                                                           cagccagcACACCACCGGCTACGGAACCGAGCAAAATcgccgc                          <  1:520011/43‑1 (MQ=255)
                                                           cagccagcACACCACCGGCTACGGAACCGAGCAAAATcgccgc                          <  1:1889120/43‑1 (MQ=255)
                                                           cagccagcACACCACCGGCTACGGAACCGAGCAAAATcgccgc                          <  1:68454/43‑1 (MQ=255)
                                                           cagccagcACACCACCGGCTACGGAACCGAGCAAAATcgccgc                          <  1:1866859/43‑1 (MQ=255)
                                                           cagccagcACACCACCGGCTACGGAACCGAGCAAAATcgccgc                          <  1:1568162/43‑1 (MQ=255)
                                                           cagccagcACACCACCGGCTACGGAACCGAGCAAAATcgccgc                          <  1:2368766/43‑1 (MQ=255)
                                                           cagccagcACACCACCGGCTACGGAACCGAGCAAAATcgccgc                          <  1:1149006/43‑1 (MQ=255)
                                                           cagccagcACACCACCGGCTACGGAACCGAGCAAAATcgccgc                          <  1:2850871/43‑1 (MQ=255)
                                                           cagccagcACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTCCATTAAAcc  >  1:2244264/1‑67 (MQ=255)
                                                            agccagcACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTc           >  1:3033605/1‑57 (MQ=255)
                                                                  |                                                           
AACGACCGCACCACCGTAGGCCAGTGCGCATGCGGCCAGGGCGACGAGGACATGCCAGTCAGCCAGCACACCACCGGCTACGGAACCGAGCAAAATCGCCGCTATGGTAGAAGCTTCCATTAAACC  >  minE/1920363‑1920488

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: