Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 268,233 T→C 52.2% E138E (GAA→GAG 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,2330TC52.2% 3.9 / 58.1 69E138E (GAA→GAGdxs1‑deoxyxylulose‑5‑phosphate synthase, thiamine‑requiring, FAD‑requiring
Reads supporting (aligned to +/- strand):  ref base T (18/15);  new base C (27/9);  total (45/24)
Fisher's exact test for biased strand distribution p-value = 8.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATGCCTGCGGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATAC  >  minE/268169‑268298
                                                                |                                                                 
ccATGCCAGCAGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTtctc                                                                 >  1:2981963/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)
ccATGCCAGCAGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTtctc                                                                 >  1:1428984/1‑67 (MQ=255)
ccATGCCAGCAGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTtctc                                                                 >  1:3412708/1‑67 (MQ=255)
ccATGCCAGCAGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTtctc                                                                 >  1:861650/1‑67 (MQ=255)
ccATGCCAGCAGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTtctc                                                                 >  1:2568986/1‑67 (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:2866782/1‑69 (MQ=255)
     ccTGCGGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGgcagca                                                          >  1:85389/1‑69 (MQ=255)
           gTAATCGCTCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGc                                                     <  1:2716174/68‑1 (MQ=255)
           gTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGc                                                     <  1:2402753/68‑1 (MQ=255)
           gTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGc                                                     <  1:41694/68‑1 (MQ=255)
            tAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCa                                                    <  1:623854/68‑1 (MQ=255)
             aaTCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAAt                                                  >  1:2944143/1‑69 (MQ=255)
                cgcgCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATAcc                                               <  1:934551/69‑1 (MQ=255)
                   gcCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCGGCAGCAACTGCGATCCCAAt                                            <  1:3163126/69‑1 (MQ=255)
                   gcCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCGGCAGCAACTGCGATCCCAAt                                            <  1:564752/69‑1 (MQ=255)
                   gcCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCGGCAGCAACTGCGATCCCAAt                                            <  1:2536042/69‑1 (MQ=255)
                    ccATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCGGCAGCAACTGCGATCCCAAt                                            <  1:2365986/68‑1 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCTATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCgg                                       >  1:204622/1‑68 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCgg                                       >  1:156228/1‑68 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCgg                                       >  1:1747242/1‑68 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCgg                                       <  1:3149840/68‑1 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCgg                                       <  1:1997009/68‑1 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCgg                                       <  1:1003877/68‑1 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCgg                                                               >  1:628101/1‑44 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCgg                                                               >  1:1581571/1‑44 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCgg                                                               >  1:3093009/1‑44 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCgg                                                               >  1:2785388/1‑44 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCgg                                                               >  1:1333985/1‑44 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCgg                                                               >  1:1218990/1‑44 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCgg                                                               >  1:2677087/1‑44 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCgg                                                               >  1:201929/1‑44 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCgg                                                               >  1:2511824/1‑44 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCgg                                                               >  1:1053719/1‑44 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCgg                                                               >  1:711874/1‑44 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCgg                                                               >  1:2367097/1‑44 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCgg                                                               >  1:2099293/1‑44 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCgg                                                               >  1:2340660/1‑44 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCgg                                                               >  1:2106411/1‑44 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCgg                                                               >  1:2247982/1‑44 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCgg                                                               >  1:682532/1‑44 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCgg                                                               >  1:1161782/1‑44 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCgg                                                               >  1:115968/1‑44 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTATTTTTCGGCAGCAACCGCAATACCAATTCCgg                                       <  1:2153394/68‑1 (MQ=255)
                         gCCAATGACACAGACGGTGCGGCGATTTTTGCCTTATTTCTCgg                                                               >  1:1281485/1‑44 (MQ=38)
                         cccAATGACACAGACGGTGCCGCGATTTTTGCCTTCTTTCTCgg                                                               >  1:2247642/2‑44 (MQ=38)
                            aaTGACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACt                                   <  1:26749/69‑1 (MQ=255)
                            aaTGACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACt                                   <  1:1350571/69‑1 (MQ=255)
                            aaTGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCGGCAGCAACTGCGATCCCAAt                                            >  1:2067511/1‑60 (MQ=255)
                            aaTGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCGGCAGCAACTGCGATCCCAATTCCGGCACt                                   >  1:859591/1‑69 (MQ=255)
                            aaTGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCGGCAGCAACTGCGATCCCAATTCCGGCACt                                   >  1:813747/1‑69 (MQ=255)
                            aaTGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCGGCAGCAACTGCGATCCCAATTCCGGCACt                                   >  1:1653720/1‑69 (MQ=255)
                            aaTGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCGGCAGCAACTGCGATCCCAATTCCGGCACt                                   >  1:2913380/1‑69 (MQ=255)
                            aaTGACACAGACGGTGCGGCGATTTTTGCCTTCTTTCTCGGCAGCAACTGCGATCCCAATTCCGGCACt                                   >  1:2348271/1‑69 (MQ=255)
                               gACACAGACGGTGCGGCGATTTTTGCCTTGTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGa                                 >  1:2904175/1‑68 (MQ=255)
                               gACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGa                                 >  1:2089793/1‑68 (MQ=255)
                               gACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGa                                 >  1:2172970/1‑68 (MQ=255)
                               gACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGa                                 >  1:2387713/1‑68 (MQ=255)
                               gACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGa                                 >  1:484887/1‑68 (MQ=255)
                               gACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGa                                 >  1:3330574/1‑68 (MQ=255)
                               gACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGa                                 >  1:3246796/1‑68 (MQ=255)
                               gACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGa                                 >  1:1901467/1‑68 (MQ=255)
                                          tGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTtgatg                     <  1:3367076/69‑1 (MQ=255)
                                          tGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTtgatg                     <  1:3368578/69‑1 (MQ=255)
                                          tGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTtgatg                     <  1:3370984/69‑1 (MQ=255)
                                              gcgATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCgg                                       <  1:2680659/47‑1 (MQ=255)
                                              gcgATTTTTGCCTTCTTTCTCGGCAGCAACTGCGATCCCAATTCCGGCACTGATGGAGGTTGATGAAt                  <  1:126922/68‑1 (MQ=255)
                                              gcgATTTTTGCCTTCTTTCTCGGCAGCAACTGCGATCCCAATTCCGGCACTGATGGAGGTTGATGAAt                  <  1:3085491/68‑1 (MQ=255)
                                              gcgATTTTTGCCTTCTTTCTCGGCAGCAACTGCGATCCCAATTCCGGCACTGATGGAGGTTGATGAAt                  <  1:2679693/68‑1 (MQ=255)
                                              gcgATTTTTGCCTTCTTTCTCGGCAGCAACTGCGATCCCAATTCCGGCACTGATGGAGGTTGATGAAt                  <  1:2473658/68‑1 (MQ=255)
                                              gcgATTTTTGCCTTCTTTCTCGGCAGCAACTGCGATCCCAATTCCGGCACTGATGGAGGTTGATGAAt                  <  1:2827641/68‑1 (MQ=255)
                                                          ttctttctCGGCAGCAACTGCGATCCCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTaa     >  1:1974350/1‑69 (MQ=255)
                                                          ttctttctCGGCAGCAACTGCGATCCCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTaa     >  1:908564/1‑69 (MQ=255)
                                                          ttctttctCGGCAGCAACTGCGATCCCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTaa     >  1:2185265/1‑69 (MQ=255)
                                                          ttctttctCGGCAGCAACTGCGATCCCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTaa     >  1:602324/1‑69 (MQ=255)
                                                          ttctttctCGGCAGCAACTGCGATCCCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTaa     >  1:3156395/1‑69 (MQ=255)
                                                              ttttCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATAc  <  1:740743/68‑1 (MQ=255)
                                                              ttttCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATAc  <  1:2422566/68‑1 (MQ=255)
                                                              ttttCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATAc  <  1:2212631/68‑1 (MQ=255)
                                                              ttttCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATAc  <  1:1253161/68‑1 (MQ=255)
                                                              ttttCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATAc  <  1:1737907/68‑1 (MQ=255)
                                                                |                                                                 
CCATGCCTGCGGTAATCGCGCCATCGCCAATGACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATAC  >  minE/268169‑268298

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