Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,190 T→A 7.1% G1216G (GGT→GGA narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,1900TA 7.1% 278.7 / 11.5 99G1216G (GGT→GGAnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (46/46);  new base A (4/3);  total (50/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

AACCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGG  >  minE/835123‑835255
                                                                   |                                                                 
aaCCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCt                                                                <  1:1006126/71‑1 (MQ=255)
  ccGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTaa                                                              <  1:3191659/71‑1 (MQ=255)
   cGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAAc                                                             <  1:2074633/71‑1 (MQ=255)
      cGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACcgc                                                          <  1:357503/71‑1 (MQ=255)
      cGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACcgc                                                          <  1:199442/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)
            tGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGtt                                                                  >  1:436609/1‑57 (MQ=255)
            tGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTaa                                                              >  1:3481785/1‑61 (MQ=255)
            tGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAg                                                    <  1:3793958/71‑1 (MQ=255)
            tGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAg                                                    >  1:3108269/1‑71 (MQ=255)
            tGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAg                                                    >  1:1762871/1‑71 (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:3170432/1‑44 (MQ=255)
                              aTCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAAc                                                             >  1:2290763/1‑44 (MQ=255)
                              aTCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGt                                     >  1:1865410/1‑68 (MQ=255)
                               tCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttg                                                <  1:1606182/56‑1 (MQ=255)
                               tCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttg                                                <  1:2799776/56‑1 (MQ=255)
                               tCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCg                                      >  1:3939148/1‑66 (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:1560120/1‑59 (MQ=255)
                                 tGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTa                                           >  1:2168586/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:2334095/52‑1 (MQ=255)
                                      tACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                             <  1:1714691/52‑1 (MQ=255)
                                      tACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                             <  1:1260439/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:601140/1‑69 (MQ=255)
                                        cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaagaa                          >  1:3651660/1‑69 (MQ=255)
                                        cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAAc                         <  1:2453992/70‑1 (MQ=255)
                                        cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAAc                         <  1:3309420/70‑1 (MQ=255)
                                        cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                        >  1:2788969/1‑71 (MQ=255)
                                        cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                        <  1:1816317/71‑1 (MQ=255)
                                        cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                        >  1:1743056/1‑71 (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)
                                               aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgt                                               >  1:3011530/1‑41 (MQ=255)
                                               aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTa                                           >  1:584510/1‑45 (MQ=255)
                                               aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTg                 >  1:1524225/1‑71 (MQ=255)
                                               aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTg                 >  1:1446072/1‑71 (MQ=255)
                                               aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTg                 >  1:3964302/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:3308921/1‑68 (MQ=255)
                                                    ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGt               >  1:436350/1‑68 (MQ=255)
                                                    ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTa             <  1:3464078/70‑1 (MQ=255)
                                                    ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTa             <  1:1400548/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)
                                                          cACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCg      <  1:834316/71‑1 (MQ=255)
                                                            ccGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                        <  1:70176/51‑1 (MQ=255)
                                                            ccGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATg         <  1:2883082/66‑1 (MQ=255)
                                                            ccGTAGGTTCTAACCGCGATGAGTTGTTGTAGTGCGTAAGATGAAGAACa                        <  1:3281758/50‑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)
                                                                   |                                                                 
AACCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGG  >  minE/835123‑835255

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: