Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,007,036 T→C 31.5% I289V (ATT→GTT) ‡ 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,0360TC31.5% 80.7 / 47.8 95I289V (ATT→GTT) ‡hybAhydrogenase 2 4Fe‑4S ferredoxin‑type component
Reads supporting (aligned to +/- strand):  ref base T (44/21);  new base C (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 = 1.00e+00

Only 100 of 101 total aligned reads displayed.

CGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAGGCACAC  >  minE/2006976‑2007104
                                                            |                                                                    
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTctg                                                              <  1:910981/69‑3 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTc                                                                >  1:3384911/1‑67 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTc                                                                >  1:1772431/1‑67 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                              >  1:1348982/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                              <  1:3315969/69‑1 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                              <  1:3320038/69‑1 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                              >  1:1645665/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                              <  1:155206/69‑1 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                              >  1:427512/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                              >  1:597354/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                              >  1:784762/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                              >  1:1982066/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                              <  1:1408643/69‑1 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                              <  1:826113/69‑1 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                              >  1:2069923/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                              >  1:865440/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                              <  1:2174470/69‑1 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                              >  1:2994149/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                              >  1:105606/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                              >  1:235211/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                              <  1:2421767/69‑1 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                              >  1:2709212/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                              >  1:2614432/1‑69 (MQ=255)
cGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGGACATGTTCgg                                                              >  1:2600987/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:1124058/1‑69 (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:424940/1‑69 (MQ=255)
    tAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAAcg                                                          >  1:1391043/1‑69 (MQ=255)
    tAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAAcg                                                          >  1:1294219/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:955446/1‑69 (MQ=255)
                              gCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGtt                                >  1:2331771/1‑69 (MQ=255)
                              gCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGtt                                >  1:79026/1‑69 (MQ=255)
                                   atGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTa                                                >  1:1150134/1‑48 (MQ=255)
                                   atGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGc                           >  1:286309/1‑69 (MQ=255)
                                                gggTGTGTTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGAtttt              <  1:299472/69‑1 (MQ=255)
                                                gggTGTGTTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGAtttt              <  1:122617/69‑1 (MQ=255)
                                                gggTGTGTTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGAtttt              <  1:1090773/69‑1 (MQ=255)
                                                gggTGTGTTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGAtttt              <  1:2609215/69‑1 (MQ=255)
                                                     tgtTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATaa         >  1:884816/1‑69 (MQ=255)
                                                     tgtTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATaa         >  1:2984889/1‑69 (MQ=255)
                                                     tgtTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATaa         >  1:2420142/1‑69 (MQ=255)
                                                     tgtTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATaa         >  1:1735054/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGt                                 >  1:302880/1‑43 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGt                                 >  1:3062999/1‑43 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGt                                 >  1:3127181/1‑43 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGt                                 >  1:925769/1‑43 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGAt                 >  1:2240079/1‑59 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATaaag       <  1:2799241/69‑3 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg       >  1:1419154/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg       >  1:1764248/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg       >  1:2533203/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg       >  1:3323726/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg       >  1:163398/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg       <  1:1540040/69‑1 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg       >  1:153832/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg       >  1:644939/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg       <  1:1475527/69‑1 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg       >  1:1718824/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg       >  1:1354870/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg       >  1:1285003/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg       >  1:1272321/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg       <  1:1204832/69‑1 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg       >  1:1113145/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg       <  1:1089899/69‑1 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg       <  1:2484651/69‑1 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg       >  1:3131671/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg       >  1:2440496/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg       <  1:2949304/69‑1 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg       >  1:2202648/1‑69 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg       <  1:1780045/69‑1 (MQ=255)
                                                       ttGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAgg       <  1:198313/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:1569817/69‑1 (MQ=255)
                                                         gAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAGGca     <  1:399343/69‑1 (MQ=255)
                                                         gAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAGGca     <  1:1680641/69‑1 (MQ=255)
                                                         gAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAGGca     <  1:1128442/69‑1 (MQ=255)
                                                            tatTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAGGcaca   <  1:609750/68‑1 (MQ=255)
                                                            tatTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAGGcaca   <  1:2094164/68‑1 (MQ=255)
                                                            tatTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAGGcaca   >  1:2215148/1‑68 (MQ=255)
                                                            tatTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAGGcaca   <  1:1939125/68‑1 (MQ=255)
                                                            tatTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAGGcaca   <  1:1730204/68‑1 (MQ=255)
                                                            tatTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAGGCacac  >  1:2676961/1‑69 (MQ=255)
                                                            |                                                                    
CGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTCCAGTTTCGGCAGGTCGAGATTTTCATAAGGCACAC  >  minE/2006976‑2007104

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: