Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,639,808 A→G 25.0% T245A (ACG→GCG) 
*147Q (TAA→CAA) 
ugpQ →
yhhA ←
glycerophosphodiester phosphodiesterase, cytosolic
conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,639,8080AG25.0% 155.6 / 56.5 84T245A (ACG→GCG) 
*147Q (TAA→CAA) 
ugpQ
yhhA
glycerophosphodiester phosphodiesterase, cytosolic
conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (41/22);  new base G (10/11);  total (51/33)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

CATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCAGCATGTCGCCGTT  >  minE/2639768‑2639874
                                        |                                                                  
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAg                                       <  1:796794/70‑1 (MQ=255)
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAg                                       <  1:787841/70‑1 (MQ=255)
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAg                                       <  1:670423/70‑1 (MQ=255)
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAg                                       <  1:1804245/70‑1 (MQ=255)
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAg                                       <  1:2674662/70‑1 (MQ=255)
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAgt                                      <  1:434703/71‑2 (MQ=255)
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAgg                                      <  1:2887985/71‑1 (MQ=255)
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAgg                                      <  1:1642300/71‑1 (MQ=255)
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAgg                                      <  1:186651/71‑1 (MQ=255)
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAgg                                      <  1:2017018/71‑1 (MQ=255)
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAgg                                      <  1:208388/71‑1 (MQ=255)
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAgg                                      <  1:2218422/71‑1 (MQ=255)
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAgg                                      <  1:2323507/71‑1 (MQ=255)
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAgg                                      <  1:2729773/71‑1 (MQ=255)
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAgg                                      <  1:2309411/71‑1 (MQ=255)
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAgg                                      <  1:378922/71‑1 (MQ=255)
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAgg                                      <  1:461468/71‑1 (MQ=255)
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAgg                                      <  1:47930/71‑1 (MQ=255)
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAgg                                      <  1:608450/71‑1 (MQ=255)
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAgg                                      <  1:987873/71‑1 (MQ=255)
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAgg                                      <  1:1312982/71‑1 (MQ=255)
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAgg                                      <  1:1587334/71‑1 (MQ=255)
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAgg                                      <  1:1553695/71‑1 (MQ=255)
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAgg                                      <  1:1401519/71‑1 (MQ=255)
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAgg                                      <  1:1218266/71‑1 (MQ=255)
cATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAgg                                      <  1:1196843/71‑1 (MQ=255)
    tGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAg                                                         >  1:311867/1‑48 (MQ=255)
    tGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAg                                                         >  1:1026200/1‑48 (MQ=255)
    tGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAg                                                         >  1:622766/1‑48 (MQ=255)
    tGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAg                                                         >  1:1062983/1‑48 (MQ=255)
             gCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAgg                                      <  1:2452711/58‑1 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGttt                                                      >  1:1983461/1‑37 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGt                             >  1:1768194/1‑62 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:2739209/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:1376482/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:1370057/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:2686989/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:1625752/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:2892287/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:265842/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:299912/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:1328770/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:494632/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:634459/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:674300/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:744040/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:1131545/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:850014/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:985445/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:1060313/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:2062688/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:1607285/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:1721900/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:1579664/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:1871682/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:1924792/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:194682/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:1507703/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:2218968/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:2271125/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:2368560/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:2488405/1‑70 (MQ=255)
                  tGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                     >  1:1507353/1‑70 (MQ=255)
                   gACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAAt                                            >  1:2670118/1‑46 (MQ=255)
                   gACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGa                   >  1:287945/1‑71 (MQ=255)
                   gACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGa                   >  1:1292465/1‑71 (MQ=255)
                   gACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGa                   >  1:974000/1‑71 (MQ=255)
                   gACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGa                   >  1:1969385/1‑71 (MQ=255)
                   gACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGa                   >  1:1934040/1‑71 (MQ=255)
                   gACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGa                   >  1:1570505/1‑71 (MQ=255)
                   gACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCAGCGGCGTGCTGGGCTg                    >  1:1999564/1‑70 (MQ=255)
                   gACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAAAGGAATGTCAGGCTGCGGCGTGCTGGGCTGa                   >  1:2400356/1‑71 (MQ=255)
                    aCGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGa                   <  1:1909322/70‑1 (MQ=255)
                    aCGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGa                   <  1:860150/70‑1 (MQ=255)
                    aCGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGa                   <  1:831282/70‑1 (MQ=255)
                    aCGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGa                   <  1:1366699/70‑1 (MQ=255)
                    aCGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGa                   <  1:234508/70‑1 (MQ=255)
                                 gAACTTTGCGGCCCAATAGTTTTCAACGGAATGGCAGGCTGCGGCGTGCTGGGCTGa                   <  1:2070592/57‑1 (MQ=255)
                                    cTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCAGCATGTCGCCGtt  >  1:799504/1‑71 (MQ=255)
                                    cTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCAGCATGTCGCCGtt  >  1:2748740/1‑71 (MQ=255)
                                    cTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCAGCATGTCGCCGtt  >  1:2404268/1‑71 (MQ=255)
                                    cTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCAGCATGTCGCCGtt  >  1:979592/1‑71 (MQ=255)
                                    cTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCAGCATGTCGCCGtt  >  1:1553420/1‑71 (MQ=255)
                                    cTTTGCGGCCCAATAGTTATCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCAGCATGTCGCCGtt  >  1:2612757/1‑71 (MQ=255)
                                        |                                                                  
CATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCAGCATGTCGCCGTT  >  minE/2639768‑2639874

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: