Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,187 A→C 31.0% 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,1870AC31.0% 124.7 / 134.0 142V1215V (GTA→GTCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (44/54);  new base C (25/19);  total (69/73)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

Only 100 of 142 total aligned reads displayed.

GCCGAAACCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGG  >  minE/835118‑835255
                                                                     |                                                                    
gCCGAAACCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAg                                                                     <  1:829241/71‑1 (MQ=255)
  cGAAACCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGt                                                                   >  1:188003/1‑71 (MQ=255)
   gAAACCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAg                                                                     >  1:1201454/1‑68 (MQ=255)
   gAAACCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGtt                                                                  <  1:2364718/71‑1 (MQ=255)
   gAAACCAACGCATATGATTGGCGGTTACGCGCAGCTGGCCTGGGGCTTTAACTACTACGGCACCGTCggat                                                                  >  1:1957960/1‑69 (MQ=255)
   gAAACCAACGCATATGATTGGCGGTTACGCGCAGCTGGCCTGGGGCTTTAACTACTACGGCACCGTCggat                                                                  >  1:1801658/1‑69 (MQ=255)
    aaaCCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTc                                                                 >  1:1058404/1‑71 (MQ=255)
     aaCCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCt                                                                <  1:5628827/71‑1 (MQ=255)
     aaCCAACGCATATGATTGGCGGTTACGCGCAGCTGGCCTGGGGCTTTAACTACTACGGCACCGTCGGATcg                                                                <  1:2752210/71‑2 (MQ=255)
     aaCCAACGCATATGATTGGCGGTTACGCGCAGCTGGCCTGGGGCTTTAACTACTACGGCACCGTCGGATc                                                                 >  1:2485470/1‑70 (MQ=255)
      aCCAACGCATATGATTGGCGGTTACGCGCAGCTGGCCTGGGGCTTTAACTACTACGGCACCGTCGGATcga                                                               <  1:1155060/71‑3 (MQ=255)
       ccGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTaa                                                              <  1:3133452/71‑1 (MQ=255)
           cGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACcgc                                                          <  1:3139153/71‑1 (MQ=255)
              atatGATTGGCGGTTACGCGCAGCTGGCCTGGGGCTTTAACTACTACGGCACCGTCGGATCGAACCGCGAt                                                       <  1:2074863/71‑1 (MQ=255)
              atatGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGAt                                                       <  1:1623867/71‑1 (MQ=255)
               tatGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCt                                                                <  1:2636473/61‑1 (MQ=255)
               tatGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCt                                                                <  1:1760513/61‑1 (MQ=255)
               tatGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATg                                                      <  1:387810/71‑1 (MQ=255)
               tatGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATg                                                      <  1:5284471/71‑1 (MQ=255)
                atGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCt                                                                <  1:3015777/60‑1 (MQ=255)
                 tGATTGGCGGTTACGCGCAGCTGGCCTGGGGCTTTAACTACTACGGCACCGTCGGATCGAACCGCGAt                                                       >  1:480016/1‑68 (MQ=255)
                 tGATTGGCGGTTACGCGCAGCTGGCCTGGGGCTTTAACTACTACGGCACCGTCGGATCGAACCGCGATGa                                                     <  1:3098752/70‑1 (MQ=255)
                 tGATTGGCGGTTACGCGCAGCTGGCCTGGGGCTTTAACTACTACGGCACCGTCGGATCGAACCGCGATGa                                                     >  1:3248356/1‑70 (MQ=255)
                 tGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAg                                                    <  1:5188796/71‑1 (MQ=255)
                  gATTGGCGGTTACGCGCAGCTGGCCTGGGGCTTTAACTACTACGGCACCGTCGGATCGAACCGCGATGAGt                                                   >  1:805842/1‑71 (MQ=255)
                  gATTGGCGGTTACGCGCAGCTGGCCTGGGGCTTTAACTACTACGGCACCGTCGGATCGAACCGCGATGAGt                                                   >  1:1260521/1‑71 (MQ=255)
                  gATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGt                                                   >  1:2967366/1‑71 (MQ=255)
                    ttggcggtTACGCGCAGCTGGCCTGGGGCTTTAACTACTACGGCACCGTCGGATCGAACCGCGATGAGttc                                                 <  1:4466787/69‑2 (MQ=255)
                     tggcggtTACGCGCAGCTGGCCTGGGGCTTTAACTACTACGGCACCGTCGGATCGAACCGCGATGAGttca                                                >  1:4320174/2‑69 (MQ=255)
                        cggcTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                             <  1:3619986/71‑1 (MQ=255)
                           cTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTa                                           >  1:3173179/1‑70 (MQ=255)
                            tATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGt                                         >  1:235901/1‑71 (MQ=255)
                             aTGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTATGTTCTAACcgcg                                                         <  1:3713206/54‑1 (MQ=255)
                             aTGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACcgcg                                                         <  1:324015/54‑1 (MQ=255)
                              tGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttg                                                <  1:4504329/62‑1 (MQ=255)
                                cccATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCt                                                                <  1:2125127/44‑1 (MQ=255)
                                cccATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCt                                                                <  1:4788922/44‑1 (MQ=255)
                                cccATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTaa                                                              >  1:725216/1‑46 (MQ=255)
                                 ccATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAg                                          >  1:4361515/1‑65 (MQ=255)
                                 ccATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAg                                          >  1:4510044/1‑65 (MQ=255)
                                      tGGCATCCGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAg                                          >  1:731687/1‑60 (MQ=255)
                                      tGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                             >  1:3993820/1‑57 (MQ=255)
                                         cATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaag                            <  1:3865259/71‑1 (MQ=255)
                                           tACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                             <  1:4254182/52‑1 (MQ=255)
                                            ggggCTTTAACTACTACGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAAc                         >  1:4667421/3‑71 (MQ=255)
                                            aCGGCTATAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaagaa                          >  1:730620/1‑70 (MQ=255)
                                             cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGttt                                                 >  1:2608722/1‑46 (MQ=255)
                                             cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGttt                                                 >  1:1959672/1‑46 (MQ=255)
                                             cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTa                                           >  1:5170029/1‑52 (MQ=255)
                                             cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTa                                           >  1:787744/1‑52 (MQ=255)
                                             cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                        >  1:5614683/1‑71 (MQ=255)
                                             cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                        <  1:442185/71‑1 (MQ=255)
                                             cGGCTTTAACTACTATGGCACCCTAGGTTCTAACCGCGa                                                        >  1:5509859/1‑39 (MQ=255)
                                              ggCTTTAACTACTACGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAAcgt                       <  1:2579900/71‑3 (MQ=255)
                                               gCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACAtt                      <  1:2893467/71‑1 (MQ=255)
                                                  ttAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGtt                                                  <  1:4115027/40‑1 (MQ=255)
                                                  ttAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGtt                                                  <  1:4224053/40‑1 (MQ=255)
                                                    aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACt                  <  1:3100500/70‑1 (MQ=255)
                                                    aaCTACTACGGGACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATgaagaa                          >  1:5171552/1‑62 (MQ=255)
                                                    aaCTACTACGGCACCGTCGGATCGAACCGCGATGAGttc                                                 >  1:1521123/1‑38 (MQ=25)
                                                     actactATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                             <  1:1108825/42‑1 (MQ=255)
                                                     actactATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACAt                       >  1:3442705/1‑64 (MQ=255)
                                                     actactATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACt                  <  1:5176581/69‑1 (MQ=255)
                                                     actactATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTgg                <  1:4961804/71‑1 (MQ=255)
                                                      ctactaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTgg                <  1:4001945/70‑1 (MQ=255)
                                                      ctactaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGt               <  1:105446/71‑1 (MQ=255)
                                                      ctactaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGt               <  1:2812972/71‑1 (MQ=255)
                                                      ctactaCGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTggc               <  1:3275169/71‑2 (MQ=255)
                                                         ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGt               >  1:5368718/1‑68 (MQ=255)
                                                           acggCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAACAACGTTAACTGGCTGGAt          >  1:359637/3‑71 (MQ=255)
                                                            tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaa                             >  1:1035745/1‑51 (MQ=255)
                                                            tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaa                             >  1:3205365/1‑51 (MQ=255)
                                                            tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaa                             >  1:1552395/1‑51 (MQ=255)
                                                            tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                        <  1:5095584/56‑1 (MQ=255)
                                                            tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGa           <  1:1103606/69‑1 (MQ=255)
                                                            tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGa           <  1:1230388/69‑1 (MQ=255)
                                                            tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGa           <  1:1782557/69‑1 (MQ=255)
                                                            tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATg         >  1:1580240/1‑71 (MQ=255)
                                                            tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATg         >  1:5242741/1‑71 (MQ=255)
                                                            tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATg         >  1:1344664/1‑71 (MQ=255)
                                                            cggCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAAcg                        >  1:1312757/2‑55 (MQ=255)
                                                            cggCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGttaa                    >  1:5650863/2‑58 (MQ=255)
                                                            cggCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGttaa                    >  1:2452135/2‑58 (MQ=255)
                                                            cggCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGAt          <  1:4956914/69‑1 (MQ=255)
                                                            cggCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGAt          <  1:3390824/69‑1 (MQ=255)
                                                            cggCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGAt          <  1:1590219/69‑1 (MQ=255)
                                                              gCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGATgat       >  1:227615/1‑69 (MQ=255)
                                                              gCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGATgat       >  1:3678124/1‑69 (MQ=255)
                                                              gCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATgg        <  1:5209487/70‑1 (MQ=255)
                                                               cACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGATgat       >  1:3810990/1‑68 (MQ=255)
                                                               cACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAg            >  1:1898385/1‑65 (MQ=255)
                                                               cACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAg            >  1:5296291/1‑65 (MQ=255)
                                                               cACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCg      <  1:2292372/71‑1 (MQ=255)
                                                               cACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCg      <  1:2319062/71‑1 (MQ=255)
                                                                aCCGTAGGTTGTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCg      <  1:4327716/70‑1 (MQ=255)
                                                                 ccGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCg      >  1:2935844/1‑69 (MQ=255)
                                                                   gTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGGAGAACATTGACTGGTTa             >  1:2039359/1‑60 (MQ=255)
                                                                   gTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTa             >  1:1116169/1‑60 (MQ=255)
                                                                   gTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAgg  >  1:932341/1‑71 (MQ=255)
                                                                    tAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                        >  1:3197425/1‑48 (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: