Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,478,347 T→G 34.5% N73H (AAT→CAT)  glpQ ← periplasmic glycerophosphodiester phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,478,3470TG34.5% 107.9 / 112.2 136N73H (AAT→CAT) glpQperiplasmic glycerophosphodiester phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (69/20);  new base G (26/21);  total (95/41)
Fisher's exact test for biased strand distribution p-value = 1.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 137 total aligned reads displayed.

CCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTGCGCATACGCCATCGCTTTTG  >  minE/1478284‑1478415
                                                               |                                                                    
ccGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTgtcgt                                                                 >  1:1056507/1‑69 (MQ=255)
ccGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTgtcgt                                                                 >  1:820304/1‑69 (MQ=255)
ccGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTgtcgt                                                                 >  1:3430518/1‑69 (MQ=255)
ccGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTgtcgt                                                                 >  1:3066284/1‑69 (MQ=255)
ccGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTgtcgt                                                                 >  1:2757899/1‑69 (MQ=255)
ccGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTgtcgt                                                                 >  1:2745468/1‑69 (MQ=255)
ccGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTgtcgt                                                                 >  1:2580908/1‑69 (MQ=255)
ccGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTgtcgt                                                                 >  1:2487432/1‑69 (MQ=255)
ccGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTgtcgt                                                                 >  1:2199872/1‑69 (MQ=255)
ccGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTgtcgt                                                                 >  1:2186867/1‑69 (MQ=255)
ccGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTgtcgt                                                                 >  1:1761670/1‑69 (MQ=255)
ccGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTgtcgt                                                                 >  1:1735514/1‑69 (MQ=255)
ccGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTgtcgt                                                                 >  1:1105758/1‑69 (MQ=255)
   tgAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                                                              <  1:959685/68‑1 (MQ=255)
   ggAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                                                              <  1:37479/69‑1 (MQ=255)
   ggAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                                                              <  1:2405065/69‑1 (MQ=255)
   ggAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                                                              <  1:1138316/69‑1 (MQ=255)
   ggAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                                                              <  1:2600779/69‑1 (MQ=255)
   ggAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                                                              <  1:132948/69‑1 (MQ=255)
   ggAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                                                              <  1:2970539/69‑1 (MQ=255)
   ggAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                                                              <  1:548923/69‑1 (MQ=255)
   ggAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                                                              <  1:1411602/69‑1 (MQ=255)
   ggAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                                                              <  1:1317201/69‑1 (MQ=255)
   ggAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                                                              <  1:1709032/69‑1 (MQ=255)
    gAAACGATCGGCAAGATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                                                              <  1:577256/68‑1 (MQ=255)
    gAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                                                              <  1:2668744/68‑1 (MQ=255)
     aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtttt                                                            >  1:2463386/1‑68 (MQ=255)
     aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCttt                                                             >  1:852224/1‑68 (MQ=255)
     aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCttt                                                             >  1:1644553/1‑68 (MQ=255)
     aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCttt                                                             >  1:2298165/1‑68 (MQ=255)
     aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                                                            >  1:1067789/1‑69 (MQ=255)
     aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                                                            >  1:926776/1‑69 (MQ=255)
     aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                                                            >  1:1669331/1‑69 (MQ=255)
     aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                                                            >  1:2856186/1‑69 (MQ=255)
     aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                                                            >  1:2887505/1‑69 (MQ=255)
     aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                                                            >  1:2918570/1‑69 (MQ=255)
     aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                                                            >  1:2130923/1‑69 (MQ=255)
     aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                                                            >  1:2080745/1‑69 (MQ=255)
     aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                                                            >  1:371344/1‑69 (MQ=255)
     aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                                                            >  1:488718/1‑69 (MQ=255)
     aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                                                            >  1:1968120/1‑69 (MQ=255)
     aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                                                            >  1:62362/1‑69 (MQ=255)
     aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                                                            >  1:683833/1‑69 (MQ=255)
                              gATCGAGGTAATGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGAt                                   <  1:272239/69‑1 (MQ=255)
                              gATCGAGGTAATGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGAt                                   <  1:1013334/69‑1 (MQ=255)
                              gATCGAGGTAATGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGAt                                   <  1:699228/69‑1 (MQ=255)
                               aTCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                                                              <  1:2000381/41‑1 (MQ=38)
                                       aaTGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCa                                      >  1:3175056/1‑57 (MQ=255)
                                       aaTGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCa                                      >  1:1356479/1‑57 (MQ=255)
                                       aaTGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCa                                      >  1:80220/1‑57 (MQ=255)
                                       aaTGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCa                                      >  1:1924303/1‑57 (MQ=255)
                                       aaTGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCa                                      >  1:2018431/1‑57 (MQ=255)
                                       aaTGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCa                                      >  1:3211643/1‑57 (MQ=255)
                                       aaTGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCa                                      >  1:2703161/1‑57 (MQ=255)
                                        aTGGTCATGCAGAACAACCAGATTGTCGTCTTTGGt                                                          >  1:2466214/1‑36 (MQ=255)
                                        aTGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCtac                         >  1:2647633/1‑67 (MQ=255)
                                        aTGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCt                           >  1:809368/1‑67 (MQ=255)
                                        aTGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTcc                         >  1:864199/1‑69 (MQ=255)
                                        aTGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTcc                         >  1:2310111/1‑69 (MQ=255)
                                        aTGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTcc                         >  1:1088323/1‑69 (MQ=255)
                                        aTGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTcc                         >  1:1235223/1‑69 (MQ=255)
                                        aTGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTcc                         >  1:1860786/1‑69 (MQ=255)
                                        aTGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTcc                         >  1:1972142/1‑69 (MQ=255)
                                        aTGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTcc                         >  1:515285/1‑69 (MQ=255)
                                        aTGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTcc                         >  1:2365545/1‑69 (MQ=255)
                                        aTGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTcc                         >  1:1217747/1‑69 (MQ=255)
                                        aTGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTcc                         >  1:2938348/1‑69 (MQ=255)
                                        aTGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTc                          >  1:770787/1‑68 (MQ=255)
                                           gTCATGCATAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      >  1:2352694/1‑69 (MQ=255)
                                           gTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      >  1:2335117/1‑69 (MQ=255)
                                           gTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      >  1:671371/1‑69 (MQ=255)
                                           gTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      >  1:2782839/1‑69 (MQ=255)
                                           gTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      >  1:2749448/1‑69 (MQ=255)
                                           gTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      >  1:2724638/1‑69 (MQ=255)
                                           gTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      >  1:1477930/1‑69 (MQ=255)
                                           gTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      >  1:206795/1‑69 (MQ=255)
                                           gTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      >  1:2229910/1‑69 (MQ=255)
                                           gTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      >  1:3309538/1‑69 (MQ=255)
                                           gTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      >  1:3218973/1‑69 (MQ=255)
                                           gTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      >  1:88099/1‑69 (MQ=255)
                                           gTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      >  1:24110/1‑69 (MQ=255)
                                           gTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      >  1:1077841/1‑69 (MQ=255)
                                           gTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      >  1:1403981/1‑69 (MQ=255)
                                           gTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      >  1:1159675/1‑69 (MQ=255)
                                             taTGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      <  1:2877329/66‑1 (MQ=255)
                                             cATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTgcg                    <  1:2670756/69‑1 (MQ=255)
                                                gCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTGCGCa                  <  1:3349535/68‑1 (MQ=255)
                                                 cAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      <  1:2806204/63‑1 (MQ=255)
                                                 cAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      <  1:2567398/63‑1 (MQ=255)
                                                 cAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      <  1:726952/63‑1 (MQ=255)
                                                 cAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      <  1:2244794/63‑1 (MQ=255)
                                                 cAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      <  1:239464/63‑1 (MQ=255)
                                                 cAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      <  1:1728880/63‑1 (MQ=255)
                                                 cAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      <  1:397536/63‑1 (MQ=255)
                                                 cAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      <  1:2211236/63‑1 (MQ=255)
                                                 cAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTGCGCATa                >  1:1147279/1‑69 (MQ=255)
                                                 cAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTGCGCATa                >  1:3050147/1‑69 (MQ=255)
                                                 cAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTGCGCATa                >  1:2969914/1‑69 (MQ=255)
                                                  aGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTg                      <  1:2848041/62‑1 (MQ=255)
                                                               tGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTGCGCATACGCCATCGCTTTTg  >  1:1752547/1‑69 (MQ=255)
                                                               |                                                                    
CCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGATAATCCGCTCCCTGCGCATACGCCATCGCTTTTG  >  minE/1478284‑1478415

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: