Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,330,278 C→T 45.9% S436S (TCG→TCA)  cpsG ← phosphomannomutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,330,2780CT45.9% 25.4 / 111.0 134S436S (TCG→TCA) cpsGphosphomannomutase
Reads supporting (aligned to +/- strand):  ref base C (38/34);  new base T (26/35);  total (65/69)
Fisher's exact test for biased strand distribution p-value = 2.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 134 total aligned reads displayed.

TTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTC  >  minE/1330215‑1330340
                                                               |                                                              
ttACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTc                                                           >  1:1073235/1‑69 (MQ=255)
ttACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTc                                                           <  1:2441627/69‑1 (MQ=255)
ttACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTc                                                           >  1:486463/1‑69 (MQ=255)
ttACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTc                                                           >  1:267097/1‑69 (MQ=255)
ttACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTc                                                           <  1:2786207/69‑1 (MQ=255)
ttACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTc                                                           >  1:1872264/1‑69 (MQ=255)
ttACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTc                                                           <  1:1609856/69‑1 (MQ=255)
 tACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTc                                                           >  1:640950/1‑68 (MQ=255)
 tACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTc                                                           >  1:2749429/1‑68 (MQ=255)
  aCTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGGGATTCca                                                         >  1:507298/1‑69 (MQ=255)
  aCTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCca                                                         >  1:2838448/1‑69 (MQ=255)
  aCTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCca                                                         >  1:2205710/1‑69 (MQ=255)
  aCTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCca                                                         >  1:565842/1‑69 (MQ=255)
  aCTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCca                                                         >  1:2991780/1‑69 (MQ=255)
  aCTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCca                                                         >  1:895333/1‑69 (MQ=255)
  aCTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCca                                                         >  1:3227270/1‑69 (MQ=255)
                     cagAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCgg                                      >  1:582077/1‑69 (MQ=255)
                     cagAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCgg                                      >  1:531817/1‑69 (MQ=255)
                     cagAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCg                                       >  1:3086243/1‑68 (MQ=255)
                     cagAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCg                                       >  1:1501596/1‑68 (MQ=255)
                         gttcgcgttcgcGCTTCCATCAGCGGCACATCACCGCGCGATTcc                                                          >  1:1267521/1‑45 (MQ=255)
                         gttcgcgttcgcGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGtt                                    >  1:42286/1‑67 (MQ=255)
                         gttcgcgttcgcGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGtt                                    >  1:372971/1‑67 (MQ=255)
                         gttcgcgttcgcGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTc                                   >  1:2824646/1‑68 (MQ=255)
                         gttcgcgttcgcGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTc                                   >  1:711438/1‑68 (MQ=255)
                         gttcgcgttcgcGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTc                                   >  1:37881/1‑68 (MQ=255)
                         gttcgcgttcgcGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTc                                   >  1:2576260/1‑68 (MQ=255)
                         gttcgcgttcgcGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTc                                   >  1:773716/1‑68 (MQ=255)
                          ttcgcgttcgcgCTTCCCTCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCgg                                 <  1:1206685/69‑1 (MQ=255)
                          ttcgcgttcgcgCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCgg                                 <  1:2416623/69‑1 (MQ=255)
                          ttcgcgttcgcgCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCgg                                 <  1:787169/69‑1 (MQ=255)
                          ttcgcgttcgcgCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCgg                                 <  1:2527699/69‑1 (MQ=255)
                          ttcgcgttcgcgCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCgg                                 <  1:2279158/69‑1 (MQ=255)
                          ttcgcgttcgcgCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCgg                                 <  1:656948/69‑1 (MQ=255)
                          ttcgcgttcgcgCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCgg                                 <  1:2254057/69‑1 (MQ=255)
                          ttcgcgttcgcgCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCgg                                 <  1:2660502/69‑1 (MQ=255)
                          ttcgcgttcgcgCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCgg                                 <  1:881552/69‑1 (MQ=255)
                          ttcgcgttcgcgCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCgg                                 <  1:1544138/69‑1 (MQ=255)
                          ttcgcgttcgcgCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCgg                                 <  1:719405/69‑1 (MQ=255)
                          ttcgcgttcgcgCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCg                                  >  1:489769/1‑68 (MQ=255)
                           tcgcgttcgcgCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGTCGCACCACCGGTTCGGt                                <  1:1583030/69‑1 (MQ=255)
                           tcgcgttcgcgCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCgg                                 <  1:2004585/68‑1 (MQ=255)
                           tcgcgttcgcgCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGt                                <  1:2359460/69‑1 (MQ=255)
                           tcgcgttcgcgCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGt                                <  1:525867/69‑1 (MQ=255)
                            cgcgttcgcgCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGt                                <  1:2074473/68‑1 (MQ=255)
                                 tcgcgCTTCCATCAGCGGCAGATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                          <  1:3253008/69‑1 (MQ=255)
                                 tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                          <  1:1890067/69‑1 (MQ=255)
                                 tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                          <  1:1267697/69‑1 (MQ=255)
                                 tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                          <  1:1393830/69‑1 (MQ=255)
                                 tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                          <  1:816503/69‑1 (MQ=255)
                                 tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                          <  1:3295237/69‑1 (MQ=255)
                                 tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                          <  1:3206035/69‑1 (MQ=255)
                                 tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                          <  1:1604972/69‑1 (MQ=255)
                                 tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                          <  1:2953898/69‑1 (MQ=255)
                                 tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                          <  1:1785917/69‑1 (MQ=255)
                                 tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                          <  1:2820640/69‑1 (MQ=255)
                                 tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                          <  1:2073922/69‑1 (MQ=255)
                                 tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                          <  1:663409/69‑1 (MQ=255)
                                 tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                          <  1:2728727/69‑1 (MQ=255)
                                 tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                          <  1:2280865/69‑1 (MQ=255)
                                 tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                          <  1:583703/69‑1 (MQ=255)
                                 tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                          <  1:227347/69‑1 (MQ=255)
                                 tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                          <  1:2186149/69‑1 (MQ=255)
                                 tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                          <  1:2697082/69‑1 (MQ=255)
                                 tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                          <  1:2095422/69‑1 (MQ=255)
                                 tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                          <  1:693399/69‑1 (MQ=255)
                                 tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                          <  1:2182274/69‑1 (MQ=255)
                                 tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACAGGTTCGGTGTTGGa                          <  1:3090356/69‑1 (MQ=255)
                                  cgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                          <  1:2184501/68‑1 (MQ=255)
                                  cgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                          <  1:2436783/68‑1 (MQ=255)
                                   gcgcTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                          <  1:1734240/67‑1 (MQ=255)
                                        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGt                                     >  1:2839560/1‑51 (MQ=255)
                                        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTTGAggag                      >  1:2504417/1‑64 (MQ=255)
                                        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                      >  1:760706/1‑64 (MQ=255)
                                        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                      >  1:790005/1‑64 (MQ=255)
                                        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                      >  1:803040/1‑64 (MQ=255)
                                        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                      >  1:979566/1‑64 (MQ=255)
                                        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                      >  1:1942829/1‑64 (MQ=255)
                                        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                      >  1:1078622/1‑64 (MQ=255)
                                        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                      >  1:1377872/1‑64 (MQ=255)
                                        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                      >  1:1547909/1‑64 (MQ=255)
                                        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                      >  1:1562179/1‑64 (MQ=255)
                                        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                      >  1:1991307/1‑64 (MQ=255)
                                        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                      >  1:2979991/1‑64 (MQ=255)
                                        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                      >  1:2707665/1‑64 (MQ=255)
                                        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                      >  1:2803474/1‑64 (MQ=255)
                                        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                      >  1:2558524/1‑64 (MQ=255)
                                        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAGGAgcg                    >  1:2480653/1‑68 (MQ=255)
                                        tCCATCAGCGGCACATCACAGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                      >  1:3002623/1‑64 (MQ=255)
                                                 ggCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAg          >  1:2710097/1‑69 (MQ=255)
                                                    acaTCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGa                          <  1:1614442/50‑1 (MQ=255)
                                                    acaTCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgcgc       <  1:1483328/69‑1 (MQ=255)
                                                    acaTCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgcgc       <  1:2551808/69‑1 (MQ=255)
                                                     caTCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgcgc       <  1:2003240/68‑1 (MQ=255)
                                                         aCCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTc  >  1:1764510/1‑69 (MQ=255)
                                                         aCCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTc  >  1:1560126/1‑69 (MQ=255)
                                                         aCCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTc  >  1:391454/1‑69 (MQ=255)
                                                         aCCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTc  >  1:831511/1‑69 (MQ=255)
                                                         aCCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTc  >  1:1156901/1‑69 (MQ=255)
                                                         aCCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTc  >  1:950078/1‑69 (MQ=255)
                                                               |                                                              
TTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTC  >  minE/1330215‑1330340

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: ‑