Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,209 G→A 23.7% V1223I (GTT→ATT) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,2090GA23.7% 206.3 / 60.5 114V1223I (GTT→ATT) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (40/47);  new base A (16/11);  total (56/58)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 130 total aligned reads displayed.

CGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG  >  minE/835142‑835278
                                                                   |                                                                     
cggcTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                                                    <  1:809850/71‑1 (MQ=255)
    tATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGt                                                                >  1:235901/1‑71 (MQ=255)
      tGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttg                                                                       <  1:4504329/62‑1 (MQ=255)
         ccATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAg                                                                 >  1:4361515/1‑65 (MQ=255)
              tGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                                                    >  1:3993820/1‑57 (MQ=255)
                 cATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaag                                                   <  1:3865259/71‑1 (MQ=255)
                   tACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                                                    <  1:4648466/52‑1 (MQ=255)
                   tACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                                                    <  1:4254182/52‑1 (MQ=255)
                    aCGGCTATAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaagaa                                                 >  1:730620/1‑70 (MQ=255)
                     cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTa                                                                  >  1:4306227/1‑52 (MQ=255)
                     cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTa                                                                  >  1:5170029/1‑52 (MQ=255)
                     cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                                               >  1:5614683/1‑71 (MQ=255)
                     cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                                               <  1:442185/71‑1 (MQ=255)
                     cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                                               <  1:5282743/71‑1 (MQ=255)
                      ggCTTTAACTACTACGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAAcgt                                              <  1:2579900/71‑3 (MQ=255)
                       gCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACAtt                                             <  1:2893467/71‑1 (MQ=255)
                        cTTTAACTACTACGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGtta                                            <  1:3280931/71‑2 (MQ=255)
                            aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACt                                         <  1:3100500/70‑1 (MQ=255)
                            aaCTACTACGGGACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATgaagaa                                                 >  1:5171552/1‑62 (MQ=255)
                             actactATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                                                    <  1:1108825/42‑1 (MQ=255)
                             actactATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                                                    <  1:3178041/42‑1 (MQ=255)
                             actactATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACAt                                              >  1:3442705/1‑64 (MQ=255)
                             actactATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTgg                                       <  1:4961804/71‑1 (MQ=255)
                              ctactaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTgg                                       <  1:4001945/70‑1 (MQ=255)
                              ctactaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGt                                      <  1:2812972/71‑1 (MQ=255)
                              ctactaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGt                                      <  1:105446/71‑1 (MQ=255)
                              ctactaCGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTggc                                      <  1:3275169/71‑2 (MQ=255)
                               tactaTGGCACCGTAGGTTCTAACCGCGATGAGTTtgttgt                                                                   <  1:4172156/41‑1 (MQ=255)
                                 ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGt                                      >  1:4961987/1‑68 (MQ=255)
                                 ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGt                                      >  1:5368718/1‑68 (MQ=255)
                                   acggCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGAt                                 >  1:256060/3‑71 (MQ=255)
                                   acggCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAACAACGTTAACTGGCTGGAt                                 >  1:359637/3‑71 (MQ=255)
                                    tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaa                                                    >  1:1035745/1‑51 (MQ=255)
                                    tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaa                                                    >  1:1552395/1‑51 (MQ=255)
                                    tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaa                                                    >  1:3205365/1‑51 (MQ=255)
                                    tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGa                                  <  1:1782557/69‑1 (MQ=255)
                                    tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGa                                  <  1:1103606/69‑1 (MQ=255)
                                    tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGa                                  <  1:1230388/69‑1 (MQ=255)
                                    tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGa                                  <  1:5171278/69‑1 (MQ=255)
                                    tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATg                                >  1:1344664/1‑71 (MQ=255)
                                    tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATg                                >  1:5242741/1‑71 (MQ=255)
                                    tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATg                                >  1:1580240/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:3390824/69‑1 (MQ=255)
                                    cggCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGAt                                 <  1:1590219/69‑1 (MQ=255)
                                    cggCACCGTCGGATCGAACCGCGATGAGCTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGAt                                 >  1:283218/2‑70 (MQ=255)
                                     ggCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGATga                               <  1:5684417/71‑2 (MQ=255)
                                      gCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGATgat                              >  1:3678124/1‑69 (MQ=255)
                                      gCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGATgat                              >  1:227615/1‑69 (MQ=255)
                                      gCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGATga                               >  1:5190126/1‑69 (MQ=255)
                                      gCACCGTCGGATCGAACCGCGATGACTTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGATgat                              <  1:733046/71‑3 (MQ=255)
                                      gCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGAc                                          >  1:3683516/1‑59 (MQ=255)
                                      gCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATgg                               <  1:5209487/70‑1 (MQ=255)
                                       cACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGATgat                              >  1:3810990/1‑68 (MQ=255)
                                       cACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGc                                                              <  1:1586383/38‑1 (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)
                                          cGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTgg                                       <  1:3136321/58‑1 (MQ=255)
                                          cGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGt                                                            <  1:2258041/37‑1 (MQ=255)
                                           gTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGa                                                        >  1:4806926/1‑40 (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)
                                               gTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACAtt                                             <  1:30916/47‑1 (MQ=255)
                                               gTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTa                                    <  1:1946457/56‑1 (MQ=255)
                                               gTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTa                                    <  1:5485123/56‑1 (MQ=255)
                                               gTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAAt                     <  1:959522/71‑1 (MQ=255)
                                                ttCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGa                            <  1:1123275/63‑1 (MQ=255)
                                                ttCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGa                            <  1:3286952/63‑1 (MQ=255)
                                                 tCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGa                                  >  1:4934101/1‑56 (MQ=255)
                                                 tCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGc                              <  1:4522556/60‑1 (MQ=255)
                                                 tCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGc                              <  1:1346366/60‑1 (MQ=255)
                                                 tCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGa                   >  1:785147/1‑71 (MQ=255)
                                                  cTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc                  >  1:2503013/1‑71 (MQ=255)
                                                   tAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTa                                    <  1:628252/52‑1 (MQ=255)
                                                   tAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTa                                    <  1:4870463/52‑1 (MQ=255)
                                                   tAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc                  >  1:4231039/1‑70 (MQ=255)
                                                    aaCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc                  >  1:5509325/1‑69 (MQ=255)
                                                     aCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGa                                  <  1:4227824/52‑1 (MQ=255)
                                                     aCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGa                   <  1:2149009/67‑1 (MQ=255)
                                                     aCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAcc                 <  1:5054707/69‑1 (MQ=255)
                                                     aCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc                  <  1:1877255/68‑1 (MQ=255)
                                                     aCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc                  <  1:1694742/68‑1 (MQ=255)
                                                     aCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc                  <  1:4592101/68‑1 (MQ=255)
                                                     aCCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAAc                                          >  1:657043/1‑44 (MQ=255)
                                                     aCCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAAc                                          >  1:3547809/1‑44 (MQ=255)
                                                      ccGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGATg                                <  1:2007834/53‑1 (MQ=255)
                                                        gcgATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTa            <  1:3591986/71‑1 (MQ=255)
                                                         cgATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTa            >  1:3728797/1‑70 (MQ=255)
                                                                tttGTTGTAGTGCGTAAGATGAAGAACATTGACTg                                        <  1:3443788/35‑1 (MQ=255)
                                                                  tgttgtAGTGCGTAAGATGAAGAACATTGACTGGTTAg                                   >  1:5434985/1‑38 (MQ=255)
                                                                  tgttgtAGTGCGTAAGATGAAGAACATTGACTGGTTAg                                   >  1:786608/1‑38 (MQ=255)
                                                                  tgttgtAGTGCGTAAGATGAAGAACATTGACTGGTTAGAt                                 >  1:4855383/1‑40 (MQ=255)
                                                                  tgttgtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCaa                      >  1:4375949/1‑51 (MQ=255)
                                                                  tgttgtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCg  >  1:3439300/1‑71 (MQ=255)
                                                                   |                                                                     
CGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG  >  minE/835142‑835278

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: