Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,325,790 C→T 35.7% E407E (GAG→GAA wcaK ← predicted pyruvyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,325,7900CT35.7% 46.6 / 36.0 67E407E (GAG→GAAwcaKpredicted pyruvyl transferase
Reads supporting (aligned to +/- strand):  ref base C (19/24);  new base T (9/15);  total (28/39)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCAT  >  minE/1325728‑1325858
                                                              |                                                                    
cTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGa                                                                >  1:1100955/1‑69 (MQ=255)
cTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGa                                                                >  1:2431438/1‑69 (MQ=255)
cTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGa                                                                >  1:3140069/1‑69 (MQ=255)
      ttCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAc                                                           <  1:1681056/68‑1 (MQ=255)
      ttCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAc                                                           <  1:1950095/68‑1 (MQ=255)
      ttCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAc                                                           <  1:3001321/68‑1 (MQ=255)
      ttCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAc                                                           <  1:3052545/68‑1 (MQ=255)
      ttCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAc                                                           <  1:1526288/68‑1 (MQ=255)
      ttCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAc                                                           <  1:3425260/68‑1 (MQ=255)
      ttCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAc                                                           <  1:637596/68‑1 (MQ=255)
       tCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCt                                                                     <  1:136255/57‑1 (MQ=255)
       tCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCt                                                                     <  1:2280132/57‑1 (MQ=255)
           cTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGACTGACAGCTTc                                                     <  1:3288939/69‑1 (MQ=255)
           cTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGACTGACAGCTTc                                                     <  1:385639/69‑1 (MQ=255)
           cTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGACTGACAGCTTc                                                     <  1:2904518/69‑1 (MQ=255)
           cTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGACTGACAGCTTc                                                     <  1:2507002/69‑1 (MQ=255)
           cTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGACTGACAGCTTc                                                     <  1:551730/69‑1 (MQ=255)
           cTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGACTGACAGCTTc                                                     <  1:2101050/69‑1 (MQ=255)
                                    tGCACCATCTGCATTCCCGTCTGACGTTCGCGACTGACAGCTTc                                                     <  1:1553656/44‑1 (MQ=255)
                                             tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                        >  1:2757351/1‑48 (MQ=255)
                                             tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                        >  1:3065589/1‑48 (MQ=255)
                                             tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                        >  1:2455628/1‑48 (MQ=255)
                                             tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                        >  1:1681624/1‑48 (MQ=255)
                                             tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                        >  1:1804963/1‑48 (MQ=255)
                                             tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                        >  1:2038858/1‑48 (MQ=255)
                                             tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                        >  1:1809017/1‑48 (MQ=255)
                                               cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa                 >  1:219926/1‑69 (MQ=255)
                                               cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa                 >  1:3279818/1‑69 (MQ=255)
                                               cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa                 >  1:2850444/1‑69 (MQ=255)
                                               cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa                 >  1:2812535/1‑69 (MQ=255)
                                               cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa                 >  1:1029790/1‑69 (MQ=255)
                                               cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa                 >  1:3408956/1‑69 (MQ=255)
                                               cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa                 >  1:1876401/1‑69 (MQ=255)
                                               cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa                 >  1:1127694/1‑69 (MQ=255)
                                               cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa                 >  1:706214/1‑69 (MQ=255)
                                                 ttCCTGTCTGACGCTCACTACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCt                   <  1:2404167/65‑1 (MQ=255)
                                                   ccTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTa             <  1:1925704/69‑1 (MQ=255)
                                                   ccTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTa             >  1:3304678/1‑69 (MQ=255)
                                                   ccTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTa             <  1:859209/69‑1 (MQ=255)
                                                   ccTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTa             >  1:3178922/1‑69 (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:929719/69‑1 (MQ=255)
                                                      gTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTAAGGTATcc          <  1:2260783/69‑1 (MQ=255)
                                                      gTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTAAGGTATcc          <  1:3016424/69‑1 (MQ=255)
                                                      gTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTAAGGTATcc          <  1:2890764/69‑1 (MQ=255)
                                                      gTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTAAGGTATcc          <  1:2504157/69‑1 (MQ=255)
                                                        cTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTTGCCTAAGGTATCCGc        >  1:1351866/1‑69 (MQ=255)
                                                        cTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc        >  1:320986/1‑69 (MQ=255)
                                                        cTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc        >  1:2383085/1‑69 (MQ=255)
                                                        cTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc        >  1:1301265/1‑69 (MQ=255)
                                                        cTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc        >  1:439813/1‑69 (MQ=255)
                                                        cTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc        >  1:2690655/1‑69 (MQ=255)
                                                         tGACGCTCACGACTTACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc        <  1:1974217/68‑1 (MQ=255)
                                                         tGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc        <  1:1177218/68‑1 (MQ=255)
                                                         tGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc        <  1:645918/68‑1 (MQ=255)
                                                         tGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc        <  1:642028/68‑1 (MQ=255)
                                                         tGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc        <  1:610278/68‑1 (MQ=255)
                                                         tGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc        <  1:304324/68‑1 (MQ=255)
                                                         tGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc        <  1:3394256/68‑1 (MQ=255)
                                                         tGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc        <  1:3383041/68‑1 (MQ=255)
                                                             gCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAAcc    <  1:343588/68‑1 (MQ=255)
                                                             gCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAAcc    <  1:896967/68‑1 (MQ=255)
                                                             gCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAAcc    <  1:227244/68‑1 (MQ=255)
                                                             gCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAAcc    <  1:937258/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)
                                                              |                                                                    
CTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCAT  >  minE/1325728‑1325858

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: