Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,193 T→G 19.0% S1217S (TCT→TCG narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,1930TG19.0% 167.0 / 51.7 100S1217S (TCT→TCGnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (42/39);  new base G (15/4);  total (57/43)
Fisher's exact test for biased strand distribution p-value = 4.01e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 101 total aligned reads displayed.

ACCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATG  >  minE/835124‑835260
                                                                     |                                                                   
aCCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTa                                                                    <  1:2975098/71‑1 (MQ=255)
 ccGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTaa                                                                   >  1:524933/1‑71 (MQ=255)
 ccGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGACTTTAACTACTATGGCACCGTAGGTTCTaa                                                                   <  1:2054333/71‑1 (MQ=255)
    aCGCATATGATTGGCGGTTACGCGCAGCTGGCCTGGGGCTTTAACTACTACGGCACCGTCGGATCGaa                                                                   >  1:821054/1‑68 (MQ=255)
    aCGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAAcc                                                                 <  1:1736336/70‑1 (MQ=255)
    aCGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACcg                                                                >  1:1386857/1‑71 (MQ=255)
     cGCATATGATTGGCGGTTACGCGCAGCTGGCCTGGGGCTTTAACTACTACGGCACCGTCGGATCGAACcgc                                                               >  1:1024265/1‑71 (MQ=255)
     cGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACcgc                                                               >  1:1413148/1‑71 (MQ=255)
     cGCATATGATAGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCAACGTAGGTTCTAACcgc                                                               >  1:2120364/1‑71 (MQ=255)
        atatGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCGCCGTAGGTTCTAACCGCGAt                                                            <  1:3111744/71‑1 (MQ=255)
         tatGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGAt                                                            >  1:2326822/1‑70 (MQ=255)
          atGATTGGCGGTTACGCGCAGCTGGCCTGGGGCTTTAACTACTACGGCACCGTCGGATCGAACCGCGAt                                                            >  1:2767866/1‑69 (MQ=255)
           tGATTGGCGGTTACGCGCAGCTGGCCTGGGGCTTTAACTACTACGGCACCGTCGGATCGAACCGCGATGAg                                                         >  1:2430583/1‑71 (MQ=255)
           tGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAg                                                         >  1:2385409/1‑71 (MQ=255)
           tGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAg                                                         >  1:3775995/1‑71 (MQ=255)
           tGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAg                                                         >  1:2004842/1‑71 (MQ=255)
           tGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAg                                                         >  1:1928400/1‑71 (MQ=255)
             attggcggtTACGCGCAGCTGGCCTGGGGCTTTAACTACTACGGCACCGTCGGCTCGAACCGCGATGAGtt                                                       >  1:917970/4‑71 (MQ=255)
              tCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGttt                                                      <  1:2406389/71‑1 (MQ=255)
                ggcggcTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGa                                                          >  1:1534694/1‑65 (MQ=255)
                ggcggcTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgt                                                    <  1:1211359/71‑1 (MQ=255)
                ggcggcTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgt                                                    <  1:467046/71‑1 (MQ=255)
                ggcggcTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttg                                                     <  1:2815226/70‑1 (MQ=255)
                 gcggcTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGAt                                                            <  1:397028/62‑1 (MQ=255)
                  cggcTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                                  <  1:3315226/71‑1 (MQ=255)
                    gcTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGTTTCTAACCGCGATGAGtt                                                       <  1:2175920/64‑1 (MQ=255)
                    gcTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                                  >  1:1368259/1‑69 (MQ=255)
                    gcTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTa                                                <  1:3072486/71‑1 (MQ=255)
                     cTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAAcc                                                                 <  1:2392088/53‑1 (MQ=255)
                     cTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACcg                                                                <  1:2689045/54‑1 (MQ=255)
                      tATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGt                                              >  1:577755/1‑71 (MQ=255)
                       aTGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACcgcg                                                              >  1:2537532/1‑54 (MQ=255)
                       aTGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTg                                             >  1:1307414/1‑71 (MQ=255)
                          cccATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGAt                                                            >  1:3937580/1‑53 (MQ=255)
                          cccATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTtgttgt                                                 <  1:3111698/64‑1 (MQ=255)
                           ccATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGc                                            >  1:1738706/1‑68 (MQ=255)
                           ccATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGc                                            >  1:642834/1‑68 (MQ=255)
                           ccATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTa                                         <  1:1495627/71‑1 (MQ=255)
                            cATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTa                                                <  1:1549713/63‑1 (MQ=255)
                             aTCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTaa                                                                   >  1:730364/1‑43 (MQ=255)
                             aTCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTtgttgt                                                 <  1:3531113/61‑1 (MQ=255)
                             aTCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGt                                          >  1:470160/1‑68 (MQ=255)
                               cTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAg                                               >  1:1887257/1‑61 (MQ=255)
                                tGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTtgttgt                                                 >  1:175682/1‑58 (MQ=255)
                                tGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGt                                              >  1:179754/1‑61 (MQ=255)
                                 ggCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATg                                                           <  1:1965582/47‑1 (MQ=255)
                                  gCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                                  <  1:2804415/55‑1 (MQ=255)
                                  gCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCg                                           >  1:1782375/1‑62 (MQ=255)
                                  gCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATg                                    <  1:3848771/69‑1 (MQ=255)
                                    aTACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAg                                               >  1:285395/1‑56 (MQ=255)
                                     tACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                                  <  1:1845858/52‑1 (MQ=255)
                                     tACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaagaa                               >  1:3157522/1‑71 (MQ=255)
                                     tACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaaga                                >  1:2413535/1‑70 (MQ=255)
                                      aCGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaagaa                               >  1:2662020/1‑70 (MQ=255)
                                       cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAAc                              <  1:1441248/70‑1 (MQ=255)
                                       cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                             >  1:3578172/1‑71 (MQ=255)
                                        ggCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                             >  1:1272169/1‑70 (MQ=255)
                                        ggCTTTAACTACTACGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAAcgt                            >  1:1695585/1‑69 (MQ=255)
                                          cTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTtgttgt                                                 <  1:2013526/48‑1 (MQ=255)
                                          cTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACAtt                           >  1:3558716/1‑70 (MQ=255)
                                          cTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACAtt                           >  1:2849592/1‑70 (MQ=255)
                                          cTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTg                          <  1:3580214/71‑1 (MQ=255)
                                          cTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTg                          <  1:1633377/71‑1 (MQ=255)
                                          cTTTAACTACTACGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGtta                          <  1:297796/71‑2 (MQ=255)
                                            ttAACTACTATGGCACCGTAGGTTCTAACCGCGATGa                                                          >  1:1594689/1‑37 (MQ=255)
                                            ttAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGt                                              >  1:1880774/1‑49 (MQ=255)
                                            ttAACTACTACGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAAc                        >  1:110190/1‑71 (MQ=255)
                                            ttAACTACTACGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAAc                        >  1:2657265/1‑71 (MQ=255)
                                            ttAACTACTACGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAAc                        >  1:2248875/1‑71 (MQ=255)
                                            atAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTa                                                >  1:314686/2‑47 (MQ=255)
                                             tAACTACTACGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACt                       >  1:3348933/1‑71 (MQ=255)
                                              aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTg                      >  1:48316/1‑71 (MQ=255)
                                               actactACGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAAcg                             <  1:1522885/63‑2 (MQ=255)
                                               actactACGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTgg                     <  1:914913/71‑1 (MQ=255)
                                               actactACGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTgg                     >  1:1627332/1‑71 (MQ=255)
                                                ctactaCGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTggc                    >  1:3641095/1‑70 (MQ=255)
                                                   ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTa                  <  1:3456589/70‑1 (MQ=255)
                                                   ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTa                  <  1:1907475/70‑1 (MQ=255)
                                                   ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTa                  <  1:2382585/70‑1 (MQ=255)
                                                   ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTa                  <  1:3220800/70‑1 (MQ=255)
                                                   ctaCGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTggctg                  <  1:1429012/70‑4 (MQ=255)
                                                   ctaCGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTggctg                  >  1:2086464/1‑67 (MQ=255)
                                                     acggCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGAt               >  1:2152391/3‑71 (MQ=255)
                                                      tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAAc                              >  1:3636905/1‑55 (MQ=255)
                                                      tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGa                         >  1:2879013/1‑60 (MQ=255)
                                                      tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGa                         >  1:1677448/1‑60 (MQ=255)
                                                      tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGa                         >  1:543449/1‑60 (MQ=255)
                                                      tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGa                         >  1:3271910/1‑60 (MQ=255)
                                                      tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGAt               <  1:1095397/70‑1 (MQ=255)
                                                        gCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATg                                    >  1:3240278/1‑47 (MQ=255)
                                                         cACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGAt                                     <  1:1714743/45‑1 (MQ=255)
                                                           ccGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAGGAACATTGACTGGTTAGATGGCGaa         <  1:487614/71‑1 (MQ=255)
                                                           ccGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                             <  1:3599427/51‑1 (MQ=255)
                                                           ccGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGaa         <  1:725893/71‑1 (MQ=255)
                                                              tAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGc      <  1:464788/71‑1 (MQ=255)
                                                                 gTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGt                    <  1:1175320/54‑1 (MQ=255)
                                                                  ttCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaagaa                               <  1:1097543/42‑1 (MQ=255)
                                                                  ttCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATg  <  1:1214518/71‑1 (MQ=255)
                                                                   tCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATg              >  1:1144809/1‑58 (MQ=255)
                                                                     tAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAAt   >  1:82640/1‑67 (MQ=255)
                                                                     |                                                                   
ACCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATG  >  minE/835124‑835260

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: