Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,833,650 A→T 48.5% G131G (GGA→GGT ygbL → predicted class II aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,833,6500AT48.5% 3.1 / 126.7 103G131G (GGA→GGTygbLpredicted class II aldolase
Reads supporting (aligned to +/- strand):  ref base A (41/12);  new base T (27/23);  total (68/35)
Fisher's exact test for biased strand distribution p-value = 1.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.21e-01

Only 100 of 103 total aligned reads displayed.

ACAGCAGCAACGTTATTCGTCCGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTGGCGGAACT  >  minE/1833598‑1833715
                                                    |                                                                 
aCAGCAATAACGTTATTCGGCCGTTCACACCATACGTGGTGATGCGGATGGGTAATGTCCCGCTGGTg                                                    >  1:2012609/1‑68 (MQ=255)
aCAGCAATAACGTTATTCGGCCGTTCACACCATACGTGGTGATGCGGATGGGTAATGTCCCGCTGGTGc                                                   >  1:1005487/1‑69 (MQ=255)
aCAGCAATAACGTTATTCGGCCGTTCACACCATACGTGGTGATGCGGATGGGTAATGTCCCGCTGGTGc                                                   >  1:1778365/1‑69 (MQ=255)
aCAGCAATAACGTTATTCGGCCGTTCACACCATACGTGGTGATGCGGATGGGTAATGTCCCGCTGGTGc                                                   >  1:2413124/1‑69 (MQ=255)
aCAGCAATAACGTTATTCGGCCGTTCACACCATACGTGGTGATGCGGATGGGTAATGTCCCGCTGGTGc                                                   >  1:2570212/1‑69 (MQ=255)
aCAGCAATAACGTTATTCGGCCGTTCACACCATACGTGGTGATGCGGATGGGTAATGTCCCGCTGGTGc                                                   >  1:3376086/1‑69 (MQ=255)
aCAGCAATAACGTTATTCGGCCGTTCACACCATACGTGGTGATGCGGATGGGTAATGTCCCGCTGGTGc                                                   >  1:1340701/1‑69 (MQ=255)
aCAGCAATAACGTTATTCGGCCGTTCACACCATACGTGGTGATGCGGATGGGTAATGTCCCGCTGGTGc                                                   >  1:1273327/1‑69 (MQ=255)
aCAGCAATAACGTTATTCGGCCGTTCACACCATACGTGGTGATGCGGATGGGTAATGTCCCGCTGGTGc                                                   >  1:2920044/1‑69 (MQ=255)
aCAGCAATAACGTTATTCGGCCGTTCACACCATACGTGGTGATGCGGATGGGTAATGTCCCGCTGGTGc                                                   >  1:1217105/1‑69 (MQ=255)
aCAGCAATAACGTTATTCGGCCGTTCACACCATACGTGGTGATGCGGATGGGTAATGTCCCGCTGGTGc                                                   >  1:1116314/1‑69 (MQ=255)
aCAGCAATAACGTTATTCGGCCGTTCACACCATACGTGGTGATGCGGATGGGTAATGTCCCGCTGGTGc                                                   >  1:3166675/1‑69 (MQ=255)
aCAGCAATAACGTTATTCGGCCGTTCACACCATACGTGGTGATGCGGATGGTAATGTCCCGCTGGTGc                                                   >  1:1245140/1‑68 (MQ=255)
 cagcagCAACGTTATTCGTCCGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGc                                                   >  1:1378116/1‑68 (MQ=255)
 cagcagCAACGTTATTCGTCCGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGc                                                   >  1:2688603/1‑68 (MQ=255)
 cagcagCAACGTTATTCGTCCGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGc                                                   >  1:391886/1‑68 (MQ=255)
 cagcagCAACGTTATTCGTCCGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGc                                                   >  1:2803514/1‑68 (MQ=255)
 cagcagCAACGTTATTCGTCCGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGc                                                   >  1:2748745/1‑68 (MQ=255)
 cagcagCAACGTTATTCGTCCGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGc                                                   >  1:3083613/1‑68 (MQ=255)
 cagcagCAACGTTATTCGTCCGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGc                                                   >  1:515412/1‑68 (MQ=255)
           gTTATTCGTCCGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTgg                                                      >  1:2288808/1‑55 (MQ=255)
           gTTATTCGTCCGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTgg                                                      >  1:635886/1‑55 (MQ=255)
           gTTATTCGTCCGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTgg                                                      >  1:2840932/1‑55 (MQ=255)
           gTTATTCGTCCGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTgg                                                      >  1:2531642/1‑55 (MQ=255)
           gTTATTCGTCCGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTgg                                                      >  1:1590473/1‑55 (MQ=255)
                     cGTTCAGACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATa                              >  1:436975/1‑69 (MQ=255)
                     cGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATa                              >  1:1628626/1‑69 (MQ=255)
                     cGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATa                              >  1:2899340/1‑69 (MQ=255)
                     cGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATa                              >  1:2689787/1‑69 (MQ=255)
                     cGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGAACGGGCGATa                              >  1:2075661/1‑69 (MQ=255)
                       ttCACACCATACGTGGTGATGCGGATGGGAAATGt                                                              >  1:2563560/1‑35 (MQ=255)
                       ttCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATaaa                            >  1:2085672/1‑69 (MQ=255)
                       ttCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATaaa                            >  1:1946851/1‑69 (MQ=255)
                       ttCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATaaa                            >  1:730949/1‑69 (MQ=255)
                       ttCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATaaa                            >  1:1606307/1‑69 (MQ=255)
                       ttCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACAGACCGGGCGATaaa                            >  1:370054/1‑69 (MQ=255)
                         cacaCCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACg                          >  1:780842/1‑69 (MQ=255)
                         cacaCCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACg                          >  1:2134047/1‑69 (MQ=255)
                             ccATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGCCt                                                 <  1:1066200/42‑1 (MQ=255)
                              cATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCg                     >  1:1044555/1‑69 (MQ=255)
                              cATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCg                     >  1:1178282/1‑69 (MQ=255)
                                    tggtgATGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGcaca                 >  1:2502022/1‑67 (MQ=255)
                                    tggtgATGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAgg               >  1:1210791/1‑69 (MQ=255)
                                    tggtgATGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAgg               >  1:1399063/1‑69 (MQ=255)
                                    tggtgATGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAgg               >  1:1398262/1‑69 (MQ=255)
                                    tggtgATGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAgg               >  1:2763480/1‑69 (MQ=255)
                                    tggtgATGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAgg               >  1:1959432/1‑69 (MQ=255)
                                    tggtgATGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAgg               >  1:3123113/1‑69 (MQ=255)
                                    tggtgATGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAgg               >  1:3065918/1‑69 (MQ=255)
                                    tggtgATGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAgg               >  1:2482495/1‑69 (MQ=255)
                                    tggtgATGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAgg               >  1:603922/1‑69 (MQ=255)
                                    tggtgATGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAgg               >  1:507263/1‑69 (MQ=255)
                                    tggtgATGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAgg               >  1:1588837/1‑69 (MQ=255)
                                    tggtgATGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAg                >  1:1949111/1‑68 (MQ=255)
                                     ggtgATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGa              <  1:1433418/69‑1 (MQ=255)
                                     ggtgATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGa              <  1:1569019/69‑1 (MQ=255)
                                     ggtgATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGa              <  1:3098522/69‑1 (MQ=255)
                                     ggtgATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGa              <  1:2998359/69‑1 (MQ=255)
                                     ggtgATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGa              <  1:2972900/69‑1 (MQ=255)
                                     ggtgATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGa              <  1:1808930/69‑1 (MQ=255)
                                     ggtgATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGa              <  1:235563/69‑1 (MQ=255)
                                     ggtgATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGa              <  1:2403099/69‑1 (MQ=255)
                                      gtgATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGa              <  1:1473528/68‑1 (MQ=255)
                                      gtgATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGa              <  1:344605/68‑1 (MQ=255)
                                         aTGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCTTCGCACAGGATCTg          <  1:2868693/69‑1 (MQ=255)
                                         aTGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTg          <  1:844914/69‑1 (MQ=255)
                                         aTGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTg          <  1:543673/69‑1 (MQ=255)
                                         aTGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTg          <  1:907964/69‑1 (MQ=255)
                                         aTGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTg          <  1:82663/69‑1 (MQ=255)
                                         aTGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTg          <  1:571715/69‑1 (MQ=255)
                                         aTGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTg          <  1:917530/69‑1 (MQ=255)
                                         aTGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTg          <  1:3170402/69‑1 (MQ=255)
                                         aTGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTg          <  1:818037/69‑1 (MQ=255)
                                         aTGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTg          <  1:2417333/69‑1 (MQ=255)
                                         aTGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTg          <  1:1023582/69‑1 (MQ=255)
                                         aTGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTg          <  1:1200039/69‑1 (MQ=255)
                                         aTGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTg          <  1:1671454/69‑1 (MQ=255)
                                         aTGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTg          <  1:1779318/69‑1 (MQ=255)
                                         aTGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTg          <  1:190998/69‑1 (MQ=255)
                                         aTGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTg          <  1:2307392/69‑1 (MQ=255)
                                         aTGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTg          <  1:2583651/69‑1 (MQ=255)
                                         aTGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTg          <  1:2684289/69‑1 (MQ=255)
                                         aTGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTg          <  1:2688070/69‑1 (MQ=255)
                                         aTGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTg          <  1:3018022/69‑1 (MQ=255)
                                         aTGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTg          <  1:2870116/69‑1 (MQ=255)
                                         aTGCGGATGGGTAATGTCCCGCTGGTGCCCTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTg          <  1:2948463/69‑1 (MQ=255)
                                           gCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTgg         >  1:2720879/1‑68 (MQ=255)
                                           gCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTgg         >  1:2302006/1‑68 (MQ=255)
                                                 gggAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTGGCGGAACt  >  1:3088517/1‑69 (MQ=255)
                                                 gggAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTGGCGGAACt  >  1:1306042/1‑69 (MQ=255)
                                                 gggAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTGGCGGAACt  >  1:1446454/1‑69 (MQ=255)
                                                 gggAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTGGCGGAACt  >  1:149303/1‑69 (MQ=255)
                                                 gggAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTGGCGGAACt  >  1:2670489/1‑69 (MQ=255)
                                                   gAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTGGCGGAACt  >  1:2133073/1‑67 (MQ=255)
                                                   gAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTGGCGGAACt  >  1:3025285/1‑67 (MQ=255)
                                                   gAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTGGCGGAACt  >  1:246616/1‑67 (MQ=255)
                                                   gAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTGGCGGAACt  >  1:3356175/1‑67 (MQ=255)
                                                   gAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTGGCGGAACt  >  1:835828/1‑67 (MQ=255)
                                                   gAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTGGCGGAACt  >  1:1270537/1‑67 (MQ=255)
                                                   gAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTGGCGGAACt  >  1:3095188/1‑67 (MQ=255)
                                                    |                                                                 
ACAGCAGCAACGTTATTCGTCCGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCGATAAACGCATCGCACAGGATCTGGCGGAACT  >  minE/1833598‑1833715

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