Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,792 C→G 8.1% S170R (AGC→AGG 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,7920CG 8.1% 173.0 / 3.9 61S170R (AGC→AGGnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base C (28/26);  new base G (3/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGT  >  minE/835722‑835860
                                                                      |                                                                    
aCCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGc                                                                      <  1:684421/71‑1 (MQ=255)
         gTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaa                                                             <  1:9350/71‑1 (MQ=255)
         gTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGACAGTTCGaa                                                             <  1:2735343/71‑1 (MQ=255)
           gAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa                                                           >  1:1350532/1‑71 (MQ=255)
            aGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa                                                           <  1:3080244/70‑1 (MQ=255)
            aGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa                                                           >  1:2475322/1‑70 (MQ=255)
            aGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa                                                           <  1:1817962/70‑1 (MQ=255)
            aGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAac                                                          <  1:828789/71‑1 (MQ=255)
                  aCAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaa                                                             <  1:1257208/62‑1 (MQ=255)
                  aCAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCtt                                                     <  1:3736150/70‑1 (MQ=255)
                   cAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCa                                                   >  1:3185386/1‑71 (MQ=255)
                       cTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCatgat                                               >  1:2773951/1‑71 (MQ=255)
                          gCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa                                                           >  1:1860346/1‑56 (MQ=255)
                            gAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTAtt                                          <  1:3874513/71‑1 (MQ=255)
                              aaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTAttt                                         <  1:1911298/70‑1 (MQ=255)
                              aaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGACAGTTCGAAAACACCTTCATGATGt                                             >  1:2040788/1‑66 (MQ=255)
                               aGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCa                                                                    >  1:482701/1‑42 (MQ=255)
                               aGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAg                                                                   <  1:3372591/43‑1 (MQ=255)
                               aGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGc                                       <  1:3845997/71‑1 (MQ=255)
                                 aCAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGtt                                                                 >  1:3337209/1‑43 (MQ=255)
                                  cAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCcgc                                    >  1:2973612/1‑71 (MQ=255)
                                   aaGAACTTCGACAACATCCAGAAGGCGATGTATAGCCa                                                                    >  1:3595282/1‑38 (MQ=255)
                                    aGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGc                                       >  1:2139478/1‑66 (MQ=255)
                                        cTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCtgt                              <  1:2042434/71‑1 (MQ=255)
                                        cTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCtgt                              >  1:2842541/1‑71 (MQ=255)
                                         ttCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTc                                                    <  1:1021256/48‑1 (MQ=255)
                                          tCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGc                            >  1:2295280/1‑71 (MQ=255)
                                           cGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCatga                                                >  1:2748015/1‑50 (MQ=255)
                                            gACAACATCCAGAAGGCGATGTATAGGCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCg                           <  1:3809266/70‑1 (MQ=255)
                                            gACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCg                           <  1:327566/70‑1 (MQ=255)
                                              caacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCat                                                  <  1:549767/45‑1 (MQ=255)
                                              caacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCat                                                  <  1:1667375/45‑1 (MQ=255)
                                               aacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGcc                                 <  1:668319/61‑1 (MQ=255)
                                                 caTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCatga                                                >  1:2481738/1‑44 (MQ=255)
                                                 caTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCatga                                                >  1:1318590/1‑44 (MQ=255)
                                                 caTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCatga                                                >  1:8274/1‑44 (MQ=255)
                                                  aTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGt                                             <  1:220207/46‑1 (MQ=255)
                                                   tCCAGAAGGCGATGTATAGCCAGTTCGAAAACAcc                                                       <  1:1186687/35‑1 (MQ=255)
                                                    ccAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAacac                      <  1:48615/67‑1 (MQ=255)
                                                     caaaaGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCt                 >  1:1609970/4‑71 (MQ=255)
                                                     caaaaGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCt                 >  1:3979858/4‑71 (MQ=255)
                                                     caaaaGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCt                 >  1:1437971/4‑71 (MQ=255)
                                                     cAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCatg                                                 >  1:448204/1‑39 (MQ=255)
                                                     cAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAaca                       >  1:1616683/1‑65 (MQ=255)
                                                        aaGGAGATGTACGGGCAGTTTGAAAACACCTTC‑‑‑ATGTACCTGCCGCGCCTGTGCGAac                        <  1:801635/58‑1 (MQ=255)
                                                         aGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAaca                       >  1:3820804/1‑61 (MQ=255)
                                                          ggCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAAcc            <  1:3816693/71‑1 (MQ=255)
                                                           gCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTTCGAac                        >  1:1861873/1‑58 (MQ=255)
                                                             gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGACCACTGCCTGAACCCg          >  1:3526373/1‑70 (MQ=255)
                                                             gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCg          >  1:2089863/1‑70 (MQ=255)
                                                              aTGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGcc                                 <  1:2474934/46‑1 (MQ=255)
                                                              aTGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGcc                                 <  1:2674594/46‑1 (MQ=255)
                                                              aTGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAac                        <  1:1813714/55‑1 (MQ=255)
                                                              aTGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCgg         >  1:2079504/1‑70 (MQ=255)
                                                                 tataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGaa                         >  1:2888282/1‑51 (MQ=255)
                                                                 tataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtg     <  1:1048690/71‑1 (MQ=255)
                                                                  ataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGc                            >  1:3322785/1‑47 (MQ=255)
                                                                  ataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtg     <  1:3191144/70‑1 (MQ=255)
                                                                   taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCtgt                              <  1:2586662/44‑1 (MQ=255)
                                                                   taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGaa                         >  1:422970/1‑49 (MQ=255)
                                                                   taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtg   >  1:1669244/1‑71 (MQ=255)
                                                                     gCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtgt  >  1:2683894/1‑70 (MQ=255)
                                                                      |                                                                    
ACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGT  >  minE/835722‑835860

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: