Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,478,321 G→A 47.4% Y81Y (TAC→TAT glpQ ← periplasmic glycerophosphodiester phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,478,3210GA47.4% 13.4 / 103.0 118Y81Y (TAC→TATglpQperiplasmic glycerophosphodiester phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (46/16);  new base A (35/21);  total (81/37)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

Only 100 of 118 total aligned reads displayed.

GTAGTAACGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGAT  >  minE/1478255‑1478382
                                                                  |                                                             
gtagtaACGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGt                                                              <  1:485365/68‑1 (MQ=255)
gtagtaACGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTa                                                             <  1:2563930/69‑1 (MQ=255)
gtagtaACGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTa                                                             <  1:2933992/69‑1 (MQ=255)
gtagtaACGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTa                                                             <  1:2206673/69‑1 (MQ=255)
 tagtaACGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTa                                                             <  1:1935810/68‑1 (MQ=255)
      aCGACCGTCTTTGCGCGCTCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTc                                                       >  1:2226249/1‑69 (MQ=255)
      aCGACCGTCTTTGCGCGCTCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTc                                                       >  1:2892917/1‑69 (MQ=255)
      aCGACCGTCTTTGCGCGCTCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTc                                                       >  1:537388/1‑69 (MQ=255)
      aCGACCGTCTTTGCGCGCTCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTc                                                       >  1:181016/1‑69 (MQ=255)
      aCGACCGTCTTTGCGCGCTCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTc                                                       >  1:3193399/1‑69 (MQ=255)
      aCGACCGTCTTTGCGCGCTCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTc                                                       >  1:1512131/1‑69 (MQ=255)
      aCGACCGTCTTTGCGCGCTCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTc                                                       >  1:3013943/1‑69 (MQ=255)
      aCGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATgg                                                         >  1:389865/1‑67 (MQ=255)
      aCGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATgg                                                         >  1:2567180/1‑67 (MQ=255)
      aCGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATgg                                                         >  1:2036325/1‑67 (MQ=255)
      aCGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATgg                                                         >  1:2075280/1‑67 (MQ=255)
      aCGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATgg                                                         >  1:3266220/1‑67 (MQ=255)
      aCGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATgg                                                         >  1:2084707/1‑67 (MQ=255)
      aCGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATgg                                                         >  1:682233/1‑67 (MQ=255)
      aCGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATgg                                                         >  1:147570/1‑67 (MQ=255)
      aCGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATgg                                                         >  1:767037/1‑67 (MQ=255)
      aCGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATgg                                                         >  1:3027977/1‑67 (MQ=255)
      aCGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATgg                                                         >  1:2825488/1‑67 (MQ=255)
      aCGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATgg                                                         >  1:244457/1‑67 (MQ=255)
      aCGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATgg                                                         >  1:2527226/1‑67 (MQ=255)
      aCGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATgg                                                         >  1:250381/1‑67 (MQ=255)
      aCGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATg                                                          >  1:1679023/1‑66 (MQ=255)
      aCGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATg                                                          >  1:901444/1‑66 (MQ=255)
      aCGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGt                                                        <  1:72400/68‑1 (MQ=255)
      aCGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTc                                                       >  1:2153165/1‑69 (MQ=255)
      aCGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTc                                                       >  1:1932754/1‑69 (MQ=255)
      aCGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTc                                                       >  1:1972852/1‑69 (MQ=255)
      aCGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTc                                                       >  1:798359/1‑69 (MQ=255)
      aCGACAGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATgg                                                         >  1:905639/1‑67 (MQ=255)
           cGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGc                                                   >  1:3268596/1‑68 (MQ=255)
           cGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGc                                                   >  1:1817310/1‑68 (MQ=255)
           cGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGc                                                   >  1:2785742/1‑68 (MQ=255)
                   cgcgcCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGaacaac                                           <  1:431547/68‑1 (MQ=255)
                   cgcgcCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGaacaac                                           <  1:792119/68‑1 (MQ=255)
                             ccGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTgtcgt                                >  1:3430518/1‑69 (MQ=255)
                             ccGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTgtcgt                                >  1:2580908/1‑69 (MQ=255)
                             ccGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTgtcgt                                >  1:2681611/1‑69 (MQ=255)
                             ccGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTgtcgt                                >  1:2745468/1‑69 (MQ=255)
                             ccGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTgtcgt                                >  1:2487432/1‑69 (MQ=255)
                             ccGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTgtcgt                                >  1:3066284/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:820304/1‑69 (MQ=255)
                             ccGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTgtcgt                                >  1:128662/1‑69 (MQ=255)
                             ccGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTgtcgt                                >  1:1105758/1‑69 (MQ=255)
                             ccGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTgtcgt                                >  1:1056507/1‑69 (MQ=255)
                                tgAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                             <  1:959685/68‑1 (MQ=255)
                                ggAAACGATTGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                             <  1:2750607/69‑1 (MQ=255)
                                ggAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                             <  1:2970539/69‑1 (MQ=255)
                                ggAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                             <  1:37479/69‑1 (MQ=255)
                                ggAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                             <  1:1998217/69‑1 (MQ=255)
                                ggAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                             <  1:1709032/69‑1 (MQ=255)
                                ggAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                             <  1:548923/69‑1 (MQ=255)
                                ggAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                             <  1:1138316/69‑1 (MQ=255)
                                ggAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                             <  1:1411602/69‑1 (MQ=255)
                                ggAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                             <  1:785413/69‑1 (MQ=255)
                                ggAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                             <  1:2600779/69‑1 (MQ=255)
                                ggAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                             <  1:132948/69‑1 (MQ=255)
                                ggAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                             <  1:1317201/69‑1 (MQ=255)
                                ggAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                             <  1:2607046/69‑1 (MQ=255)
                                ggAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                             <  1:2723099/69‑1 (MQ=255)
                                 gAAACGATCGGCAACATCAGTAACACGATGGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                             <  1:625990/68‑1 (MQ=255)
                                 gAAACGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                             <  1:2668744/68‑1 (MQ=255)
                                  aaaGGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                           >  1:1918503/1‑69 (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:3021050/1‑68 (MQ=255)
                                  aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCttt                            >  1:2298165/1‑68 (MQ=255)
                                  aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCttt                            >  1:2246555/1‑68 (MQ=255)
                                  aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                             >  1:3100614/1‑67 (MQ=255)
                                  aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                           >  1:3432326/1‑69 (MQ=255)
                                  aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                           >  1:2856186/1‑69 (MQ=255)
                                  aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                           >  1:2887505/1‑69 (MQ=255)
                                  aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                           >  1:1368777/1‑69 (MQ=255)
                                  aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                           >  1:2918570/1‑69 (MQ=255)
                                  aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                           >  1:2534616/1‑69 (MQ=255)
                                  aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                           >  1:683833/1‑69 (MQ=255)
                                  aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                           >  1:1669331/1‑69 (MQ=255)
                                  aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                           >  1:62362/1‑69 (MQ=255)
                                  aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                           >  1:2212891/1‑69 (MQ=255)
                                  aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                           >  1:488718/1‑69 (MQ=255)
                                  aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                           >  1:1067789/1‑69 (MQ=255)
                                  aaaCGATCGGCAACATCAGTAACACGATCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCTTTg                           >  1:1968120/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:926776/1‑69 (MQ=255)
                                                           gATCTAGGTAATGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGAt  <  1:3254477/69‑1 (MQ=255)
                                                           gATCGAGGTAATGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGAt  <  1:1013334/69‑1 (MQ=255)
                                                           gATCGAGGTAATGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGAt  <  1:699228/69‑1 (MQ=255)
                                                           gATCGAGGTAATGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGAt  <  1:4487/69‑1 (MQ=255)
                                                           gATCGAGGTAATGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGAt  <  1:272239/69‑1 (MQ=255)
                                                            aTCGAGATAATGGTCATGCAGAACAACCAGATGGTCGTCtt                             <  1:2000381/41‑1 (MQ=38)
                                                                  |                                                             
GTAGTAACGACCGTCTTTGCGCGCCCGATCCGGGAAACGATCGGCAACATCAGTAACACGATCGAGGTAATGGTCATGCAGAACAACCAGATTGTCGTCTTTGGTCATCACCAAATCCTGTTCCAGAT  >  minE/1478255‑1478382

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: