Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,622,837 T→C 35.6% E603E (GAA→GAG maeB ← fused malic enzyme predicted oxidoreductase and predicted phosphotransacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,622,8370TC35.6% 81.1 / 79.6 129E603E (GAA→GAGmaeBfused malic enzyme predicted oxidoreductase and predicted phosphotransacetylase
Reads supporting (aligned to +/- strand):  ref base T (30/53);  new base C (40/6);  total (70/59)
Fisher's exact test for biased strand distribution p-value = 1.46e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-02

Only 100 of 129 total aligned reads displayed.

CCAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGTATCGGCAATAAAGGTGTTACCACTCGGCAGCAGCAGCGCGTTCATGGC  >  minE/1622775‑1622902
                                                              |                                                                 
ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcatc                                                              >  1:1066055/1‑68 (MQ=255)
ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcatc                                                              >  1:1481388/1‑68 (MQ=255)
ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcatc                                                              >  1:2225783/1‑68 (MQ=255)
ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcatc                                                              >  1:1349014/1‑68 (MQ=255)
ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcatc                                                              >  1:2816531/1‑68 (MQ=255)
ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcatc                                                              >  1:2839139/1‑68 (MQ=255)
ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcatc                                                              >  1:3389829/1‑68 (MQ=255)
ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcatc                                                              >  1:739402/1‑68 (MQ=255)
ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcatc                                                              >  1:1140680/1‑68 (MQ=255)
ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcat                                                               >  1:1621865/1‑67 (MQ=255)
ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcat                                                               >  1:1814415/1‑67 (MQ=255)
ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcat                                                               >  1:2349717/1‑67 (MQ=255)
ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcat                                                               >  1:2503416/1‑67 (MQ=255)
  aaaaCGACGAACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCAGGCTCATCAtt                                                           <  1:2396684/69‑1 (MQ=255)
  aaaaCGACGAACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCAGGCTCATCAtt                                                           <  1:2011391/69‑1 (MQ=255)
  aaaaCGACGAACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCAGGCTCATCAtt                                                           <  1:2039153/69‑1 (MQ=255)
  aaaaCGACGAACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCAGGCTCATCAtt                                                           <  1:882526/69‑1 (MQ=255)
  aaaaCGACGAACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCAGGCTCATCAtt                                                           <  1:1318082/69‑1 (MQ=255)
      cgacgGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACa                                                       <  1:2706509/69‑1 (MQ=255)
      cgacgGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACa                                                       <  1:871230/69‑1 (MQ=255)
      cgacgGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACa                                                       <  1:1019936/69‑1 (MQ=255)
      cgacgGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACa                                                       <  1:2173916/69‑1 (MQ=255)
      cgacgGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACa                                                       <  1:287371/69‑1 (MQ=255)
      cgacgGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACa                                                       <  1:980800/69‑1 (MQ=255)
      cgacgGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACa                                                       <  1:3025529/69‑1 (MQ=255)
      cgacgGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACa                                                       <  1:3302927/69‑1 (MQ=255)
      cgacgGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACa                                                       <  1:3317568/69‑1 (MQ=255)
      cgacgGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACa                                                       <  1:1327998/69‑1 (MQ=255)
        acgGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACata                                                     <  1:602429/69‑1 (MQ=255)
           gACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGt                                                  <  1:1124038/69‑1 (MQ=255)
           gACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGt                                                  <  1:1333329/69‑1 (MQ=255)
           gACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGt                                                  <  1:3181425/69‑1 (MQ=255)
           gACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGt                                                  <  1:3385806/69‑1 (MQ=255)
           gACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGt                                                  <  1:391619/69‑1 (MQ=255)
           gACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGt                                                  <  1:439527/69‑1 (MQ=255)
            aCAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGt                                                  <  1:2280226/68‑1 (MQ=255)
                  tCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGTATCGGc                                            >  1:115445/1‑68 (MQ=255)
                  tCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGTATCGGc                                            >  1:1114487/1‑68 (MQ=255)
                        tccATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGTATCGGCAATaaa                                      <  1:790176/67‑1 (MQ=255)
                        gccATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGTATCGGCAATaaa                                      <  1:591290/68‑1 (MQ=255)
                        gccATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGTATCGGCAATaaa                                      <  1:287652/68‑1 (MQ=255)
                        gccATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGTATCGGCAATaaa                                      <  1:1084464/68‑1 (MQ=255)
                        gccATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGTATCGGCAATaaa                                      <  1:2828949/68‑1 (MQ=255)
                        gccATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGTATCGGCAATaaa                                      <  1:2490570/68‑1 (MQ=255)
                        gccATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGTATCGGCAATaaa                                      <  1:2470840/68‑1 (MQ=255)
                        gccATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGTATCGGCAATaaa                                      <  1:978578/68‑1 (MQ=255)
                        gccATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGTATCGGCAATaaa                                      <  1:1721937/68‑1 (MQ=255)
                        gacATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGTATCGGCAATaaa                                      <  1:2730452/66‑1 (MQ=255)
                              aaGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACa                                                       <  1:1349927/45‑1 (MQ=255)
                              aaGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGTATCGGCAATAAAGGTGtt                                >  1:1708511/1‑68 (MQ=255)
                              aaGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGTATCGGCAATAAAGGTGTTa                               >  1:2512554/1‑69 (MQ=255)
                                       tcCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGt                                                  <  1:1394466/41‑1 (MQ=255)
                                          gCCAGCTCTTCTGCATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGgctac                   >  1:1588938/1‑66 (MQ=255)
                                          gCCAGCTCTTCTGCATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGgcaac                   >  1:3076881/1‑67 (MQ=255)
                                          gCCAGCTCTTCTGCATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGgcaac                   >  1:3104571/1‑67 (MQ=255)
                                          gCCAGCTCTTCTGCATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGgcaac                   >  1:1176476/1‑67 (MQ=255)
                                          gCCAGCTCTTCTGCATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGgcaac                   >  1:3362321/1‑67 (MQ=255)
                                          gCCAGCTCTTCTGCATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGgcaac                   >  1:3329914/1‑67 (MQ=255)
                                          gCCAGCTCTTCTGCATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGgcaac                   >  1:1272121/1‑67 (MQ=255)
                                          gCCAGCTCTTCTGCATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGgcaac                   >  1:1889509/1‑67 (MQ=255)
                                          gCCAGCTCTTCTGCATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGgcaac                   >  1:1153807/1‑67 (MQ=255)
                                          gCCAGCTCTTCTGCATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGgcaac                   >  1:1098911/1‑67 (MQ=255)
                                          gCCAGCTCTTCTGCATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGgcaac                   >  1:2034656/1‑67 (MQ=255)
                                          gCCAGCTCTTCTGCATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGgcaac                   >  1:759708/1‑67 (MQ=255)
                                          gCCAGCTCTTCTGCATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGgcaac                   >  1:766845/1‑67 (MQ=255)
                                          gCCAGCTCTTCTGCATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGgcaac                   >  1:2145147/1‑67 (MQ=255)
                                          gCCAGCTCTTCTGCATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGgcaac                   >  1:879063/1‑67 (MQ=255)
                                          gCCAGCTCTTCTGCATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGgcaac                   >  1:247770/1‑67 (MQ=255)
                                          gCCAGCTCTTCTGCATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGgcaac                   >  1:225068/1‑67 (MQ=255)
                                          gCCAGCTCTTCTGCATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGgcaac                   >  1:2281199/1‑67 (MQ=255)
                                          gCCAGCTCTTCTGCATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGgcaac                   >  1:2383159/1‑67 (MQ=255)
                                          gCCAGCTCTTCTGCATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGgcaac                   >  1:236669/1‑67 (MQ=255)
                                          gCCAGCTCTTCTGCATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGgca                     >  1:144459/1‑67 (MQ=255)
                                          gCCAGCTCTTCTGCATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGgc                      >  1:1173903/1‑66 (MQ=255)
                                                    ctGCATCCGGTTCATCATTAACATATGTATCGGCAATAAAGGTGTTACCACTCGGCAGCAGCAGcgcg          >  1:471572/1‑68 (MQ=255)
                                                    ctGCATCCGGTTCATCATTAACATATGTATCGGCAATAAAGGTGTTACCACTCGGCAGCAGCAGcgcg          >  1:3161774/1‑68 (MQ=255)
                                                    ctGCATCCGGTTCATCATTAACATATGTATCGGCAATAAAGGTGTTACCACTCGGCAGCAGCAGcgcg          >  1:314644/1‑68 (MQ=255)
                                                    ctGCATCCGGTTCATCATTAACATATGTATCGGCAATAAAGGTGTTACCACTCGGCAGCAGCAGcgcg          >  1:2988808/1‑68 (MQ=255)
                                                    ctGCATCCGGTTCATCATTAACATATGTATCGGCAATAAAGGTGTTACCACTCGGCAGCAGCAGcgcg          >  1:960554/1‑68 (MQ=255)
                                                    ctGCATCCGGTTCATCATTAACATATGTATCGGCAATAAAGGTGTTACCACTCGGCAGCAGCAGcgcg          >  1:771615/1‑68 (MQ=255)
                                                      gCATCCGGTTCATCATTAACATATGTATCGGCAATGAAGGTGTTACCACTCGGCAGCAGCAGCGCGtt        <  1:1421549/68‑1 (MQ=255)
                                                      gCATCCGGTTCATCATTAACATATGTATCGGCAATAAAGGTGTTACCACTCGGCAGCAGCAGCGCGtt        <  1:1255160/68‑1 (MQ=255)
                                                      gCATCCGGTTCATCATTAACATATGTATCGGCAATAAAGGTGTTACCACTCGGCAGCAGCAGCGCGtt        <  1:570666/68‑1 (MQ=255)
                                                      gCATCCGGTTCATCATTAACATATGTATCGGCAATAAAGGTGTTACCACTCGGCAGCAGCAGCGCGtt        <  1:1868568/68‑1 (MQ=255)
                                                      gCATCCGGTTCATCATTAACATATGTATCGGCAATAAAGGTGTTACCACTCGGCAGCAGCAGCGCGtt        <  1:1517257/68‑1 (MQ=255)
                                                      gCATCCGGTTCATCATTAACATATGTATCGGCAATAAAGGTGTTACCACTCGGCAGCAGCAGCGCGtt        <  1:2189861/68‑1 (MQ=255)
                                                      gCATCCGGTTCATCATTAACATATGTATCGGCAATAAAGGTGTTACCACTCGGCAGCAGCAGCGCGtt        <  1:2716995/68‑1 (MQ=255)
                                                      gCATCCGGTTCATCATTAACATATGTATCGGCAATAAAGGTGTTACCACTCGGCAGCAGCAGCGCGtt        <  1:2576106/68‑1 (MQ=255)
                                                       cATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGGCAACAGCAGCGCGTTc       >  1:2961017/1‑68 (MQ=255)
                                                       cATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGGCAACAGCAGCGCGTTCa      >  1:1641911/1‑69 (MQ=255)
                                                       cATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGGCAACAGCAGCGCGTTCa      >  1:449800/1‑69 (MQ=255)
                                                       cATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGGCAACAGCAGCGCGTTCa      >  1:1906935/1‑69 (MQ=255)
                                                       cATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGGCAACAGCAGCGCGTTCa      >  1:806027/1‑69 (MQ=255)
                                                       cATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGGCAACAGCAGCGCGTTCa      >  1:263565/1‑69 (MQ=255)
                                                       cATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGGCAACAGCAGCGCGTTCa      >  1:1265427/1‑69 (MQ=255)
                                                       cATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGGCAACAGCAGCGCGTTCa      >  1:654662/1‑69 (MQ=255)
                                                       cATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGGCAACAGCAGCGCGTTCa      >  1:624280/1‑69 (MQ=255)
                                                       cATCAGGCTCATCATTGACGTAGGTATCGGCAATAAAGGTGTTACCACTCGGCAACAGCAGCGCGTTCa      >  1:1614387/1‑69 (MQ=255)
                                                           cGGTTCATCATTAACATATGTATCGGCAATAAAGgtgt                                 >  1:2807273/1‑38 (MQ=255)
                                                           cGGTTCATCATTAACATATGTATCGGCAATAAAGGTGTTACCACTCGGCAGCAGCAGCGCGTTCATGGc  <  1:1284407/69‑1 (MQ=255)
                                                              |                                                                 
CCAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGTATCGGCAATAAAGGTGTTACCACTCGGCAGCAGCAGCGCGTTCATGGC  >  minE/1622775‑1622902

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: