Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,217 G→C 22.7% V1225V (GTG→GTC) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,2170GC22.7% 226.9 / 63.1 119V1225V (GTG→GTC) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (41/51);  new base C (16/11);  total (57/62)
Fisher's exact test for biased strand distribution p-value = 1.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

CCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGA  >  minE/835150‑835286
                                                                   |                                                                     
cccATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGt                                                                    <  1:140015/71‑1 (MQ=255)
         cATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaag                                                           <  1:3865259/71‑1 (MQ=255)
         cATACGGCTTTAACTACTATGGCACCGTAGGTCCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaag                                                           <  1:3211259/71‑1 (MQ=255)
            ggggCTTTAACTACTACGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAAc                                                        >  1:4667421/3‑71 (MQ=255)
            aCGGCTATAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaagaa                                                         >  1:730620/1‑70 (MQ=255)
             cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                                                       <  1:442185/71‑1 (MQ=255)
             cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                                                       >  1:5614683/1‑71 (MQ=255)
             cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                                                       <  1:5282743/71‑1 (MQ=255)
              ggCTTTAACTACTACGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAAcgt                                                      <  1:2579900/71‑3 (MQ=255)
               gCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACAtt                                                     <  1:2893467/71‑1 (MQ=255)
                cTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACAtt                                                     >  1:3862488/1‑70 (MQ=255)
                 tttAACTACTACGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGtt                                                     <  1:2956660/69‑1 (MQ=255)
                    aaCTACTACGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATgaagaa                                                         >  1:3456343/1‑62 (MQ=255)
                     actactATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACAt                                                      >  1:3442705/1‑64 (MQ=255)
                     actactATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACt                                                 <  1:5176581/69‑1 (MQ=255)
                     actactATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTgg                                               <  1:4961804/71‑1 (MQ=255)
                      ctactaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGt                                              <  1:105446/71‑1 (MQ=255)
                      ctactaCGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTggc                                              <  1:3275169/71‑2 (MQ=255)
                         ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGt                                              >  1:4961987/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)
                            tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                                                       <  1:5095584/56‑1 (MQ=255)
                            tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGa                                          <  1:1103606/69‑1 (MQ=255)
                            tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGa                                          <  1:1230388/69‑1 (MQ=255)
                            tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGa                                          <  1:5171278/69‑1 (MQ=255)
                            tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGa                                          <  1:1782557/69‑1 (MQ=255)
                            tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATg                                        >  1:1344664/1‑71 (MQ=255)
                            tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATg                                        >  1:5242741/1‑71 (MQ=255)
                            cggCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAAcg                                                       >  1:1312757/2‑55 (MQ=255)
                            cggCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGttaa                                                   >  1:2452135/2‑58 (MQ=255)
                            cggCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGAt                                         <  1:1590219/69‑1 (MQ=255)
                            cggCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGAt                                         <  1:3390824/69‑1 (MQ=255)
                            cggCACCGTCGGATCGAACCGCGATGAGCTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGAt                                         >  1:283218/2‑70 (MQ=255)
                             ggCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGATga                                       <  1:5684417/71‑2 (MQ=255)
                              gCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGATgat                                      >  1:2619988/1‑69 (MQ=255)
                              gCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGATgat                                      >  1:227615/1‑69 (MQ=255)
                              gCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGATga                                       >  1:5190126/1‑69 (MQ=255)
                              gCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGAc                                                  >  1:3683516/1‑59 (MQ=255)
                              gCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATgg                                       <  1:5495953/70‑1 (MQ=255)
                              gCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATgg                                       <  1:5209487/70‑1 (MQ=255)
                               cACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGATgat                                      >  1:3810990/1‑68 (MQ=255)
                               cACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGc                                                                      <  1:1586383/38‑1 (MQ=255)
                               cACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGTAGAACATTGACTGGTTAGATGGCg                                     <  1:2297451/71‑1 (MQ=255)
                               cACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAg                                           >  1:1898385/1‑65 (MQ=255)
                               cACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAg                                           >  1:5296291/1‑65 (MQ=255)
                               cACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCg                                     <  1:2319062/71‑1 (MQ=255)
                               cACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCg                                     <  1:2292372/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)
                                   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:5485123/56‑1 (MQ=255)
                                       gTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTa                                            <  1:1946457/56‑1 (MQ=255)
                                       gTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAAt                             <  1:959522/71‑1 (MQ=255)
                                       gTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAAt                             <  1:4239680/71‑1 (MQ=255)
                                        ttCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGa                                    <  1:3286952/63‑1 (MQ=255)
                                        ttCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGa                                    <  1:1123275/63‑1 (MQ=255)
                                         tCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGa                                          >  1:4039788/1‑56 (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)
                                          cTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc                          >  1:340032/1‑71 (MQ=255)
                                           tAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTa                                            <  1:4870463/52‑1 (MQ=255)
                                           tAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTa                                            <  1:628252/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: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)
                                                gcgATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTa                    >  1:2143009/1‑71 (MQ=255)
                                                 cgATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTa                    >  1:3728797/1‑70 (MQ=255)
                                                  gATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTa                    >  1:5679238/1‑69 (MQ=255)
                                                        tttGTTGTAGTGCGTAAGATGAAGAACATTGACTg                                                <  1:4366415/35‑1 (MQ=255)
                                                         ttgttgTAGTGCGTAAGATGAAGAACATTGACTGGTTAGa                                          <  1:1614334/40‑1 (MQ=255)
                                                          tgttgtAGTGCGTAAGATGAAGAACATTGACTGGTTAg                                           >  1:5434985/1‑38 (MQ=255)
                                                          tgttgtAGTGCGTAAGATGAAGAACATTGACTGGTTAg                                           >  1:786608/1‑38 (MQ=255)
                                                          tgttgtAGTGCGTAAGATGAAGAACATTGACTGGTTAGAt                                         >  1:4855383/1‑40 (MQ=255)
                                                          tgttgtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCg          >  1:3439300/1‑71 (MQ=255)
                                                              gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGAt                                         >  1:3267547/1‑36 (MQ=255)
                                                              gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGAt                                         >  1:4466483/1‑36 (MQ=255)
                                                               tAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTaaaa     <  1:3577390/71‑1 (MQ=255)
                                                                 gTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGaa                                   >  1:4810196/1‑39 (MQ=255)
                                                                  tGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAcc                         <  1:795173/48‑1 (MQ=255)
                                                                  tGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGa  <  1:2664627/71‑1 (MQ=255)
                                                                   |                                                                     
CCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGA  >  minE/835150‑835286

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: