Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,178 T→C 6.2% Y1212Y (TAT→TAC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,1780TC 6.2% 279.3 / 8.2 96Y1212Y (TAT→TACnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (47/43);  new base C (3/3);  total (50/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

CCCGTATCACGCCGAAACCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTA  >  minE/835108‑835244
                                                                      |                                                                  
cccGTATCACGCCGAAACCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTAt                                                                    >  1:680478/1‑71 (MQ=255)
   gTATCACGCCGAAACCGGCGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATgg                                                                  >  1:1359461/1‑70 (MQ=255)
   gTATCACGCCGAAACCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATgg                                                                  >  1:2527176/1‑70 (MQ=255)
    tATCACGCCGAAACCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCa                                                                >  1:717722/1‑71 (MQ=255)
      tCACGCCGAAACCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCAcc                                                              >  1:1245001/1‑71 (MQ=255)
        aCGCCGAAACCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGt                                                            <  1:3322285/71‑1 (MQ=255)
          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:1762871/1‑71 (MQ=255)
                           tGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAg                                         <  1:3793958/71‑1 (MQ=255)
                           tGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAg                                         >  1:3108269/1‑71 (MQ=255)
                             aTCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGt                                                            <  1:1800052/50‑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)
                                       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: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)
                                              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: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:601140/1‑69 (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:1816317/71‑1 (MQ=255)
                                                       cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa             >  1:2788969/1‑71 (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:1446072/1‑71 (MQ=255)
                                                              aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTg      >  1:3964302/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: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)
                                                                      |                                                                  
CCCGTATCACGCCGAAACCGACGCATATGATCGGCGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTA  >  minE/835108‑835244

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: