Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,325,820 A→G 25.3% N397N (AAT→AAC wcaK ← predicted pyruvyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,325,8200AG25.3% 76.5 / 23.0 67N397N (AAT→AACwcaKpredicted pyruvyl transferase
Reads supporting (aligned to +/- strand):  ref base A (23/27);  new base G (9/8);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 7.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCTGCCGTCTAATAAATGACG  >  minE/1325773‑1325888
                                               |                                                                    
tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                                                      >  1:2038858/1‑48 (MQ=255)
tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                                                      >  1:3065589/1‑48 (MQ=255)
tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                                                      >  1:2757351/1‑48 (MQ=255)
tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                                                      >  1:2455628/1‑48 (MQ=255)
tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                                                      >  1:1809017/1‑48 (MQ=255)
tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                                                      >  1:1804963/1‑48 (MQ=255)
tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                                                      >  1:1681624/1‑48 (MQ=255)
  cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa                                               >  1:1029790/1‑69 (MQ=255)
  cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa                                               >  1:706214/1‑69 (MQ=255)
  cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa                                               >  1:2850444/1‑69 (MQ=255)
  cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa                                               >  1:2812535/1‑69 (MQ=255)
  cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa                                               >  1:1127694/1‑69 (MQ=255)
  cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa                                               >  1:3279818/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:859209/69‑1 (MQ=255)
      ccTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTa                                           >  1:3178922/1‑69 (MQ=255)
      ccTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTa                                           >  1:3304678/1‑69 (MQ=255)
      ccTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTa                                           <  1:1925704/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:2260783/69‑1 (MQ=255)
         gTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTAAGGTATcc                                        <  1:929719/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)
         gTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTAAGGTATcc                                        <  1:58148/69‑1 (MQ=255)
           cTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTTGCCTAAGGTATCCGc                                      >  1:1351866/1‑69 (MQ=255)
           cTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc                                      >  1:320986/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)
           cTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc                                      >  1:439813/1‑69 (MQ=255)
            tGACGCTCACGACTTACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc                                      <  1:1974217/68‑1 (MQ=255)
            tGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc                                      <  1:304324/68‑1 (MQ=255)
            tGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc                                      <  1:645918/68‑1 (MQ=255)
            tGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc                                      <  1:3394256/68‑1 (MQ=255)
            tGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc                                      <  1:642028/68‑1 (MQ=255)
            tGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc                                      <  1:3383041/68‑1 (MQ=255)
            tGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc                                      <  1:1177218/68‑1 (MQ=255)
            tGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc                                      <  1:610278/68‑1 (MQ=255)
                gCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAAcc                                  <  1:343588/68‑1 (MQ=255)
                gCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAAcc                                  <  1:937258/68‑1 (MQ=255)
                gCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAAcc                                  <  1:896967/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:3179805/1‑69 (MQ=255)
                   cACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCg                              >  1:1379696/1‑69 (MQ=255)
                   cACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCg                              >  1:1663880/1‑69 (MQ=255)
                   cACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCg                              >  1:69348/1‑69 (MQ=255)
                   cACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCg                              >  1:1777696/1‑69 (MQ=255)
                   cACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCg                              >  1:1777824/1‑69 (MQ=255)
                   cACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCg                              >  1:634031/1‑69 (MQ=255)
                   cACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCg                              >  1:1987622/1‑69 (MQ=255)
                   cACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCg                              >  1:2478419/1‑69 (MQ=255)
                   cACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCg                              >  1:2534953/1‑69 (MQ=255)
                   cACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCg                              >  1:2607373/1‑69 (MQ=255)
                   cACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCg                              >  1:2792434/1‑69 (MQ=255)
                       aCTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTg                          <  1:882277/69‑1 (MQ=255)
                        cTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTg                          <  1:2702791/68‑1 (MQ=255)
                         tGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCa                        <  1:1944178/69‑1 (MQ=255)
                         tGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCa                        <  1:3290615/69‑1 (MQ=255)
                         tGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCa                        <  1:3184871/69‑1 (MQ=255)
                         tGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCa                        <  1:136864/69‑1 (MQ=255)
                         tGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCa                        <  1:932546/69‑1 (MQ=255)
                         tGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCa                        <  1:1002307/69‑1 (MQ=255)
                                          cgcgCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCTGCCGTCtaataa        >  1:3403532/1‑68 (MQ=255)
                                          cgcgCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCTGCCGTCtaataa        >  1:1921627/1‑68 (MQ=255)
                                             gcATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCa                        <  1:3403639/49‑1 (MQ=255)
                                               aTTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCTGCCGTCTAATAAATGACg  <  1:3050463/69‑1 (MQ=255)
                                               |                                                                    
TGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCTGCCGTCTAATAAATGACG  >  minE/1325773‑1325888

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: