Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,817,164 T→C 20.7% A95A (GCT→GCC nrfA → nitrite reductase, formate‑dependent, cytochrome

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,817,1640TC20.7% 128.3 / 24.5 91A95A (GCT→GCCnrfAnitrite reductase, formate‑dependent, cytochrome
Reads supporting (aligned to +/- strand):  ref base T (45/27);  new base C (5/14);  total (50/41)
Fisher's exact test for biased strand distribution p-value = 8.42e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGAAGATG  >  minE/2817096‑2817225
                                                                    |                                                             
ccACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCt                                                               >  1:3098426/1‑69 (MQ=255)
ccACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCt                                                               >  1:1610387/1‑69 (MQ=255)
          tGATCCTGTGGGCAGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCCTTTGCtgtg                                                      >  1:1181342/1‑68 (MQ=255)
          tGATCCTGTGGGCAGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCCTTTGCTGTGa                                                     >  1:1047984/1‑69 (MQ=255)
          tGATCCTGTGGGCAGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCCTTTGCTGTGa                                                     >  1:2031712/1‑69 (MQ=255)
          tGATCCTGTGGGCAGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCCTTTGCTGTGa                                                     >  1:298522/1‑69 (MQ=255)
          tGATCCTGTGGGCAGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCCTTTGCTGTGa                                                     >  1:1503669/1‑69 (MQ=255)
           gATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGa                                                     >  1:2654239/1‑68 (MQ=255)
           gATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGa                                                     >  1:939794/1‑68 (MQ=255)
           gATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGa                                                     >  1:2954516/1‑68 (MQ=255)
           gATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGa                                                     >  1:191563/1‑68 (MQ=255)
           gATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGa                                                     >  1:423950/1‑68 (MQ=255)
           gATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGa                                                     >  1:17145/1‑68 (MQ=255)
           gATCCTGTGGGCGGAGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGa                                                     >  1:585867/1‑68 (MQ=255)
              ccTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGc                                                          >  1:940331/1‑60 (MQ=255)
              ccTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGc                                                          >  1:2729603/1‑60 (MQ=255)
              ccTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGc                                                          >  1:2785136/1‑60 (MQ=255)
              ccTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGc                                                          >  1:2814507/1‑60 (MQ=255)
              ccTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGc                                                          >  1:2122044/1‑60 (MQ=255)
              ccTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGc                                                          >  1:1104785/1‑60 (MQ=255)
              ccTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGc                                                          >  1:1260360/1‑60 (MQ=255)
              ccTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGc                                                          >  1:1609463/1‑60 (MQ=255)
              ccTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGc                                                          >  1:720145/1‑60 (MQ=255)
              ccTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGa                                                 >  1:1108795/1‑69 (MQ=255)
              ccTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGa                                                 >  1:70089/1‑69 (MQ=255)
                    ggCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGAtgtg                                             <  1:2594198/67‑1 (MQ=255)
                    ggCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGAtgtg                                             <  1:2960120/67‑1 (MQ=255)
                      cGGTGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTg                                         <  1:1852959/69‑1 (MQ=255)
                      cGGGGTATCCCTTCTCGCGCTCTTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTg                                         <  1:776532/69‑1 (MQ=255)
                      cGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTg                                         <  1:3435940/69‑1 (MQ=255)
                      cGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTg                                         <  1:2595165/69‑1 (MQ=255)
                      cGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTg                                         <  1:2716346/69‑1 (MQ=255)
                      cGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTg                                         <  1:22899/69‑1 (MQ=255)
                      cGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTg                                         <  1:2076481/69‑1 (MQ=255)
                      cGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTg                                         <  1:947027/69‑1 (MQ=255)
                      cGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTg                                         <  1:1192915/69‑1 (MQ=255)
                      cGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTg                                         <  1:1144187/69‑1 (MQ=255)
                      cGGGCTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTg                                         <  1:2940203/69‑1 (MQ=255)
                          gTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAAc                                     >  1:2887545/1‑69 (MQ=255)
                          gTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAAc                                     >  1:1392214/1‑69 (MQ=255)
                          gTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAAc                                     >  1:1959433/1‑69 (MQ=255)
                          gTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAAc                                     >  1:37528/1‑69 (MQ=255)
                              cccTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTg                                 >  1:2814583/1‑69 (MQ=255)
                              cccTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTg                                 >  1:2750730/1‑69 (MQ=255)
                              cccTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTg                                 >  1:254032/1‑69 (MQ=255)
                              cccTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTg                                 >  1:237234/1‑69 (MQ=255)
                              cccTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTg                                 >  1:2289441/1‑69 (MQ=255)
                               ccTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCtgt                                                       <  1:2661843/46‑1 (MQ=255)
                               ccTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCtgt                                                       <  1:76980/46‑1 (MQ=255)
                               ccTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGa                                                 >  1:1817767/1‑52 (MQ=255)
                               ccTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGa                                                 >  1:3015859/1‑52 (MQ=255)
                               ccTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGa                                                 >  1:3407624/1‑52 (MQ=255)
                               ccTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGa                                                 >  1:2155815/1‑52 (MQ=255)
                               ccTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGa                                                 >  1:1464517/1‑52 (MQ=255)
                               ccTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGa                                                 >  1:2244585/1‑52 (MQ=255)
                               ccTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGa                                                 >  1:942423/1‑52 (MQ=255)
                               ccTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGa                                                 >  1:120767/1‑52 (MQ=255)
                               ccTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGa                                                 >  1:2530678/1‑52 (MQ=255)
                               ccTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGa                                                 >  1:2467706/1‑52 (MQ=255)
                               ccTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGa                                                 >  1:2494895/1‑52 (MQ=255)
                                   ctcGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTg                                         <  1:324735/56‑1 (MQ=255)
                                                     gCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTg                                         <  1:1689283/38‑1 (MQ=255)
                                                       cgcgTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTg        <  1:461454/69‑1 (MQ=255)
                                                       cgcgTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTg        <  1:997420/69‑1 (MQ=255)
                                                       cgcgTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTg        <  1:2766531/69‑1 (MQ=255)
                                                       cgcgTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTg        <  1:1875975/69‑1 (MQ=255)
                                                       cgcgTGGACATGCCTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTg        <  1:572190/69‑1 (MQ=255)
                                                       cgcgTGGACATGCCTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTg        <  1:1314768/69‑1 (MQ=255)
                                                       cgcgTGGACATGCCTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTg        <  1:890392/69‑1 (MQ=255)
                                                       cgcgTGGACATGCCTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTg        <  1:13042/69‑1 (MQ=255)
                                                       cgcgTGGACATGCCTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTg        <  1:75940/69‑1 (MQ=255)
                                                       cgcgTGGACATGCCTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTg        <  1:734653/69‑1 (MQ=255)
                                                       cgcgTGGACATGCCTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTg        <  1:234316/69‑1 (MQ=255)
                                                       cgcgTGGACATGCCTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTg        <  1:703081/69‑1 (MQ=255)
                                                       cgcgTGGACATGCCTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTg        <  1:971057/69‑1 (MQ=255)
                                                       cgcgTGGACATGCCTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTg        <  1:2069187/69‑1 (MQ=255)
                                                       cgcgTGGACATGCCTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTg        <  1:3377591/69‑1 (MQ=255)
                                                       cgcgTGGACATGCCTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTg        <  1:3162075/69‑1 (MQ=255)
                                                       cgcgTGGACATGCCTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTg        <  1:1779190/69‑1 (MQ=255)
                                                       cgcgTGGACATGCCTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTg        <  1:1980208/69‑1 (MQ=255)
                                                           tGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGAAGa    <  1:2242730/69‑1 (MQ=255)
                                                           tGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGAAGa    <  1:2426842/69‑1 (MQ=255)
                                                           tGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGAAGa    <  1:1700132/69‑1 (MQ=255)
                                                           tGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGAAGa    <  1:1983869/69‑1 (MQ=255)
                                                            ggACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGAAGa    <  1:667670/68‑1 (MQ=255)
                                                            ggACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGAAGa    <  1:976150/68‑1 (MQ=255)
                                                             gACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGAAGATg  >  1:2436588/1‑69 (MQ=255)
                                                             gACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGAAGATg  >  1:1319453/1‑69 (MQ=255)
                                                             gACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGAAGATg  >  1:142812/1‑69 (MQ=255)
                                                             gACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGAAGATg  >  1:2799865/1‑69 (MQ=255)
                                                             gACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGAAGATg  >  1:1478096/1‑69 (MQ=255)
                                                             gACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGAAGATg  >  1:2816026/1‑69 (MQ=255)
                                                             gACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGAAGATg  >  1:1733889/1‑69 (MQ=255)
                                                                    |                                                             
CCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGAAGATG  >  minE/2817096‑2817225

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: