Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 552,639 C→T 58.4% P295P (CCG→CCA ybiW ← predicted pyruvate formate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE552,6390CT58.4% 26.5 / 62.2 84P295P (CCG→CCAybiWpredicted pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base C (15/20);  new base T (24/25);  total (39/45)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCA  >  minE/552571‑552706
                                                                    |                                                                   
cTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTAc                                                                     <  1:3292406/69‑1 (MQ=255)
cTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTAc                                                                     <  1:519068/69‑1 (MQ=255)
cTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTAc                                                                     <  1:1967663/69‑1 (MQ=255)
cTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTAc                                                                     <  1:52555/69‑1 (MQ=255)
cTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTAc                                                                     <  1:2343482/69‑1 (MQ=255)
cTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTAc                                                                     <  1:2384340/69‑1 (MQ=255)
cTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTAc                                                                     <  1:2034262/69‑1 (MQ=255)
cTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTAc                                                                     <  1:65499/69‑1 (MQ=255)
cTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTAc                                                                     <  1:65848/69‑1 (MQ=255)
cTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTAc                                                                     <  1:1667705/69‑1 (MQ=255)
cTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTAc                                                                     <  1:1310733/69‑1 (MQ=255)
cTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTAc                                                                     <  1:3013567/69‑1 (MQ=255)
 tATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTAc                                                                     <  1:104371/68‑1 (MQ=255)
   tGCAGCATCTCGATGGCGTGTTCGCGATCGAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTAc                                                                     <  1:3424823/66‑1 (MQ=255)
     cagcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACggg                                                                  >  1:3068446/1‑67 (MQ=255)
     cagcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACggg                                                                  >  1:2772418/1‑67 (MQ=255)
     cagcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACggg                                                                  >  1:360069/1‑67 (MQ=255)
     cagcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACggg                                                                  >  1:3043837/1‑67 (MQ=255)
     cagcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACggg                                                                  >  1:2919035/1‑67 (MQ=255)
     cagcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACggg                                                                  >  1:1858812/1‑67 (MQ=255)
     cagcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACggg                                                                  >  1:1884032/1‑67 (MQ=255)
     cagcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACggg                                                                  >  1:796820/1‑67 (MQ=255)
     cagcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTa                                                                <  1:2039974/69‑1 (MQ=255)
     cagcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTa                                                                <  1:2874359/69‑1 (MQ=255)
     cagcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTa                                                                <  1:336475/69‑1 (MQ=255)
     cagcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTa                                                                <  1:732180/69‑1 (MQ=255)
      agcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTa                                                                <  1:2187334/68‑1 (MQ=255)
         aTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAgaga                                                            <  1:2266969/69‑1 (MQ=255)
         aTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAgaga                                                            <  1:311646/69‑1 (MQ=255)
                  gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATagag                                                             >  1:1814285/1‑59 (MQ=255)
                  gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGAt                                                           >  1:3182202/1‑61 (MQ=255)
                  gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTgtt                                                     >  1:2258843/1‑65 (MQ=255)
                  gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:3270680/1‑68 (MQ=255)
                  gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:3215557/1‑68 (MQ=255)
                  gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:3051598/1‑68 (MQ=255)
                  gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:2672397/1‑68 (MQ=255)
                  gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:2765119/1‑68 (MQ=255)
                  gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:285360/1‑68 (MQ=255)
                  gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:521025/1‑68 (MQ=255)
                  gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:987987/1‑68 (MQ=255)
                  gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:111101/1‑68 (MQ=255)
                  gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:771867/1‑68 (MQ=255)
                  gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:753528/1‑68 (MQ=255)
                  gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:1306149/1‑68 (MQ=255)
                  gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:1404161/1‑68 (MQ=255)
                  gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:1412836/1‑68 (MQ=255)
                  gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:1431766/1‑68 (MQ=255)
                  gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:1458542/1‑68 (MQ=255)
                  gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:731751/1‑68 (MQ=255)
                  gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:677928/1‑68 (MQ=255)
                  gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:1859705/1‑68 (MQ=255)
                  gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:1878207/1‑68 (MQ=255)
                  gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:2176948/1‑68 (MQ=255)
                      gtTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGATATACTGGTCCATAc                                               <  1:557936/69‑1 (MQ=255)
                      gtTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGTTACTGGTCCATAc                                               <  1:1059807/69‑1 (MQ=255)
                      gtTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTCCATAc                                               <  1:2104811/69‑1 (MQ=255)
                      gtTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTCCATAc                                               <  1:92417/69‑1 (MQ=255)
                      gtTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTCCATAc                                               <  1:2251906/69‑1 (MQ=255)
                      gtTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTCCATAc                                               <  1:1761593/69‑1 (MQ=255)
                      gtTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTCCATAc                                               <  1:3220422/69‑1 (MQ=255)
                      gtTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTCCATAc                                               <  1:25562/69‑1 (MQ=255)
                       ttCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTCCATAc                                               <  1:570488/68‑1 (MQ=255)
                       ttCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTCCATAc                                               <  1:680714/68‑1 (MQ=255)
                       ttCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTCCATAc                                               <  1:2468129/68‑1 (MQ=255)
                       ttCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTCCATAc                                               <  1:2540282/68‑1 (MQ=255)
                       ttCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACg                                              <  1:1301486/69‑1 (MQ=255)
                       ttCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACg                                              <  1:1737838/69‑1 (MQ=255)
                        tCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACg                                              <  1:1080531/68‑1 (MQ=255)
                         cgcgATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGa                                             >  1:3048753/1‑68 (MQ=255)
                         cgcgATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGa                                             >  1:1570941/1‑68 (MQ=255)
                         cgcgATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGa                                             >  1:2320135/1‑68 (MQ=255)
                         cgcgATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGa                                             >  1:1191639/1‑68 (MQ=255)
                         cgcgATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGa                                             >  1:1857697/1‑68 (MQ=255)
                                  gCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGAt                                    <  1:825260/68‑1 (MQ=255)
                                                 cAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGa                                     <  1:3181928/52‑1 (MQ=255)
                                                 cAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGa                                     <  1:1553478/52‑1 (MQ=255)
                                                 cAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGt                             >  1:2256545/1‑60 (MQ=255)
                                                 cAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGAc                          >  1:2423510/1‑63 (MQ=255)
                                                      tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg               <  1:3369700/66‑1 (MQ=255)
                                                      tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg               <  1:2002041/66‑1 (MQ=255)
                                                      tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg               <  1:2015891/66‑1 (MQ=255)
                                                      tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg               <  1:268360/66‑1 (MQ=255)
                                                      tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg               <  1:281923/66‑1 (MQ=255)
                                                      tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg               <  1:3362800/66‑1 (MQ=255)
                                                      tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg               <  1:2955335/66‑1 (MQ=255)
                                                      tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg               <  1:2578753/66‑1 (MQ=255)
                                                      tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg               <  1:3044464/66‑1 (MQ=255)
                                                      tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg               <  1:2676572/66‑1 (MQ=255)
                                                      tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg               <  1:3029292/66‑1 (MQ=255)
                                                       cgcgGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGa              <  1:406181/69‑1 (MQ=255)
                                                       cgcgGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGa              <  1:2949840/69‑1 (MQ=255)
                                                       cgcgGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGa              <  1:2306229/69‑1 (MQ=255)
                                                       cgcgGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGa              <  1:1763854/69‑1 (MQ=255)
                                                       cacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg               <  1:2948374/66‑1 (MQ=255)
                                                        gcggcgATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGa              <  1:1196037/68‑1 (MQ=255)
                                                        acggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg               <  1:2591301/66‑1 (MQ=255)
                                                                   aCGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCa  <  1:950838/69‑1 (MQ=255)
                                                                   aCGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCa  <  1:1888020/69‑1 (MQ=255)
                                                                    |                                                                   
CTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCA  >  minE/552571‑552706

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: