Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,275 Δ1 bp 7.6% coding (3733/3744 nt) narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,2750A. 7.6% 301.5 / 21.4 105coding (3733/3744 nt)narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (35/62);  new base . (6/2);  total (41/64)
Fisher's exact test for biased strand distribution p-value = 5.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 114 total aligned reads displayed.

TTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACA  >  minE/835206‑835343
                                                                     |                                                                     
tttGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGagag                                                                      <  1:1979568/71‑1 (MQ=255)
  tgttgtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCg                                                                    >  1:888426/1‑71 (MQ=255)
      gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑Taa                                                                <  1:4080190/70‑1 (MQ=255)
      gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑Taa                                                                <  1:4134512/70‑1 (MQ=255)
      gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑Taa                                                                <  1:1597182/70‑1 (MQ=255)
       tAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑Taaaa                                                              <  1:804020/71‑1 (MQ=255)
       tAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑Taa                                                                <  1:3671209/69‑1 (MQ=255)
         gTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATg                                                            <  1:323748/71‑1 (MQ=255)
              tAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAAt                                                       <  1:1762831/71‑1 (MQ=255)
              tAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAAt                                                       <  1:2971032/71‑1 (MQ=255)
                aGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑Taaaa                                                              >  1:1501658/1‑62 (MQ=255)
                aGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTc                                                     >  1:1499077/1‑71 (MQ=255)
                  aTGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGa                                                           <  1:1567267/63‑1 (MQ=255)
                   ttaagaaCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGtt                                                  >  1:2154045/3‑71 (MQ=255)
                   tGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGtt                                                  <  1:319193/71‑1 (MQ=255)
                     aagaaCGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAG‑GCGAAAAAATGAAAATCCGTTca                                                >  1:858897/1‑71 (MQ=255)
                     aagaaCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTc                                                 <  1:2883549/70‑1 (MQ=255)
                     aagaaCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTc                                                 <  1:1757746/70‑1 (MQ=255)
                     aagaaCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTc                                                 <  1:2921813/70‑1 (MQ=255)
                     aagaaCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTc                                                 <  1:4019999/70‑1 (MQ=255)
                      agaaCGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAG‑GCGAAAAAATGAAAATCCGTTcac                                               <  1:1250350/71‑1 (MQ=255)
                      agaaCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTcac                                               <  1:3708720/71‑1 (MQ=255)
                      agaaCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTcac                                               <  1:2396118/71‑1 (MQ=255)
                       gaaCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTcac                                               <  1:2240454/70‑1 (MQ=255)
                        aaCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTc                                                     <  1:3636421/63‑1 (MQ=255)
                         aCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAg                                            <  1:1157365/71‑1 (MQ=255)
                         aCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAg                                            <  1:126830/71‑1 (MQ=255)
                          cATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAAt                                                       <  1:2554304/59‑1 (MQ=255)
                          cATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGt                                           <  1:1872410/71‑1 (MQ=255)
                           aTTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTcac                                               <  1:426998/66‑1 (MQ=255)
                           aTTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTcac                                               <  1:1300697/66‑1 (MQ=255)
                           aTTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTc                                          <  1:3359358/71‑1 (MQ=255)
                            ttGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTc                                          <  1:2556883/70‑1 (MQ=255)
                             tGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCgg                                        >  1:1462467/1‑71 (MQ=255)
                              gACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑Taaaa                                                              >  1:1341279/1‑48 (MQ=255)
                              gACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCgg                                        <  1:489407/70‑1 (MQ=255)
                                cTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTcaca                                              >  1:932227/1‑62 (MQ=255)
                                    ttAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTg                                 >  1:104270/1‑71 (MQ=255)
                                      aGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAg                                            <  1:326096/58‑1 (MQ=255)
                                      aGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGc                                <  1:3769537/70‑1 (MQ=255)
                                      aGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGc                                <  1:1468284/70‑1 (MQ=255)
                                       gATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTg                              <  1:1145921/71‑1 (MQ=255)
                                        aTGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTc                                          >  1:2464507/1‑58 (MQ=255)
                                        aTGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTc                                          >  1:3995761/1‑58 (MQ=255)
                                        aTGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTc                                          >  1:864810/1‑58 (MQ=255)
                                        aTGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCgg                                        <  1:2099645/60‑1 (MQ=255)
                                        aTGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCgg                                        <  1:2427586/60‑1 (MQ=255)
                                         tGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTg                                 <  1:3949191/66‑1 (MQ=255)
                                          ggCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAt                           >  1:1741163/1‑71 (MQ=255)
                                          ggCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAt                           >  1:1561820/1‑71 (MQ=255)
                                            cGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGaa                                                          <  1:1622925/38‑1 (MQ=255)
                                            cGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTca                                                >  1:2752775/1‑48 (MQ=255)
                                            cGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtct                         <  1:1533646/71‑1 (MQ=255)
                                              aaGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtc                          <  1:2869373/68‑1 (MQ=255)
                                                ggCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGAt                     <  1:538969/71‑1 (MQ=255)
                                                ggCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGAt                     <  1:3187346/71‑1 (MQ=255)
                                                 gCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTcac                                               <  1:2928468/44‑1 (MQ=255)
                                                 gCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATg                                    <  1:315717/55‑1 (MQ=255)
                                                  cAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAAt                                                       <  1:972866/35‑1 (MQ=255)
                                                   aaTGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAt                           <  1:3007371/62‑1 (MQ=255)
                                                    aTGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGt                 <  1:3008300/71‑1 (MQ=255)
                                                     tGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTcac                                               <  1:3702460/40‑1 (MQ=255)
                                                      gACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTca                                                >  1:2146996/1‑38 (MQ=255)
                                                      gACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTg                >  1:2091815/1‑70 (MQ=255)
                                                       aCCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa              <  1:271334/71‑1 (MQ=255)
                                                       aCCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa              >  1:3291216/1‑71 (MQ=255)
                                                       aCCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa              >  1:3325566/1‑71 (MQ=255)
                                                       aCCAGGTAAAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa              >  1:3377333/1‑71 (MQ=255)
                                                        tcAGGTACAGGAG‑GCGAAAAAATGAAAATCCGTTCACAAGTc                                          >  1:146600/2‑42 (MQ=12)
                                                        tcAGGTACAGGAG‑GCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTaa                            >  1:1778164/2‑56 (MQ=255)
                                                        tcAGGTACAGGAG‑GCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTaa                            >  1:2259667/2‑56 (MQ=255)
                                                         cAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCg                       <  1:1992091/60‑1 (MQ=255)
                                                         cAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATc            <  1:3986898/71‑1 (MQ=255)
                                                         cAGGTACAGGAG‑GCGAAAAAATGAAAATCCGTTcaca                                              >  1:3893664/1‑37 (MQ=9)
                                                          aGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATgg                                   <  1:4096275/47‑1 (MQ=255)
                                                          aGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCg           <  1:3659142/71‑1 (MQ=255)
                                                          aGGTACAGGAG‑GCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATc            <  1:3340659/70‑1 (MQ=255)
                                                            gTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGc                                       <  1:1710886/41‑1 (MQ=255)
                                                            gTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCAt                                     <  1:926956/43‑1 (MQ=255)
                                                            gTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCAt                                     <  1:1538958/43‑1 (MQ=255)
                                                            gTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTg                              <  1:1424238/50‑1 (MQ=255)
                                                            gTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTg                >  1:1295622/1‑64 (MQ=255)
                                                              aCAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCa                                      >  1:2820936/1‑40 (MQ=255)
                                                              aCAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCAt                                     >  1:3472985/1‑41 (MQ=255)
                                                              aCAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCAt                                     >  1:1701981/1‑41 (MQ=255)
                                                              aCAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCAt             >  1:3697308/1‑65 (MQ=255)
                                                              aCAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTg       >  1:1720348/1‑71 (MQ=255)
                                                              aCAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTg       >  1:650800/1‑71 (MQ=255)
                                                               cAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGaa                            >  1:3334287/1‑49 (MQ=255)
                                                               cAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAt                           >  1:3197154/1‑50 (MQ=255)
                                                               cAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGc      <  1:1878940/71‑1 (MQ=255)
                                                               cAGGAGAGCG‑TAAAATGAAAATTCGATCACACGTCGTCATGGTGCTGAAtc                          >  1:2101013/1‑51 (MQ=255)
                                                                aGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtctc                        <  1:2356822/52‑1 (MQ=255)
                                                                aGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGAt                     <  1:2451042/55‑1 (MQ=255)
                                                                  gagaGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCAt                                     <  1:3417621/37‑1 (MQ=255)
                                                                  gagaGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCAt                                     >  1:2335326/1‑37 (MQ=255)
                                                                  gagaGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATaa                   >  1:56236/1‑55 (MQ=255)
                                                                  gagaGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa              >  1:2403855/1‑60 (MQ=255)
                                                                   agagCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCcaca  <  1:2470985/71‑1 (MQ=255)
                                                                    gagCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGa                      >  1:53865/1‑50 (MQ=255)
                                                                     |                                                                     
TTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACA  >  minE/835206‑835343

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: ‑