Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,790,341 A→G 22.1% V462V (GTA→GTG pgi → glucosephosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,790,3410AG22.1% 108.1 / 23.1 81V462V (GTA→GTGpgiglucosephosphate isomerase
Reads supporting (aligned to +/- strand):  ref base A (26/37);  new base G (11/7);  total (37/44)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGAGCAGGAATATCGTGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAGCCTGGGTGCG  >  minE/2790273‑2790407
                                                                    |                                                                  
gTTGAGCAGGAATATCGTGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTa                                                                    <  1:1901944/69‑1 (MQ=255)
       agtaaTATCGTCATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAg                                                             <  1:3114461/66‑1 (MQ=255)
       aGGAATATCGTGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAg                                                             <  1:597233/69‑1 (MQ=255)
         gAATATCGCGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTGTTCGAAGGt                                                           <  1:852567/69‑1 (MQ=255)
         gAATATCGCGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTGTTCGAAGGt                                                           <  1:3012227/69‑1 (MQ=255)
         gAATATCGCGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTGTTCGAAGGt                                                           <  1:733645/69‑1 (MQ=255)
         gAATATCGCGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTGTTCGAAGGt                                                           <  1:613915/69‑1 (MQ=255)
         gAATATCGCGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTGTTCGAAGGt                                                           <  1:1161077/69‑1 (MQ=255)
         gAATATCGCGATCAGGGTAAAGATCCGGCAACGCTTGACTACGGGGTGCCGTTCAAAGTGTTCGAAGGt                                                           <  1:1673132/69‑1 (MQ=255)
          aaTATCGCGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTGTTCGAAGGt                                                           <  1:1720583/68‑1 (MQ=255)
              tCGTGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAAcc                                                       >  1:1128170/1‑68 (MQ=255)
              tCGTGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAAcc                                                       >  1:1032796/1‑68 (MQ=255)
              tCGTGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCg                                                      >  1:327533/1‑69 (MQ=255)
              tCGTGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCg                                                      >  1:3381342/1‑69 (MQ=255)
              tCGTGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCg                                                      >  1:2728826/1‑69 (MQ=255)
              tCGTGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCg                                                      >  1:2480934/1‑69 (MQ=255)
              tCGTGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCg                                                      >  1:2384150/1‑69 (MQ=255)
              tCGTGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCg                                                      >  1:2316196/1‑69 (MQ=255)
              tCGTGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCg                                                      >  1:2266440/1‑69 (MQ=255)
              tCGTGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCg                                                      >  1:2071623/1‑69 (MQ=255)
              tCGTGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCg                                                      >  1:903228/1‑69 (MQ=255)
              tCGTGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCg                                                      >  1:1835479/1‑69 (MQ=255)
              tCGTGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCg                                                      >  1:945602/1‑69 (MQ=255)
              tCGTGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCg                                                      >  1:1033661/1‑69 (MQ=255)
              tCGTGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCg                                                      >  1:1048577/1‑69 (MQ=255)
              tCGTGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCg                                                      >  1:1205204/1‑69 (MQ=255)
              tCGTGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCg                                                      >  1:1346305/1‑69 (MQ=255)
              tCGTGATAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCGc                                                     >  1:279815/1‑69 (MQ=255)
                         gTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAAc                                            >  1:1811460/1‑68 (MQ=255)
                         gTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAAc                                            >  1:2411582/1‑68 (MQ=255)
                          tAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTGTTCGAAGGTAACCGCCCGACTAACt                                           >  1:491264/1‑68 (MQ=255)
                          tAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTGTTCGAAGGTAACCGCCCGACTAACTc                                          >  1:2945126/1‑69 (MQ=255)
                          tAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTGTTCGAAGGTAACCGCCCGACTAACTc                                          >  1:1241072/1‑69 (MQ=255)
                          tAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTGTTCGAAGGTAACCGCCCGACTAACTc                                          >  1:3162159/1‑69 (MQ=255)
                          tAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTGTTCGAAGGTAACCGCCCGACTAACTc                                          >  1:2832866/1‑69 (MQ=255)
                          tAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTGTTCGAAGGTAACCGCCCGACTAACTc                                          >  1:1392192/1‑69 (MQ=255)
                          tAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTGTTCGAAGGTAACCGCCCGACTAACTc                                          >  1:2273335/1‑69 (MQ=255)
                          tAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTGTTCGAAGGTAACCGCCCGACTAACTc                                          >  1:1143218/1‑69 (MQ=255)
                          tAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTGTTCGAAGGTAACCGCCCGACTAACTc                                          >  1:2230343/1‑69 (MQ=255)
                          tAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTGTTCGAAGGTAACCGCCCGACTAACTc                                          >  1:1631713/1‑69 (MQ=255)
                          tAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTGTTCGAAGGTAACCGCCCGACTAACTc                                          >  1:1138935/1‑69 (MQ=255)
                          tAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTc                                          <  1:1556962/69‑1 (MQ=255)
                          tAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTc                                          <  1:2985445/69‑1 (MQ=255)
                          tAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTc                                          <  1:3239682/69‑1 (MQ=255)
                          tAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTc                                          <  1:2844644/69‑1 (MQ=255)
                          tAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTc                                          <  1:800049/69‑1 (MQ=255)
                          tAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTc                                          <  1:400420/69‑1 (MQ=255)
                          tAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTc                                          <  1:2386668/69‑1 (MQ=255)
                          tAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTc                                          <  1:2382691/69‑1 (MQ=255)
                           aaaGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTc                                          <  1:2800227/68‑1 (MQ=255)
                            aaGATCCGGCAACGGTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCa                                        <  1:2611067/69‑1 (MQ=255)
                            aaGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCa                                        <  1:1906265/69‑1 (MQ=255)
                              gATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCAt                                       <  1:876045/68‑1 (MQ=255)
                                      aaCGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCtgctgc                               >  1:1187133/1‑68 (MQ=255)
                                               cTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATc                       <  1:1937384/67‑1 (MQ=255)
                                                        gCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCTCTCCGTTCAg            <  1:2949598/69‑1 (MQ=255)
                                                        gCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAg            <  1:695260/69‑1 (MQ=255)
                                                        gCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAg            <  1:673495/69‑1 (MQ=255)
                                                        gCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAg            <  1:1021994/69‑1 (MQ=255)
                                                        gCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAg            <  1:2594002/69‑1 (MQ=255)
                                                        gCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAg            <  1:1686105/69‑1 (MQ=255)
                                                        gCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAg            <  1:1847772/69‑1 (MQ=255)
                                                        gCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAg            <  1:1884267/69‑1 (MQ=255)
                                                        gCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAg            <  1:1926546/69‑1 (MQ=255)
                                                        gCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAg            <  1:2643113/69‑1 (MQ=255)
                                                        gCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAg            <  1:1286199/69‑1 (MQ=255)
                                                        gCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAg            <  1:2865596/69‑1 (MQ=255)
                                                        gCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAg            <  1:2914337/69‑1 (MQ=255)
                                                        gCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAg            <  1:2918942/69‑1 (MQ=255)
                                                         ccGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAg            <  1:1140291/68‑1 (MQ=255)
                                                            ttCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAGCCTg        <  1:3170899/69‑1 (MQ=255)
                                                            ttCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAGCCTg        <  1:2889180/69‑1 (MQ=255)
                                                            ttCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAGCCTg        <  1:2567775/69‑1 (MQ=255)
                                                            ttCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAGCCTg        <  1:1575916/69‑1 (MQ=255)
                                                            ttCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAGCCTg        <  1:1609013/69‑1 (MQ=255)
                                                            ttCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAGCCTg        <  1:125492/69‑1 (MQ=255)
                                                            ttCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAGCCTg        <  1:2974171/69‑1 (MQ=255)
                                                              cAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCg                 >  1:3284749/1‑58 (MQ=255)
                                                              cAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAGCCTggg      >  1:2922884/1‑69 (MQ=255)
                                                              cAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAGCCTggg      >  1:2159632/1‑69 (MQ=255)
                                                              cAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAGCCTgg       >  1:2100312/1‑68 (MQ=255)
                                                                   tATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAGCCTGGGTGCg  >  1:2028366/1‑68 (MQ=255)
                                                                    |                                                                  
GTTGAGCAGGAATATCGTGATCAGGGTAAAGATCCGGCAACGCTTGACTACGTGGTGCCGTTCAAAGTATTCGAAGGTAACCGCCCGACCAACTCCATCCTGCTGCGTGAAATCACTCCGTTCAGCCTGGGTGCG  >  minE/2790273‑2790407

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: