Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,472,750 G→A 8.7% A9T (GCG→ACG)  glmS → L‑glutamine:D‑fructose‑6‑phosphate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,472,7500GA 8.7% 328.3 / 10.5 156A9T (GCG→ACG) glmSL‑glutamine:D‑fructose‑6‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (62/80);  new base A (3/11);  total (65/91)
Fisher's exact test for biased strand distribution p-value = 1.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 159 total aligned reads displayed.

CAGGTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTGGAATACCGCGGA  >  minE/2472687‑2472809
                                                               |                                                           
cAGGTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCaa                                                        >  1:1023883/1‑69 (MQ=255)
cAGGTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCaa                                                        >  1:983167/1‑69 (MQ=255)
cAGGTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCaa                                                        >  1:776071/1‑69 (MQ=255)
cAGGTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCaa                                                        >  1:684663/1‑69 (MQ=255)
cAGGTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCaa                                                        >  1:3292323/1‑69 (MQ=255)
cAGGTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCaa                                                        >  1:3078604/1‑69 (MQ=255)
cAGGTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCaa                                                        >  1:2823269/1‑69 (MQ=255)
cAGGTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCaa                                                        >  1:2813729/1‑69 (MQ=255)
cAGGTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCaa                                                        >  1:2608703/1‑69 (MQ=255)
cAGGTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCaa                                                        >  1:2575130/1‑69 (MQ=255)
cAGGTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCaa                                                        >  1:2451012/1‑69 (MQ=255)
cAGGTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCaa                                                        >  1:2362243/1‑69 (MQ=255)
cAGGTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCaa                                                        >  1:2316941/1‑69 (MQ=255)
cAGGTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCaa                                                        >  1:1801611/1‑69 (MQ=255)
cAGGTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCaa                                                        >  1:1792186/1‑69 (MQ=255)
cAGGTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCaa                                                        >  1:1544704/1‑69 (MQ=255)
cAGGTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCaa                                                        >  1:1088437/1‑69 (MQ=255)
cAGGTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCa                                                         >  1:829970/1‑68 (MQ=255)
cAGGTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCa                                                         >  1:517281/1‑68 (MQ=255)
cAGGTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCa                                                         >  1:427897/1‑68 (MQ=255)
cAGGTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCa                                                         >  1:350377/1‑68 (MQ=255)
cAGGTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATAGTTGGCGCGATCGCGCaa                                                        >  1:2145694/1‑69 (MQ=255)
   gTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGt                                                     <  1:956509/69‑1 (MQ=255)
   gTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGt                                                     <  1:2337005/69‑1 (MQ=255)
   gTTGACCGACAACGATATAAATCGGAATAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTg                                                    <  1:1615095/69‑1 (MQ=255)
    ttGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGt                                                     <  1:1478897/68‑1 (MQ=255)
    ttGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTg                                                    <  1:3254931/69‑1 (MQ=255)
    ttGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTg                                                    <  1:2649380/69‑1 (MQ=255)
    ttGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTg                                                    <  1:2701563/69‑1 (MQ=255)
    ttGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTg                                                    <  1:2058006/69‑1 (MQ=255)
    ttGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTg                                                    <  1:2928724/69‑1 (MQ=255)
    ttGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTg                                                    <  1:3208436/69‑1 (MQ=255)
    ttGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTg                                                    <  1:1848631/69‑1 (MQ=255)
    ttGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTg                                                    <  1:3274227/69‑1 (MQ=255)
    ttGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTg                                                    <  1:3347040/69‑1 (MQ=255)
    ttGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTg                                                    <  1:711097/69‑1 (MQ=255)
    ttGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTg                                                    <  1:1569904/69‑1 (MQ=255)
     tGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTg                                                    <  1:2628570/68‑1 (MQ=255)
     tGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGa                                                   >  1:2929126/1‑69 (MQ=255)
         cGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAATGTGATGt                                                <  1:1344568/68‑1 (MQ=255)
         cGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGt                                                <  1:1022325/68‑1 (MQ=255)
         cGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGt                                                <  1:754503/68‑1 (MQ=255)
         cGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGt                                                <  1:2661054/68‑1 (MQ=255)
         cGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGt                                                <  1:914667/68‑1 (MQ=255)
                 atatAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCACGCAACGTGATGTAGCAGAAAt                                       >  1:1298578/1‑69 (MQ=255)
                 atatAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCACGCAACGTGATGTAGCAGAAAt                                       >  1:3282753/1‑69 (MQ=255)
                 atatAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCACGCAACGTGATGTAGCAGAAAt                                       >  1:2626117/1‑69 (MQ=255)
                   ataAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATcc                                     <  1:288196/69‑1 (MQ=255)
                   ataAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATcc                                     <  1:474362/69‑1 (MQ=255)
                   ataAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATcc                                     <  1:2798318/69‑1 (MQ=255)
                   ataAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATcc                                     <  1:1045716/69‑1 (MQ=255)
                   ataAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATcc                                     <  1:1245466/69‑1 (MQ=255)
                    taAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATcc                                     <  1:1698213/68‑1 (MQ=255)
                                aaaaaCTATGTGTGGAATTGTTGGCGCGATCGCGCa                                                         >  1:1265404/1‑36 (MQ=255)
                                aaaaaCTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGttt                        >  1:2359208/1‑69 (MQ=255)
                                aaaaaCTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGttt                        >  1:1638447/1‑69 (MQ=255)
                                aaaaaCTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGttt                        >  1:1471749/1‑69 (MQ=255)
                                aaaaaCTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGttt                        >  1:1752297/1‑69 (MQ=255)
                                aaaaaCTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGttt                        >  1:2786852/1‑69 (MQ=255)
                                aaaaaCTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGtt                         >  1:310983/1‑68 (MQ=255)
                                aaaaaCTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGtt                         >  1:370278/1‑68 (MQ=255)
                                aaaaaCTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGtt                         <  1:1389252/68‑1 (MQ=255)
                                aaaaaCTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGtt                         <  1:2209209/68‑1 (MQ=255)
                                aaaaaCTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGt                          >  1:3171556/1‑67 (MQ=255)
                                aaaaaCTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCCTCTTGAAGGttt                        >  1:1046636/1‑69 (MQ=255)
                                     cTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTAcgtc                   >  1:848219/1‑69 (MQ=255)
                                     cTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTAcgtc                   >  1:2736400/1‑69 (MQ=255)
                                     cTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTAcgt                    >  1:194932/1‑68 (MQ=255)
                                          tgtgGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACgtcgtc                >  1:71512/1‑67 (MQ=255)
                                          tgtgGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCt               >  1:488302/1‑68 (MQ=255)
                                          tgtgGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCt               >  1:143952/1‑68 (MQ=255)
                                            tgGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTgg             <  1:2375999/68‑1 (MQ=255)
                                            tgGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTgg             <  1:2501099/68‑1 (MQ=255)
                                            tgGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTgg             <  1:2505196/68‑1 (MQ=255)
                                            tgGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTgg             <  1:2553770/68‑1 (MQ=255)
                                            tgGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTgg             <  1:2633561/68‑1 (MQ=255)
                                            tgGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTgg             <  1:1376520/68‑1 (MQ=255)
                                              gAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTGGaa           <  1:1195896/68‑1 (MQ=255)
                                                     tgGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTGGAATACcgcg    <  1:3005023/68‑1 (MQ=255)
                                                     tgGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTGGAATACcgcg    <  1:1646042/68‑1 (MQ=255)
                                                     tgGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTGGAATACcgcg    <  1:1959588/68‑1 (MQ=255)
                                                      ggCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTGGAATACCGCgg   <  1:2437177/68‑1 (MQ=255)
                                                      ggCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTGGAATACCGCgg   <  1:1357141/68‑1 (MQ=255)
                                                      ggCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTGGAATACCGCgg   <  1:1188944/68‑1 (MQ=255)
                                                      ggCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTGGAATACCGCgg   <  1:3365718/68‑1 (MQ=255)
                                                      ggCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTGGAATACCGCgg   <  1:2329178/68‑1 (MQ=255)
                                                      ggCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTGGAATACCGCGGa  <  1:1937459/69‑1 (MQ=255)
                                                      ggCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTGGAATACCGCGGa  <  1:927613/69‑1 (MQ=255)
                                                      ggCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTGGAATACCGCGGa  <  1:1612121/69‑1 (MQ=255)
                                                      ggCGCGATCACGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTGGAATACCGCGGa  <  1:628313/69‑1 (MQ=255)
                                                      ggCGCGATCACGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTGGAATACCGCGGa  <  1:610712/69‑1 (MQ=255)
                                                      ggCGCGATCACGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTGGAATACCGCGGa  <  1:3001854/69‑1 (MQ=255)
                                                      ggCGCGATCACGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTGGAATACCGCGGa  <  1:1800763/69‑1 (MQ=255)
                                                      ggCGCGATCACGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTGGAATACCGCGGa  <  1:422338/69‑1 (MQ=255)
                                                      ggCGCGATCACGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTGGAATACCGCGGa  <  1:2957802/69‑1 (MQ=255)
                                                      ggCGCGAGCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTGGAATACCGCGGa  <  1:1136616/69‑1 (MQ=255)
                                                       gcgcGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTGGAATACCGCGGa  <  1:1849832/68‑1 (MQ=255)
                                                            aTCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTGGAATACCGCGGa  <  1:1078760/63‑1 (MQ=255)
                                                              cgcgCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTGGAATACCGCGGa  >  1:2240081/1‑61 (MQ=255)
                                                              cgcgCAACGTGATCTAGCAGAAATCCTTCTTGAAGGTTTAcgtc                   >  1:1752935/1‑44 (MQ=255)
                                                               |                                                           
CAGGTTGACCGACAACGATATAAATCGGAATCAAAAACTATGTGTGGAATTGTTGGCGCGATCGCGCAACGTGATGTAGCAGAAATCCTTCTTGAAGGTTTACGTCGTCTGGAATACCGCGGA  >  minE/2472687‑2472809

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: