Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,193 T→G 6.8% 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,1930TG 6.8% 258.5 / 15.6 103S1217S (TCT→TCGnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (52/44);  new base G (4/3);  total (56/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.36e-01

Only 100 of 111 total aligned reads displayed.

AACCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAC  >  minE/835123‑835262
                                                                      |                                                                     
aaCCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCt                                                                       <  1:1006126/71‑1 (MQ=255)
  ccGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTaa                                                                     <  1:3191659/71‑1 (MQ=255)
   cGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAAc                                                                    <  1:2074633/71‑1 (MQ=255)
       gCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACcgcg                                                                <  1:355622/71‑1 (MQ=255)
          tatGATTGGCGGTTACGCGCAGCTGGCCTGGGGCTTTAACTACTACGGCACCGTCGGATCGAACcgc                                                                 <  1:2481663/67‑1 (MQ=255)
          tatGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCt                                                                       <  1:1229878/61‑1 (MQ=255)
          tatGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACGACTATGGCACCGTAGGTTCTAACCGCGAt                                                              >  1:3786174/1‑70 (MQ=255)
            tGATTGGAGGTTACGCGCAGCTGGCCTGGGGCTTTAACTACTACGGCACCGTCGGATCGAACCGCGATGAg                                                           >  1:640141/1‑71 (MQ=255)
            tGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTaa                                                                     >  1:3481785/1‑61 (MQ=255)
            tGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAg                                                           >  1:3108269/1‑71 (MQ=255)
            tGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAg                                                           <  1:3793958/71‑1 (MQ=255)
              aTCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGt                                                          <  1:691408/70‑1 (MQ=255)
                   cggcTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGTTTCTAACCGCGATGAGTttgttg                                                    <  1:245151/71‑1 (MQ=255)
                   cggcTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                                    <  1:2732991/71‑1 (MQ=255)
                   cggcTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGGttgttg                                                    <  1:6775/71‑1 (MQ=255)
                    ggcTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                                    >  1:3279475/1‑70 (MQ=255)
                    ggcTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                                    >  1:671076/1‑70 (MQ=255)
                    ggcTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgt                                                      >  1:1233875/1‑68 (MQ=255)
                    ggcTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTtgttgt                                                   <  1:3688752/71‑1 (MQ=255)
                       tATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGt                                                >  1:1902457/1‑71 (MQ=255)
                       tATGCCCACCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGt                                                >  1:3935296/1‑71 (MQ=255)
                        aTGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAAGCGCGATGAGTTTGTTGTAGTg                                               >  1:71732/1‑71 (MQ=255)
                        aTGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTg                                               >  1:3152808/1‑71 (MQ=255)
                         tGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCt                                                                       <  1:2876800/46‑1 (MQ=255)
                            ccATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTa                                           <  1:230824/71‑1 (MQ=255)
                             cATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAg                                                           <  1:1089940/54‑1 (MQ=255)
                              aTCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAAc                                                                    >  1:2290763/1‑44 (MQ=255)
                              aTCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAAc                                                                    >  1:3170432/1‑44 (MQ=255)
                              aTCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGt                                            >  1:1865410/1‑68 (MQ=255)
                               tCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttg                                                       <  1:2799776/56‑1 (MQ=255)
                               tCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttg                                                       <  1:1606182/56‑1 (MQ=255)
                                cTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGa                                        >  1:1912363/1‑70 (MQ=255)
                                 tGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                                    <  1:1387017/57‑1 (MQ=255)
                                 tGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTtgttgt                                                   <  1:2129653/58‑1 (MQ=255)
                                 tGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTa                                                  >  1:2168586/1‑59 (MQ=255)
                                 tGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTa                                                  >  1:1560120/1‑59 (MQ=255)
                                 tGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTaa                                          <  1:1844840/67‑1 (MQ=255)
                                 tGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGAt                                       >  1:192605/1‑70 (MQ=255)
                                 tGGCATACGGCTTTAACTACTATGGCACCGAAGGTTCTAACCGCGATGa                                                            >  1:3531384/1‑49 (MQ=255)
                                  ggCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                                    <  1:3085191/56‑1 (MQ=255)
                                   gCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGttt                                                        >  1:2605526/1‑51 (MQ=255)
                                   gCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGAt                                       <  1:2851865/68‑1 (MQ=255)
                                   gCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaa                                    <  1:604381/71‑1 (MQ=255)
                                    cATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaag                                   <  1:4072356/71‑1 (MQ=255)
                                     aTACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaaga                                  <  1:182280/71‑1 (MQ=255)
                                      tACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGa                                                               >  1:3629713/1‑41 (MQ=255)
                                      tACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                                    <  1:1260439/52‑1 (MQ=255)
                                      tACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                                    <  1:1714691/52‑1 (MQ=255)
                                        gggCTTTAACTACTACGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAAc                                <  1:2679682/69‑1 (MQ=255)
                                        cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATg                                                             >  1:3656587/1‑41 (MQ=255)
                                        cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGa                                                            >  1:558235/1‑42 (MQ=255)
                                        cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaagaa                                 >  1:3651660/1‑69 (MQ=255)
                                        cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAAc                                <  1:3309420/70‑1 (MQ=255)
                                        cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAAc                                <  1:2453992/70‑1 (MQ=255)
                                        cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                               >  1:1743056/1‑71 (MQ=255)
                                        cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                               <  1:1816317/71‑1 (MQ=255)
                                         ggCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaagaa                                 <  1:232103/68‑1 (MQ=255)
                                          gCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttg                                                       >  1:1984950/1‑45 (MQ=255)
                                          gCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACAtt                             <  1:674326/71‑1 (MQ=255)
                                          gCTTTAACTACTACGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGtt                             <  1:1072415/71‑1 (MQ=255)
                                           cTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                               >  1:2875070/1‑68 (MQ=255)
                                           cTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTg                            <  1:1492233/71‑1 (MQ=255)
                                            tttAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTtgttgt                                                   <  1:2800791/47‑1 (MQ=255)
                                               aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgt                                                      >  1:157987/1‑41 (MQ=255)
                                               aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTa                                                  >  1:584510/1‑45 (MQ=255)
                                               aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTg                        >  1:3964302/1‑71 (MQ=255)
                                               aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTg                        >  1:1446072/1‑71 (MQ=255)
                                               aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTg                        >  1:1524225/1‑71 (MQ=255)
                                                actactATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTgg                       >  1:3977686/1‑71 (MQ=255)
                                                 ctactaTGGCACCGTAGGTTCTAACCGCGATGAGTTtgttgt                                                   <  1:1107993/42‑1 (MQ=255)
                                                   actaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGt                      >  1:2601759/1‑69 (MQ=255)
                                                   actaCGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAAcg                               >  1:2101860/1‑59 (MQ=255)
                                                   actaCGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAAcg                               >  1:761130/1‑59 (MQ=255)
                                                    ctaTGGCACCGTAGGTTCTAACCGCGATGAGTttgt                                                      >  1:3940512/1‑36 (MQ=255)
                                                    ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGt                      >  1:436350/1‑68 (MQ=255)
                                                    ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGt                      >  1:3308921/1‑68 (MQ=255)
                                                    ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTa                    <  1:3464078/70‑1 (MQ=255)
                                                       tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATga                                     <  1:3062234/50‑1 (MQ=255)
                                                       tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAAc                                >  1:3293592/1‑55 (MQ=255)
                                                       cggCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATgaagaa                                 >  1:1825414/2‑54 (MQ=255)
                                                         gCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaagaa                                 >  1:3418545/1‑52 (MQ=255)
                                                          cACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCg             <  1:3299478/71‑1 (MQ=255)
                                                            ccGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                               <  1:70176/51‑1 (MQ=255)
                                                            ccGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATg                <  1:2883082/66‑1 (MQ=255)
                                                               tAGGTTCTAACCGCGATGAGTTTGTTGTATTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAgg         >  1:4022729/1‑70 (MQ=255)
                                                               tAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACt                         >  1:2933549/1‑54 (MQ=255)
                                                               tAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAgg         >  1:233002/1‑70 (MQ=255)
                                                                  gTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGa            <  1:2801162/64‑1 (MQ=255)
                                                                    tCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGa   >  1:2495360/1‑71 (MQ=255)
                                                                    tCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGa   >  1:2742322/1‑71 (MQ=255)
                                                                    tCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGa   >  1:3131654/1‑71 (MQ=255)
                                                                     cTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGa                  >  1:3203561/1‑55 (MQ=255)
                                                                     cTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGa   <  1:4099664/70‑1 (MQ=255)
                                                                     cTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc  >  1:1286771/1‑71 (MQ=255)
                                                                     cTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc  >  1:3929115/1‑71 (MQ=255)
                                                                      tAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACAtt                             <  1:2065101/43‑1 (MQ=255)
                                                                      tAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACt                         <  1:2998402/47‑1 (MQ=255)
                                                                      tAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGc        <  1:1604771/64‑1 (MQ=255)
                                                                      tAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGc        <  1:1360076/64‑1 (MQ=255)
                                                                      tAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATg    <  1:3862420/68‑1 (MQ=255)
                                                                      |                                                                     
AACCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAC  >  minE/835123‑835262

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: