Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 551,567 C→T 30.6% V653I (GTT→ATT)  ybiW ← predicted pyruvate formate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE551,5670CT30.6% 114.1 / 55.1 117V653I (GTT→ATT) ybiWpredicted pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base C (36/45);  new base T (18/18);  total (54/63)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 117 total aligned reads displayed.

GTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGC  >  minE/551503‑551630
                                                                |                                                               
gTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCgg                                                              >  1:2676591/1‑68 (MQ=255)
gTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCgg                                                              >  1:900834/1‑68 (MQ=255)
gTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCgg                                                              >  1:2883100/1‑68 (MQ=255)
gTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCgg                                                              >  1:437097/1‑68 (MQ=255)
gTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCgg                                                              >  1:3420962/1‑68 (MQ=255)
gTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCgg                                                              >  1:3305839/1‑68 (MQ=255)
   tGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGAcc                                                           >  1:2364032/1‑68 (MQ=255)
   tGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGAcc                                                           >  1:2446158/1‑68 (MQ=255)
   tGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGAcc                                                           >  1:261067/1‑68 (MQ=255)
   tGCGCGCCAAACGGTACGTTAGCGGAGATTGAAGACGTACCCGCGTAATAGTTGCCGCCAACCGGAc                                                            >  1:3254541/1‑67 (MQ=255)
       cgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGAc                                                      <  1:1665900/69‑1 (MQ=255)
       cgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGAc                                                      <  1:1917496/69‑1 (MQ=255)
       cgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGAc                                                      <  1:20953/69‑1 (MQ=255)
       cgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGAc                                                      <  1:2785347/69‑1 (MQ=255)
       cgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGAc                                                      <  1:3354407/69‑1 (MQ=255)
             aCGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAgc                                                >  1:1621764/1‑69 (MQ=255)
             aCGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAgc                                                >  1:428656/1‑69 (MQ=255)
             aCGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACAGTAgc                                                >  1:2149288/1‑69 (MQ=255)
                    gTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGAc                                                      <  1:1140016/56‑1 (MQ=255)
                     ttAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTa                                        >  1:856210/1‑69 (MQ=255)
                     ttAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTa                                        >  1:55803/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:1754245/69‑1 (MQ=255)
                                         aCCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTtcgtc                    <  1:2879422/69‑1 (MQ=255)
                                            cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATg                                      >  1:2817336/1‑48 (MQ=255)
                                            cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                  <  1:3251052/68‑1 (MQ=255)
                                            cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                  >  1:2630074/1‑68 (MQ=255)
                                            cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                  <  1:308712/68‑1 (MQ=255)
                                            cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                  <  1:308617/68‑1 (MQ=255)
                                            cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                  <  1:2274804/68‑1 (MQ=255)
                                            cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                  <  1:2727031/68‑1 (MQ=255)
                                            cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                  <  1:2968801/68‑1 (MQ=255)
                                            cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                  <  1:2756524/68‑1 (MQ=255)
                                            cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                  >  1:1040270/1‑68 (MQ=255)
                                            cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                  >  1:468607/1‑68 (MQ=255)
                                            cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                  <  1:889548/68‑1 (MQ=255)
                                            cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                  <  1:877861/68‑1 (MQ=255)
                                            cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                  >  1:746327/1‑68 (MQ=255)
                                            cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                  >  1:1263146/1‑68 (MQ=255)
                                            cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                  >  1:1380783/1‑68 (MQ=255)
                                            cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                  >  1:1381863/1‑68 (MQ=255)
                                            cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                  <  1:1447333/68‑1 (MQ=255)
                                            cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                  <  1:632148/68‑1 (MQ=255)
                                            cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                  >  1:1823702/1‑68 (MQ=255)
                                             gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGatt                <  1:1618396/69‑2 (MQ=255)
                                             gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                <  1:1612113/69‑1 (MQ=255)
                                             gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                <  1:2647858/69‑1 (MQ=255)
                                             gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                <  1:2675630/69‑1 (MQ=255)
                                             gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                <  1:1433513/69‑1 (MQ=255)
                                             gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                <  1:1667197/69‑1 (MQ=255)
                                             gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                <  1:1811461/69‑1 (MQ=255)
                                             gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                <  1:656589/69‑1 (MQ=255)
                                             gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                <  1:239740/69‑1 (MQ=255)
                                             gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                <  1:3275337/69‑1 (MQ=255)
                                             gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                <  1:1139952/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:962830/1‑43 (MQ=255)
                                                  aTAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATgg                                     >  1:2868709/1‑43 (MQ=255)
                                                  aTAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATgg                                     >  1:1591219/1‑43 (MQ=255)
                                                  aTAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATgg                                     >  1:1358796/1‑43 (MQ=255)
                                                  aTAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATgg                                     >  1:2248421/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:3133001/69‑1 (MQ=255)
                                                     gTTGCCGCCAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTg        <  1:788999/69‑1 (MQ=255)
                                                     gTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCt         >  1:1882143/1‑68 (MQ=255)
                                                     gTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCt         >  1:3076916/1‑68 (MQ=255)
                                                           gccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc  <  1:1190083/69‑1 (MQ=255)
                                                           gccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc  <  1:1154880/69‑1 (MQ=255)
                                                           gccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc  <  1:2940394/69‑1 (MQ=255)
                                                           gccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc  <  1:1358155/69‑1 (MQ=255)
                                                           gccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc  <  1:3401075/69‑1 (MQ=255)
                                                           gccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc  <  1:2627800/69‑1 (MQ=255)
                                                           gccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc  <  1:3156883/69‑1 (MQ=255)
                                                            ccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc  <  1:2086929/68‑1 (MQ=255)
                                                            ccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc  <  1:1243987/68‑1 (MQ=255)
                                                            ccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc  <  1:1859989/68‑1 (MQ=255)
                                                             cAATCGGACCACGCCCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc  >  1:373391/1‑67 (MQ=255)
                                                             cAATCGGACCACGCCCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc  >  1:639406/1‑67 (MQ=255)
                                                             cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc  >  1:1256445/1‑67 (MQ=255)
                                                             cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc  >  1:1107120/1‑67 (MQ=255)
                                                             cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc  >  1:820418/1‑67 (MQ=255)
                                                             cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc  >  1:1972055/1‑67 (MQ=255)
                                                             cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc  >  1:2095337/1‑67 (MQ=255)
                                                             cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc  >  1:2328452/1‑67 (MQ=255)
                                                             cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc  >  1:3182747/1‑67 (MQ=255)
                                                             cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc  >  1:2659977/1‑67 (MQ=255)
                                                             cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc  >  1:2678056/1‑67 (MQ=255)
                                                             cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc  >  1:2692347/1‑67 (MQ=255)
                                                             cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc  >  1:2751404/1‑67 (MQ=255)
                                                             cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc  >  1:2787839/1‑67 (MQ=255)
                                                             cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc  >  1:2795191/1‑67 (MQ=255)
                                                             cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAg   >  1:2354645/1‑66 (MQ=255)
                                                                |                                                               
GTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGC  >  minE/551503‑551630

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: