Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,639,777 C→A 9.0% T234T (ACC→ACA ugpQ → glycerophosphodiester phosphodiesterase, cytosolic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,639,7770CA 9.0% 289.1 / 8.6 88T234T (ACC→ACAugpQglycerophosphodiester phosphodiesterase, cytosolic
Reads supporting (aligned to +/- strand):  ref base C (40/40);  new base A (6/2);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

GTTTATACCGTCAATAAACCCCAGCGCGCAGCAGAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGC  >  minE/2639709‑2639845
                                                                    |                                                                    
gTTTATACCGTCAATAAACCCCAGCGCGCAGCAGAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCg                                                                     <  1:3151088/70‑1 (MQ=255)
gTTTATACCGTCAATAAACCCCAGCGCGCAGCAGAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGa                                                                    >  1:1893877/1‑71 (MQ=255)
gTTTATACCGTCAATAAACCCCAGCGCGCAGCAGAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGa                                                                    >  1:2054930/1‑71 (MQ=255)
  ttATACCGTCAATAAACCCCAGCGCGCAGCAGAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGAt                                                                   >  1:1424794/1‑70 (MQ=255)
   tataCCGTCAATAAACCCCAGCGCGCAGCAGAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATg                                                                  <  1:1763242/70‑1 (MQ=255)
     taCCGTCAATAAACCCCAGCGCTCAGCAGACTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGa                                                               <  1:2156367/71‑1 (MQ=255)
          tCAATAAACCCCAGCGCGCAGCAGAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGAc                                                          >  1:956289/1‑71 (MQ=255)
          tCAATAAACCCCAGCGCGCAGCAGAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGAc                                                          >  1:1529102/1‑71 (MQ=255)
          tCAATAAACCCCAGCGCGCAGCAGAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGAc                                                          >  1:2152719/1‑71 (MQ=255)
           cAATAAACCCCAGCGCGCAGCAGGGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACg                                                         >  1:2408062/1‑71 (MQ=255)
             aTAAACCCCAGCGCGCAGCAGAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTg                                                       <  1:1719584/71‑1 (MQ=255)
                aaCCCCAGCGCGCAGCAGAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGAtt                                                    >  1:731269/1‑71 (MQ=255)
                 aCCCCAGCGCGCAGCAGAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTg                                                   <  1:38378/71‑1 (MQ=255)
                         gcgcAGCAGAGTTGCTGCGCTGTGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAAc                                           >  1:2358279/1‑71 (MQ=255)
                         gcgcAGCAGAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAAc                                           <  1:1525095/71‑1 (MQ=255)
                            cagcagAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACtt                                         <  1:996165/70‑1 (MQ=255)
                             agcagAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTa                                       <  1:1141894/71‑1 (MQ=255)
                             agcagAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTa                                       <  1:3342294/71‑1 (MQ=255)
                             agcagAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTa                                       >  1:1934246/1‑71 (MQ=255)
                             agcagAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTa                                       <  1:1089598/71‑1 (MQ=255)
                             agcagAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTa                                       >  1:181197/1‑71 (MQ=255)
                             agcagAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTa                                       >  1:1257570/1‑71 (MQ=255)
                             agcagAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTa                                       <  1:1462127/71‑1 (MQ=255)
                             agcagAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTa                                       <  1:1594859/71‑1 (MQ=255)
                             agcagAGTTGCTGCGCTGGGGTGTGGATCGCATCGGTACAGATg                                                                  >  1:59702/1‑44 (MQ=25)
                             agcagAGTTGCTGCGCTGGGGTGTGGATCGCATCGGTACAGATg                                                                  >  1:2491028/1‑44 (MQ=25)
                             agcagAGTTGCTGCGCTGGGGTGTGGATCGCATCGGTACAGATg                                                                  >  1:2377748/1‑44 (MQ=25)
                             agcagAGTTGCTGCGCTGGGGTGTGGATCGCATCGGTACAGATg                                                                  >  1:2205999/1‑44 (MQ=25)
                              gcagAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGt                                                 <  1:2762977/60‑1 (MQ=255)
                              gcagAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGt                                                 <  1:343533/60‑1 (MQ=255)
                              gcagAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGa                                             <  1:3591616/64‑1 (MQ=255)
                               cagAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGa                                                                    >  1:2541757/1‑40 (MQ=255)
                               cagAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGa                                                                    >  1:259333/1‑40 (MQ=255)
                               cagAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGa                                                                    >  1:1305446/1‑40 (MQ=255)
                               cagAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGa                                                                    >  1:3475367/1‑40 (MQ=255)
                               cagAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGa                                                                    >  1:2061878/1‑40 (MQ=255)
                               cagAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACg                                     <  1:3598848/71‑1 (MQ=255)
                                    ttGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCa                                <  1:3610186/71‑1 (MQ=255)
                                        tgcGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGAc                                                          <  1:684557/41‑1 (MQ=255)
                                        tgcGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCg                                              <  1:3249729/53‑1 (MQ=255)
                                        tgcGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCg                                              <  1:2946055/53‑1 (MQ=255)
                                        tgcGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCg                                              <  1:446929/53‑1 (MQ=255)
                                        tgcGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGa                                             <  1:3107011/54‑1 (MQ=255)
                                        tgcGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGa                                             <  1:3044968/54‑1 (MQ=255)
                                        tgcGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGa                                             <  1:1745254/54‑1 (MQ=255)
                                        tgcGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAg                            >  1:1211886/1‑71 (MQ=255)
                                         gcgcTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACg                                     >  1:437439/1‑61 (MQ=255)
                                         gcgcTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGt                           >  1:3374671/1‑71 (MQ=255)
                                         gcgcTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGt                           >  1:2922585/1‑71 (MQ=255)
                                          cgcTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCTCAATAg                            <  1:4002/69‑1 (MQ=255)
                                           gcTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAAc                                           >  1:1911871/1‑53 (MQ=255)
                                           gcTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGttt                         >  1:2135097/1‑71 (MQ=255)
                                               gggTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTc                       <  1:16025/69‑1 (MQ=255)
                                               gggTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCaa                     <  1:995802/71‑1 (MQ=255)
                                                ggTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGtt                          >  1:3635684/1‑65 (MQ=255)
                                                  tgtgGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACg                                     <  1:68800/52‑1 (MQ=255)
                                                  tgtgGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACg                                     <  1:827376/52‑1 (MQ=255)
                                                  tgtgGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACg                                     <  1:2755108/52‑1 (MQ=255)
                                                  tgtgGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACg                                     <  1:2739948/52‑1 (MQ=255)
                                                  tgtgGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACg                                     <  1:1901731/52‑1 (MQ=255)
                                                      gATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATg              >  1:592000/1‑71 (MQ=255)
                                                      gATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATg              >  1:386677/1‑71 (MQ=255)
                                                      gATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATg              >  1:3323047/1‑71 (MQ=255)
                                                      gATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATg              >  1:2627419/1‑71 (MQ=255)
                                                       aTTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGt             <  1:188041/71‑1 (MQ=255)
                                                       aTTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGt             <  1:2496893/71‑1 (MQ=255)
                                                        ttGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGt             >  1:716917/1‑70 (MQ=255)
                                                        ttGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGt             >  1:2217813/1‑70 (MQ=255)
                                                        tcgCATCGGTACAGATGCGATTGACGTGATTGGTCCGAACtt                                         <  1:1240878/40‑1 (MQ=25)
                                                        tcgCATCGGTACAGATGCGATTGACGTGATTGGTCCGAACtt                                         <  1:2179467/40‑1 (MQ=25)
                                                        tcgCATCGGTACAGATGCAATTGACGTGATTGGTCCGAACTTTACGGCCCa                                >  1:3424455/3‑51 (MQ=25)
                                                        tcgCATCGGTACAGATGCAATTGACGTGATTGGTCCGAACTTTACGGCCCa                                >  1:1341274/3‑51 (MQ=25)
                                                          gCATCTGTACCGATGCGTTTGACGTGATTGGTCCGAACTTTACGGCCCAATTGttt                         >  1:3105722/1‑56 (MQ=255)
                                                          gCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCa           <  1:2011086/70‑1 (MQ=255)
                                                          gCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCa           <  1:881114/70‑1 (MQ=255)
                                                          gCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCa           <  1:3257172/70‑1 (MQ=255)
                                                          gCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCa           >  1:2015295/1‑70 (MQ=255)
                                                          gCATCTGTACCGATGCGATTGACGTGATTGGGCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCa           <  1:2870965/70‑1 (MQ=255)
                                                            aTCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACgg                                    <  1:2495285/43‑1 (MQ=255)
                                                              cTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCa                                >  1:3182988/1‑45 (MQ=255)
                                                              cTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCa                                >  1:2231048/1‑45 (MQ=255)
                                                              cTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCa                                >  1:3024099/1‑45 (MQ=255)
                                                              cTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCa                                >  1:2483508/1‑45 (MQ=255)
                                                              cTGTACCGATGCGATTGACGTAATTGGTCCGAACTTTACGGCCCa                                >  1:2418452/1‑45 (MQ=255)
                                                                gTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTgcg    >  1:1307415/1‑71 (MQ=255)
                                                                gTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTgcg    >  1:1936529/1‑71 (MQ=255)
                                                                 tACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTgcgg   <  1:2599914/71‑1 (MQ=255)
                                                                  aCCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTgcggc  >  1:581557/1‑71 (MQ=255)
                                                                   ccGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTg      <  1:3411868/66‑1 (MQ=255)
                                                                    |                                                                    
GTTTATACCGTCAATAAACCCCAGCGCGCAGCAGAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGC  >  minE/2639709‑2639845

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: