Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 552,642 G→A 64.5% Y294Y (TAC→TAT ybiW ← predicted pyruvate formate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE552,6420GA64.5% 47.9 / 44.3 76Y294Y (TAC→TATybiWpredicted pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base G (7/20);  new base A (24/25);  total (31/45)
Fisher's exact test for biased strand distribution p-value = 5.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

CAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCA  >  minE/552576‑552706
                                                                  |                                                                
cagcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACggg                                                                  >  1:3043837/1‑67 (MQ=255)
cagcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACggg                                                                  >  1:1884032/1‑67 (MQ=255)
cagcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACggg                                                                  >  1:3068446/1‑67 (MQ=255)
cagcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACggg                                                                  >  1:360069/1‑67 (MQ=255)
cagcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACggg                                                                  >  1:2919035/1‑67 (MQ=255)
cagcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACggg                                                                  >  1:1858812/1‑67 (MQ=255)
cagcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACggg                                                                  >  1:2772418/1‑67 (MQ=255)
cagcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACggg                                                                  >  1:796820/1‑67 (MQ=255)
cagcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTa                                                                <  1:2039974/69‑1 (MQ=255)
cagcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTa                                                                <  1:732180/69‑1 (MQ=255)
cagcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTa                                                                <  1:336475/69‑1 (MQ=255)
cagcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTa                                                                <  1:2874359/69‑1 (MQ=255)
 agcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTa                                                                <  1:2187334/68‑1 (MQ=255)
    aTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAgaga                                                            <  1:311646/69‑1 (MQ=255)
    aTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAgaga                                                            <  1:2266969/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:2765119/1‑68 (MQ=255)
             gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:3215557/1‑68 (MQ=255)
             gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:2672397/1‑68 (MQ=255)
             gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:3270680/1‑68 (MQ=255)
             gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:285360/1‑68 (MQ=255)
             gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:521025/1‑68 (MQ=255)
             gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:3051598/1‑68 (MQ=255)
             gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:2176948/1‑68 (MQ=255)
             gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:1878207/1‑68 (MQ=255)
             gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:771867/1‑68 (MQ=255)
             gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:1458542/1‑68 (MQ=255)
             gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:1431766/1‑68 (MQ=255)
             gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:753528/1‑68 (MQ=255)
             gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:1412836/1‑68 (MQ=255)
             gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:1404161/1‑68 (MQ=255)
             gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:1306149/1‑68 (MQ=255)
             gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:1859705/1‑68 (MQ=255)
             gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:731751/1‑68 (MQ=255)
             gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:987987/1‑68 (MQ=255)
             gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:111101/1‑68 (MQ=255)
             gCGTGTTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTc                                                    >  1:677928/1‑68 (MQ=255)
                 gtTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGATATACTGGTCCATAc                                               <  1:557936/69‑1 (MQ=255)
                 gtTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGTTACTGGTCCATAc                                               <  1:1059807/69‑1 (MQ=255)
                 gtTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTCCATAc                                               <  1:1761593/69‑1 (MQ=255)
                 gtTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTCCATAc                                               <  1:92417/69‑1 (MQ=255)
                 gtTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTCCATAc                                               <  1:3220422/69‑1 (MQ=255)
                 gtTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTCCATAc                                               <  1:2104811/69‑1 (MQ=255)
                 gtTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTCCATAc                                               <  1:25562/69‑1 (MQ=255)
                 gtTCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTCCATAc                                               <  1:2251906/69‑1 (MQ=255)
                  ttCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTCCATAc                                               <  1:2540282/68‑1 (MQ=255)
                  ttCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTCCATAc                                               <  1:2468129/68‑1 (MQ=255)
                  ttCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTCCATAc                                               <  1:680714/68‑1 (MQ=255)
                  ttCGCGATCCAGTGTCTGGTTCAGTTCAACGTCACGGCGATAGTATGGATAGAGATACTGGTCCATAc                                               <  1:570488/68‑1 (MQ=255)
                  ttCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACg                                              <  1:1737838/69‑1 (MQ=255)
                  ttCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACg                                              <  1:1301486/69‑1 (MQ=255)
                   tCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACg                                              <  1:1080531/68‑1 (MQ=255)
                    cgcgATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGa                                             >  1:1191639/1‑68 (MQ=255)
                    cgcgATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGa                                             >  1:3048753/1‑68 (MQ=255)
                    cgcgATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGa                                             >  1:1570941/1‑68 (MQ=255)
                    cgcgATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGa                                             >  1:1857697/1‑68 (MQ=255)
                    cgcgATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGa                                             >  1:2320135/1‑68 (MQ=255)
                             gCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGAt                                    <  1:825260/68‑1 (MQ=255)
                                            cAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGa                                     <  1:1553478/52‑1 (MQ=255)
                                            cAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGa                                     <  1:3181928/52‑1 (MQ=255)
                                            cAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGt                             >  1:2256545/1‑60 (MQ=255)
                                            cAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGAc                          >  1:2423510/1‑63 (MQ=255)
                                                 tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg               <  1:3044464/66‑1 (MQ=255)
                                                 tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg               <  1:2002041/66‑1 (MQ=255)
                                                 tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg               <  1:2015891/66‑1 (MQ=255)
                                                 tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg               <  1:2578753/66‑1 (MQ=255)
                                                 tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg               <  1:3029292/66‑1 (MQ=255)
                                                 tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg               <  1:2955335/66‑1 (MQ=255)
                                                 tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg               <  1:281923/66‑1 (MQ=255)
                                                 tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg               <  1:268360/66‑1 (MQ=255)
                                                 tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg               <  1:2676572/66‑1 (MQ=255)
                                                 tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg               <  1:3362800/66‑1 (MQ=255)
                                                 tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg               <  1:3369700/66‑1 (MQ=255)
                                                  cgcgGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGa              <  1:2306229/69‑1 (MQ=255)
                                                  cgcgGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGa              <  1:1763854/69‑1 (MQ=255)
                                                  cgcgGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGa              <  1:2949840/69‑1 (MQ=255)
                                                  cgcgGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGa              <  1:406181/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:1888020/69‑1 (MQ=255)
                                                              aCGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCa  <  1:950838/69‑1 (MQ=255)
                                                                  |                                                                
CAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCA  >  minE/552576‑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: