Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,930,802 G→A 9.1% T120T (ACC→ACT kduI ← predicted 5‑keto 4‑deoxyuronate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,930,8020GA 9.1% 200.4 / 8.7 66T120T (ACC→ACTkduIpredicted 5‑keto 4‑deoxyuronate isomerase
Reads supporting (aligned to +/- strand):  ref base G (27/33);  new base A (4/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

TGTGACTTTTTTGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGTGACCGATT  >  minE/1930733‑1930871
                                                                     |                                                                     
tgtgACTTTTTTGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGGCGGAGTGCCgg                                                                      <  1:2095783/71‑1 (MQ=255)
tgtgACTTTTTTGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCgg                                                                      <  1:907491/71‑1 (MQ=255)
tgtgACTTTTTTGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCgg                                                                      >  1:2825677/1‑71 (MQ=255)
tgtgACTTTTTTGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCgg                                                                      >  1:1168393/1‑71 (MQ=255)
 gtgACTTTTTTGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCgg                                                                      >  1:2970096/1‑70 (MQ=255)
  tgACTTTTTTGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCgg                                                                      <  1:2387075/69‑1 (MQ=255)
  tgACTTTTTTGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTa                                                                    <  1:3065485/71‑1 (MQ=255)
  tgACTTTTTTGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTa                                                                    <  1:3634700/71‑1 (MQ=255)
    aCTTTTTTGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTa                                                                    >  1:686461/1‑69 (MQ=255)
    aCTTTTTTGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTTTc                                                                  <  1:2685976/71‑1 (MQ=255)
     cTTTTTTGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCa                                                                 <  1:3178908/71‑1 (MQ=255)
     cTTTTTTGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCa                                                                 <  1:403974/71‑1 (MQ=255)
      ttttttGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCaa                                                                >  1:513165/1‑71 (MQ=255)
      ttttttGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCaa                                                                >  1:437467/1‑71 (MQ=255)
        ttttGGTGGGATACGTCGTATGCTCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCaa                                                                <  1:2277962/69‑1 (MQ=255)
        ttttGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCaa                                                                <  1:3387307/69‑1 (MQ=255)
        ttttGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCaa                                                                <  1:142619/69‑1 (MQ=255)
            ggtggGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTg                                                           <  1:81791/70‑1 (MQ=255)
            ggtggGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTg                                                           <  1:3431682/70‑1 (MQ=255)
                ggATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCa                                                        >  1:3319498/1‑69 (MQ=255)
                ggATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCa                                                        >  1:3469807/1‑69 (MQ=255)
                  aTACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGGCGGTATCAATACTGGCAAAGa                                                    >  1:1744217/1‑71 (MQ=255)
                   tACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGAc                                                   <  1:3226122/71‑1 (MQ=255)
                    aCGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCAGTATTGATAc                                                             >  1:449108/1‑60 (MQ=255)
                    aCGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCAGTATTGATAc                                                             >  1:1164403/1‑60 (MQ=255)
                     cgtcgtATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCaaa                                                      >  1:3605944/1‑66 (MQ=255)
                     cgtcgtATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACa                                                  <  1:2988316/70‑1 (MQ=255)
                        cgtATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAActt                                              >  1:1895446/1‑71 (MQ=255)
                         gtATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCa                                                        <  1:2423012/60‑1 (MQ=255)
                         gtATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAActtc                                             >  1:612928/1‑71 (MQ=255)
                          tATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAActtct                                            <  1:2671428/71‑1 (MQ=255)
                          tATGCGCGGGAGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAActtct                                            <  1:947642/71‑1 (MQ=255)
                            tGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAActtctt                                           <  1:565002/70‑1 (MQ=255)
                             gcgcGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGAc                                                   <  1:1352139/61‑1 (MQ=255)
                               gcgGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGc                                       <  1:2424189/71‑1 (MQ=255)
                                cgGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAAt                                                               <  1:3001819/46‑1 (MQ=255)
                                   gTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCt                                   <  1:2701017/71‑1 (MQ=255)
                                    tGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCtt                                  <  1:2224317/71‑1 (MQ=255)
                                      cgcAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCtttt                                         >  1:885397/1‑62 (MQ=255)
                                       gcAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTAc                               >  1:137173/1‑71 (MQ=255)
                                       gcAATTGTAATAAAACTTCGCCGGAGTGCCAGTATTGATAcc                                                            <  1:2577009/42‑2 (MQ=25)
                                       gcAATTGTAATAAAACTTCGCCGGAGTGCCAGTATTGATAcc                                                            <  1:2250560/42‑2 (MQ=25)
                                         aaTTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCttt                                 >  1:3481515/1‑67 (MQ=255)
                                         aaTTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTAcc                              >  1:3707324/1‑70 (MQ=255)
                                           ttGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCaa                            <  1:2236018/70‑1 (MQ=255)
                                           ttGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCaa                            >  1:2268791/1‑70 (MQ=255)
                                           ttGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCaa                            >  1:2807280/1‑70 (MQ=255)
                                                ataaAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATa                        >  1:2245411/1‑69 (MQ=255)
                                                ataaAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATa                        >  1:3545943/1‑69 (MQ=255)
                                                ataaAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACa                      <  1:2289020/71‑1 (MQ=255)
                                                ataaAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACa                      <  1:3182695/71‑1 (MQ=255)
                                                    aaCTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCttt                                          >  1:1703522/1‑47 (MQ=255)
                                                      cTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGt                >  1:1703172/1‑71 (MQ=255)
                                                        tCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTcg              >  1:3147437/1‑71 (MQ=255)
                                                           ccGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCttt                                 >  1:2257988/1‑49 (MQ=255)
                                                           ccGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCgg           >  1:608095/1‑71 (MQ=255)
                                                            cGGAGTGCCGGTATCAATACTGGCGAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGt          >  1:3392912/1‑71 (MQ=255)
                                                            cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCttt                                 >  1:3645956/1‑48 (MQ=255)
                                                            cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCttt                                 >  1:88924/1‑48 (MQ=255)
                                                            cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCttt                                 >  1:776939/1‑48 (MQ=255)
                                                            cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCTACATACAGGGCGTCGCGGt          <  1:2523111/71‑1 (MQ=255)
                                                            cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGt          <  1:458369/71‑1 (MQ=255)
                                                            cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGt          <  1:237062/71‑1 (MQ=255)
                                                            cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGt          <  1:2336836/71‑1 (MQ=255)
                                                            cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGt          <  1:2407159/71‑1 (MQ=255)
                                                            cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGt          <  1:3167692/71‑1 (MQ=255)
                                                            cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGt          <  1:3684194/71‑1 (MQ=255)
                                                              gAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTAcc                              <  1:3445148/49‑1 (MQ=255)
                                                              gAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGTGa        >  1:3327075/1‑71 (MQ=255)
                                                                 tGCCAGTATTGATACCGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGTGAcc      >  1:3200235/1‑70 (MQ=38)
                                                                 tGCCAGTATTGATACCGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGTGAcc      >  1:1006146/1‑70 (MQ=38)
                                                                    cGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGTGACCGAtt  <  1:1066341/71‑1 (MQ=255)
                                                                     |                                                                     
TGTGACTTTTTTGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGTGACCGATT  >  minE/1930733‑1930871

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: