Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 551,601 T→C 37.0% K641K (AAA→AAG ybiW ← predicted pyruvate formate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE551,6010TC37.0% 60.1 / 69.1 108K641K (AAA→AAGybiWpredicted pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base T (15/53);  new base C (18/22);  total (33/75)
Fisher's exact test for biased strand distribution p-value = 1.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

GACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTA  >  minE/551539‑551646
                                                              |                                             
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)
        cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:2756524/68‑1 (MQ=255)
        cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  >  1:1823702/1‑68 (MQ=255)
        cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:2727031/68‑1 (MQ=255)
        cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:3251052/68‑1 (MQ=255)
        cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:1447333/68‑1 (MQ=255)
        cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:632148/68‑1 (MQ=255)
        cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  >  1:746327/1‑68 (MQ=255)
        cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  >  1:1381863/1‑68 (MQ=255)
        cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  >  1:1380783/1‑68 (MQ=255)
        cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  >  1:468607/1‑68 (MQ=255)
        cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:308712/68‑1 (MQ=255)
        cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:2274804/68‑1 (MQ=255)
        cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:308617/68‑1 (MQ=255)
        cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  >  1:2630074/1‑68 (MQ=255)
        cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:877861/68‑1 (MQ=255)
        cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:2968801/68‑1 (MQ=255)
        cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:889548/68‑1 (MQ=255)
        cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:241519/68‑1 (MQ=255)
        cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  >  1:1040270/1‑68 (MQ=255)
         gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:1520178/69‑1 (MQ=255)
         gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:3275337/69‑1 (MQ=255)
         gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:2647858/69‑1 (MQ=255)
         gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:239740/69‑1 (MQ=255)
         gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:2675630/69‑1 (MQ=255)
         gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:656589/69‑1 (MQ=255)
         gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:1667197/69‑1 (MQ=255)
         gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:1612113/69‑1 (MQ=255)
         gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:1811461/69‑1 (MQ=255)
         gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:1139952/69‑1 (MQ=255)
         gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:1433513/69‑1 (MQ=255)
         gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:384010/67‑1 (MQ=255)
          cgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGatat                               <  1:852247/69‑1 (MQ=255)
           gTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:3172423/65‑1 (MQ=255)
             aaTAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:2842351/63‑1 (MQ=255)
              aTAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:2505662/64‑1 (MQ=255)
                aGTTCCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:2034816/60‑1 (MQ=255)
                 gTTGCCGCCAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTg                        <  1:3133001/69‑1 (MQ=255)
                 gTTGCCGCCAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTg                        <  1:1240856/69‑1 (MQ=255)
                 gTTGCCGCCAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTg                        <  1:1218173/69‑1 (MQ=255)
                 gTTGCCGCCAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTg                        <  1:788999/69‑1 (MQ=255)
                 gTTGCCGCCAATCGGACCACGACCGTAGCGCGCGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTg                        <  1:425944/69‑1 (MQ=255)
                 gTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCt                         >  1:3076916/1‑68 (MQ=255)
                 gTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCt                         >  1:1882143/1‑68 (MQ=255)
                   tGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                  <  1:644013/57‑1 (MQ=255)
                     ccgccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                <  1:1263993/57‑1 (MQ=255)
                       gccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:3401075/69‑1 (MQ=255)
                       gccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:1154880/69‑1 (MQ=255)
                       gccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:2940394/69‑1 (MQ=255)
                       gccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:1171927/69‑1 (MQ=255)
                       gccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:1190083/69‑1 (MQ=255)
                       gccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:1358155/69‑1 (MQ=255)
                       gccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:3156883/69‑1 (MQ=255)
                       gccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:2682453/69‑1 (MQ=255)
                       gccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:2627800/69‑1 (MQ=255)
                        ccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:1859989/68‑1 (MQ=255)
                        ccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:1919200/68‑1 (MQ=255)
                        ccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:1243987/68‑1 (MQ=255)
                        ccAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  <  1:2086929/68‑1 (MQ=255)
                         cAATCGGACCACGCCCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  >  1:639406/1‑67 (MQ=255)
                         cAATCGGACCACGCCCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  >  1:373391/1‑67 (MQ=255)
                         cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGat                                 >  1:2926382/1‑52 (MQ=255)
                         cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  >  1:820418/1‑67 (MQ=255)
                         cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  >  1:1107120/1‑67 (MQ=255)
                         cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  >  1:2795191/1‑67 (MQ=255)
                         cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  >  1:1256445/1‑67 (MQ=255)
                         cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  >  1:1972055/1‑67 (MQ=255)
                         cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  >  1:2189440/1‑67 (MQ=255)
                         cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  >  1:2328452/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:2787839/1‑67 (MQ=255)
                         cAATCGGACCACGACCGTAGCGCGGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTGATAAgc                  >  1:2751404/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:627229/69‑1 (MQ=255)
                                      aCCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGt   <  1:489835/69‑1 (MQ=255)
                                      aCCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGt   <  1:719720/69‑1 (MQ=255)
                                      aCCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGt   <  1:2728597/69‑1 (MQ=255)
                                       ccGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTa  <  1:1191192/69‑1 (MQ=255)
                                       ccGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTa  <  1:905594/69‑1 (MQ=255)
                                       ccGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTa  <  1:1405549/69‑1 (MQ=255)
                                       ccGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTa  <  1:780962/69‑1 (MQ=255)
                                       ccGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTa  <  1:2242289/69‑1 (MQ=255)
                                       ccGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTATGATAAGCGCGAGCCAGCAGCGTa  <  1:785086/69‑1 (MQ=255)
                                               cgGGTTATGGTACTGCTTCAGTTCGTCGATATAGGTCTg                        <  1:2059990/39‑1 (MQ=38)
                                                              |                                             
GACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTA  >  minE/551539‑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: