Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,769 A→G 12.7% N163D (AAC→GAC) ‡ 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,7690AG12.7% 118.5 / 8.9 47N163D (AAC→GAC) ‡narHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base A (18/23);  new base G (4/2);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

GGCCGAACTGGGAAGATGACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAAC  >  minE/835704‑835838
                                                                 |                                                                     
ggCCGAACTGGGAAGATGACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACAt                                                                   <  1:1907633/70‑1 (MQ=255)
             aGATGACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATc                                                                  >  1:3223780/1‑58 (MQ=255)
                 gACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAg                                                            >  1:3178455/1‑60 (MQ=255)
                 gACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAg                                                 >  1:2417978/1‑71 (MQ=255)
                      gggCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCg                                                         >  1:2033963/1‑58 (MQ=255)
                      gggCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCg                                                         >  1:159614/1‑58 (MQ=255)
                       ggCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGc                                                          <  1:3930149/56‑1 (MQ=255)
                          ggTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGa                                                              <  1:2630266/49‑1 (MQ=255)
                          ggTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGa                                        <  1:3582048/71‑1 (MQ=255)
                           gTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATcc                                                                 <  1:1136780/45‑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)
                                     cAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGta                                                    <  1:3410205/48‑1 (MQ=255)
                                     cAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCa                             <  1:1209796/71‑1 (MQ=255)
                                      aaaCTGGCGAAAGACAAGAACTTCGCCAACATCCAGAAGGCGATGtat                                                   >  1:3596208/1‑48 (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:1396351/1‑67 (MQ=255)
                                                      aGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGcc                >  1:2299302/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:1546133/1‑71 (MQ=255)
                                                         aCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCtg         <  1:2015105/71‑1 (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)
                                                                caacaTCCAGAAGGCGATGTATAGCCAGTTCGaaa                                      >  1:1676077/1‑35 (MQ=255)
                                                                caacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAac  <  1:3166509/71‑1 (MQ=255)
                                                                 |                                                                     
GGCCGAACTGGGAAGATGACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAAC  >  minE/835704‑835838

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: