Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,820 T→C 19.6% L180L (TTG→CTG)  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,8200TC19.6% 211.6 / 46.7 102L180L (TTG→CTG) 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.87e-01

Only 100 of 102 total aligned reads displayed.

AAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATT  >  minE/835752‑835885
                                                                    |                                                                  
 aGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGc                                                                 >  1:1473338/1‑71 (MQ=255)
    cAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCcgc                                                              <  1:3433131/71‑1 (MQ=255)
      aGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGcc                                                                >  1:2299302/1‑67 (MQ=255)
      aGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGcc                                                                >  1:1396351/1‑67 (MQ=255)
      aGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCgcgc                                                            >  1:3482164/1‑71 (MQ=255)
      aGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCgcgc                                                            >  1:1582637/1‑71 (MQ=255)
       gAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTg                                                                  <  1:3467917/64‑1 (MQ=255)
       caaCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGcc                                                           >  1:2392205/2‑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:614351/1‑71 (MQ=255)
         aCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCtg                                                         <  1:2015105/71‑1 (MQ=255)
         aCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCtg                                                         >  1:1546133/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:1714013/1‑71 (MQ=255)
                     tCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGc                                             >  1:3526955/1‑71 (MQ=255)
                     tCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGc                                             >  1:2269706/1‑71 (MQ=255)
                      ccAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGcc                                            >  1:304097/1‑71 (MQ=255)
                      ccAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGcc                                            >  1:903867/1‑71 (MQ=255)
                       cAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCt                                           >  1:608740/1‑71 (MQ=255)
                        aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg                                          <  1:945051/71‑1 (MQ=255)
                        aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg                                          <  1:940361/71‑1 (MQ=255)
                        aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg                                          <  1:489827/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:2271441/1‑71 (MQ=255)
                                 tGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCa                                >  1:2081803/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:3181309/71‑1 (MQ=255)
                                     taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtg                            <  1:1472040/71‑1 (MQ=255)
                                     taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtg                            >  1:3336679/1‑71 (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:3018493/70‑1 (MQ=255)
                                       gCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtgt                           <  1:3867308/70‑1 (MQ=255)
                                        ccAGTTCGAAAACACCTTCATGATGTATTTGCCgcgc                                                            <  1:3935445/37‑1 (MQ=255)
                                        ccAGTTCGAAAACACCTTCATGATGTATTTGCCgcgc                                                            <  1:2860980/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:2028352/70‑1 (MQ=255)
                                              cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGAcc                    <  1:3935453/70‑1 (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)
                                                                    |                                                                  
AAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATT  >  minE/835752‑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: