Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,790 A→G 14.3% S170G (AGC→GGC) ‡ narH → nitrate reductase 1, beta (Fe‑S) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,7900AG14.3% 156.6 / 15.8 63S170G (AGC→GGC) ‡narHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base A (27/27);  new base G (7/2);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTG  >  minE/835721‑835859
                                                                     |                                                                     
gACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAg                                                                      >  1:2417978/1‑71 (MQ=255)
         ggTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGa                                                             <  1:3582048/71‑1 (MQ=255)
          gTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaa                                                            <  1:456674/71‑1 (MQ=255)
          gTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaa                                                            <  1:2308116/71‑1 (MQ=255)
             aGTTTGACAAACTGGCGAAAGACAAGAGCTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa                                                          >  1:2126997/1‑70 (MQ=255)
             aGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa                                                          >  1:2025126/1‑70 (MQ=255)
             aGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaa                                                            >  1:1919247/1‑68 (MQ=255)
               tttGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAacac                                                       >  1:2167966/1‑71 (MQ=255)
                 tGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGt                                                                 <  1:375901/59‑1 (MQ=255)
                    cAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCa                                                  <  1:1209796/71‑1 (MQ=255)
                          ggCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCatgatg                                             >  1:27783/1‑70 (MQ=255)
                           gCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaa                                                           <  1:2872116/55‑1 (MQ=255)
                           gCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaa                                                           <  1:981039/55‑1 (MQ=255)
                           gCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTa                                           <  1:1794791/71‑1 (MQ=255)
                            cGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTAt                                          <  1:2822654/71‑1 (MQ=255)
                              aaaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCa                                                  <  1:2760474/61‑1 (MQ=255)
                               aaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaa                                                            <  1:724718/50‑1 (MQ=255)
                               aaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTg                                       >  1:99782/1‑71 (MQ=255)
                                aGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGc                                      >  1:1473338/1‑71 (MQ=255)
                                   cAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCcgc                                   <  1:3433131/71‑1 (MQ=255)
                                     aGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGcc                                     >  1:2299302/1‑67 (MQ=255)
                                     aGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGcc                                     >  1:1396351/1‑67 (MQ=255)
                                     aGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCgcgc                                 >  1:1582637/1‑71 (MQ=255)
                                     aGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCgcgc                                 >  1:3482164/1‑71 (MQ=255)
                                      gAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGt                                                                 <  1:1399271/38‑1 (MQ=255)
                                      gAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTg                                       <  1:3467917/64‑1 (MQ=255)
                                      caaCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGcc                                >  1:2392205/2‑71 (MQ=255)
                                      caaCTTCGAGGCCATGCAAAAGGAGATGTACGCGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGcc                                >  1:2376532/2‑71 (MQ=255)
                                       aaCTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCt                               <  1:1312012/71‑1 (MQ=255)
                                       aaCTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCt                               <  1:2457881/71‑1 (MQ=255)
                                        aCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGACAACACCTTCATGATGTACCTGCCGCGCCtg                              <  1:3120492/71‑1 (MQ=255)
                                        aCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCtg                              >  1:614351/1‑71 (MQ=255)
                                        aCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCtg                              <  1:2015105/71‑1 (MQ=255)
                                        aCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCtg                              >  1:1546133/1‑71 (MQ=255)
                                          ttCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCcgc                                   <  1:1348976/64‑1 (MQ=255)
                                            cGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGc                           >  1:960591/1‑70 (MQ=255)
                                             gACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCg                          <  1:3804234/70‑1 (MQ=255)
                                              acaacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGc                           <  1:3075465/68‑1 (MQ=255)
                                               ggccaTGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGaa                        >  1:403939/4‑70 (MQ=255)
                                               caacaTCCAGAAGGCGATGTATAGCCAGTTCGaaa                                                           >  1:1676077/1‑35 (MQ=255)
                                               caacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAac                       <  1:3166509/71‑1 (MQ=255)
                                                  caTGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCatga                                               >  1:1306037/1‑44 (MQ=12)
                                                    tCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCgcgc                                 <  1:1159264/56‑1 (MQ=255)
                                                    tCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGc                  >  1:1714013/1‑71 (MQ=255)
                                                    tCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGc                  >  1:3526955/1‑71 (MQ=255)
                                                    tCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGc                  >  1:2269706/1‑71 (MQ=255)
                                                     ccAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGcc                 >  1:304097/1‑71 (MQ=255)
                                                     ccAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGcc                 >  1:903867/1‑71 (MQ=255)
                                                      cAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCt                >  1:608740/1‑71 (MQ=255)
                                                       aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg               <  1:945051/71‑1 (MQ=255)
                                                       aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg               <  1:489827/71‑1 (MQ=255)
                                                       aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg               <  1:940361/71‑1 (MQ=255)
                                                       aGAAGGCGAGGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCt                >  1:3215519/1‑70 (MQ=255)
                                                        gAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGa              <  1:374529/71‑1 (MQ=255)
                                                          aGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAAc            >  1:1848602/1‑71 (MQ=255)
                                                          aGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAAc            <  1:135266/71‑1 (MQ=255)
                                                            gCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAAccc          <  1:3876236/71‑1 (MQ=255)
                                                             cGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCg         >  1:3273268/1‑71 (MQ=255)
                                                              gATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAaca                      >  1:3508944/1‑57 (MQ=255)
                                                                tGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGc                  >  1:2309667/1‑59 (MQ=255)
                                                                tGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGc                  >  1:937686/1‑59 (MQ=255)
                                                                tGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCa      >  1:2271441/1‑71 (MQ=255)
                                                                tGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCa      >  1:2081803/1‑71 (MQ=255)
                                                                  tataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAac                       >  1:2661842/1‑52 (MQ=255)
                                                                  tataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtg    >  1:2041892/1‑71 (MQ=255)
                                                                  tataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAt     <  1:1593906/70‑1 (MQ=255)
                                                                    taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtg  <  1:1472040/71‑1 (MQ=255)
                                                                    taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtg  >  1:3336679/1‑71 (MQ=255)
                                                                    taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtg  <  1:3181309/71‑1 (MQ=255)
                                                                    taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgt   <  1:3621570/70‑1 (MQ=255)
                                                                     aGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtg  <  1:3689200/70‑1 (MQ=255)
                                                                     |                                                                     
GACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTG  >  minE/835721‑835859

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: