Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,930,808 A→G 9.8% I118I (ATT→ATC kduI ← predicted 5‑keto 4‑deoxyuronate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,930,8080AG 9.8% 171.5 / 9.1 61I118I (ATT→ATCkduIpredicted 5‑keto 4‑deoxyuronate isomerase
Reads supporting (aligned to +/- strand):  ref base A (27/28);  new base G (4/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-01

CTTTTTTGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGTGACCGATTTCATAG  >  Exported/1930738‑1930877
                                                                      |                                                                     
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:3431682/70‑1 (MQ=255)
       ggtggGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTg                                                                 <  1:81791/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:2250560/42‑2 (MQ=25)
                                  gcAATTGTAATAAAACTTCGCCGGAGTGCCAGTATTGATAcc                                                                  <  1:2577009/42‑2 (MQ=25)
                                    aaTTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCttt                                       >  1:3481515/1‑67 (MQ=255)
                                    aaTTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTAcc                                    >  1:3707324/1‑70 (MQ=255)
                                      ttGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCaa                                  >  1:2268791/1‑70 (MQ=255)
                                      ttGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCaa                                  >  1:2807280/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:3182695/71‑1 (MQ=255)
                                           ataaAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACa                            <  1:2289020/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:88924/1‑48 (MQ=255)
                                                       cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCttt                                       >  1:776939/1‑48 (MQ=255)
                                                       cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCttt                                       >  1:3645956/1‑48 (MQ=255)
                                                       cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCTACATACAGGGCGTCGCGGt                <  1:2523111/71‑1 (MQ=255)
                                                       cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGt                <  1:2336836/71‑1 (MQ=255)
                                                       cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGt                <  1:458369/71‑1 (MQ=255)
                                                       cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGt                <  1:237062/71‑1 (MQ=255)
                                                       cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGt                <  1:2407159/71‑1 (MQ=255)
                                                       cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGt                <  1:3684194/71‑1 (MQ=255)
                                                       cGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGt                <  1:3167692/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:212928/1‑71 (MQ=255)
                                                                     cAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGTGACCGATTTCATAg  >  1:1565157/1‑71 (MQ=255)
                                                                     cAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGTGACCGATTTCATAg  >  1:142991/1‑71 (MQ=255)
                                                                      aaTACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACa                            >  1:2719368/1‑44 (MQ=255)
                                                                      aaTACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGTGACCGATTTCATAg  <  1:3266815/70‑1 (MQ=255)
                                                                      |                                                                     
CTTTTTTGGTGGGATACGTCGTATGCGCGGGTGCGCAATTGTAATAAAACTTCGCCGGAGTGCCGGTATCAATACTGGCAAAGACAACTTCTTTTGCACCTTTACCAACATACAGGGCGTCGCGGTGACCGATTTCATAG  >  Exported/1930738‑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: