Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,768 C→G 12.7% D162E (GAC→GAG 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,7680CG12.7% 118.2 / 8.5 47D162E (GAC→GAGnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base C (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 = 7.99e-01

AAAAGGGCCGAACTGGGAAGATGACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAAC  >  minE/835699‑835838
                                                                     |                                                                      
aaaaGGGCCGAACTGGGAAGATGACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGac                                                                        >  1:1376093/1‑70 (MQ=255)
aaaaGGGCCGAACTGGGAAGATGACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGac                                                                        <  1:1709799/70‑1 (MQ=255)
     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:2308116/71‑1 (MQ=255)
                                gTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaa                                       <  1:456674/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:3482164/1‑71 (MQ=255)
                                                           aGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCgcgc            >  1:1582637/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:1546133/1‑71 (MQ=255)
                                                              aCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCtg         <  1:2015105/71‑1 (MQ=255)
                                                              aCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCtg         >  1:614351/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)
                                                                     caacaTCCAGAAGGCGATGTATAGCCAGTTCGaaa                                      >  1:1676077/1‑35 (MQ=255)
                                                                     caacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAac  <  1:3166509/71‑1 (MQ=255)
                                                                     |                                                                      
AAAAGGGCCGAACTGGGAAGATGACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAAC  >  minE/835699‑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: