Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,287,830 A→G 41.6% T97T (ACA→ACG argE → acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,287,8300AG41.6% 41.8 / 90.6 120T97T (ACA→ACGargEacetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base A (35/35);  new base G (30/20);  total (65/55)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 120 total aligned reads displayed.

GGGCATACCGATACGGTGCCATTTGATGACGGTCGCTGGACGCGCGATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATGAAAGGCTTTTTTGCGTT  >  minE/2287774‑2287895
                                                        |                                                                 
gggCATACCGATACGGTGCCATTTGATGACGGTCGCTGGACGCGCGATCCGTTTACACTGACGGAGCAt                                                       <  1:2116588/69‑1 (MQ=255)
gggCATACCGATACGGTGCCATTTGATGACGGTCGCTGGACGCGCGATCCGTTTACACTGACGGAGCAt                                                       <  1:3257407/69‑1 (MQ=255)
    aTACCGATACGGTGCCATTTGATGACGGTCGCTGGACGCGCGATCCGTTTACACTGACGGAg                                                          <  1:96267/62‑1 (MQ=255)
       ccGATACGGTGCCATTTGATGACGGTCGCTGGACGCGCGATCCGTTTACACTGACGGAGCATGACGGCa                                                >  1:1939200/1‑69 (MQ=255)
       ccGATACGGTGCCATTTGATGACGGTCGCTGGACGCGCGATCCGTTTACACTGACGGAGCATGACGGCa                                                >  1:389125/1‑69 (MQ=255)
       ccGATACGGTGCCATTTGATGACGGTCGCTGGACGCGCGATCCGTTTACACTGACGGAGCATGACGGCa                                                >  1:3394808/1‑69 (MQ=255)
       ccGATACGGTGCCATTTGATGACGGTCGCTGGACGCGCGATCCGTTTACACTGACGGAGCATGACGGCa                                                >  1:2314431/1‑69 (MQ=255)
       ccGATACGGTGCCATTTGATGACGGTCGCTGGACGCGCGATCCGTTTACACTGACGGAGCATGACGGCa                                                >  1:795270/1‑69 (MQ=255)
       ccGATACGGTGCCATTTGATGACGGTCGCTGGACGCGCGATCCGTTTACACTGACGGAGCATGACGGCa                                                >  1:2657080/1‑69 (MQ=255)
       ccGATACGGTGCCATTTGATGACGGTCGCTGGACGCGCGATCCGTTTACACTGACGGAGCATGACGGCa                                                >  1:770817/1‑69 (MQ=255)
       ccGATACGGTGCCATTTGATGACGGTCGCTGGACGCGCGATCCGTTTACACTGACGGAGCATGACGGCa                                                >  1:290533/1‑69 (MQ=255)
       ccGATACGGTGCCATTTGATGACGGTCGCTGGACGCGCGATCCGTTTACACTGACGGAGCATGACGGCa                                                >  1:2786214/1‑69 (MQ=255)
            aCGGTGCCATTTGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCtt                                           <  1:2995223/69‑1 (MQ=255)
            aCGGTGCCATTTGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCtt                                           <  1:2009935/69‑1 (MQ=255)
            aCGGTGCCATTTGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCtt                                           <  1:967607/69‑1 (MQ=255)
            aCGGTGCCATTTGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCtt                                           <  1:2162786/69‑1 (MQ=255)
            aCGGTGCCATTTGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCtt                                           <  1:271395/69‑1 (MQ=255)
            aCGGTGCCATTTGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCtt                                           <  1:3042362/69‑1 (MQ=255)
            aCGGTGCCATTTGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCtt                                           <  1:3298319/69‑1 (MQ=255)
            aCGGTGCCATTTGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCtt                                           <  1:368109/69‑1 (MQ=255)
            aCGGTGCCATTTGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCtt                                           <  1:559167/69‑1 (MQ=255)
            aCGGTGCCATTTGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCtt                                           <  1:658918/69‑1 (MQ=255)
            aCGGTGCCATTTGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCtt                                           <  1:710851/69‑1 (MQ=255)
            aCGGTGCCATTTGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCtt                                           <  1:1100472/69‑1 (MQ=255)
            aCGGTGCCATTTGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCtt                                           <  1:1210523/69‑1 (MQ=255)
            aCGGTGCCATTTGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCtt                                           <  1:1297888/69‑1 (MQ=255)
            aCGGTGCCATTTGATGACGGACGCTGGACGCGCGATCCATTTACGTTGAAGGAGCATGACGGCAAGCtt                                           <  1:2093092/69‑1 (MQ=255)
              ggTGCCATTTGATGACGGTCGCTGGACGCGCGATCCGTTTACACTGACGGAGCAt                                                       <  1:3288302/55‑1 (MQ=255)
                 gCCATTTGATGACGGTCGCTGGACGCGCGATCCGTTTACACTGACGGAg                                                          >  1:1234201/1‑49 (MQ=255)
                 gCCATTTGATGACGGTCGCTGGACGCGCGATCCGTTTACACTGACGGAg                                                          >  1:2486725/1‑49 (MQ=255)
                 gCCATTTGATGACGGTCGCTGGACGCGCGATCCGTTTACACTGACGGAg                                                          >  1:2776493/1‑49 (MQ=255)
                  ccATTTGATGACGGTCGCTGGACGCGCGATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGc                                     <  1:2464681/69‑1 (MQ=255)
                      ttGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAg                                                          >  1:816020/1‑44 (MQ=255)
                      ttGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTa                                  >  1:2829201/1‑68 (MQ=255)
                      ttGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTa                                  >  1:507753/1‑68 (MQ=255)
                      ttGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTa                                  >  1:499361/1‑68 (MQ=255)
                      ttGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTa                                  >  1:3196043/1‑68 (MQ=255)
                      ttGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTa                                  >  1:507934/1‑68 (MQ=255)
                      ttGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTa                                  >  1:2688668/1‑68 (MQ=255)
                      ttGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTa                                  >  1:533846/1‑68 (MQ=255)
                      ttGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTa                                  >  1:2501953/1‑68 (MQ=255)
                      ttGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTa                                  >  1:1235976/1‑68 (MQ=255)
                      ttGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTa                                  >  1:3319377/1‑68 (MQ=255)
                      ttGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTa                                  >  1:2476565/1‑68 (MQ=255)
                      ttGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTa                                  >  1:1022743/1‑68 (MQ=255)
                      ttGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTa                                  >  1:635958/1‑68 (MQ=255)
                      ttGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTa                                  >  1:1709977/1‑68 (MQ=255)
                      ttGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTa                                  >  1:1049989/1‑68 (MQ=255)
                      ttGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTa                                  >  1:1710545/1‑68 (MQ=255)
                      ttGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTa                                  >  1:1778256/1‑68 (MQ=255)
                      ttGATGACGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTa                                  >  1:2102945/1‑68 (MQ=255)
                       tgatgaCGGTCGCTGGCCGCGCGATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAgg                                >  1:337147/1‑69 (MQ=255)
                       tgatgaCGGTCGCTGGACGCGCGATCGGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAgg                                >  1:1644540/1‑69 (MQ=255)
                       tgatgaCGGTCGCTGGACGCGCGATCCGTTTACACTGACGGAGCATGACGGCAGGCTTTACGGCTTAgg                                >  1:1960267/1‑69 (MQ=255)
                       tgatgaCGGTCGCTGGACGCGCGATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAgg                                >  1:3033238/1‑69 (MQ=255)
                       tgatgaCGGTCGCTGGACGCGCGATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAgg                                >  1:862529/1‑69 (MQ=255)
                       tgatgaCGGTCGCTGGACGCGCGATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAgg                                >  1:1388528/1‑69 (MQ=255)
                        gatgaCGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCtt                                           <  1:1251860/57‑1 (MQ=255)
                           gaCGGTCGCTGGACGCGCGATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCAcc                            <  1:3336398/69‑1 (MQ=255)
                           gaCGGTCGCTGGACGCGCGATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCAcc                            <  1:2224640/69‑1 (MQ=255)
                             cGGACGCTGGCCGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCAccg                           >  1:602110/1‑68 (MQ=255)
                             cGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCAccg                           >  1:1107842/1‑68 (MQ=255)
                             cGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCAccg                           >  1:3077656/1‑68 (MQ=255)
                             cGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCAccg                           >  1:1816432/1‑68 (MQ=255)
                             cGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCAccg                           >  1:1664908/1‑68 (MQ=255)
                             cGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCAccg                           >  1:1271417/1‑68 (MQ=255)
                             cGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCAccg                           >  1:2367339/1‑68 (MQ=255)
                             cGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCAccg                           >  1:2585773/1‑68 (MQ=255)
                             cGGACGCTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCAccg                           >  1:703527/1‑68 (MQ=255)
                                tCGCTGGACGCGCGATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGa                       <  1:1804441/69‑1 (MQ=255)
                                   cTGGACGCGCGATCCATTTACGTTGACGGAGCATGACGGCAAGCtt                                           <  1:1070510/46‑1 (MQ=255)
                                    tGGACGCGCGATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATg                   <  1:1372866/69‑1 (MQ=255)
                                    tGGACGCGCGATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATg                   <  1:773036/69‑1 (MQ=255)
                                    tGGACGCGCGATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATg                   <  1:1583236/69‑1 (MQ=255)
                                    tGGACGCGCGATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATg                   <  1:1104770/69‑1 (MQ=255)
                                          cgcgATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATTAAAGGc             <  1:1450628/69‑1 (MQ=255)
                                          cgcgATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATGAAAGGc             <  1:2889732/69‑1 (MQ=255)
                                          cgcgATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATGAAAGGc             <  1:1239639/69‑1 (MQ=255)
                                          cgcgATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATGAAAGGc             <  1:2074862/69‑1 (MQ=255)
                                          cgcgATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATGAAAGGc             <  1:232092/69‑1 (MQ=255)
                                          cgcgATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATGAAAGGc             <  1:2765179/69‑1 (MQ=255)
                                          cgcgATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATGAAAGGc             <  1:2977992/69‑1 (MQ=255)
                                           gcgATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATGAAAGGc             <  1:2677160/68‑1 (MQ=255)
                                           gcgATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATGAAAGGc             <  1:2673158/68‑1 (MQ=255)
                                           gcgATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATGAAAGGc             <  1:2276782/68‑1 (MQ=255)
                                           gcgATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATGAAAGGc             <  1:932949/68‑1 (MQ=255)
                                           gcgATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATGAAAGGc             <  1:900709/68‑1 (MQ=255)
                                           gcgATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATGAAAGGCt            <  1:3121518/69‑1 (MQ=255)
                                           gcgATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATGAAAGGCt            <  1:1983589/69‑1 (MQ=255)
                                           gcgATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATGAAAGGCt            <  1:1992074/69‑1 (MQ=255)
                                           gcgATCCGTTTACACTGACGGAGCATGACGGCAAGCTGTACGGCTTAGGCACCGCCGa                       <  1:1088301/58‑1 (MQ=255)
                                             gATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATGAAAGGCttt          <  1:1510214/69‑1 (MQ=255)
                                             gATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATGAAAGGCttt          <  1:1622201/69‑1 (MQ=255)
                                                     tACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATGAAAGGCttt          >  1:2675911/1‑61 (MQ=255)
                                                     tACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATGAAAGGCTTTTTTGCgct  >  1:2801907/1‑67 (MQ=255)
                                                     tACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATGAAAGGCTTTTTTGCGtt  >  1:1685836/1‑69 (MQ=255)
                                                     tACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATGAAAGGCTTTTTTGCGtt  >  1:1059002/1‑69 (MQ=255)
                                                     tACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATGAAAGGCTTTTTTGCGtt  >  1:746843/1‑69 (MQ=255)
                                                     tACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATGAAAGGCTTTTTTGCGtt  >  1:1475408/1‑69 (MQ=255)
                                                     tACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATGAAAGGCTTTTTTGCGtt  >  1:380282/1‑69 (MQ=255)
                                                        |                                                                 
GGGCATACCGATACGGTGCCATTTGATGACGGTCGCTGGACGCGCGATCCGTTTACACTGACGGAGCATGACGGCAAGCTTTACGGCTTAGGCACCGCCGACATGAAAGGCTTTTTTGCGTT  >  minE/2287774‑2287895

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: