Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,185,040 T→G 8.4% I335M (ATT→ATG) ‡ 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,0400TG 8.4% 333.1 / 25.2 131I335M (ATT→ATG) ‡yeaGconserved hypothetical protein with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base T (50/70);  new base G (9/2);  total (59/72)
Fisher's exact test for biased strand distribution p-value = 2.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 141 total aligned reads displayed.

TAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTGAATTGACTCACGCCCCATGCGCCCCTG  >  minE/1184974‑1185110
                                                                  |                                                                      
taaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                                                                    <  1:4328090/71‑1 (MQ=255)
taaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                                                                    <  1:2230309/71‑1 (MQ=255)
taaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                                                                    <  1:4966717/71‑1 (MQ=255)
taaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                                                                    <  1:1500764/71‑1 (MQ=255)
taaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                                                                    <  1:1347661/71‑1 (MQ=255)
  aCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                                                                    >  1:2724458/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)
  aCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                                                                    >  1:308727/1‑69 (MQ=255)
  aCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGACGTATTGCTTGCGCATTTCCg                                                                    >  1:4213012/1‑69 (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)
     aaaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                                                                  >  1:3909642/1‑68 (MQ=255)
     aaaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                                                                  >  1:1947465/1‑68 (MQ=255)
     aaaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                                                                  >  1:2166856/1‑68 (MQ=255)
     aaaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                                                                  >  1:4838604/1‑68 (MQ=255)
     aaaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                                                                  >  1:3188188/1‑68 (MQ=255)
      aaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                                                                  >  1:1081638/1‑67 (MQ=255)
         caacGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATc                                                            >  1:2834938/1‑70 (MQ=255)
           acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGATGAAATc                                                            >  1:408721/1‑68 (MQ=255)
           acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGATGAAATc                                                            >  1:3615779/1‑68 (MQ=255)
           acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGATGAAATc                                                            >  1:658818/1‑68 (MQ=255)
           acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGATGAAATc                                                            >  1:4423473/1‑68 (MQ=255)
           acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGATGAAATc                                                            >  1:1061672/1‑68 (MQ=255)
           acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGATGAAATc                                                            >  1:1582473/1‑68 (MQ=255)
           acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGATGAAATc                                                            >  1:531556/1‑68 (MQ=255)
               aGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCt                                                     <  1:4073208/71‑1 (MQ=255)
                   ttCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGAtga                                                                <  1:4543288/56‑1 (MQ=255)
                      cTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACgaga                                               <  1:2976790/70‑1 (MQ=255)
                      cTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACgaga                                               <  1:4069448/70‑1 (MQ=255)
                      cTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGaa                                              <  1:736434/71‑1 (MQ=255)
                       tCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACgaga                                               <  1:2323156/69‑1 (MQ=255)
                       tCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACgaga                                               <  1:2563254/69‑1 (MQ=255)
                       tCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACgaga                                               <  1:3825626/69‑1 (MQ=255)
                       tCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGaa                                              <  1:2424165/70‑1 (MQ=255)
                        cGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGaa                                              <  1:3375686/69‑1 (MQ=255)
                             gtgtTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCtt                                       <  1:1779497/71‑1 (MQ=255)
                             gtgtTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCtt                                       <  1:3686387/71‑1 (MQ=255)
                                tttACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCaaa                                                         <  1:1879344/50‑1 (MQ=255)
                                tttACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCaaa                                                         <  1:3965795/50‑1 (MQ=255)
                                tttACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCaaa                                                         <  1:4740196/50‑1 (MQ=255)
                                tttACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCaaa                                                         <  1:2447266/50‑1 (MQ=255)
                                tttACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCaaa                                                         <  1:1884815/50‑1 (MQ=255)
                                tttACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAAtt                                           <  1:2767203/64‑1 (MQ=255)
                                tttACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAAt                                    >  1:3043214/1‑71 (MQ=255)
                                  tACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGc                                         >  1:1378579/1‑64 (MQ=255)
                                  tACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGc                                         >  1:127299/1‑64 (MQ=255)
                                  tACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGc                                         >  1:2727891/1‑64 (MQ=255)
                                  tACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGc                                         >  1:321991/1‑64 (MQ=255)
                                  tACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGc                                         >  1:4222585/1‑64 (MQ=255)
                                  tACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATca                                  >  1:728910/1‑71 (MQ=255)
                                  tACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATca                                  >  1:874782/1‑71 (MQ=255)
                                   aCATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                                                                  >  1:2762159/1‑38 (MQ=255)
                                    cATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAACCTACGAGAAATTGc                                         >  1:4261839/1‑62 (MQ=255)
                                        gTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGGAATTGCTTAATCACAGTg                             <  1:2356672/70‑1 (MQ=255)
                                        gTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATca                                  <  1:5196985/65‑1 (MQ=255)
                                        gTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATca                                  <  1:2437766/65‑1 (MQ=255)
                                        gTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTg                             <  1:1840053/70‑1 (MQ=255)
                                        gTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTg                             <  1:3073442/70‑1 (MQ=255)
                                        gTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTg                             <  1:1781544/70‑1 (MQ=255)
                                         tGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTGaa                           <  1:430561/71‑1 (MQ=255)
                                         tGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTGaa                           <  1:2881223/71‑1 (MQ=255)
                                            aGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCtt                                       <  1:1816405/56‑1 (MQ=255)
                                            aGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCtt                                       <  1:5008081/56‑1 (MQ=255)
                                            aGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCtt                                       <  1:174743/56‑1 (MQ=255)
                                            aGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCtt                                       <  1:2472414/56‑1 (MQ=255)
                                            aGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTGAATTg                        <  1:1719619/71‑1 (MQ=255)
                                            aGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTGAATTg                        <  1:4285152/71‑1 (MQ=255)
                                            aGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTGAATTg                        <  1:5046679/71‑1 (MQ=255)
                                            aGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTGAATTg                        <  1:1040723/71‑1 (MQ=255)
                                             ggTGCCGTATTGCTTGCGCATTTCCGAAGAGATCa                                                           <  1:3574465/35‑1 (MQ=255)
                                             ggTGCCGTATTGCTTGCGCATTTCCGAAGAGATCa                                                           <  1:4930530/35‑1 (MQ=255)
                                             ggTGCCGTATTGCTTGCGCATTTCCGAAGAGATCa                                                           <  1:1076060/35‑1 (MQ=255)
                                             ggTGCCGTATTGCTTGCGCATTTCCGAAGAGATCa                                                           <  1:341536/35‑1 (MQ=255)
                                             ggTGCCGTATTGCTTGCGCATTTCCGAAGAGATCa                                                           <  1:2553684/35‑1 (MQ=255)
                                                      ttgcttgcGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTGAATTGACTCACGccc              >  1:2527608/1‑71 (MQ=255)
                                                      ttgcttgcGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTGAATTGACTCACGccc              >  1:1520917/1‑71 (MQ=255)
                                                        gcttgcGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTGaa                           >  1:2820106/1‑56 (MQ=255)
                                                        gcttgcGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTGaa                           >  1:625134/1‑56 (MQ=255)
                                                        gcttgcGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTGaa                           >  1:5121104/1‑56 (MQ=255)
                                                        gcttgcGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTGAATTGACTCACGCCCCa            >  1:4295035/1‑71 (MQ=255)
                                                          ttgcGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTg                             <  1:5069451/52‑1 (MQ=255)
                                                           tgcGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTGAAtt                         <  1:3288585/55‑1 (MQ=255)
                                                           tgcGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGAGAAtt                         <  1:148272/55‑1 (MQ=255)
                                                              gcATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATca                                  <  1:850367/43‑1 (MQ=255)
                                                              gcATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTGaa                           <  1:4564777/50‑1 (MQ=255)
                                                                 tttCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTGAATTGACTCACGCCCCATGCGCCCCt   <  1:3419402/71‑1 (MQ=255)
                                                                 tttCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTGAATTGACTCACGCCCCATGCGCCCCt   <  1:4211307/71‑1 (MQ=255)
                                                                 tttCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTGAATTGACTCACGCCCCATGCGCCCCt   <  1:2337622/71‑1 (MQ=255)
                                                                  ttCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTGAATTGACTCACGCCCCATGCGCCCCTg  >  1:3367641/1‑71 (MQ=255)
                                                                  ttCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTGAATTGACTCACGCCCCATGCGCCCCTg  >  1:2926631/1‑71 (MQ=255)
                                                                  ttCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTGAATTGACTCACGCCCCATGCGCCCCTg  >  1:1671366/1‑71 (MQ=255)
                                                                  ttCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTGAATTGACTCACGCCCCATGCGCCCCTg  >  1:1202735/1‑71 (MQ=255)
                                                                  ttCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTGAATTGACTCACGCCCCATGCGCCCCTg  >  1:1983160/1‑71 (MQ=255)
                                                                  |                                                                      
TAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCACAGTGAATTGACTCACGCCCCATGCGCCCCTG  >  minE/1184974‑1185110

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: