Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,930,807 C→T 9.2% D119N (GAT→AAT)  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,8070CT 9.2% 197.5 / 8.6 65D119N (GAT→AAT) kduIpredicted 5‑keto 4‑deoxyuronate isomerase
Reads supporting (aligned to +/- strand):  ref base C (26/33);  new base T (4/2);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

ACTTTTTTGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGTGACCGATTTCATAG  >  minE/1930737‑1930877
                                                                      |                                                                      
aCTTTTTTGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTTTc                                                                        <  1:2685976/71‑1 (MQ=255)
 cTTTTTTGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCa                                                                       <  1:3178908/71‑1 (MQ=255)
 cTTTTTTGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCa                                                                       <  1:403974/71‑1 (MQ=255)
  ttttttGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCaa                                                                      >  1:437467/1‑71 (MQ=255)
  ttttttGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCaa                                                                      >  1:513165/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:2807280/1‑70 (MQ=255)
                                       ttGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCaa                                  >  1:2268791/1‑70 (MQ=255)
                                       ttGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCaa                                  <  1:2236018/70‑1 (MQ=255)
                                            ataaAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATa                              >  1:3545943/1‑69 (MQ=255)
                                            ataaAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATa                              >  1:2245411/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:776939/1‑48 (MQ=255)
                                                        cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCttt                                       >  1:88924/1‑48 (MQ=255)
                                                        cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCttt                                       >  1:3645956/1‑48 (MQ=255)
                                                        cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCTACATACAGGGCGTCGCGGt                <  1:2523111/71‑1 (MQ=255)
                                                        cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGt                <  1:3684194/71‑1 (MQ=255)
                                                        cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGt                <  1:3167692/71‑1 (MQ=255)
                                                        cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGt                <  1:2336836/71‑1 (MQ=255)
                                                        cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGt                <  1:458369/71‑1 (MQ=255)
                                                        cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGt                <  1:2407159/71‑1 (MQ=255)
                                                        cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGt                <  1:237062/71‑1 (MQ=255)
                                                          gAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTAcc                                    <  1:3445148/49‑1 (MQ=255)
                                                          gAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGTGa              >  1:3327075/1‑71 (MQ=255)
                                                             tGCCAGTATTGATACCGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGTGAcc            >  1:1006146/1‑70 (MQ=38)
                                                             tGCCAGTATTGATACCGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGTGAcc            >  1:3200235/1‑70 (MQ=38)
                                                                cGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGTGACCGAtt        <  1:1066341/71‑1 (MQ=255)
                                                                   tATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACTGGGCGTCGCGGTGACCGAttt       <  1:1482438/69‑1 (MQ=255)
                                                                   tATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGTGACCGATTTCa     <  1:914265/71‑1 (MQ=255)
                                                                     tCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCa                                   <  1:701104/39‑1 (MQ=255)
                                                                      cAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGTGACCGATTTCATAg  >  1:142991/1‑71 (MQ=255)
                                                                      cAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGTGACCGATTTCATAg  >  1:1565157/1‑71 (MQ=255)
                                                                      cAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGTGACCGATTTCATAg  >  1:212928/1‑71 (MQ=255)
                                                                      |                                                                      
ACTTTTTTGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGTGACCGATTTCATAG  >  minE/1930737‑1930877

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: