Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,325,793 A→G 29.9% R406R (CGT→CGC wcaK ← predicted pyruvyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,325,7930AG29.9% 73.7 / 37.1 80R406R (CGT→CGCwcaKpredicted pyruvyl transferase
Reads supporting (aligned to +/- strand):  ref base A (31/25);  new base G (9/15);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

CTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCG  >  minE/1325728‑1325860
                                                                 |                                                                   
cTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGa                                                                  >  1:1100955/1‑69 (MQ=255)
cTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGa                                                                  >  1:3140069/1‑69 (MQ=255)
cTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGa                                                                  >  1:2431438/1‑69 (MQ=255)
      ttCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAc                                                             <  1:1681056/68‑1 (MQ=255)
      ttCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAc                                                             <  1:1950095/68‑1 (MQ=255)
      ttCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAc                                                             <  1:637596/68‑1 (MQ=255)
      ttCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAc                                                             <  1:3001321/68‑1 (MQ=255)
      ttCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAc                                                             <  1:1526288/68‑1 (MQ=255)
      ttCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAc                                                             <  1:3052545/68‑1 (MQ=255)
      ttCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAc                                                             <  1:3425260/68‑1 (MQ=255)
           cTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGACTGACAGCTTc                                                       <  1:385639/69‑1 (MQ=255)
           cTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGACTGACAGCTTc                                                       <  1:2904518/69‑1 (MQ=255)
           cTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGACTGACAGCTTc                                                       <  1:551730/69‑1 (MQ=255)
           cTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGACTGACAGCTTc                                                       <  1:2507002/69‑1 (MQ=255)
           cTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGACTGACAGCTTc                                                       <  1:3288939/69‑1 (MQ=255)
           cTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGACTGACAGCTTc                                                       <  1:2101050/69‑1 (MQ=255)
                                    tGCACCATCTGCATTCCCGTCTGACGTTCGCGACTGACAGCTTc                                                       <  1:1553656/44‑1 (MQ=255)
                                             tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                          >  1:3065589/1‑48 (MQ=255)
                                             tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                          >  1:2455628/1‑48 (MQ=255)
                                             tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                          >  1:2757351/1‑48 (MQ=255)
                                             tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                          >  1:2038858/1‑48 (MQ=255)
                                             tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                          >  1:1681624/1‑48 (MQ=255)
                                             tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                          >  1:1804963/1‑48 (MQ=255)
                                             tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                          >  1:1809017/1‑48 (MQ=255)
                                               cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa                   >  1:1029790/1‑69 (MQ=255)
                                               cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa                   >  1:2850444/1‑69 (MQ=255)
                                               cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa                   >  1:2812535/1‑69 (MQ=255)
                                               cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa                   >  1:706214/1‑69 (MQ=255)
                                               cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa                   >  1:3279818/1‑69 (MQ=255)
                                               cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa                   >  1:1127694/1‑69 (MQ=255)
                                               cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa                   >  1:3408956/1‑69 (MQ=255)
                                               cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa                   >  1:1876401/1‑69 (MQ=255)
                                               cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa                   >  1:219926/1‑69 (MQ=255)
                                                 ttCCTGTCTGACGCTCACTACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCt                     <  1:2404167/65‑1 (MQ=255)
                                                   ccTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTa               >  1:3304678/1‑69 (MQ=255)
                                                   ccTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTa               >  1:3178922/1‑69 (MQ=255)
                                                   ccTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTa               <  1:1925704/69‑1 (MQ=255)
                                                   ccTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTa               <  1:859209/69‑1 (MQ=255)
                                                   ccTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGACGGAAGCTGGCCTAAGGTa               >  1:1520830/1‑69 (MQ=255)
                                                      gTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTAAGGTAccc            <  1:182495/69‑1 (MQ=255)
                                                      gTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTAAGGTATcc            <  1:571624/69‑1 (MQ=255)
                                                      gTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTAAGGTATcc            <  1:58148/69‑1 (MQ=255)
                                                      gTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTAAGGTATcc            <  1:2504157/69‑1 (MQ=255)
                                                      gTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTAAGGTATcc            <  1:2890764/69‑1 (MQ=255)
                                                      gTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTAAGGTATcc            <  1:2260783/69‑1 (MQ=255)
                                                      gTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTAAGGTATcc            <  1:3016424/69‑1 (MQ=255)
                                                      gTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTAAGGTATcc            <  1:929719/69‑1 (MQ=255)
                                                        cTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTTGCCTAAGGTATCCGc          >  1:1351866/1‑69 (MQ=255)
                                                        cTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc          >  1:320986/1‑69 (MQ=255)
                                                        cTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc          >  1:439813/1‑69 (MQ=255)
                                                        cTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc          >  1:1301265/1‑69 (MQ=255)
                                                        cTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc          >  1:2690655/1‑69 (MQ=255)
                                                        cTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc          >  1:2383085/1‑69 (MQ=255)
                                                         tGACGCTCACGACTTACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc          <  1:1974217/68‑1 (MQ=255)
                                                         tGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc          <  1:304324/68‑1 (MQ=255)
                                                         tGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc          <  1:610278/68‑1 (MQ=255)
                                                         tGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc          <  1:1177218/68‑1 (MQ=255)
                                                         tGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc          <  1:642028/68‑1 (MQ=255)
                                                         tGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc          <  1:3383041/68‑1 (MQ=255)
                                                         tGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc          <  1:3394256/68‑1 (MQ=255)
                                                         tGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc          <  1:645918/68‑1 (MQ=255)
                                                             gCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAAcc      <  1:937258/68‑1 (MQ=255)
                                                             gCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAAcc      <  1:896967/68‑1 (MQ=255)
                                                             gCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAAcc      <  1:343588/68‑1 (MQ=255)
                                                             gCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAAcc      <  1:227244/68‑1 (MQ=255)
                                                             ccTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAAcc      <  1:245203/67‑1 (MQ=255)
                                                              cTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCAt    <  1:1250490/69‑1 (MQ=255)
                                                              cTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCAt    <  1:2702942/69‑1 (MQ=255)
                                                                cACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCg  >  1:2607373/1‑69 (MQ=255)
                                                                cACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCg  >  1:2534953/1‑69 (MQ=255)
                                                                cACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCg  >  1:2478419/1‑69 (MQ=255)
                                                                cACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCg  >  1:634031/1‑69 (MQ=255)
                                                                cACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCg  >  1:2792434/1‑69 (MQ=255)
                                                                cACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCg  >  1:1987622/1‑69 (MQ=255)
                                                                cACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCg  >  1:69348/1‑69 (MQ=255)
                                                                cACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCg  >  1:1777824/1‑69 (MQ=255)
                                                                cACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCg  >  1:1777696/1‑69 (MQ=255)
                                                                cACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCg  >  1:1663880/1‑69 (MQ=255)
                                                                cACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCg  >  1:1379696/1‑69 (MQ=255)
                                                                cACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCg  >  1:3179805/1‑69 (MQ=255)
                                                                 |                                                                   
CTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCG  >  minE/1325728‑1325860

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: