Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,982,115 C→T 10.1% Q222Q (CAG→CAA 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,1150CT10.1% 162.6 / 7.2 59Q222Q (CAG→CAAkduIpredicted 5‑keto 4‑deoxyuronate isomerase
Reads supporting (aligned to +/- strand):  ref base C (22/31);  new base T (2/4);  total (24/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.96e-01

GATCGACCAGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGG  >  W3110S.gb/2982049‑2982181
                                                                  |                                                                  
caTCGACCAGCTCGGGGAGATCACCGCCTGCTCCTTATGCATCACAATATGACGCGTTTCTTGCGGCTGcc                                                                >  1:1922557/2‑71 (MQ=255)
   cGACCAGCTCGGGGAGATCACCGCCTTCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCca                                                             >  1:700590/1‑71 (MQ=255)
    gACCAGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCcat                                                            <  1:2135335/71‑1 (MQ=255)
    gACCAGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCcat                                                            <  1:1467355/71‑1 (MQ=255)
     aCCAGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATGTGACGCGTTTCTTGCGGTTGCCCcatc                                                           <  1:3126847/71‑1 (MQ=255)
     aCCAGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATGTGACGCGTTTCTTGCGGTTGCCCcatc                                                           <  1:3435807/71‑1 (MQ=255)
     aCCAGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATGTGACGCGTTTCTTGCGGTTGCCCcatc                                                           <  1:754096/71‑1 (MQ=255)
     aCCAGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATGTGACGCGTTTCTTGCGGTTGCCCcatc                                                           <  1:1250621/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:2554145/61‑1 (MQ=255)
             ggggAGATCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCca                                                             <  1:1869622/61‑1 (MQ=255)
               ggAGATCACCGCCTGCTCGTTATGCATCACAATGTGACGCGTTTCTTGCGGTTGCCCCATCAtg                                                        >  1:267016/1‑64 (MQ=255)
               ggAGATCACCGCCTGCTCGTTATGCATCACAATGTGACGCGTTTCTTGCGGTTGCCCCATCAtg                                                        >  1:2670373/1‑64 (MQ=255)
                 agaTCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACg                                               <  1:923364/71‑1 (MQ=255)
                 agaTCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACg                                               <  1:511130/71‑1 (MQ=255)
                 agaTCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACg                                               <  1:3627914/71‑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:3472425/1‑71 (MQ=255)
                       ccGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCg                                         <  1:97642/71‑1 (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:1970514/1‑46 (MQ=255)
                                           acaATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGc                                              >  1:3708728/1‑46 (MQ=255)
                                           acaATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTgaaa                       >  1:1215307/1‑69 (MQ=255)
                                           acaATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTgaaa                       >  1:1841632/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)
                                                   aCGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAACaa              <  1:1652502/70‑1 (MQ=255)
                                                     gcgTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTga                         <  1:2203942/57‑1 (MQ=255)
                                                        tttCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTgaa                        <  1:3449934/55‑1 (MQ=255)
                                                        tttCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTgaa                        <  1:267610/55‑1 (MQ=255)
                                                        tttCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTcc        <  1:3392173/71‑1 (MQ=255)
                                                          tCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCAt      >  1:187068/1‑71 (MQ=255)
                                                          tCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCAt      >  1:3338887/1‑71 (MQ=255)
                                                           cTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATaa              >  1:254352/1‑62 (MQ=255)
                                                              gCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATc     <  1:1033328/68‑1 (MQ=255)
                                                               cGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCgg  <  1:1114614/70‑1 (MQ=255)
                                                                 gCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATaa              <  1:717318/56‑1 (MQ=255)
                                                                 gCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATaa              <  1:2147326/56‑1 (MQ=255)
                                                                  |                                                                  
GATCGACCAGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGG  >  W3110S.gb/2982049‑2982181

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: