Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 481,239 C→T 49.3% R11R (CGC→CGT ybgC → predicted acyl‑CoA thioesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE481,2390CT49.3% ‑2.3 / 63.9 71R11R (CGC→CGTybgCpredicted acyl‑CoA thioesterase
Reads supporting (aligned to +/- strand):  ref base C (23/13);  new base T (19/16);  total (42/29)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

AGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAGATACCGATGCCGGTGGTGTGGTGTACCACGCCAGT  >  minE/481173‑481287
                                                                  |                                                
aGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTcgcg                                                 <  1:2083001/68‑1 (MQ=255)
aGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTcgcg                                                 <  1:401867/68‑1 (MQ=255)
aGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTcgcg                                                 <  1:683506/68‑1 (MQ=255)
aGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTcgcg                                                 <  1:289030/68‑1 (MQ=255)
aGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTcgcg                                                 <  1:1325125/68‑1 (MQ=255)
aGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTcgcg                                                 <  1:133833/68‑1 (MQ=255)
aGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTcgcg                                                 <  1:2355473/68‑1 (MQ=255)
aGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTcgcg                                                 <  1:2384162/68‑1 (MQ=255)
aGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTcgcg                                                 <  1:1808583/68‑1 (MQ=255)
aGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTcgcg                                                 <  1:334090/68‑1 (MQ=255)
aGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTcgcg                                                 <  1:2527540/68‑1 (MQ=255)
 gCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTc                                               <  1:3358982/69‑1 (MQ=255)
 gCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTc                                               <  1:1732708/69‑1 (MQ=255)
 gCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTc                                               <  1:2586013/69‑1 (MQ=255)
 gCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTc                                               <  1:1655450/69‑1 (MQ=255)
 gCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTc                                               <  1:2663154/69‑1 (MQ=255)
  ctctAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTct                                              >  1:141688/1‑69 (MQ=255)
  ctctAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTct                                              >  1:1073342/1‑69 (MQ=255)
    ctAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTctac                                            <  1:2202252/69‑1 (MQ=255)
        cTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTctacta                                          >  1:411454/1‑67 (MQ=255)
        cTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTctacta                                          >  1:3418275/1‑67 (MQ=255)
        cTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTctacta                                          >  1:267743/1‑67 (MQ=255)
        cTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTctacta                                          >  1:2842014/1‑67 (MQ=255)
        cTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTctacta                                          >  1:3003091/1‑67 (MQ=255)
        cTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTctacta                                          >  1:2335721/1‑67 (MQ=255)
        cTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTctacta                                          >  1:3009975/1‑67 (MQ=255)
        cTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTctacta                                          >  1:2195908/1‑67 (MQ=255)
        cTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTctacta                                          >  1:1286777/1‑67 (MQ=255)
        cTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTctacta                                          >  1:1479490/1‑67 (MQ=255)
        cTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTctacta                                          >  1:86392/1‑67 (MQ=255)
        cTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTctacta                                          >  1:1451701/1‑67 (MQ=255)
        cTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTctacta                                          >  1:1784488/1‑67 (MQ=255)
        cTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTctacta                                          >  1:181908/1‑67 (MQ=255)
        cTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTctacta                                          >  1:963752/1‑67 (MQ=255)
        cTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTctacta                                          >  1:1928906/1‑67 (MQ=255)
        cTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTctacta                                          >  1:581900/1‑67 (MQ=255)
        cTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTctacta                                          >  1:1114416/1‑67 (MQ=255)
        cTTTTGTTGCATTAACGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTctacta                                          >  1:821319/1‑67 (MQ=255)
                gCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTCTACTATGAAGATAcc                                <  1:3198680/69‑1 (MQ=255)
                gCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTCTACTATGAAGATAcc                                <  1:331888/69‑1 (MQ=255)
                gCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTCTACTATGAAGATAcc                                <  1:307736/69‑1 (MQ=255)
                gCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTCTACTATGAAGATAcc                                <  1:3320154/69‑1 (MQ=255)
                gCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTCTACTATGAAGATAcc                                <  1:883393/69‑1 (MQ=255)
                gCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTCTACTATGAAGATAcc                                <  1:2813776/69‑1 (MQ=255)
                gCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTCTACTATGAAGATAcc                                <  1:506976/69‑1 (MQ=255)
                gCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTCTACTATGAAGATAcc                                <  1:684086/69‑1 (MQ=255)
                gCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTCTACTATGAAGATAcc                                <  1:2456069/69‑1 (MQ=255)
                gCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTCTACTATGAAGATAcc                                <  1:371471/69‑1 (MQ=255)
                gCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTCTACTATGAAGATAcc                                <  1:1814880/69‑1 (MQ=255)
                gCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTCTACTATGAAGATAcc                                <  1:2353880/69‑1 (MQ=255)
                gCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTCTACTATGAAGATAcc                                <  1:2001322/69‑1 (MQ=255)
                gCATTACCGGGATGTAAAGTGAATACAACGCTGTTACGATGGCCGGTTCGTGTCTACTATGAAGATAcc                                <  1:2940478/69‑1 (MQ=255)
                ccATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGTGTCTACTATGAAGATAc                                 <  1:1676937/67‑1 (MQ=255)
                        gggATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGaa                                      <  1:1878986/55‑1 (MQ=255)
                             gTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAGa                                    >  1:351785/1‑52 (MQ=255)
                             gTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAGa                                    >  1:1103074/1‑52 (MQ=255)
                             gTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAGa                                    >  1:2990011/1‑52 (MQ=255)
                             gTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAGa                                    >  1:1489258/1‑52 (MQ=255)
                             gTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAGa                                    >  1:2400993/1‑52 (MQ=255)
                             gTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAGATACCGATGCCGgtggtg                   >  1:2059508/1‑69 (MQ=255)
                             gTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAGATACCGATGCCGgtggtg                   >  1:979490/1‑69 (MQ=255)
                             gTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAGATACCGATGCCGgtggtg                   >  1:3123482/1‑69 (MQ=255)
                             gTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAGATACCGATGCCGgtggtg                   >  1:3410355/1‑69 (MQ=255)
                             gTAAAGTGAATACAACGCTGTTTCGATGGCGGTTCGCGTCTACTATGAAGATACCGATGCCGgtggtgt                  >  1:577123/1‑69 (MQ=255)
                              tAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAGATACCGATGCCggt                       >  1:3180053/1‑64 (MQ=255)
                              tAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAGATACCGATGCCggt                       >  1:1255991/1‑64 (MQ=255)
                              tAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAGATACCGATGCCggt                       >  1:1413599/1‑64 (MQ=255)
                              tAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAGATACCGATGCCggt                       >  1:2175599/1‑64 (MQ=255)
                              tAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAGATACCGATGCCggt                       >  1:1803184/1‑64 (MQ=255)
                              tAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAGATACCGATGCCggt                       >  1:678735/1‑64 (MQ=255)
                              tAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAGATACCGATGCCgg                        >  1:1714515/1‑63 (MQ=255)
                              tAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAGATACCGATGCCGgtggtgt                  >  1:1843937/1‑69 (MQ=255)
                                   tGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAGATACCGATGCCGGtggtatggtg              >  1:2627506/1‑68 (MQ=255)
                                   tGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAGATACCGATGCCGGTggtgtggtgt             >  1:2472436/1‑69 (MQ=255)
                                   tGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAGATACCGATGCCGGTggtgtggtgt             >  1:2257689/1‑69 (MQ=255)
                                      aTACAACGCTGTTTCGATGGCCGGTTCGTGTCTACTATGAAGATAcc                                <  1:2862211/47‑1 (MQ=255)
                                              cTGTTTCGATGGCCGGTTCGCGTCTACTATGAAGATACCGATGCCGGTGGTGTGGTGTACCACGCCAGt  <  1:718966/69‑1 (MQ=255)
                                              cTGTTTCGATGGCCGGTTCGCGTCTACTATGAAGATACCGATGCCGGTGGTGTGGTGTACCACGCCAGt  <  1:2668247/69‑1 (MQ=255)
                                              cTGTTTCGATGGCCGGTTCGCGTCTACTATGAAGATACCGATGCCGGTGGTGTGGTGCACCACGCCAGt  <  1:2402266/69‑1 (MQ=255)
                                                        ggCCGGTTCGCGTCTACTATGAAGATACCGATGCCGgtggtgtgg                <  1:1990000/45‑1 (MQ=255)
                                                        ggCCGGTTCGCGTCTACTATGAAGATACCGATGCCGgtggtgtgg                <  1:536604/45‑1 (MQ=255)
                                                        ggCCGGTTCGCGTCTACTATGAAGATACCGATGCCGgtggtgtgg                <  1:298018/45‑1 (MQ=255)
                                                                  |                                                
AGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAGATACCGATGCCGGTGGTGTGGTGTACCACGCCAGT  >  minE/481173‑481287

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: