Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,007,034 A→G 31.5% I289I (ATT→ATC) ‡ hybA ← hydrogenase 2 4Fe‑4S ferredoxin‑type component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,007,0340AG31.5% 81.9 / 48.1 95I289I (ATT→ATC) ‡hybAhydrogenase 2 4Fe‑4S ferredoxin‑type component
Reads supporting (aligned to +/- strand):  ref base A (44/21);  new base G (19/11);  total (63/32)
Fisher's exact test for biased strand distribution p-value = 8.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.96e-01

CGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAGGCA  >  minE/2006976‑2007101
                                                          |                                                                   
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTctg                                                           <  1:910981/69‑3 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTc                                                             >  1:3384911/1‑67 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTc                                                             >  1:1772431/1‑67 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                           >  1:105606/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                           >  1:1982066/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                           <  1:3320038/69‑1 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                           >  1:427512/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                           >  1:1645665/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                           >  1:597354/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                           >  1:2069923/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                           <  1:155206/69‑1 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                           >  1:784762/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                           <  1:2174470/69‑1 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                           <  1:826113/69‑1 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                           >  1:865440/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                           <  1:1408643/69‑1 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                           <  1:3315969/69‑1 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                           >  1:1348982/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                           >  1:235211/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                           <  1:2421767/69‑1 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                           >  1:2600987/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                           >  1:2614432/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                           >  1:2709212/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                           >  1:2994149/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                           >  1:975574/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTAAACAGGGTGTGTTGGACATGTTCgg                                                           >  1:109026/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCGGa                                                          >  1:1966919/1‑69 (MQ=255)
 ggTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                           <  1:324552/68‑1 (MQ=255)
 ggTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                           <  1:3042716/68‑1 (MQ=255)
 ggTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                           <  1:2365528/68‑1 (MQ=255)
    tAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAAcg                                                       >  1:3175634/1‑69 (MQ=255)
    tAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAAcg                                                       >  1:1890557/1‑69 (MQ=255)
    tAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAAcg                                                       >  1:1892838/1‑69 (MQ=255)
    tAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAAcg                                                       >  1:889957/1‑69 (MQ=255)
    tAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAAcg                                                       >  1:1294219/1‑69 (MQ=255)
    tAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAAcg                                                       >  1:1391043/1‑69 (MQ=255)
    tAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAAcg                                                       >  1:424940/1‑69 (MQ=255)
    tAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAAcg                                                       >  1:1124058/1‑69 (MQ=255)
    tAAGCCCGCCAGCACAGCCAGTGGAAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAAcg                                                       >  1:2214784/1‑69 (MQ=255)
                cacaGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTc                                                             >  1:404893/1‑51 (MQ=255)
                cacaGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTa                                             >  1:1856500/1‑67 (MQ=255)
                cacaGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTa                                             >  1:2180853/1‑67 (MQ=255)
                   aGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTAGAAAg                                        <  1:1217913/69‑1 (MQ=255)
                        gTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTAGAAAg                                        <  1:1953281/64‑1 (MQ=255)
                              gCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGtt                             >  1:79026/1‑69 (MQ=255)
                              gCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGtt                             >  1:955446/1‑69 (MQ=255)
                              gCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGtt                             >  1:2331771/1‑69 (MQ=255)
                                   atGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTa                                             >  1:1150134/1‑48 (MQ=255)
                                   atGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGc                        >  1:286309/1‑69 (MQ=255)
                                                gggTGTGTTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGAtttt           <  1:1090773/69‑1 (MQ=255)
                                                gggTGTGTTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGAtttt           <  1:299472/69‑1 (MQ=255)
                                                gggTGTGTTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGAtttt           <  1:122617/69‑1 (MQ=255)
                                                gggTGTGTTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGAtttt           <  1:2609215/69‑1 (MQ=255)
                                                     tgtTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATaa      >  1:2420142/1‑69 (MQ=255)
                                                     tgtTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATaa      >  1:884816/1‑69 (MQ=255)
                                                     tgtTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATaa      >  1:2984889/1‑69 (MQ=255)
                                                     tgtTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATaa      >  1:1735054/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGt                              >  1:3127181/1‑43 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGt                              >  1:3062999/1‑43 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGt                              >  1:925769/1‑43 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGt                              >  1:302880/1‑43 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGAt              >  1:2240079/1‑59 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATaaag    <  1:2799241/69‑3 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg    >  1:3131671/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg    <  1:1089899/69‑1 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg    >  1:644939/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg    >  1:3323726/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg    >  1:2533203/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg    >  1:1113145/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg    <  1:1204832/69‑1 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg    >  1:1272321/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg    >  1:1285003/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg    >  1:1354870/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg    >  1:1419154/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg    <  1:1475527/69‑1 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg    >  1:153832/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg    <  1:1540040/69‑1 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg    >  1:163398/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg    >  1:1718824/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg    >  1:1764248/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg    <  1:1780045/69‑1 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg    <  1:198313/69‑1 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg    >  1:2202648/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg    >  1:2440496/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg    <  1:2484651/69‑1 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg    <  1:2949304/69‑1 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCAAGTTTCGGCAGGTCGAGATTTTCATAAgg    >  1:3124878/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGACCAGTTTCGGCAGGTCGAGATTTTCATAAgg    <  1:640484/69‑1 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCAGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg    >  1:2221693/1‑69 (MQ=255)
                                                       ttGAATATTTTCCGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGTCAGGTCGAGATTTTCATAAgg    >  1:2256112/1‑69 (MQ=255)
                                                         gAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAGTca  <  1:2053138/69‑1 (MQ=255)
                                                         gAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAGGca  <  1:399343/69‑1 (MQ=255)
                                                         gAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAGGca  <  1:1680641/69‑1 (MQ=255)
                                                         gAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAGGca  <  1:1569817/69‑1 (MQ=255)
                                                         gAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAGGca  <  1:1128442/69‑1 (MQ=255)
                                                          |                                                                   
CGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAGGCA  >  minE/2006976‑2007101

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: