Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 67,547 C→A 34.6% L241L (CTG→CTT leuA ← 2‑isopropylmalate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE67,5470CA34.6% 132.0 / 87.9 124L241L (CTG→CTTleuA2‑isopropylmalate synthase
Reads supporting (aligned to +/- strand):  ref base C (48/33);  new base A (10/33);  total (58/66)
Fisher's exact test for biased strand distribution p-value = 1.45e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 124 total aligned reads displayed.

TGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGC  >  minE/67482‑67606
                                                                 |                                                           
tgattgatGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCAAg                                                           >  1:405529/1‑68 (MQ=255)
tgattgatGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCAAg                                                           >  1:1806816/1‑68 (MQ=255)
tgattgatGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCAAg                                                           >  1:2497304/1‑68 (MQ=255)
tgattgatGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCAAg                                                           >  1:618998/1‑68 (MQ=255)
tgattgatGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCAAg                                                           >  1:950702/1‑68 (MQ=255)
tgattgatGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCAAg                                                           >  1:1078893/1‑68 (MQ=255)
tgattgatGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCCTCAAg                                                           >  1:2697892/1‑68 (MQ=255)
          cGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATTATGACTTCTTCAAggg                                                         <  1:3283545/60‑1 (MQ=255)
          cGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCAAggg                                                         <  1:3305621/60‑1 (MQ=255)
          cGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCAAggg                                                         <  1:1234374/60‑1 (MQ=255)
          cGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCAAggg                                                         <  1:3235026/60‑1 (MQ=255)
          cGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCAAggg                                                         <  1:3005506/60‑1 (MQ=255)
          cGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCAAggg                                                         <  1:296133/60‑1 (MQ=255)
          cGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCAAggg                                                         <  1:894389/60‑1 (MQ=255)
          cGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCAAggg                                                         <  1:1534184/60‑1 (MQ=255)
          cGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCAAggg                                                         <  1:1581705/60‑1 (MQ=255)
          cGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCAAggg                                                         <  1:878539/60‑1 (MQ=255)
          cGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCAAggg                                                         <  1:2540149/60‑1 (MQ=255)
          cGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCAAggg                                                         <  1:591509/60‑1 (MQ=255)
            gtgtGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTcc                                               >  1:367693/1‑68 (MQ=255)
            gtgtGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTcc                                               >  1:2293618/1‑68 (MQ=255)
            gtgtGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTcc                                               >  1:2628740/1‑68 (MQ=255)
            gtgtGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTcc                                               >  1:1002746/1‑68 (MQ=255)
            gtgtGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTcc                                               >  1:3229904/1‑68 (MQ=255)
            gtgtGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTcc                                               >  1:3377294/1‑68 (MQ=255)
            gtgtGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTcc                                               >  1:129941/1‑68 (MQ=255)
            gtgtGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTcc                                               >  1:1491021/1‑68 (MQ=255)
              gtgGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGc                                            >  1:412612/1‑69 (MQ=255)
              gtgGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGc                                            >  1:3019332/1‑69 (MQ=255)
              gtgGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGc                                            >  1:704983/1‑69 (MQ=255)
              gtgGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGc                                            >  1:2771185/1‑69 (MQ=255)
              gtgGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGc                                            >  1:2690841/1‑69 (MQ=255)
              gtgGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGc                                            >  1:2661112/1‑69 (MQ=255)
               tgGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCa                                           >  1:432736/1‑69 (MQ=255)
                 gACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACg                                         >  1:542293/1‑69 (MQ=255)
                 gACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACg                                         >  1:3299078/1‑69 (MQ=255)
                 gACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACg                                         >  1:3206876/1‑69 (MQ=255)
                 gACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACg                                         >  1:3052791/1‑69 (MQ=255)
                 gACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACg                                         >  1:1750295/1‑69 (MQ=255)
                 gACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACg                                         >  1:1880855/1‑69 (MQ=255)
                 gACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACg                                         <  1:1932499/69‑1 (MQ=255)
                  aCGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACg                                         <  1:601452/68‑1 (MQ=255)
                          aaTATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGa                                 <  1:968827/68‑1 (MQ=255)
                          aaTATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGa                                 <  1:3128733/68‑1 (MQ=255)
                          aaTATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGa                                 <  1:3059014/68‑1 (MQ=255)
                          aaTATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGa                                 <  1:2922694/68‑1 (MQ=255)
                          aaTATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGa                                 <  1:3330948/68‑1 (MQ=255)
                           atatCCTTACGAACTTTGATCGCCATGATGACTTCTTCAAggg                                                         <  1:3298593/43‑1 (MQ=255)
                             atCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATccc                             >  1:739724/1‑69 (MQ=255)
                             atCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATccc                             >  1:1497132/1‑69 (MQ=255)
                             atCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATccc                             >  1:2429024/1‑69 (MQ=255)
                             atCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATccc                             >  1:83014/1‑69 (MQ=255)
                              tCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCgg                                             >  1:2866572/1‑52 (MQ=255)
                              tCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCgg                                             >  1:12082/1‑52 (MQ=255)
                              tCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCgg                                             >  1:2220433/1‑52 (MQ=255)
                              tCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCg                            >  1:2113900/1‑69 (MQ=255)
                              tCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCg                            >  1:3396930/1‑69 (MQ=255)
                              tCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCg                            >  1:179998/1‑69 (MQ=255)
                              tCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCg                            >  1:2667056/1‑69 (MQ=255)
                              tCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCg                            >  1:2892937/1‑69 (MQ=255)
                                cTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGtt                          <  1:2500462/69‑1 (MQ=255)
                                cTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGtt                          <  1:126302/69‑1 (MQ=255)
                                cTTACGAACTTTGAACGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGtt                          <  1:2327423/69‑1 (MQ=255)
                                      aaCTTTGATCGCCATGATGACTTCTTCCAGGGATCAGTTTCCGGCACGCTCGCCGATCCCGTTc                         <  1:2774063/64‑1 (MQ=255)
                                      aaCTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTc                         <  1:2042153/64‑1 (MQ=255)
                                      aaCTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTc                         <  1:528073/64‑1 (MQ=255)
                                      aaCTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTc                         <  1:173919/64‑1 (MQ=255)
                                        cTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTgcgc                  <  1:1152235/69‑1 (MQ=255)
                                        cTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTgcgc                  <  1:2140627/69‑1 (MQ=255)
                                        cTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTgcgc                  <  1:909538/69‑1 (MQ=255)
                                        cTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTgcgc                  <  1:549249/69‑1 (MQ=255)
                                        cTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTgcgc                  <  1:553250/69‑1 (MQ=255)
                                         tttGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTgcgc                  <  1:2142219/68‑1 (MQ=255)
                                         tttGATCGCCATGATGACTTCTTCAAGGGAACAGTTACCGGCACGCTCGCCGATCCCGTTCATTGCGcc                 <  1:950812/69‑1 (MQ=255)
                                         tttGATCGCCATGATGACTTCTTCAAGGGAACAGTTACCGGCACGCTCGCCGATCCCGTTCATTGCGcc                 <  1:389390/69‑1 (MQ=255)
                                         tttGATCGCCATGATGACTTCTTCAAGGGAACAGTTACCGGCACGCTCGCCGATCCCGTTCATTGCGcc                 <  1:2815759/69‑1 (MQ=255)
                                         tttGATCGCCATGATGACTTCTTCAAGGGAACAGTTACCGGCACGCTCGCCGATCCCGTTCATTGCGcc                 <  1:3384217/69‑1 (MQ=255)
                                         tttGATCGCCATGATGACTTCTTCAAGGGAACAGTTACCGGCACGCTCGCCGATCCCGTTCATTGCGcc                 <  1:1020664/69‑1 (MQ=255)
                                         tttGATCGCCATGATGACTTCTTCAAGGGAACAGTTACCGGCACGCTCGCCGATCCCGTTCATTGCGcc                 <  1:1128346/69‑1 (MQ=255)
                                         tttGATCGCCATGATGACTTCTTCAAGGGAACAGTTACCGGCACGCTCGCCGATCCCGTTCATTGCGcc                 <  1:1287469/69‑1 (MQ=255)
                                         tttGATCGCCATGATGACTTCTTCAAGGGAACAGTTACCGGCACGCTCGCCGATCCCGTTCATTGCGcc                 <  1:2065017/69‑1 (MQ=255)
                                         tttGATCGCCATGATGACTTCTTCAAGGGAACAGTTACCGGCACGCTCGCCGATCCCGTTCATTGCGcc                 <  1:2175254/69‑1 (MQ=255)
                                         tttGATCGCCATGATGACTTCTTCAAGGGAACAGTTACCGGCACGCTCGCCGATCCCGTTCATTGCGcc                 <  1:2317967/69‑1 (MQ=255)
                                         tttGATCGCCATGATGACTTCTTCAAGGGAACAGTTACCGGCACGCTCGCCGATCCCGTTCATTGCGcc                 <  1:255200/69‑1 (MQ=255)
                                         tttGATCGCCATGATGACTTCTTCAAGGGAACAGTTACCGGCACGCTCGCCGATCCCGTTCATTGCGcc                 <  1:262866/69‑1 (MQ=255)
                                         tttGATCGCCATGATGACTTCTTCAAGGGAACAGTTACCGGCACGCTCGCCGATCCCGTTCATTGCGcc                 <  1:2721420/69‑1 (MQ=255)
                                         tttGATCGCCATGATGACTTCTTCAAGGGAACAGTTACCGGCACGCTCGCCGATCCCGTTCATTGCGcc                 <  1:3167728/69‑1 (MQ=255)
                                           tGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCt                >  1:2304658/1‑68 (MQ=255)
                                             aTCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTc              <  1:739951/68‑1 (MQ=255)
                                             aTCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTc              <  1:895818/68‑1 (MQ=255)
                                             aTCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTc              <  1:1219043/68‑1 (MQ=255)
                                               cGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCa            <  1:1689084/68‑1 (MQ=255)
                                               cGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCa            <  1:963073/68‑1 (MQ=255)
                                                    tgatgaCTTCTTCCAGTGAACAGTTTCCGGCACGCTCGCCGa                                 <  1:1517413/42‑1 (MQ=255)
                                                     gatgaCTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCg     <  1:613282/69‑1 (MQ=255)
                                                      atgaCTTCTTCCAGGGAACAGTTTCCTGCACGCTCGCCGATCCCGTTCATTGCGCCTTcc             >  1:2487183/1‑60 (MQ=255)
                                                      atgaCTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCg     >  1:2767027/1‑68 (MQ=255)
                                                      atgaCTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCg     >  1:2908553/1‑68 (MQ=255)
                                                      atgaCTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCg     >  1:3410340/1‑68 (MQ=255)
                                                        gaCTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGc  >  1:2957797/1‑69 (MQ=255)
                                                                 |                                                           
TGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGC  >  minE/67482‑67606

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: