Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,982,097 A→G 9.8% H228H (CAT→CAC kduI ← predicted 5‑keto 4‑deoxyuronate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,982,0970AG 9.8% 160.3 / 7.8 61H228H (CAT→CACkduIpredicted 5‑keto 4‑deoxyuronate isomerase
Reads supporting (aligned to +/- strand):  ref base A (24/31);  new base G (2/4);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

GGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAG  >  W3110S.gb/2982031‑2982163
                                                                  |                                                                  
ggTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATATGa                                                                 >  1:2950898/1‑70 (MQ=255)
ggTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATATGAc                                                                >  1:712206/1‑71 (MQ=255)
    ccGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATATGACGCGt                                                            >  1:3197565/1‑71 (MQ=255)
          ccGGAATGGATCGACCAGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTg                                                      <  1:1957455/71‑1 (MQ=255)
            ggAATGGATCGACCAGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCg                                                    >  1:3217133/1‑71 (MQ=255)
             gAATGGATCGACCAGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCg                                                    <  1:3390002/70‑1 (MQ=255)
                  caTCGACCAGCTCGGGGAGATCACCGCCTGCTCCTTATGCATCACAATATGACGCGTTTCTTGCGGCTGcc                                              >  1:1922557/2‑71 (MQ=255)
                     cGACCAGCTCGGGGAGATCACCGCCTTCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCca                                           >  1:700590/1‑71 (MQ=255)
                      gACCAGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCcat                                          <  1:1467355/71‑1 (MQ=255)
                      gACCAGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCcat                                          <  1:2135335/71‑1 (MQ=255)
                       aCCAGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATGTGACGCGTTTCTTGCGGTTGCCCcatc                                         <  1:1250621/71‑1 (MQ=255)
                       aCCAGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATGTGACGCGTTTCTTGCGGTTGCCCcatc                                         <  1:754096/71‑1 (MQ=255)
                       aCCAGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATGTGACGCGTTTCTTGCGGTTGCCCcatc                                         <  1:3126847/71‑1 (MQ=255)
                       aCCAGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATGTGACGCGTTTCTTGCGGTTGCCCcatc                                         <  1:3435807/71‑1 (MQ=255)
                         cAGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCcatcat                                       <  1:1917973/71‑1 (MQ=255)
                          aGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCAtg                                      <  1:247520/71‑1 (MQ=255)
                               ggggAGATCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCca                                           <  1:1758451/61‑1 (MQ=255)
                               ggggAGATCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCca                                           <  1:1869622/61‑1 (MQ=255)
                               ggggAGATCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCca                                           <  1:2554145/61‑1 (MQ=255)
                                gggAGATCACCGCCTGCTCGTTATGCATCACAATATGACGCGttt                                                          <  1:2886719/45‑1 (MQ=255)
                                gggAGATCACCGCCTGCTCGTTATGCATCACAATATGACGCGttt                                                          <  1:2446488/45‑1 (MQ=255)
                                gggAGATCACCGCCTGCTCGTTATGCATCACAATATGACGCGttt                                                          <  1:2593172/45‑1 (MQ=255)
                                 ggAGATCACCGCCTGCTCGTTATGCATCACAATGTGACGCGTTTCTTGCGGTTGCCCCATCAtg                                      >  1:267016/1‑64 (MQ=255)
                                 ggAGATCACCGCCTGCTCGTTATGCATCACAATGTGACGCGTTTCTTGCGGTTGCCCCATCAtg                                      >  1:2670373/1‑64 (MQ=255)
                                   agaTCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACg                             <  1:3627914/71‑1 (MQ=255)
                                   agaTCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACg                             <  1:511130/71‑1 (MQ=255)
                                   agaTCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACg                             <  1:923364/71‑1 (MQ=255)
                                     aTCACCGCCTGCTCGTTATGCATCACAATATGACGCGttt                                                          <  1:2126587/40‑1 (MQ=255)
                                     aTCACCGCCTGCTCGTTATGCATCACAATATGACGCGttt                                                          <  1:2556396/40‑1 (MQ=255)
                                     aTCACCGCCTGCTCGTTATGCATCACAATATGACGCGttt                                                          <  1:121500/40‑1 (MQ=255)
                                      tCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCtt                                                       <  1:731640/42‑1 (MQ=255)
                                       cACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGc                                                  >  1:1340508/1‑46 (MQ=255)
                                       cACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGaaaa                               >  1:1438109/1‑65 (MQ=255)
                                       cACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATAATGTGGaaaa                               >  1:2660888/1‑65 (MQ=255)
                                        aCCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCAtgtga                                   <  1:977138/60‑2 (MQ=255)
                                        aCCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGa                                  <  1:1813915/61‑1 (MQ=255)
                                         ccGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCg                       <  1:97642/71‑1 (MQ=255)
                                         ccGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCg                       >  1:3472425/1‑71 (MQ=255)
                                           gCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGt                      >  1:3825520/1‑70 (MQ=255)
                                             cTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGtcat                   <  1:2163268/71‑1 (MQ=255)
                                                   gTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCatcatc               <  1:444529/69‑1 (MQ=255)
                                                   gTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCatcatc               <  1:143508/69‑1 (MQ=255)
                                                    ttATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCAtg                                      <  1:790316/45‑1 (MQ=255)
                                                    ttATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACg                             <  1:2377613/54‑1 (MQ=255)
                                                      aTGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACg                             >  1:1790018/1‑52 (MQ=255)
                                                      aTGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGt                      >  1:1318880/1‑59 (MQ=255)
                                                        gCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAgg                         >  1:1844444/1‑54 (MQ=255)
                                                        gCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAgg                         >  1:1894237/1‑54 (MQ=255)
                                                         cATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTga       <  1:2673878/71‑1 (MQ=255)
                                                           tCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTg        <  1:3235839/68‑1 (MQ=255)
                                                            cacaATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTgaaa     >  1:179857/1‑70 (MQ=255)
                                                             acaATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGaaaa                               >  1:1911417/1‑43 (MQ=255)
                                                             acaATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGc                            >  1:3708728/1‑46 (MQ=255)
                                                             acaATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGc                            >  1:1970514/1‑46 (MQ=255)
                                                             acaATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTgaaa     >  1:1841632/1‑69 (MQ=255)
                                                             acaATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTgaaa     >  1:1215307/1‑69 (MQ=255)
                                                              caATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGaaaa                               >  1:1630301/1‑42 (MQ=255)
                                                              caATATGACGCGTTTCTTGCGGCTGCCCAATCATGTGGaaaa                               >  1:2683962/1‑42 (MQ=255)
                                                               aaTATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTgaaatag  >  1:1999982/1‑70 (MQ=255)
                                                               aaTATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTgaaa     <  1:644686/67‑1 (MQ=255)
                                                               aaTATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTgaaa     <  1:1513197/67‑1 (MQ=255)
                                                                  |                                                                  
GGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAG  >  W3110S.gb/2982031‑2982163

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: