Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,187 A→C 6.9% V1215V (GTA→GTC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,1870AC 6.9% 265.6 / 15.7 102V1215V (GTA→GTCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (48/47);  new base C (4/3);  total (52/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

Only 100 of 102 total aligned reads displayed.

GCCGAAACCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGG  >  minE/835118‑835255
                                                                     |                                                                    
gCCGAAACCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAg                                                                     <  1:1039256/71‑1 (MQ=255)
gCCGAAACCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAg                                                                     <  1:675757/71‑1 (MQ=255)
     aaCCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCt                                                                <  1:1006126/71‑1 (MQ=255)
       ccGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTaa                                                              <  1:3191659/71‑1 (MQ=255)
        cGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAAc                                                             <  1:2074633/71‑1 (MQ=255)
           cGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACcgc                                                          <  1:199442/71‑1 (MQ=255)
           cGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACcgc                                                          <  1:357503/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:1762871/1‑71 (MQ=255)
                 tGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAg                                                    >  1:3108269/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:671076/1‑70 (MQ=255)
                         ggcTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                             >  1:3279475/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)
                             aTGCCCATCTGGCATACGGCTTAAACTACTATGGCACCGTAg                                                                     >  1:13495/1‑42 (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: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:2334095/52‑1 (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)
                                             cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaagaa                          >  1:601140/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)
                                             cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                        >  1:2788969/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)
                                                cTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTg                     <  1:1492233/71‑1 (MQ=255)
                                                 tttAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTtgttgt                                            <  1:2800791/47‑1 (MQ=255)
                                                    aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgt                                               >  1:3011530/1‑41 (MQ=255)
                                                    aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgt                                               >  1:157987/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:761130/1‑59 (MQ=255)
                                                        actaCGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAAcg                        >  1:2101860/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:1400548/70‑1 (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)
                                                                 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)
                                                                     |                                                                    
GCCGAAACCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGG  >  minE/835118‑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: