Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,819 T→C 19.6% Y179Y (TAT→TAC)  narH → nitrate reductase 1, beta (Fe‑S) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,8190TC19.6% 211.8 / 46.8 102Y179Y (TAT→TAC) narHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base T (38/44);  new base C (13/7);  total (51/51)
Fisher's exact test for biased strand distribution p-value = 2.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

Only 100 of 103 total aligned reads displayed.

CGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATT  >  minE/835749‑835885
                                                                      |                                                                   
cGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTAt                                                                     <  1:2822654/71‑1 (MQ=255)
    aGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGc                                                                 >  1:1473338/1‑71 (MQ=255)
       cAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCcgc                                                              <  1:3433131/71‑1 (MQ=255)
         aGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGcc                                                                >  1:1396351/1‑67 (MQ=255)
         aGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGcc                                                                >  1:2299302/1‑67 (MQ=255)
         aGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCgcgc                                                            >  1:3482164/1‑71 (MQ=255)
         aGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCgcgc                                                            >  1:1582637/1‑71 (MQ=255)
          caaCTTCGAGGCCATGCAAAAGGAGATGTACGCGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGcc                                                           >  1:2376532/2‑71 (MQ=255)
           aaCTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCt                                                          <  1:1312012/71‑1 (MQ=255)
           aaCTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCt                                                          <  1:2457881/71‑1 (MQ=255)
            aCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGACAACACCTTCATGATGTACCTGCCGCGCCtg                                                         <  1:3120492/71‑1 (MQ=255)
            aCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCtg                                                         <  1:2015105/71‑1 (MQ=255)
            aCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCtg                                                         >  1:1546133/1‑71 (MQ=255)
            aCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCtg                                                         >  1:614351/1‑71 (MQ=255)
              ttCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCcgc                                                              <  1:1348976/64‑1 (MQ=255)
                cGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGc                                                      >  1:960591/1‑70 (MQ=255)
                 gACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCg                                                     <  1:3804234/70‑1 (MQ=255)
                  acaacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGc                                                      <  1:3075465/68‑1 (MQ=255)
                   ggccaTGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGaa                                                   >  1:403939/4‑70 (MQ=255)
                   caacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAac                                                  <  1:3166509/71‑1 (MQ=255)
                        tCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCgcgc                                                            <  1:1159264/56‑1 (MQ=255)
                        tCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGc                                             >  1:3526955/1‑71 (MQ=255)
                        tCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGc                                             >  1:2269706/1‑71 (MQ=255)
                        tCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGc                                             >  1:1714013/1‑71 (MQ=255)
                         ccAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGcc                                            >  1:903867/1‑71 (MQ=255)
                         ccAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGcc                                            >  1:304097/1‑71 (MQ=255)
                          cAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCt                                           >  1:608740/1‑71 (MQ=255)
                           aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg                                          <  1:489827/71‑1 (MQ=255)
                           aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg                                          <  1:940361/71‑1 (MQ=255)
                           aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg                                          <  1:945051/71‑1 (MQ=255)
                           aGAAGGCGAGGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCt                                           >  1:3215519/1‑70 (MQ=255)
                            gAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGa                                         <  1:374529/71‑1 (MQ=255)
                              aGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑c                                      >  1:1848602/1‑71 (MQ=255)
                              aGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑c                                      <  1:135266/71‑1 (MQ=255)
                                gCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑ccc                                    <  1:3876236/71‑1 (MQ=255)
                                 cGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCg                                   >  1:3273268/1‑71 (MQ=255)
                                  gATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAaca                                                 >  1:3508944/1‑57 (MQ=255)
                                    tGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGc                                             >  1:937686/1‑59 (MQ=255)
                                    tGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGc                                             >  1:2309667/1‑59 (MQ=255)
                                    tGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCa                                >  1:2081803/1‑71 (MQ=255)
                                    tGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCa                                >  1:2271441/1‑71 (MQ=255)
                                      tataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAac                                                  >  1:2661842/1‑52 (MQ=255)
                                      tataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtg                              >  1:2041892/1‑71 (MQ=255)
                                      tataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAt                               <  1:1593906/70‑1 (MQ=255)
                                        taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgtg                            >  1:3336679/1‑71 (MQ=255)
                                        taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgtg                            <  1:3181309/71‑1 (MQ=255)
                                        taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgtg                            <  1:1472040/71‑1 (MQ=255)
                                        taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgt                             <  1:3621570/70‑1 (MQ=255)
                                         aGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgtg                            <  1:3689200/70‑1 (MQ=255)
                                          gCCAGTTCGAAAACACCTTCATGATGTATTTGCcg                                                               >  1:1766361/1‑35 (MQ=255)
                                          gCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgtgt                           <  1:3867308/70‑1 (MQ=255)
                                          gCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgtgt                           <  1:3018493/70‑1 (MQ=255)
                                           ccAGTTCGAAAACACCTTCATGATGTATTTGCCgcgc                                                            <  1:2860980/37‑1 (MQ=255)
                                           ccAGTTCGAAAACACCTTCATGATGTATTTGCCgcgc                                                            <  1:3935445/37‑1 (MQ=255)
                                           ccAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTgg                         <  1:3039361/71‑1 (MQ=255)
                                            cAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGc                        <  1:3913462/71‑1 (MQ=255)
                                             aGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAA‑tccca                                  >  1:1090451/1‑58 (MQ=255)
                                             aGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTTCGAACACTGCCTGAA‑CCCGGCATGTGTGGc                        >  1:349795/1‑70 (MQ=255)
                                             aGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑cc                                     >  1:3218393/1‑57 (MQ=255)
                                             aGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGc                        >  1:2971718/1‑70 (MQ=255)
                                             aGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGc                        >  1:289658/1‑70 (MQ=255)
                                              gTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACt                                               <  1:2882346/47‑1 (MQ=38)
                                              gTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACt                                               <  1:1907587/47‑1 (MQ=38)
                                               ttCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGAc                     <  1:3779713/71‑1 (MQ=255)
                                               ttCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGAc                     <  1:1829041/71‑1 (MQ=255)
                                                tCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCg                                                     <  1:3163158/39‑1 (MQ=255)
                                                tCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCgg                                  >  1:2069830/1‑57 (MQ=255)
                                                 tgAAAACACCTTCATTATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCt                   >  1:741533/2‑71 (MQ=255)
                                                 tgAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCt                   >  1:2829823/2‑71 (MQ=255)
                                                 cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGAcc                    <  1:3935453/70‑1 (MQ=255)
                                                 cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGAcc                    <  1:2028352/70‑1 (MQ=255)
                                                 cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCt                   >  1:2973485/1‑71 (MQ=255)
                                                  gAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGaa                                                   <  1:3476178/39‑1 (MQ=255)
                                                  gAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCg                       <  1:333273/66‑1 (MQ=255)
                                                  gAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCt                   <  1:834442/70‑1 (MQ=255)
                                                   aaaaCACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCgg                                  <  1:1960819/54‑1 (MQ=255)
                                                   aaaaCACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGc                 <  1:2907425/71‑1 (MQ=255)
                                                   aaaaCACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGc                 >  1:1628771/1‑71 (MQ=255)
                                                    aaaCACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGAc                     >  1:3534493/1‑66 (MQ=255)
                                                     aaCACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGc                                 >  1:2851102/1‑53 (MQ=255)
                                                      acacCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGccc               >  1:1145245/1‑70 (MQ=255)
                                                        acCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAg            >  1:1576109/1‑71 (MQ=255)
                                                          cTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCg          <  1:3489727/71‑1 (MQ=255)
                                                          cTTCATGATGTACCTGCCTCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCg          <  1:3854788/71‑1 (MQ=255)
                                                           ttCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCg                       <  1:961350/57‑1 (MQ=255)
                                                            tCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGt        >  1:740119/1‑71 (MQ=255)
                                                            tCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGt        >  1:2237282/1‑71 (MQ=255)
                                                             cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCggcg       >  1:1591014/1‑69 (MQ=255)
                                                              atgatgTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCgg                                  >  1:875595/1‑43 (MQ=255)
                                                              atgatgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCggcg       <  1:231505/70‑3 (MQ=255)
                                                               tgatgTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCg                       <  1:2988497/53‑1 (MQ=255)
                                                               tgatgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGGGGCGACCTGCCCAAGCGgcgc      <  1:2500398/70‑1 (MQ=255)
                                                                 atgTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGc                                 <  1:2487162/41‑1 (MQ=255)
                                                                 atgTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGAc                     >  1:3899069/1‑53 (MQ=255)
                                                                 atgTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAt   <  1:2439608/71‑1 (MQ=255)
                                                                 atgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCgg         >  1:2599389/1‑65 (MQ=255)
                                                                 atgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTAt   >  1:2411484/1‑71 (MQ=255)
                                                                  tgTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAtt  <  1:2838530/71‑1 (MQ=255)
                                                                    tATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgtg                            <  1:181752/43‑1 (MQ=255)
                                                                    tATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCg          <  1:3267660/61‑1 (MQ=255)
                                                                      |                                                                   
CGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATT  >  minE/835749‑835885

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑