Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 178,645 G→T 23.7% G631G (GGG→GGT yaeT → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE178,6450GT23.7% 240.3 / 54.1 134G631G (GGG→GGTyaeTconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (19/83);  new base T (8/24);  total (27/107)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 134 total aligned reads displayed.

AAGTGACGTTAGACACGGCGACTTATGTGCCGATCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAG  >  minE/178581‑178702
                                                                |                                                         
aaGTGACGTTAGACACGGCGACTTATGTGCCGATCGATGACGATCACAAATGGGTTGTTCTGGGGCGTa                                                       >  1:2919323/1‑69 (MQ=255)
aaGTGACGTTAGACACGGCGACTTATGTGCCGATCGATGACGATCACAAATGGGTTGTTCTGGGGCGTa                                                       >  1:2077874/1‑69 (MQ=255)
aaGTGACGTTAGACACGGCGACTTATGTGCCGATCGATGACGATCACAAATGGGTTGTTCTGGGGCGTa                                                       >  1:67409/1‑69 (MQ=255)
                 gCGACTTATGTGCCTATCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGt                                      <  1:2495498/69‑1 (MQ=255)
                 gCGACTTATGTGCCGATCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGt                                      <  1:9537/69‑1 (MQ=255)
                 gCGACTTATGTGCCGATCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGt                                      >  1:255002/1‑69 (MQ=255)
                 gCGACTTATGTGCCGATCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGt                                      <  1:2539873/69‑1 (MQ=255)
                 gCGACTTATGTGCCGATCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGt                                      >  1:3150315/1‑69 (MQ=255)
                 gCGACTTATGTGCCGATCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGt                                      <  1:2239247/69‑1 (MQ=255)
                 gCGACTTATGTGCCGATCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGt                                      <  1:2161486/69‑1 (MQ=255)
                 gCGACTTATGTGCCGATCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGt                                      <  1:2157276/69‑1 (MQ=255)
                 gCGACTTATGTGCCGATCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGt                                      <  1:1733704/69‑1 (MQ=255)
                 gCGACTTATGTGCCGATCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGt                                      <  1:1645186/69‑1 (MQ=255)
                 gCGACTTATGTGCCGATCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGt                                      <  1:1543441/69‑1 (MQ=255)
                 gCGACTTATGTGCCGATCGATGACGATCACAAATGGGTAGTTCTGGGGCGTACCCGCTGGGGTTATGGt                                      <  1:2847618/69‑1 (MQ=255)
                      ttATGTGCCGATCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATgg                                 <  1:1701659/69‑1 (MQ=255)
                      ttATGTGCCGATCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATgg                                 <  1:753920/69‑1 (MQ=255)
                      ttATGTGCCGATCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATgg                                 <  1:2974172/69‑1 (MQ=255)
                           tgCCGATCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTa                             <  1:2378432/68‑1 (MQ=255)
                           tgCCGATCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTa                             <  1:2042730/68‑1 (MQ=255)
                             ccGATTGATGATGATCACAAATGGGTTGTTCTGGGTCGTACCCGCTGGGGTTATGGTGATGGTTTAgg                           <  1:474953/68‑1 (MQ=255)
                             ccGATTGATGATGATCACAAATGGGTTGTTCTGGGTCGTACCCGCTGGGGTTATGGTGATGGTTTAgg                           <  1:2232647/68‑1 (MQ=255)
                             ccGATTGATGATGATCACAAATGGGTTGTTCTGGGTCGTACCCGCTGGGGTTATGGTGATGGTTTAgg                           <  1:606236/68‑1 (MQ=255)
                             ccGATTGATGATGATCACAAATGGGTTGTTCTGGGTCGTACCCGCTGGGGTTATGGTGATGGTTTAgg                           <  1:808309/68‑1 (MQ=255)
                             ccGATTGATGATGATCACAAATGGGTTGTTCTGGGTCGTACCCGCTGGGGTTATGGTGATGGTTTAgg                           <  1:3222169/68‑1 (MQ=255)
                             ccGATTGATGATGATCACAAATGGGTTGTTCTGGGTCGTACCCGCTGGGGTTATGGTGATGGTTTAgg                           <  1:3258513/68‑1 (MQ=255)
                             ccGATTGATGATGATCACAAATGGGTTGTTCTGGGTCGTACCCGCTGGGGTTATGGTGATGGTTTAgg                           <  1:332408/68‑1 (MQ=255)
                             ccGATTGATGATGATCACAAATGGGTTGTTCTGGGTCGTACCCGCTGGGGTTATGGTGATGGTTTAgg                           <  1:3154290/68‑1 (MQ=255)
                             ccGATTGATGATGATCACAAATGGGTTGTTCTGGGTCGTACCCGCTGGGGTTATGGTGATGGTTTAgg                           <  1:1997107/68‑1 (MQ=255)
                             ccGATTGATGATGATCACAAATGGGTTGTTCTGGGTCGTACCCGCTGGGGTTATGGTGATGGTTTAgg                           <  1:369145/68‑1 (MQ=255)
                             ccGATTGATGATGATCACAAATGGGTTGTTCTGGGTCGTACCCGCTGGGGTTATGGTGATGGTTTAgg                           <  1:1832578/68‑1 (MQ=255)
                             ccGATTGATGATGATCACAAATGGGTTGTTCTGGGTCGTACCCGCTGGGGTTATGGTGATGGTTTAgg                           <  1:29357/68‑1 (MQ=255)
                             ccGATTGATGATGATCACAAATGGGTTGTTCTGGGTCGTACCCGCTGGGGTTATGGTGATGGTTTAgg                           <  1:487427/68‑1 (MQ=255)
                             ccGATTGATGATGATCACAAATGGGTTGTTCTGGGTCGTACCCGCTGGGGTTATGGTGATGGTTTAgg                           <  1:2997263/68‑1 (MQ=255)
                             ccGATTGATGATGATCACAAATGGGTTGTTCTGGGTCGTACCCGCTGGGGTTATGGTGATGGTTTAgg                           <  1:623143/68‑1 (MQ=255)
                             ccGATTGATGATGATCACAAATGGGTTGTTCTGGGTCGTACCCGCTGGGGTTATGGTGATGGTTTAgg                           <  1:2850389/68‑1 (MQ=255)
                             ccGATTGATGATGATCACAAATGGGTTGTTCTGGGTCGTACCCGCTGGGGTTATGGTGATGGTTTAgg                           <  1:1618546/68‑1 (MQ=255)
                             ccGATTGATGATGATCACAAATGGGTTGTTCTGGGTCGTACCCGCTGGGGTTATGGTGATGGTTTAgg                           <  1:2695172/68‑1 (MQ=255)
                             ccGATTGATGATGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAgg                           <  1:2671299/68‑1 (MQ=255)
                               gattgatGATGATCACTAATGGGTTGTTCTGGGTCGTACCCGCTGGGGTTATGGTGATGGTTTAgg                           <  1:1896408/66‑1 (MQ=255)
                                 tcgatGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCa                      <  1:1331369/69‑1 (MQ=255)
                                 tcgatGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCa                      <  1:3109427/69‑1 (MQ=255)
                                 tcgatGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCa                      <  1:1756271/69‑1 (MQ=255)
                                 tcgatGACGATCACAAATGGGGTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCa                      <  1:2530098/69‑1 (MQ=255)
                                  tgatgatgatCACAAATGGGTTGTTCTGGGTCGTACCCGCTGGGGTTATGGTGATGGTTTAggcgg                        >  1:2342827/2‑66 (MQ=255)
                                  tgatgatgatCACAAATGGGTTGTTCTGGGTCGTACCCGCTGGGGTTATGGTGATGGTTTAggcgg                        >  1:489878/2‑66 (MQ=255)
                                  tgatgatgatCACAAATGGGTTGTTCTGGGTCGTACCCGCTGGGGTTATGGTGATGGTTTAggcgg                        >  1:2526944/2‑66 (MQ=255)
                                  tgatgatgatCACAAATGGGTTGTTCTGGGTCGTACCCGCTGGGGTTATGGTGATGGTTTAggcgg                        >  1:1157530/2‑66 (MQ=255)
                                  cgatGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAggcgg                        >  1:1781247/1‑66 (MQ=255)
                                     tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:2636644/69‑1 (MQ=255)
                                     tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:3193691/69‑1 (MQ=255)
                                     tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:3076345/69‑1 (MQ=255)
                                     tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:942016/69‑1 (MQ=255)
                                     tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:3089500/69‑1 (MQ=255)
                                     tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:930860/69‑1 (MQ=255)
                                     tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:923078/69‑1 (MQ=255)
                                     tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:858145/69‑1 (MQ=255)
                                     tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:507687/69‑1 (MQ=255)
                                     tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:521178/69‑1 (MQ=255)
                                     tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:1622337/69‑1 (MQ=255)
                                     tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:1269899/69‑1 (MQ=255)
                                     tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:2885678/69‑1 (MQ=255)
                                     tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:2350701/69‑1 (MQ=255)
                                     tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:2442050/69‑1 (MQ=255)
                                     tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:1551066/69‑1 (MQ=255)
                                     tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:1310360/69‑1 (MQ=255)
                                     tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:101523/69‑1 (MQ=255)
                                     tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:1761745/69‑1 (MQ=255)
                                     tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:1137448/69‑1 (MQ=255)
                                     tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:1744816/69‑1 (MQ=255)
                                     tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:2838431/69‑1 (MQ=255)
                                     tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGCCAAaga                  <  1:4521/69‑1 (MQ=255)
                                     tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCCGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:2153431/69‑1 (MQ=255)
                                     tGACGATAACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:810735/69‑1 (MQ=255)
                                      gACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:1901823/68‑1 (MQ=255)
                                            cacaAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATgg                                 <  1:3216333/47‑1 (MQ=255)
                                            cacaAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCg           >  1:2809271/1‑69 (MQ=255)
                                                  tGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:2521248/56‑1 (MQ=255)
                                                   gggTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAaga                  <  1:2086860/55‑1 (MQ=255)
                                                     gttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACgag  >  1:197598/1‑69 (MQ=255)
                                                     gttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACgag  >  1:933611/1‑69 (MQ=255)
                                                     gttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACga   >  1:2217586/1‑68 (MQ=255)
                                                      ttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACgag  <  1:792647/68‑1 (MQ=255)
                                                      ttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACgag  <  1:1367099/68‑1 (MQ=255)
                                                      ttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACgag  <  1:42601/68‑1 (MQ=255)
                                                      ttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACgag  <  1:1293157/68‑1 (MQ=255)
                                                      ttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACgag  <  1:1219801/68‑1 (MQ=255)
                                                      ttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACgag  <  1:1216342/68‑1 (MQ=255)
                                                      ttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACgag  <  1:1101228/68‑1 (MQ=255)
                                                      ttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACgag  <  1:1402864/68‑1 (MQ=255)
                                                      ttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACgag  <  1:62264/68‑1 (MQ=255)
                                                      ttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACgag  <  1:1569560/68‑1 (MQ=255)
                                                      ttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACgag  <  1:415286/68‑1 (MQ=255)
                                                      ttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACgag  <  1:3398870/68‑1 (MQ=255)
                                                      ttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACgag  <  1:2121974/68‑1 (MQ=255)
                                                      ttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACgag  <  1:2482901/68‑1 (MQ=255)
                                                      ttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACgag  <  1:3100709/68‑1 (MQ=255)
                                                      ttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACgag  <  1:2779358/68‑1 (MQ=255)
                                                      ttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACgag  <  1:2995365/68‑1 (MQ=255)
                                                      ttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACgag  <  1:2927782/68‑1 (MQ=255)
                                                                |                                                         
AAGTGACGTTAGACACGGCGACTTATGTGCCGATCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAG  >  minE/178581‑178702

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: