Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 801,126 A→G 7.5% M68T (ATG→ACG)  dhaH ← fused predicted dihydroxyacetone‑specific PTS enzyme HPr component and EI component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE801,1260AG 7.5% 347.2 / 14.7 119M68T (ATG→ACG) dhaHfused predicted dihydroxyacetone‑specific PTS enzyme HPr component and EI component
Reads supporting (aligned to +/- strand):  ref base A (61/49);  new base G (7/2);  total (68/51)
Fisher's exact test for biased strand distribution p-value = 2.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 129 total aligned reads displayed.

GGCGATCTCGGGAGCCAGCAATTCCAGCGCAGTTTCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACTTTGACGGCATCGGTACCAAT  >  minE/801059‑801196
                                                                   |                                                                      
  cgATCTCGGGAGCCAGCAATTCCAGCGCAGTTTCAGCACTCAATAATGCGCTACCCATATCCATCGTGAc                                                                    >  1:2160866/1‑70 (MQ=255)
   gATCTCGGGAGCCAGCAATTCCAGCGCAGTTTCAGCACTCAATAATGCGCTACCCATATCCATCGTGAc                                                                    >  1:3551180/1‑69 (MQ=255)
   gATCTCGGGAGCCAGCAATTCCAGCGCAGTTTCAGCACTCAATAATGCGCTACCCATATCCATCATGAc                                                                    >  1:1159457/1‑69 (MQ=255)
         gggAGCCAGCAATTCCAGCGCAGTTTCAGCACTCAATAATGCGCTACCCATATCCATCATGa                                                                     >  1:2871206/1‑62 (MQ=255)
         gggAGCCAGCAATTCCAGCGCAGTTTCAGCACTCAATAATGCGCTACCCATATCCATCATGTCCAGca                                                               <  1:1209561/68‑1 (MQ=255)
         gggAGCCAGCAATTCCAGCGCAGTTTCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGca                                                               <  1:3481863/68‑1 (MQ=255)
         gggAGCCAGCAATTCCAGCGCAGTTTCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGca                                                               <  1:1278213/68‑1 (MQ=255)
               cAGCAATTCCAGCGCAGTTTCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCgg                                                      <  1:1463193/71‑1 (MQ=255)
               cAGCAATTCCAGCGCAGTTTCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCgg                                                      <  1:1650071/71‑1 (MQ=255)
                aGCAATTCCAGCGCAGTTTCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCGGc                                                     <  1:445978/71‑1 (MQ=255)
                 gCAATTCCAGCGCAGTTTCAGCACTCAATAATGCGCTACCCATATCCATCATGa                                                                     >  1:1582205/1‑54 (MQ=255)
                 gCAATTCCAGCGCAGTTTCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCa                                                    >  1:673941/1‑71 (MQ=255)
                 gCAATTCCAGCGCAGTTTCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCa                                                    >  1:685978/1‑71 (MQ=255)
                   aaTTCCAGCGCAGTTTCAGCACTCAATAATGCGCTACCCATATCCATCGTGACCAGCACATGGTCGGCATc                                                  <  1:2184228/71‑1 (MQ=255)
                    aTTCCAGCGCAGTTTCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGca                                                               <  1:674217/57‑1 (MQ=255)
                    aTTCCAGCGCAGTTTCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCAt                                                   <  1:1462151/69‑1 (MQ=255)
                    aTTCCAGCGCAGTTTCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCa                                                 <  1:1501541/71‑1 (MQ=255)
                     ttCCAGCGCAGTTTCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCa                                                 <  1:1708141/70‑1 (MQ=255)
                     ttCCAGCGCAGTTTCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAg                                                <  1:3586419/71‑1 (MQ=255)
                     ttCCAGCGAAGTTTCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAg                                                >  1:445638/1‑71 (MQ=255)
                       ccAGCGCAGTTTCAGCACTCAATAATGCGCTACCCATATCCATCGTGACCAGCACATGGTCGGCATCAGc                                               >  1:2179112/1‑70 (MQ=255)
                        cAGCGCAGTTTCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCa                                              <  1:887023/70‑1 (MQ=255)
                               gTTTCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTc                                      <  1:1507776/71‑1 (MQ=255)
                                tttCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACa                                           <  1:2649150/65‑1 (MQ=255)
                                tttCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACa                                           <  1:3506391/65‑1 (MQ=255)
                                tttCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCg                                     <  1:2527558/71‑1 (MQ=255)
                                tttCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCg                                     <  1:106642/71‑1 (MQ=255)
                                  tCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGa                                         >  1:3468168/1‑65 (MQ=255)
                                  tCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGAtt                                       >  1:1655723/1‑67 (MQ=255)
                                  tCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGAtt                                       >  1:3165526/1‑67 (MQ=255)
                                  tCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGa                                    <  1:2859539/70‑1 (MQ=255)
                                  tCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGAt                                   >  1:1227506/1‑71 (MQ=255)
                                  tCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGAt                                   >  1:381224/1‑71 (MQ=255)
                                    aGCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATg                                  <  1:78793/70‑1 (MQ=255)
                                     gCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATgg                                                          <  1:2220263/45‑1 (MQ=255)
                                       aCTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCt                              <  1:591740/71‑1 (MQ=255)
                                           aataatGCGCTACCCATATCCATCGTGACCAGCACATGGTCGGCATCAGCa                                              >  1:3731530/1‑51 (MQ=255)
                                           aataatGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAAc                                            >  1:2143206/1‑53 (MQ=255)
                                           aataatGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGAtt                                       >  1:1818074/1‑58 (MQ=255)
                                           aataatGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGAtt                                       >  1:3485376/1‑58 (MQ=255)
                                           aataatGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGAt                                        >  1:774653/1‑57 (MQ=255)
                                           aataatGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTcc                            <  1:2047712/69‑1 (MQ=255)
                                           aataatGCGCTACCCATATCCATAATGACCAGCACATGGTCGGCATCAGCAACAGAt                                        >  1:2415612/1‑57 (MQ=255)
                                            ataatGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGAt                                        >  1:2042317/1‑56 (MQ=255)
                                            ataatGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCt                              >  1:2044881/1‑66 (MQ=255)
                                            ataatGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCt                              >  1:3599945/1‑66 (MQ=255)
                                            ataatGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCt                              >  1:1271943/1‑66 (MQ=255)
                                             taatGCGCTACCCATATCCATCATGACCAGCATATGGTCGGCa                                                    >  1:3410049/1‑43 (MQ=255)
                                             taatGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCa                                                    >  1:1266126/1‑43 (MQ=255)
                                             taatGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCa                                                    >  1:459933/1‑43 (MQ=255)
                                             taatGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCa                                                    >  1:185309/1‑43 (MQ=255)
                                             taatGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCa                                                    >  1:2717473/1‑43 (MQ=255)
                                             taatGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCAt                                                   <  1:133033/44‑1 (MQ=255)
                                             taatGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCAt                                                   <  1:2591845/44‑1 (MQ=255)
                                             taatGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCAt                                                   <  1:3448858/44‑1 (MQ=255)
                                             taatGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCAt                                                   <  1:3249726/44‑1 (MQ=255)
                                             taatGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTc                                      >  1:1142396/1‑57 (MQ=255)
                                             taatGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTc                                      >  1:2420308/1‑57 (MQ=255)
                                             taatGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTc                                      >  1:1070191/1‑57 (MQ=255)
                                               atGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACt                      <  1:2916908/71‑1 (MQ=255)
                                                tGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGc                                               <  1:738155/45‑1 (MQ=255)
                                                tGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTcc                            <  1:3488176/64‑1 (MQ=255)
                                                tGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCAt                          <  1:1110628/66‑1 (MQ=255)
                                                tGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCAt                          <  1:1435284/66‑1 (MQ=255)
                                                 gcgcTACCCATATCCATCGTGACCAGCACATGGTCg                                                       <  1:1469044/36‑1 (MQ=255)
                                                 gcgcTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGc                                <  1:3022894/59‑1 (MQ=255)
                                                 gcgcTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGc                                <  1:3646633/59‑1 (MQ=255)
                                                  cgcTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGc                                >  1:805568/1‑58 (MQ=255)
                                                  cgcTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCAt                          >  1:3543485/1‑64 (MQ=255)
                                                  cgcTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCAt                          >  1:1728766/1‑64 (MQ=255)
                                                  cgcTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCAt                          >  1:2117476/1‑64 (MQ=255)
                                                  cgcTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCAt                          >  1:683084/1‑64 (MQ=255)
                                                  cgcTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACTTTg                   <  1:3201915/71‑1 (MQ=255)
                                                     tACCCATATCCATCGTGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACTTTGACg                >  1:1715883/1‑71 (MQ=255)
                                                       cccATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACTTTGACgg               >  1:3695992/1‑70 (MQ=255)
                                                       cccATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACTTTGACgg               >  1:3347478/1‑70 (MQ=255)
                                                         cATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACTTTGACGGCAt            >  1:2516019/1‑71 (MQ=255)
                                                         cATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACTTTGACGGCAt            <  1:424459/71‑1 (MQ=255)
                                                         cATATCCATCATGACCAG‑‑CATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACTTTGACGGCATCg          >  1:2252211/1‑71 (MQ=255)
                                                          atatCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACTTTGACGGCATc           >  1:1887152/1‑71 (MQ=255)
                                                          atatCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACTTTGACGGCATc           >  1:1678394/1‑71 (MQ=255)
                                                          atatCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACTTTGACGGCATc           >  1:555215/1‑71 (MQ=255)
                                                           tatCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCa                           <  1:2755465/54‑1 (MQ=255)
                                                           tatCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACTTTGACGGCATCg          >  1:1672378/1‑71 (MQ=255)
                                                            atCCATCATGACCAGCACATGGTCGGCATCAGCAAc                                            >  1:1052530/1‑36 (MQ=255)
                                                            atCCATCATGACCAGCACATGGTCGGCATCAGCAAc                                            >  1:2476923/1‑36 (MQ=255)
                                                            atCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACTTTGACGGCATCgg         >  1:3007656/1‑71 (MQ=255)
                                                            atCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACTTTGACGGCATCgg         >  1:3595371/1‑71 (MQ=255)
                                                             tCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATTGCCTCCATCACTTTGACGGCATCGGt        <  1:3536697/71‑1 (MQ=255)
                                                             tCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACttt                    <  1:2725532/59‑1 (MQ=255)
                                                             tCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACttt                    <  1:2551865/59‑1 (MQ=255)
                                                             tCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACTTTGACGGCATCGGt        <  1:3506395/71‑1 (MQ=255)
                                                             tCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACTTTGACGGCATCGGt        <  1:3509256/71‑1 (MQ=255)
                                                             tCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACTTTGACGGCATCGGt        <  1:826802/71‑1 (MQ=255)
                                                              ccATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACTTTGACGGCATCGGTa       >  1:151769/1‑71 (MQ=255)
                                                               catcatGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACTTTGACGGCATCGGTAc      >  1:3042573/1‑71 (MQ=255)
                                                               catcatGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCAATCACTTTGACGGCATCGGTAc      >  1:1378928/1‑71 (MQ=255)
                                                                 tcatGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACTTTGACGGCATCGGTAcc     >  1:183895/1‑70 (MQ=255)
                                                                   atGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACTTTGACGGCATCGGTACCAAt  >  1:3233709/1‑71 (MQ=255)
                                                                   atGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACTTTGACGGCAACGGTACCAAt  >  1:540175/1‑71 (MQ=255)
                                                                   |                                                                      
GGCGATCTCGGGAGCCAGCAATTCCAGCGCAGTTTCAGCACTCAATAATGCGCTACCCATATCCATCATGACCAGCACATGGTCGGCATCAGCAACAGATTCGATGGCCTCCATCACTTTGACGGCATCGGTACCAAT  >  minE/801059‑801196

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: