Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 268,179 G→A 33.6% T156T (ACC→ACT dxs ← 1‑deoxyxylulose‑5‑phosphate synthase, thiamine‑requiring, FAD‑requiring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE268,1790GA33.6% 70.7 / 46.7 89T156T (ACC→ACTdxs1‑deoxyxylulose‑5‑phosphate synthase, thiamine‑requiring, FAD‑requiring
Reads supporting (aligned to +/- strand):  ref base G (31/28);  new base A (13/17);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

GAATCACCAGCATATCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCA  >  minE/268115‑268242
                                                                |                                                               
gaATCACCAGCATATCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTa                                                              <  1:1065007/68‑1 (MQ=255)
gaATCACCAGCATATCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTa                                                              <  1:1007956/68‑1 (MQ=255)
gaATCACCAGCATATCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTa                                                              <  1:2387181/68‑1 (MQ=255)
gaATCACCAGCATATCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTa                                                              <  1:2648224/68‑1 (MQ=255)
gaATCACCAGCATATCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTa                                                              <  1:200133/68‑1 (MQ=255)
gaATCACCAGCATATCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTa                                                              <  1:1181711/68‑1 (MQ=255)
 aaTCACCAGCATATCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAAt                                                            <  1:2078168/69‑1 (MQ=255)
 aaTCACCAGCATATCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAAt                                                            <  1:1322261/69‑1 (MQ=255)
 aaTCACCAGCATATCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAAt                                                            <  1:696622/69‑1 (MQ=255)
 aaTCACCAGCATATCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAAt                                                            <  1:1945556/69‑1 (MQ=255)
 aaTCACCAGCATATCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAAt                                                            <  1:3369225/69‑1 (MQ=255)
 aaTCACCAGCATATCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAAt                                                            <  1:3030488/69‑1 (MQ=255)
 aaTCACCAGCATATCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGAGGTAAt                                                            <  1:3114/69‑1 (MQ=255)
 aaTCACCAGCATATCAGGACGGATATCGCCCGCGTGATTCATAGCTTCAAACGCCATGCCTGCGGTAAt                                                            <  1:2394490/69‑1 (MQ=255)
   tCACCAGCATATCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATcg                                                          >  1:1209742/1‑69 (MQ=255)
   tCACCAGCATATCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATcg                                                          >  1:2696494/1‑69 (MQ=255)
   tCACCAGCATATCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATcg                                                          >  1:680847/1‑69 (MQ=255)
   tCACCAGCATATCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATcg                                                          >  1:2880488/1‑69 (MQ=255)
   tCACAAGCATATCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATcg                                                          >  1:2562591/1‑69 (MQ=255)
           atatCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATc                                                   >  1:2935466/1‑68 (MQ=255)
             atCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCAGCAGTAATCGCGCCATCGcc                                                <  1:2777123/69‑1 (MQ=255)
             atCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCAGCAGTAATCGCGCCATCGcc                                                <  1:2144143/69‑1 (MQ=255)
             atCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCAGCAGTAATCGCGCCATCGcc                                                <  1:2658668/69‑1 (MQ=255)
             atCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCAGCAGTAATCGCGCCATCGcc                                                <  1:1294855/69‑1 (MQ=255)
             atCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCAGCAGTAATCGCGCCATCGcc                                                <  1:1455683/69‑1 (MQ=255)
             atCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCAGCAGTAATCGCGCCATCGcc                                                <  1:1971506/69‑1 (MQ=255)
             atCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCAGCAGTAATCGCGCCATCGcc                                                <  1:3119972/69‑1 (MQ=255)
             atCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCAGCAGTAATCGCGCCATCGcc                                                <  1:339708/69‑1 (MQ=255)
             atCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCAGCAGTAATCGCGCCATCGcc                                                <  1:1500059/69‑1 (MQ=255)
             atCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCAGCAGTAATCGCGCCATCGcc                                                <  1:1931147/69‑1 (MQ=255)
             atCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCAGCAGTAATCGCGCCATCGcc                                                <  1:2230558/69‑1 (MQ=255)
             atCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCAGCAGTAATCGCGCCATCGcc                                                <  1:2269778/69‑1 (MQ=255)
             atCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCAGCAGTAATCGCGCCATCGcc                                                <  1:157472/69‑1 (MQ=255)
             atCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCAGCAGTAATCGCGCCATCGcc                                                <  1:2558311/69‑1 (MQ=255)
              tCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCa                                               >  1:411637/1‑69 (MQ=255)
              tCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCa                                               >  1:1658781/1‑69 (MQ=255)
              tCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCAGCAGTAATCGCGCCATCGcc                                                <  1:2214173/68‑1 (MQ=255)
              tCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCAGCAGTAATCGCGCAATCGcc                                                <  1:5835/68‑1 (MQ=255)
              acAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTCCGGTAATCGCGCCATCGcc                                                >  1:1579847/2‑68 (MQ=255)
                  gACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCAGCAGTAATCGCGCCATCGcc                                                <  1:1857179/64‑1 (MQ=255)
                      gATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGAcaca                                       <  1:2718183/69‑1 (MQ=255)
                      gATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGAcaca                                       <  1:2961441/69‑1 (MQ=255)
                      gATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGAcaca                                       <  1:491224/69‑1 (MQ=255)
                      gATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGAcaca                                       <  1:563405/69‑1 (MQ=255)
                      gATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGAcaca                                       <  1:2389247/69‑1 (MQ=255)
                      gATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGAcaca                                       <  1:2382045/69‑1 (MQ=255)
                      gATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGAcaca                                       <  1:1757211/69‑1 (MQ=255)
                      gATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGAcaca                                       <  1:1370075/69‑1 (MQ=255)
                      gATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGAcaca                                       <  1:2015426/69‑1 (MQ=255)
                      gATATCGCCAGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGAcaca                                       <  1:3159115/69‑1 (MQ=255)
                       atatCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGAcaca                                       <  1:3029656/68‑1 (MQ=255)
                        tatCGCCCTCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATTAcaca                                       <  1:939121/67‑1 (MQ=255)
                        tatCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGACACAGa                                     >  1:2997104/1‑69 (MQ=255)
                        tatCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGACACAGa                                     >  1:1951564/1‑69 (MQ=255)
                        tatCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGACACAGa                                     >  1:2142824/1‑69 (MQ=255)
                        tatCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGACACAGa                                     >  1:1194318/1‑69 (MQ=255)
                        tatCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGACACAGa                                     >  1:1155800/1‑69 (MQ=255)
                        tatCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGACACAGa                                     >  1:728919/1‑69 (MQ=255)
                            gCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGACACAGACgg                                  <  1:259835/68‑1 (MQ=255)
                                       tCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGcc                                                <  1:2187294/43‑1 (MQ=255)
                                       tCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGAtttt                      >  1:1284006/1‑69 (MQ=255)
                                       tCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGAtttt                      >  1:1383728/1‑69 (MQ=255)
                                       tCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGAtttt                      >  1:1572015/1‑69 (MQ=255)
                                       tCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGAtttt                      >  1:3147898/1‑69 (MQ=255)
                                       tCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGAtttt                      >  1:960921/1‑69 (MQ=255)
                                       tCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGAtttt                      >  1:2541784/1‑69 (MQ=255)
                                       tCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGAtttt                      >  1:1894896/1‑69 (MQ=255)
                                       tCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGAtttt                      >  1:1936789/1‑69 (MQ=255)
                                       tCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGAtttt                      >  1:2373594/1‑69 (MQ=255)
                                                      ccATGCCAGCAGTAATCGCGCCATCGCCAATGACACAGAc                                    >  1:1497113/1‑40 (MQ=38)
                                                      ccATGCCAGCAGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCtttct          >  1:2551747/1‑64 (MQ=255)
                                                      ccATGCCAGCAGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTtctc         >  1:2981963/1‑67 (MQ=255)
                                                      ccATGCCAGCAGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTtctc         >  1:861650/1‑67 (MQ=255)
                                                      ccATGCCAGCAGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTtctc         >  1:2568986/1‑67 (MQ=255)
                                                      ccATGCCAGCAGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTtctc         >  1:1428984/1‑67 (MQ=255)
                                                      ccATGCCAGCAGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTtctc         >  1:3412708/1‑67 (MQ=255)
                                                      ccATGCCAGCAGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTtctc         >  1:2829357/1‑67 (MQ=255)
                                                      ccATGCCAGCAGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTtctc         >  1:286157/1‑67 (MQ=255)
                                                      ccATGCCAGCAGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTtctc         >  1:1688146/1‑67 (MQ=255)
                                                      ccATGCCAGCAGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTtctc         >  1:2928912/1‑67 (MQ=255)
                                                      ccATGCCAGCAGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTtctc         >  1:1852844/1‑67 (MQ=255)
                                                      ccATGCCAGCAGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTAGCCTTCtttct          >  1:2133944/1‑64 (MQ=255)
                                                           ccTGCGGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGAttttt                     >  1:2549868/1‑50 (MQ=255)
                                                           ccTGCGGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGgcagca  >  1:2866782/1‑69 (MQ=255)
                                                           ccTGCGGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGgcagca  >  1:2753088/1‑69 (MQ=255)
                                                           ccTGCGGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGgcagca  >  1:1994897/1‑69 (MQ=255)
                                                           ccTGCGGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGgcagca  >  1:2613191/1‑69 (MQ=255)
                                                           ccTGCGGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGgcagca  >  1:1639507/1‑69 (MQ=255)
                                                           ccTGCGGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGgcagca  >  1:85389/1‑69 (MQ=255)
                                                                |                                                               
GAATCACCAGCATATCAGGACGGATATCGCCCGCGTGATTCATCGCTTCAAACGCCATGCCTGCGGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCA  >  minE/268115‑268242

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: