Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 551,589 A→G 34.1% N645N (AAT→AAC ybiW ← predicted pyruvate formate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE551,5890AG34.1% 84.5 / 67.3 117N645N (AAT→AACybiWpredicted pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base A (24/53);  new base G (18/22);  total (42/75)
Fisher's exact test for biased strand distribution p-value = 1.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 117 total aligned reads displayed.

TTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTA  >  minE/551524‑551646
                                                                 |                                                         
ttAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTa                                                        >  1:55803/1‑69 (MQ=255)
ttAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTa                                                        >  1:856210/1‑69 (MQ=255)
               gaCGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGt                                         <  1:697927/69‑1 (MQ=255)
               gaCGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGt                                         <  1:1130763/69‑1 (MQ=255)
                aCGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGt                                         <  1:801262/68‑1 (MQ=255)
                   tACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTtcg                                      >  1:2664738/1‑68 (MQ=255)
                    aCCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTtcgtc                                    <  1:2879422/69‑1 (MQ=255)
                    aCCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTtcgtc                                    <  1:1754245/69‑1 (MQ=255)
                       cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATg                                                      >  1:2817336/1‑48 (MQ=255)
                       cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  >  1:1040270/1‑68 (MQ=255)
                       cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:1447333/68‑1 (MQ=255)
                       cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:2727031/68‑1 (MQ=255)
                       cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:308712/68‑1 (MQ=255)
                       cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  >  1:2650064/1‑68 (MQ=255)
                       cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:241519/68‑1 (MQ=255)
                       cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:632148/68‑1 (MQ=255)
                       cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  >  1:468607/1‑68 (MQ=255)
                       cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:2274804/68‑1 (MQ=255)
                       cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  >  1:1380783/1‑68 (MQ=255)
                       cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:3251052/68‑1 (MQ=255)
                       cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:889548/68‑1 (MQ=255)
                       cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:2756524/68‑1 (MQ=255)
                       cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:308617/68‑1 (MQ=255)
                       cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:2968801/68‑1 (MQ=255)
                       cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:877861/68‑1 (MQ=255)
                        gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGatt                                <  1:1618396/69‑2 (MQ=255)
                        gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:1811461/69‑1 (MQ=255)
                        gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:1139952/69‑1 (MQ=255)
                        gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:239740/69‑1 (MQ=255)
                        gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:2647858/69‑1 (MQ=255)
                        gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:3275337/69‑1 (MQ=255)
                        gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:656589/69‑1 (MQ=255)
                        gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:2675630/69‑1 (MQ=255)
                        gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:1667197/69‑1 (MQ=255)
                        gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:1520178/69‑1 (MQ=255)
                        gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:1433513/69‑1 (MQ=255)
                        gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:384010/67‑1 (MQ=255)
                          gTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:3172423/65‑1 (MQ=255)
                            aaTAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:2842351/63‑1 (MQ=255)
                             aTAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATgg                                                     >  1:2868709/1‑43 (MQ=255)
                             aTAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATgg                                                     >  1:962830/1‑43 (MQ=255)
                             aTAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATgg                                                     >  1:2248421/1‑43 (MQ=255)
                             aTAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATgg                                                     >  1:1358796/1‑43 (MQ=255)
                             aTAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATgg                                                     >  1:1591219/1‑43 (MQ=255)
                             aTAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATgg                                                     >  1:1726306/1‑43 (MQ=255)
                             aTAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:2505662/64‑1 (MQ=255)
                               aGTTCCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:2034816/60‑1 (MQ=255)
                                gTTGCCGCCAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTg                        <  1:1218173/69‑1 (MQ=255)
                                gTTGCCGCCAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTg                        <  1:1240856/69‑1 (MQ=255)
                                gTTGCCGCCAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTg                        <  1:788999/69‑1 (MQ=255)
                                gTTGCCGCCAATCGGACCACGACCGTAGCGCGCGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTg                        <  1:425944/69‑1 (MQ=255)
                                gTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCt                         >  1:1882143/1‑68 (MQ=255)
                                gTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCt                         >  1:3076916/1‑68 (MQ=255)
                                  tGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:644013/57‑1 (MQ=255)
                                    ccgccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:1263993/57‑1 (MQ=255)
                                      gccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:2627800/69‑1 (MQ=255)
                                      gccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:1358155/69‑1 (MQ=255)
                                      gccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:2940394/69‑1 (MQ=255)
                                      gccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:1190083/69‑1 (MQ=255)
                                      gccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:3156883/69‑1 (MQ=255)
                                      gccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:2682453/69‑1 (MQ=255)
                                      gccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:1154880/69‑1 (MQ=255)
                                       ccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:1919200/68‑1 (MQ=255)
                                       ccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:1859989/68‑1 (MQ=255)
                                       ccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:1243987/68‑1 (MQ=255)
                                       ccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:2086929/68‑1 (MQ=255)
                                        cAATCGGACCACGCCCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  >  1:373391/1‑67 (MQ=255)
                                        cAATCGGACCACGCCCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  >  1:639406/1‑67 (MQ=255)
                                        cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  >  1:3182747/1‑67 (MQ=255)
                                        cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  >  1:820418/1‑67 (MQ=255)
                                        cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  >  1:2659977/1‑67 (MQ=255)
                                        cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  >  1:1107120/1‑67 (MQ=255)
                                        cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  >  1:2795191/1‑67 (MQ=255)
                                        cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  >  1:2751404/1‑67 (MQ=255)
                                        cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  >  1:1256445/1‑67 (MQ=255)
                                        cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  >  1:2692347/1‑67 (MQ=255)
                                        cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  >  1:1972055/1‑67 (MQ=255)
                                        cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  >  1:2095337/1‑67 (MQ=255)
                                        cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  >  1:2189440/1‑67 (MQ=255)
                                        cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAg                   >  1:2354645/1‑66 (MQ=255)
                                            cGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTg                        <  1:3175431/57‑1 (MQ=255)
                                             ggACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGcc           >  1:1891343/1‑69 (MQ=255)
                                             ggACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGcc           >  1:1276452/1‑69 (MQ=255)
                                             ggACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGcc           >  1:2342072/1‑69 (MQ=255)
                                              gACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:1432224/61‑1 (MQ=255)
                                                  aCGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:3062475/57‑1 (MQ=255)
                                                  aCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCcagcag      <  1:76777/69‑1 (MQ=255)
                                                  aCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGAcagcag      <  1:1265130/69‑1 (MQ=255)
                                                   cGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:1768009/42‑1 (MQ=255)
                                                     aCCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCg  <  1:1142452/69‑2 (MQ=255)
                                                     aCCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGt   <  1:719720/69‑1 (MQ=255)
                                                     aCCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGt   <  1:489835/69‑1 (MQ=255)
                                                     aCCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGt   <  1:2728597/69‑1 (MQ=255)
                                                      ccGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTa  <  1:2242289/69‑1 (MQ=255)
                                                      ccGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTa  <  1:780962/69‑1 (MQ=255)
                                                      ccGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTa  <  1:1405549/69‑1 (MQ=255)
                                                      ccGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTa  <  1:1191192/69‑1 (MQ=255)
                                                      ccGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTa  <  1:905594/69‑1 (MQ=255)
                                                      ccGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTATGATAAGCGCGAGCCAGCAGCGTa  <  1:785086/69‑1 (MQ=255)
                                                       cGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGtta                                <  1:2861693/38‑1 (MQ=255)
                                                                 |                                                         
TTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTA  >  minE/551524‑551646

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: