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→TACnarHnitrate 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.

CGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATT  >  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)
                              aGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAAc                                      >  1:1848602/1‑71 (MQ=255)
                              aGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAAc                                      <  1:135266/71‑1 (MQ=255)
                                gCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAAccc                                    <  1:3876236/71‑1 (MQ=255)
                                 cGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCg                                   >  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)
                                    tGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCa                                >  1:2081803/1‑71 (MQ=255)
                                    tGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCa                                >  1:2271441/1‑71 (MQ=255)
                                      tataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAac                                                  >  1:2661842/1‑52 (MQ=255)
                                      tataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtg                              >  1:2041892/1‑71 (MQ=255)
                                      tataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAt                               <  1:1593906/70‑1 (MQ=255)
                                        taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtg                            >  1:3336679/1‑71 (MQ=255)
                                        taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtg                            <  1:3181309/71‑1 (MQ=255)
                                        taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtg                            <  1:1472040/71‑1 (MQ=255)
                                        taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgt                             <  1:3621570/70‑1 (MQ=255)
                                         aGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtg                            <  1:3689200/70‑1 (MQ=255)
                                          gCCAGTTCGAAAACACCTTCATGATGTATTTGCcg                                                               >  1:1766361/1‑35 (MQ=255)
                                          gCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtgt                           <  1:3867308/70‑1 (MQ=255)
                                          gCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtgt                           <  1:3018493/70‑1 (MQ=255)
                                           ccAGTTCGAAAACACCTTCATGATGTATTTGCCgcgc                                                            <  1:2860980/37‑1 (MQ=255)
                                           ccAGTTCGAAAACACCTTCATGATGTATTTGCCgcgc                                                            <  1:3935445/37‑1 (MQ=255)
                                           ccAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTgg                         <  1:3039361/71‑1 (MQ=255)
                                            cAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGc                        <  1:3913462/71‑1 (MQ=255)
                                             aGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAAtccca                                  >  1:1090451/1‑58 (MQ=255)
                                             aGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTTCGAACACTGCCTGAACCCGGCATGTGTGGc                        >  1:349795/1‑70 (MQ=255)
                                             aGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAAcc                                     >  1:3218393/1‑57 (MQ=255)
                                             aGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGc                        >  1:2971718/1‑70 (MQ=255)
                                             aGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGc                        >  1:289658/1‑70 (MQ=255)
                                              gTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACt                                               <  1:2882346/47‑1 (MQ=38)
                                              gTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACt                                               <  1:1907587/47‑1 (MQ=38)
                                               ttCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGAc                     <  1:3779713/71‑1 (MQ=255)
                                               ttCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGAc                     <  1:1829041/71‑1 (MQ=255)
                                                tCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCg                                                     <  1:3163158/39‑1 (MQ=255)
                                                tCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCgg                                  >  1:2069830/1‑57 (MQ=255)
                                                 tgAAAACACCTTCATTATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCt                   >  1:741533/2‑71 (MQ=255)
                                                 tgAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCt                   >  1:2829823/2‑71 (MQ=255)
                                                 cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGAcc                    <  1:3935453/70‑1 (MQ=255)
                                                 cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGAcc                    <  1:2028352/70‑1 (MQ=255)
                                                 cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCt                   >  1:2973485/1‑71 (MQ=255)
                                                  gAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGaa                                                   <  1:3476178/39‑1 (MQ=255)
                                                  gAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCg                       <  1:333273/66‑1 (MQ=255)
                                                  gAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCt                   <  1:834442/70‑1 (MQ=255)
                                                   aaaaCACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCgg                                  <  1:1960819/54‑1 (MQ=255)
                                                   aaaaCACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGc                 <  1:2907425/71‑1 (MQ=255)
                                                   aaaaCACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGc                 >  1:1628771/1‑71 (MQ=255)
                                                    aaaCACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGAc                     >  1:3534493/1‑66 (MQ=255)
                                                     aaCACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGc                                 >  1:2851102/1‑53 (MQ=255)
                                                      acacCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGccc               >  1:1145245/1‑70 (MQ=255)
                                                        acCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAg            >  1:1576109/1‑71 (MQ=255)
                                                          cTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCg          <  1:3489727/71‑1 (MQ=255)
                                                          cTTCATGATGTACCTGCCTCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCg          <  1:3854788/71‑1 (MQ=255)
                                                           ttCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCg                       <  1:961350/57‑1 (MQ=255)
                                                            tCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGt        >  1:740119/1‑71 (MQ=255)
                                                            tCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGt        >  1:2237282/1‑71 (MQ=255)
                                                             cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCggcg       >  1:1591014/1‑69 (MQ=255)
                                                              atgatgTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCgg                                  >  1:875595/1‑43 (MQ=255)
                                                              atgatgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCggcg       <  1:231505/70‑3 (MQ=255)
                                                               tgatgTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCg                       <  1:2988497/53‑1 (MQ=255)
                                                               tgatgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGGGGCGACCTGCCCAAGCGgcgc      <  1:2500398/70‑1 (MQ=255)
                                                                 atgTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGc                                 <  1:2487162/41‑1 (MQ=255)
                                                                 atgTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGAc                     >  1:3899069/1‑53 (MQ=255)
                                                                 atgTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAt   <  1:2439608/71‑1 (MQ=255)
                                                                 atgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCgg         >  1:2599389/1‑65 (MQ=255)
                                                                 atgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTAt   >  1:2411484/1‑71 (MQ=255)
                                                                  tgTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAtt  <  1:2838530/71‑1 (MQ=255)
                                                                    tATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtg                            <  1:181752/43‑1 (MQ=255)
                                                                    tATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCg          <  1:3267660/61‑1 (MQ=255)
                                                                      |                                                                   
CGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATT  >  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: