Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,393,956 A→G 14.2% G324G (GGT→GGC fruB ← fused fructose‑specific PTS enzyme IIA component and HPr component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,393,9560AG14.2% 195.4 / 16.9 111G324G (GGT→GGCfruBfused fructose‑specific PTS enzyme IIA component and HPr component
Reads supporting (aligned to +/- strand):  ref base A (27/68);  new base G (9/7);  total (36/75)
Fisher's exact test for biased strand distribution p-value = 4.18e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

Only 100 of 117 total aligned reads displayed.

ACCTTTCTTAACGCCAAGTGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGT  >  minE/1393893‑1394018
                                                               |                                                              
aCCTTTCTTAACGCCAAGTGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTGCCATc                                                           <  1:2173976/69‑1 (MQ=255)
aCCTTTCTTAACGCCAAGTGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTGCCATc                                                           <  1:2150981/69‑1 (MQ=255)
aCCTTTCTTAACGCCAAGTGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTGCCATc                                                           <  1:1099189/69‑1 (MQ=255)
aCCTTTCTTAACGCCAAGTGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTGCCATc                                                           <  1:3293314/69‑1 (MQ=255)
aCCTTTCTTAACGCCAAGTGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTGCCATc                                                           <  1:1008495/69‑1 (MQ=255)
aCCTTTCTTAACGCCAAGTGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCa                                                             >  1:2856045/1‑67 (MQ=255)
aCCTTTCTTAACGCCAAGTGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCa                                                             >  1:232483/1‑67 (MQ=255)
aCCTTTCTTAACGCCAAGTGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCa                                                             >  1:2616389/1‑67 (MQ=255)
aCCTTTCTTAACGCCAAGTGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCa                                                             >  1:2007292/1‑67 (MQ=255)
aCCTTTCTTAACGCCAAGTGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCa                                                             >  1:391655/1‑67 (MQ=255)
aCCTTTCTTAACGCCAAGTGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCa                                                             >  1:87138/1‑67 (MQ=255)
aCCTTTCTTAACGCCAAGTGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCa                                                             >  1:2930128/1‑67 (MQ=255)
aCCTTTCTTAACGCCAAGTGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCa                                                             >  1:1265400/1‑67 (MQ=255)
aCCTTTCTTAACGCCAAGTGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCa                                                             >  1:3114738/1‑67 (MQ=255)
aCCTTTCTTAACGCCAAGTGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCa                                                             >  1:1187653/1‑67 (MQ=255)
 ccTTTCTTAACGCCAAGTGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTGCCATc                                                           <  1:1364979/68‑1 (MQ=255)
                  tGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTGCCATCAAGGTTTGTCACGGTaa                                          >  1:3146331/1‑68 (MQ=255)
                  tGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTGCCATCAAGGTTTGTCACGGTaa                                          >  1:1698639/1‑68 (MQ=255)
                  tGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTGCCATCAAGGTTTGTCACGGTa                                           >  1:3213773/1‑67 (MQ=255)
                  tGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTGCCATCAAGGTTTGTCACGGTAat                                         >  1:1210322/1‑69 (MQ=255)
                  tGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTGCCATCAAGGTTTGTCACGGTAat                                         >  1:2435822/1‑69 (MQ=255)
                  tGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTGCCATCAAGGTTTGTCACGGTAat                                         >  1:1822763/1‑69 (MQ=255)
                  tGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTGCCATCAAGGTTTGTCACGGTAat                                         >  1:3436239/1‑69 (MQ=255)
                  tGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTGCCATCAAGGTTTGTCACGGTAat                                         >  1:140077/1‑69 (MQ=255)
                   gCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACgg                                             >  1:2836758/1‑64 (MQ=255)
                   gCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGt                                            >  1:3272748/1‑65 (MQ=255)
                   gCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGt                                            >  1:1088042/1‑65 (MQ=255)
                   gCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGt                                            >  1:1201144/1‑65 (MQ=255)
                      acaACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCa                                     <  1:1117237/69‑1 (MQ=255)
                      acaACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCa                                     <  1:7921/69‑1 (MQ=255)
                      acaACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCa                                     <  1:25253/69‑1 (MQ=255)
                      acaACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCa                                     <  1:2960815/69‑1 (MQ=255)
                      acaACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCa                                     <  1:771865/69‑1 (MQ=255)
                      acaACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCa                                     <  1:202005/69‑1 (MQ=255)
                      acaACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCa                                     <  1:2649851/69‑1 (MQ=255)
                      acaACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCa                                     <  1:432434/69‑1 (MQ=255)
                      acaACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCa                                     <  1:1801841/69‑1 (MQ=255)
                      acaACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCa                                     <  1:3165693/69‑1 (MQ=255)
                      acaACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCa                                     <  1:2419822/69‑1 (MQ=255)
                                    aCTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGttt                       <  1:2378930/69‑1 (MQ=255)
                                    aCTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGttt                       <  1:2689776/69‑1 (MQ=255)
                                    aCTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGttt                       <  1:2489472/69‑1 (MQ=255)
                                    aCTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGttt                       <  1:2493463/69‑1 (MQ=255)
                                    aCTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGttt                       <  1:2558882/69‑1 (MQ=255)
                                    aCTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGttt                       <  1:200308/69‑1 (MQ=255)
                                    aCTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGttt                       <  1:1598096/69‑1 (MQ=255)
                                    aCTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGttt                       <  1:2905458/69‑1 (MQ=255)
                                    aCTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGttt                       <  1:1058058/69‑1 (MQ=255)
                                    aCTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGttt                       <  1:3063973/69‑1 (MQ=255)
                                      tACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTaa                     >  1:255976/1‑69 (MQ=255)
                                      tACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTaa                     >  1:2216680/1‑69 (MQ=255)
                                      tACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTaa                     >  1:2342690/1‑69 (MQ=255)
                                      tACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTaa                     >  1:561351/1‑69 (MQ=255)
                                      tACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTaa                     >  1:2256017/1‑69 (MQ=255)
                                      tACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTaa                     >  1:2263313/1‑69 (MQ=255)
                                          tCCGTTTGCCGTTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATgt                  <  1:610598/68‑2 (MQ=255)
                                          tCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATgg                  <  1:536153/68‑1 (MQ=255)
                                          tCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATgg                  <  1:2688030/68‑1 (MQ=255)
                                          tCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATgg                  <  1:378372/68‑1 (MQ=255)
                                          tCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATgg                  <  1:277118/68‑1 (MQ=255)
                                          tCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATgg                  <  1:2855715/68‑1 (MQ=255)
                                            cGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGttt                       <  1:2429066/61‑1 (MQ=255)
                                            cGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGttt                       <  1:2259325/61‑1 (MQ=255)
                                            cGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGttt                       <  1:1168559/61‑1 (MQ=255)
                                            cGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGttt                       <  1:124710/61‑1 (MQ=255)
                                            cGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGttt                       <  1:1275859/61‑1 (MQ=255)
                                            cGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGttt                       <  1:1398335/61‑1 (MQ=255)
                                            cGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGttt                       <  1:1583240/61‑1 (MQ=255)
                                            cGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGttt                       <  1:1600297/61‑1 (MQ=255)
                                            cGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGttt                       <  1:278784/61‑1 (MQ=255)
                                            cGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGttt                       <  1:2568582/61‑1 (MQ=255)
                                            cGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTCttt                       <  1:2253676/61‑1 (MQ=255)
                                             gtttgccggtttgccgGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGttt                       <  1:1023566/60‑1 (MQ=255)
                                                  ccggtttgccggtACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGAcc          <  1:2407811/68‑1 (MQ=255)
                                                  ccggtttgccggtACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGAcc          <  1:3178613/68‑1 (MQ=255)
                                                  ccggtttgccggtACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGAcc          <  1:2139143/68‑1 (MQ=255)
                                                  ccggtttgccggtACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGAcc          <  1:1982114/68‑1 (MQ=255)
                                                    ggtttgccggtACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAg        <  1:109806/68‑1 (MQ=255)
                                                    ggtttgccggtACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAg        <  1:1512694/68‑1 (MQ=255)
                                                    ggtttgccggtACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAg        <  1:2223250/68‑1 (MQ=255)
                                                    ggtttgccggtACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAg        <  1:3321462/68‑1 (MQ=255)
                                                    ggtttgccggtACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAg        <  1:3418115/68‑1 (MQ=255)
                                                    ggtttgccggtACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAg        <  1:173001/68‑1 (MQ=255)
                                                      tttgccggtACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCAt     <  1:3051139/69‑1 (MQ=255)
                                                      tttgccggtACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCAt     <  1:416458/69‑1 (MQ=255)
                                                      tttgccggtACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCAt     <  1:1516962/69‑1 (MQ=255)
                                                      tttgccggtACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCAt     <  1:3328257/69‑1 (MQ=255)
                                                      tttgccggtACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCAt     <  1:234508/69‑1 (MQ=255)
                                                      tttgccggtACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCAt     <  1:765896/69‑1 (MQ=255)
                                                      tttgccggtACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCAt     <  1:3089527/69‑1 (MQ=255)
                                                         ggcggtACCATCAAGGTTTGTCACGTTAATATCACTGTTAAATTGttt                       <  1:1821914/46‑1 (MQ=255)
                                                          ccggtACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGt  <  1:604908/68‑1 (MQ=255)
                                                          ccggtACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGt  <  1:894808/68‑1 (MQ=255)
                                                          ccggtACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGt  <  1:506802/68‑1 (MQ=255)
                                                               aCCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGTTATTg             <  1:1233754/52‑1 (MQ=255)
                                                               aCCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTg             <  1:3431198/52‑1 (MQ=255)
                                                               aCCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTg             <  1:2962961/52‑1 (MQ=255)
                                                               aCCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTg             <  1:2638136/52‑1 (MQ=255)
                                                               aCCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTg             <  1:2483929/52‑1 (MQ=255)
                                                               aCCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTg             <  1:1450631/52‑1 (MQ=255)
                                                               |                                                              
ACCTTTCTTAACGCCAAGTGCCACAACTTTCATCAGACTACGTCCGTTTGCCGGTTTGCCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGT  >  minE/1393893‑1394018

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: