Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 552,675 T→C 42.6% S283S (TCA→TCG ybiW ← predicted pyruvate formate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE552,6750TC42.6% 31.4 / 83.6 108S283S (TCA→TCGybiWpredicted pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base T (40/22);  new base C (16/30);  total (56/52)
Fisher's exact test for biased strand distribution p-value = 3.34e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

CAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGC  >  minE/552620‑552729
                                                       |                                                      
cAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGt                                                    >  1:2256545/1‑60 (MQ=255)
cAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGAc                                                 >  1:2423510/1‑63 (MQ=255)
     tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg                                      <  1:268360/66‑1 (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:3362800/66‑1 (MQ=255)
     tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg                                      <  1:3044464/66‑1 (MQ=255)
     tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg                                      <  1:3029292/66‑1 (MQ=255)
     tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg                                      <  1:2955335/66‑1 (MQ=255)
     tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg                                      <  1:2676572/66‑1 (MQ=255)
     tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg                                      <  1:281923/66‑1 (MQ=255)
     tcacggcgATAGTATGGATAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg                                      <  1:2578753/66‑1 (MQ=255)
      cgcgGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGa                                     <  1:2949840/69‑1 (MQ=255)
      cgcgGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGa                                     <  1:2306229/69‑1 (MQ=255)
      cgcgGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGa                                     <  1:406181/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)
                  aCGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCa                         <  1:950838/69‑1 (MQ=255)
                  aCGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCa                         <  1:1888020/69‑1 (MQ=255)
                     gatAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCg                                      >  1:2208319/3‑53 (MQ=255)
                     gatAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAAc                       >  1:2513221/3‑68 (MQ=255)
                     gatAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAAc                       >  1:3219208/3‑68 (MQ=255)
                     gatAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAAc                       >  1:3012177/3‑68 (MQ=255)
                     gatAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAAc                       >  1:296870/3‑68 (MQ=255)
                     gatAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAAc                       >  1:2420731/3‑68 (MQ=255)
                     gatAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAAc                       >  1:958896/3‑68 (MQ=255)
                     gatAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAAc                       >  1:1155044/3‑68 (MQ=255)
                     gatAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAAc                       >  1:1804718/3‑68 (MQ=255)
                     gatAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAAc                       >  1:1745328/3‑68 (MQ=255)
                     gatAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAAc                       >  1:3431826/3‑68 (MQ=255)
                     gatAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAAc                       >  1:389559/3‑68 (MQ=255)
                     gatAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAAc                       >  1:1582924/3‑68 (MQ=255)
                     gatAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAAc                       >  1:1390394/3‑68 (MQ=255)
                     gatAGAGATACTGGTCCATACGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAAc                       >  1:561734/3‑68 (MQ=255)
                       tAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGa                                     <  1:2958023/52‑1 (MQ=255)
                              aCTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGt             >  1:268799/1‑69 (MQ=255)
                              aCTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGt             >  1:596647/1‑69 (MQ=255)
                              aCTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGt             >  1:1623187/1‑69 (MQ=255)
                              aCTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGt             >  1:3430254/1‑69 (MQ=255)
                                    ccATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAacac        <  1:1479989/68‑1 (MQ=255)
                                    ccATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAacac        <  1:126678/68‑1 (MQ=255)
                                    ccATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAacac        <  1:1243625/68‑1 (MQ=255)
                                    ccATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAacac        <  1:2350015/68‑1 (MQ=255)
                                    ccATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAacac        <  1:2774433/68‑1 (MQ=255)
                                    ccATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAacac        <  1:1155881/68‑1 (MQ=255)
                                        aCGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACATAGTTg   <  1:814707/69‑1 (MQ=255)
                                        aCGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACATAGTTg   <  1:3372594/69‑1 (MQ=255)
                                        aCGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACATAGTTg   <  1:1139838/69‑1 (MQ=255)
                                        aCGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACATAGTTg   <  1:2737876/69‑1 (MQ=255)
                                        aCGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACATAGTTg   <  1:873739/69‑1 (MQ=255)
                                        aCGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACATAGTTg   <  1:2853311/69‑1 (MQ=255)
                                        aCGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACATAGTTg   <  1:64477/69‑1 (MQ=255)
                                        aCGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACATAGTTg   <  1:426258/69‑1 (MQ=255)
                                        aCGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACATAGTTg   <  1:2958077/69‑1 (MQ=255)
                                        aCGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACATAGTTg   <  1:262152/69‑1 (MQ=255)
                                        aCGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACATAGTTg   <  1:20421/69‑1 (MQ=255)
                                        aCGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACATAGTTg   <  1:2528551/69‑1 (MQ=255)
                                        aCGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACATAGTTg   <  1:1260063/69‑1 (MQ=255)
                                        aCGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACATAGTTg   <  1:2211186/69‑1 (MQ=255)
                                         cGACCACACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:561/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGc                                >  1:3250831/1‑39 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGGAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:402252/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACt                      >  1:1399323/1‑49 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATTAAGTAACACAGTTGc  >  1:1141903/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:1639856/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:1999193/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  <  1:1603225/69‑1 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:1551017/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  <  1:152304/69‑1 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:475243/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:1431530/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:805144/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:1225227/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:868577/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:1191110/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:1127211/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:976847/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:987803/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:1867285/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  <  1:2747210/69‑1 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:2634458/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:2548227/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:2535605/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:2475490/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:2438399/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:238163/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:2324609/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  <  1:2972266/69‑1 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  <  1:3029120/69‑1 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:2205789/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:3039106/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:2203189/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  <  1:2025790/69‑1 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:2716244/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:1830233/1‑69 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGTTACATTCGATCTGCAAAATCAACTGGATCaa               >  1:2887397/1‑56 (MQ=255)
                                         cGACCAAACGATACTGAGTGACCGATAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGc  >  1:2645767/1‑69 (MQ=255)
                                         cGACCAAACGATACCGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACATAGTTg   <  1:44431/68‑1 (MQ=255)
                                                       |                                                      
CAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGC  >  minE/552620‑552729

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: