Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,651,936 C→T 37.4% A57V (GCG→GTG)  glgC → glucose‑1‑phosphate adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,651,9360CT37.4% 95.2 / 107.0 99A57V (GCG→GTG) glgCglucose‑1‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (37/25);  new base T (18/19);  total (55/44)
Fisher's exact test for biased strand distribution p-value = 3.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 105 total aligned reads displayed.

AAGGATTTAACCAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAGTACCAGTCC  >  minE/2651866‑2651997
                                                                      |                                                             
aaGGATTTAACCAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTgc                                                               >  1:3557886/1‑71 (MQ=255)
aaGGATTTAACCAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTgc                                                               >  1:4232496/1‑71 (MQ=255)
aaGGATTTAACCAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTgc                                                               >  1:152418/1‑71 (MQ=255)
aaGGATATAACCAATAAGCGAGCAAAACCGGCCGTACACTTAGGCGGTAAGTTCCGCATTATCGACTTTgc                                                               >  1:542098/1‑71 (MQ=255)
 aGGATTTAACCAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTGTCGACTTTgcg                                                              >  1:2014908/1‑71 (MQ=255)
 aGGATTTAACCAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTgcg                                                              >  1:1909391/1‑71 (MQ=255)
   gATTTAGCCAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCt                                                            >  1:545008/1‑71 (MQ=255)
   gATTTAACCAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCt                                                            >  1:3846376/1‑71 (MQ=255)
    aTTTAACCAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTGCGCATTATCGACTTTGTGCTg                                                           >  1:1889679/1‑71 (MQ=255)
    aTTTAACCAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGTGCTg                                                           >  1:1975712/1‑71 (MQ=255)
    aTTTAACCAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTg                                                           <  1:4072046/71‑1 (MQ=255)
       tAACCAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTc                                                         >  1:2055179/1‑70 (MQ=255)
        aaCCAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCt                                                        <  1:1466415/70‑1 (MQ=255)
         aCCAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTaa                                                      <  1:3572290/71‑1 (MQ=255)
         aCCAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTa                                                       <  1:1959813/70‑1 (MQ=255)
         aCCAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTa                                                       <  1:3833718/70‑1 (MQ=255)
          ccAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGTGCTGTCTAAc                                                     >  1:1168955/1‑71 (MQ=255)
              tAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCa                                                 >  1:3130160/1‑71 (MQ=255)
               aaGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTaa                                                      >  1:3311663/1‑65 (MQ=255)
                aGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCAt                                                >  1:303797/1‑70 (MQ=255)
                aGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCAt                                                >  1:859315/1‑70 (MQ=255)
                aGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCAt                                                >  1:111773/1‑70 (MQ=255)
                aGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCAt                                                >  1:501775/1‑70 (MQ=255)
                           ccGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGCCTTTGCGCTGTCTAACTGCATCAACTCCGGGAt                                    >  1:1007479/1‑71 (MQ=255)
                           ccGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGAt                                    >  1:1261681/1‑71 (MQ=255)
                            cGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATc                                   >  1:2127378/1‑71 (MQ=255)
                            cGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATc                                   >  1:1213321/1‑71 (MQ=255)
                            cGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTAGCGCTg                                                           <  1:1516904/47‑1 (MQ=255)
                             ggCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATcc                                  <  1:2469320/71‑1 (MQ=255)
                             ggCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATcc                                  <  1:2338135/71‑1 (MQ=255)
                              gCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATCcg                                 >  1:4103458/1‑71 (MQ=255)
                              gCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATCcg                                 >  1:1224954/1‑71 (MQ=255)
                                 gTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCcgtcg                              <  1:1242708/71‑1 (MQ=255)
                                  tACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCcgt                                <  1:3703639/68‑1 (MQ=255)
                                   acacTTCGGCGGTAAGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGTGATCcgtcgt                             <  1:1373832/70‑1 (MQ=255)
                                   acacTTCGGCGGTAAGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATTcgtcgt                             <  1:540774/70‑1 (MQ=255)
                                   acacTTCGGCGGTAAGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATCcgtcgt                             <  1:2803903/70‑1 (MQ=255)
                                   acacTTCGGCGGTAAGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATCcgtcgt                             <  1:865739/70‑1 (MQ=255)
                                   acacTTCGGCGGTAAGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATCcgtcgt                             <  1:1981682/70‑1 (MQ=255)
                                   acacTTCGGCGGTAAGTTCCGCATTATCGACTTTGTGCTGTCTAACTCCATCAACTCCGGGATCcgtcgt                             <  1:3189851/70‑1 (MQ=255)
                                    cacTTCGGCGGTAGGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCcgtcg                              >  1:1427373/1‑68 (MQ=255)
                                    cacTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCcgtcg                              >  1:2383347/1‑68 (MQ=255)
                                    cacTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCcgtcg                              >  1:2082696/1‑68 (MQ=255)
                                    cacTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCcgtcg                              >  1:4203130/1‑68 (MQ=255)
                                    cacTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTAt                           >  1:1869204/1‑71 (MQ=255)
                                    cacTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTAt                           <  1:1942415/71‑1 (MQ=255)
                                    cacTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTAt                           <  1:1028355/71‑1 (MQ=255)
                                    cacTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTATAAATGCATCAACTCCGGGAt                                    >  1:3285883/1‑62 (MQ=255)
                                      cTTCGGCGGTAAGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATgg                         <  1:3904353/71‑1 (MQ=255)
                                           gcggTAAGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGTCGTg                    <  1:2608155/71‑1 (MQ=255)
                                           gcggTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTg                    <  1:2211116/71‑1 (MQ=255)
                                            cggTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCggg                                      <  1:2159896/52‑1 (MQ=255)
                                             ggTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGAt                  >  1:2545779/1‑71 (MQ=255)
                                              gTAAGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTg                    <  1:1337115/68‑1 (MQ=255)
                                              gTAAGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTg                    <  1:101730/68‑1 (MQ=255)
                                              gTAAGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTg                    <  1:3516179/68‑1 (MQ=255)
                                              gTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGAt                                    >  1:1098628/1‑52 (MQ=255)
                                              gTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGAt                                    >  1:1054158/1‑52 (MQ=255)
                                              gTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGa                   <  1:3344063/69‑1 (MQ=255)
                                              gTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGa                   <  1:2189392/69‑1 (MQ=255)
                                                aaGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCg                      >  1:2761467/1‑64 (MQ=255)
                                                aaGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCg                      >  1:3109977/1‑64 (MQ=255)
                                                aaGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGAt                  >  1:1845091/1‑68 (MQ=255)
                                                aaGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCAc               >  1:317250/1‑71 (MQ=255)
                                                   ttCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGc                       >  1:2639470/1‑60 (MQ=255)
                                                   ttCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCg                      >  1:1332126/1‑61 (MQ=255)
                                                   ttCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCg                      >  1:3948710/1‑61 (MQ=255)
                                                   ttCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCg                      >  1:3573887/1‑61 (MQ=255)
                                                     ccGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAg           >  1:1014864/1‑70 (MQ=255)
                                                     ccGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAg           >  1:3590597/1‑70 (MQ=255)
                                                     ccGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAg           >  1:429622/1‑70 (MQ=255)
                                                     ccGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAg           >  1:713725/1‑70 (MQ=255)
                                                     ccGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAg           >  1:1031448/1‑70 (MQ=255)
                                                     ccGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAg           >  1:2755070/1‑70 (MQ=255)
                                                     ccGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAg           <  1:1810317/70‑1 (MQ=255)
                                                     ccGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAg           <  1:3553927/70‑1 (MQ=255)
                                                      cGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAGt          <  1:2833250/70‑1 (MQ=255)
                                                      cGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAGt          <  1:4217314/70‑1 (MQ=255)
                                                      cGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAGt          <  1:1815388/70‑1 (MQ=255)
                                                       gCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAGTAc        <  1:1092403/71‑1 (MQ=255)
                                                        cATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAGTAcc       <  1:2034800/71‑1 (MQ=255)
                                                        cATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAGtgc        >  1:3897425/1‑68 (MQ=255)
                                                          ttATCGACTTTGCTCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAGTACCAg     <  1:156626/71‑1 (MQ=255)
                                                           tATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGa                                     <  1:868546/38‑1 (MQ=255)
                                                           tATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGa                                     <  1:2027834/38‑1 (MQ=255)
                                                           tATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTa                            >  1:4181041/1‑47 (MQ=255)
                                                           tATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTa                            >  1:665668/1‑47 (MQ=255)
                                                           tATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTa                            >  1:4151220/1‑47 (MQ=255)
                                                           tATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTa                            >  1:2860371/1‑47 (MQ=255)
                                                           tATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTa                            >  1:3550232/1‑47 (MQ=255)
                                                            aTCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAGTACCAGTc   >  1:1785577/1‑71 (MQ=255)
                                                            aTCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAGTACCAGTc   >  1:1780425/1‑71 (MQ=255)
                                                             tCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAGTACCAGTcc  <  1:2760881/71‑1 (MQ=255)
                                                              cGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGt                     >  1:3822990/1‑51 (MQ=255)
                                                              cGACTTTGCGCTGTCTAACTGCATCAACTACGGGATCCGTCGTATGGGCGTGATCACCCAGTACCa      >  1:189449/1‑66 (MQ=255)
                                                               gACTTTGCGCTGTCTAACTGCATCAACTCCGGGATc                                   >  1:3527142/1‑36 (MQ=255)
                                                               gACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGAt                  >  1:1863667/1‑53 (MQ=255)
                                                                  tttGCGCTGTCTAACTGCATCAACTCCGGGATCcgtcg                              <  1:1803006/38‑1 (MQ=255)
                                                                  tttGCGCTGTCTAACTGCATCAACTCCGGGATCcgtcg                              <  1:1777828/38‑1 (MQ=255)
                                                                  tttGCGCTGTCTAACTGCATCAACTCCGGGATCcgtcg                              <  1:1220124/38‑1 (MQ=255)
                                                                      |                                                             
AAGGATTTAACCAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAGTACCAGTCC  >  minE/2651866‑2651997

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