Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,185,005 G→A 13.4% V324I (GTT→ATT) ‡ yeaG → conserved hypothetical protein with nucleoside triphosphate hydrolase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,185,0050GA13.4% 259.7 / 25.2 97V324I (GTT→ATT) ‡yeaGconserved hypothetical protein with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base G (28/56);  new base A (11/2);  total (39/58)
Fisher's exact test for biased strand distribution p-value = 6.45e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACACTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTT  >  minE/1184938‑1185073
                                                                   |                                                                    
acacTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGttt                                                                   <  1:4643252/71‑1 (MQ=255)
      gAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCg                                                             >  1:1649587/1‑71 (MQ=255)
         cGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                          <  1:899725/71‑1 (MQ=255)
         cGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                          <  1:218075/71‑1 (MQ=255)
         cGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                          <  1:5153118/71‑1 (MQ=255)
         cGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                          <  1:2610481/71‑1 (MQ=255)
         cGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                          <  1:3669232/71‑1 (MQ=255)
          gAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGaa                                                         <  1:4863717/71‑1 (MQ=255)
                   tGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGt                                                 >  1:841094/1‑70 (MQ=255)
                    gggTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGaa                                                         <  1:2938971/61‑1 (MQ=255)
                     tgTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTAtt                                              <  1:3752969/70‑1 (MQ=255)
                     ggTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTAtt                                              <  1:1239636/71‑1 (MQ=255)
                     ggTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTAtt                                              <  1:1176660/71‑1 (MQ=255)
                          cTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTAttgcttg                                         <  1:3882260/71‑1 (MQ=255)
                          cTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTAttgcttg                                         <  1:3220125/71‑1 (MQ=255)
                          cTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTAttgcttg                                         <  1:3491505/71‑1 (MQ=255)
                          cTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTAttgcttg                                         <  1:4523425/71‑1 (MQ=255)
                          cTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTAttgcttg                                         <  1:3346598/71‑1 (MQ=255)
                           tttCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                          <  1:2367878/53‑1 (MQ=255)
                           tttCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                          <  1:5108633/53‑1 (MQ=255)
                              ccGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCa                                     >  1:280797/1‑71 (MQ=255)
                                gTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                          <  1:4089734/48‑1 (MQ=255)
                                gTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                          <  1:919603/48‑1 (MQ=255)
                                gTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                          <  1:3250381/48‑1 (MQ=255)
                                gTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                          <  1:3869008/48‑1 (MQ=255)
                                gTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                          <  1:2128243/48‑1 (MQ=255)
                                gTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                          <  1:637328/48‑1 (MQ=255)
                                gTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                          <  1:2545427/48‑1 (MQ=255)
                                gTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                          <  1:2074062/48‑1 (MQ=255)
                                gTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                          <  1:4964384/48‑1 (MQ=255)
                                gTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                          <  1:4829017/48‑1 (MQ=255)
                                 taataaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                          <  1:4254376/47‑1 (MQ=255)
                                 taataaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                          <  1:4002593/47‑1 (MQ=255)
                                    taaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                               <  1:1500764/71‑1 (MQ=255)
                                    taaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                               <  1:1346692/71‑1 (MQ=255)
                                    taaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                               <  1:1347661/71‑1 (MQ=255)
                                    taaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                               <  1:2230309/71‑1 (MQ=255)
                                    taaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                               <  1:2328440/71‑1 (MQ=255)
                                    taaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                               <  1:4966717/71‑1 (MQ=255)
                                    taaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                               <  1:4328090/71‑1 (MQ=255)
                                     aaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTgcgc                                      <  1:1580750/63‑1 (MQ=255)
                                     aaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTgcgc                                      <  1:4457892/63‑1 (MQ=255)
                                     aaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTgcgc                                      <  1:3078019/63‑1 (MQ=255)
                                     aaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTgcgc                                      <  1:3312983/63‑1 (MQ=255)
                                      aCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                               >  1:2724458/1‑69 (MQ=255)
                                      aCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                               >  1:308727/1‑69 (MQ=255)
                                      aCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                               >  1:1417489/1‑69 (MQ=255)
                                      aCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                               >  1:1273432/1‑69 (MQ=255)
                                      aCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                               >  1:106081/1‑69 (MQ=255)
                                      aCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGACGTATTGCTTGCGCATTTCCg                               >  1:4213012/1‑69 (MQ=255)
                                       cAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGt                                                 >  1:3336316/1‑50 (MQ=255)
                                       cAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTcag                               >  1:3904249/1‑66 (MQ=255)
                                       cAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTc                                 >  1:4035432/1‑66 (MQ=255)
                                       cAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                               >  1:1907883/1‑68 (MQ=255)
                                       cAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                               >  1:2544656/1‑68 (MQ=255)
                                       cAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                               >  1:3690936/1‑68 (MQ=255)
                                       cAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                               >  1:4716675/1‑68 (MQ=255)
                                       cAAAAACAAAGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                               >  1:4307666/1‑68 (MQ=255)
                                        aaaaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                             >  1:4189999/1‑69 (MQ=255)
                                        aaaaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                             >  1:2582564/1‑69 (MQ=255)
                                         aaaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCa                                     >  1:3467636/1‑60 (MQ=255)
                                         aaaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                             >  1:4838604/1‑68 (MQ=255)
                                         aaaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                             >  1:4081742/1‑68 (MQ=255)
                                         aaaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                             >  1:3909642/1‑68 (MQ=255)
                                         aaaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                             >  1:2166856/1‑68 (MQ=255)
                                         aaaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                             >  1:1947465/1‑68 (MQ=255)
                                         aaaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                             >  1:3188188/1‑68 (MQ=255)
                                          aaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTg                                                           <  1:2529728/37‑1 (MQ=255)
                                          aaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTg                                                           <  1:3783730/37‑1 (MQ=255)
                                          aaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTg                                                           <  1:1701134/37‑1 (MQ=255)
                                          aaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                             >  1:1081638/1‑67 (MQ=255)
                                             caacGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATc                       >  1:2834938/1‑70 (MQ=255)
                                               acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTAt                                               >  1:3568936/1‑44 (MQ=16)
                                               acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTgcg                                       >  1:3466443/1‑52 (MQ=19)
                                               acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGATGAAATc                       >  1:408721/1‑68 (MQ=255)
                                               acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGATGAAATc                       >  1:1061672/1‑68 (MQ=255)
                                               acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGATGAAATc                       >  1:1582473/1‑68 (MQ=255)
                                               acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGATGAAATc                       >  1:658818/1‑68 (MQ=255)
                                               acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGATGAAATc                       >  1:531556/1‑68 (MQ=255)
                                               acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGATGAAATc                       >  1:3615779/1‑68 (MQ=255)
                                               acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGATGAAATc                       >  1:3124410/1‑68 (MQ=255)
                                               acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGATGAAATc                       >  1:2730294/1‑68 (MQ=255)
                                               acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGATGAAATc                       >  1:4423473/1‑68 (MQ=255)
                                               acGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCaaa                    <  1:289350/71‑1 (MQ=255)
                                                 gAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAAt                  <  1:2536838/71‑1 (MQ=255)
                                                   aGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCt                <  1:4073208/71‑1 (MQ=255)
                                                       ttCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGAtga                           <  1:4543288/56‑1 (MQ=255)
                                                         ccTCGATCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGAtga                           <  1:1753935/54‑1 (MQ=255)
                                                          cTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACgaga          <  1:4069448/70‑1 (MQ=255)
                                                          cTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACgaga          <  1:3221779/70‑1 (MQ=255)
                                                          cTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACgaga          <  1:2976790/70‑1 (MQ=255)
                                                          cTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGaa         <  1:736434/71‑1 (MQ=255)
                                                           tCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACgaga          <  1:3825626/69‑1 (MQ=255)
                                                           tCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACgaga          <  1:2323156/69‑1 (MQ=255)
                                                           tCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACgaga          <  1:2563254/69‑1 (MQ=255)
                                                           tCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGaa         <  1:2424165/70‑1 (MQ=255)
                                                            cGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGaa         <  1:3375686/69‑1 (MQ=255)
                                                                 gtgtTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCtt  <  1:3686387/71‑1 (MQ=255)
                                                                 gtgtTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCtt  <  1:1779497/71‑1 (MQ=255)
                                                                   |                                                                    
ACACTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTT  >  minE/1184938‑1185073

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: