Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,279 T→A 7.1% V1246E (GTA→GAA)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,2790TA 7.1% 319.3 / 13.8 113V1246E (GTA→GAA) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (40/65);  new base A (6/2);  total (46/67)
Fisher's exact test for biased strand distribution p-value = 6.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 115 total aligned reads displayed.

GTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTT  >  minE/835212‑835349
                                                                    |                                                                      
gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTaa                                                                      <  1:4080190/70‑1 (MQ=255)
gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTaa                                                                      <  1:1597182/70‑1 (MQ=255)
 tAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTaaaa                                                                    <  1:804020/71‑1 (MQ=255)
 tAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTaa                                                                      <  1:3671209/69‑1 (MQ=255)
   gTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATg                                                                  <  1:323748/71‑1 (MQ=255)
        tAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAAt                                                             <  1:2971032/71‑1 (MQ=255)
          aGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTc                                                           >  1:1499077/1‑71 (MQ=255)
            aTGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGa                                                                 <  1:1567267/63‑1 (MQ=255)
             ttaagaaCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGtt                                                        >  1:2154045/3‑71 (MQ=255)
             tGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGtt                                                        <  1:319193/71‑1 (MQ=255)
               aagaaCGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAGGCGAAAAAATGAAAATCCGTTca                                                      >  1:858897/1‑71 (MQ=255)
               aagaaCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTc                                                       <  1:2883549/70‑1 (MQ=255)
               aagaaCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTc                                                       <  1:1757746/70‑1 (MQ=255)
               aagaaCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTc                                                       <  1:4019999/70‑1 (MQ=255)
                agaaCGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAGGCGAAAAAATGAAAATCCGTTcac                                                     <  1:1250350/71‑1 (MQ=255)
                agaaCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTcac                                                     <  1:2396118/71‑1 (MQ=255)
                agaaCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTcac                                                     <  1:3708720/71‑1 (MQ=255)
                agaaCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTcac                                                     <  1:3108867/71‑1 (MQ=255)
                 gaaCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTcac                                                     <  1:2240454/70‑1 (MQ=255)
                  aaCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTc                                                           <  1:3636421/63‑1 (MQ=255)
                   aCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAg                                                  <  1:126830/71‑1 (MQ=255)
                   aCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAg                                                  <  1:1157365/71‑1 (MQ=255)
                    cATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGt                                                 <  1:1872410/71‑1 (MQ=255)
                     aTTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTcac                                                     <  1:426998/66‑1 (MQ=255)
                     aTTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTcac                                                     <  1:1300697/66‑1 (MQ=255)
                     aTTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTc                                                <  1:3359358/71‑1 (MQ=255)
                      ttGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTc                                                <  1:2556883/70‑1 (MQ=255)
                       tGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCgg                                              >  1:1462467/1‑71 (MQ=255)
                        gACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTaaaa                                                                    >  1:1341279/1‑48 (MQ=255)
                        gACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCgg                                              <  1:1879892/70‑1 (MQ=255)
                        gACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCgg                                              <  1:489407/70‑1 (MQ=255)
                          cTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTcaca                                                    >  1:932227/1‑62 (MQ=255)
                              ttAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTg                                       >  1:104270/1‑71 (MQ=255)
                                aGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAg                                                  <  1:326096/58‑1 (MQ=255)
                                aGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGc                                      <  1:1468284/70‑1 (MQ=255)
                                aGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGc                                      <  1:3769537/70‑1 (MQ=255)
                                 gATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTg                                    <  1:1145921/71‑1 (MQ=255)
                                  aTGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTc                                                >  1:3995761/1‑58 (MQ=255)
                                  aTGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTc                                                >  1:864810/1‑58 (MQ=255)
                                  aTGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTc                                                >  1:2464507/1‑58 (MQ=255)
                                  aTGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCgg                                              <  1:2099645/60‑1 (MQ=255)
                                  aTGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCgg                                              <  1:2427586/60‑1 (MQ=255)
                                  aTGGCGAAGGCAATGACAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTc                                                >  1:1343662/1‑57 (MQ=255)
                                   tGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTg                                       <  1:3949191/66‑1 (MQ=255)
                                    ggCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAt                                 >  1:1741163/1‑71 (MQ=255)
                                    ggCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAt                                 >  1:1561820/1‑71 (MQ=255)
                                     gCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGaa                                                                <  1:2820440/39‑1 (MQ=255)
                                     gCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtc                                <  1:2442591/71‑1 (MQ=255)
                                      cGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTca                                                      >  1:2752775/1‑48 (MQ=255)
                                      cGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtct                               <  1:1533646/71‑1 (MQ=255)
                                        aaGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtc                                <  1:2869373/68‑1 (MQ=255)
                                          ggCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGAt                           <  1:538969/71‑1 (MQ=255)
                                           gCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTcac                                                     <  1:2928468/44‑1 (MQ=255)
                                           gCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATg                                          <  1:315717/55‑1 (MQ=255)
                                            cAATGACCAGGTACAGGAGAGCGTAAAATGAAAAt                                                             <  1:972866/35‑1 (MQ=255)
                                             aaTGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAt                                 <  1:3007371/62‑1 (MQ=255)
                                              aTGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGt                       <  1:3008300/71‑1 (MQ=255)
                                               tGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTcac                                                     <  1:3702460/40‑1 (MQ=255)
                                                gACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTca                                                      >  1:2146996/1‑38 (MQ=255)
                                                gACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTg                      >  1:2091815/1‑70 (MQ=255)
                                                 aCCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa                    <  1:271334/71‑1 (MQ=255)
                                                 aCCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa                    >  1:3325566/1‑71 (MQ=255)
                                                 aCCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa                    >  1:3291216/1‑71 (MQ=255)
                                                 aCCAGGTAAAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa                    >  1:3377333/1‑71 (MQ=255)
                                                  tcAGGTACAGGAGGCGAAAAAATGAAAATCCGTTCACAAGTc                                                >  1:146600/2‑42 (MQ=12)
                                                  tcAGGTACAGGAGGCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTaa                                  >  1:1715360/2‑56 (MQ=255)
                                                  tcAGGTACAGGAGGCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTaa                                  >  1:1778164/2‑56 (MQ=255)
                                                  tcAGGTACAGGAGGCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTaa                                  >  1:2259667/2‑56 (MQ=255)
                                                   cAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCCCAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATc                  <  1:412281/71‑1 (MQ=255)
                                                   cAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCg                             <  1:1992091/60‑1 (MQ=255)
                                                   cAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATc                  <  1:3986898/71‑1 (MQ=255)
                                                    aGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATgg                                         <  1:4096275/47‑1 (MQ=255)
                                                    aGGTACAGGAGGCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATc                  <  1:3340659/70‑1 (MQ=255)
                                                     ggTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCg                 >  1:3315964/1‑70 (MQ=255)
                                                      gTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGc                                             <  1:1710886/41‑1 (MQ=255)
                                                      gTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCAt                                           <  1:926956/43‑1 (MQ=255)
                                                      gTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCAt                                           <  1:1538958/43‑1 (MQ=255)
                                                      gTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTg                                    <  1:1424238/50‑1 (MQ=255)
                                                      gTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTg                      >  1:1295622/1‑64 (MQ=255)
                                                       tACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtct                               >  1:2213868/1‑54 (MQ=255)
                                                        aCAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCa                                            >  1:2820936/1‑40 (MQ=255)
                                                        aCAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCAt                                           >  1:1701981/1‑41 (MQ=255)
                                                        aCAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCAt                                           >  1:3472985/1‑41 (MQ=255)
                                                        aCAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTg                                       >  1:37773/1‑45 (MQ=255)
                                                        aCAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCAt                   >  1:3697308/1‑65 (MQ=255)
                                                        aCAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTg             >  1:1720348/1‑71 (MQ=255)
                                                        aCAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTg             >  1:650800/1‑71 (MQ=255)
                                                        aCAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCCGc                                             >  1:358216/1‑39 (MQ=255)
                                                         cAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAt                                 >  1:3197154/1‑50 (MQ=255)
                                                         cAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGc            <  1:1878940/71‑1 (MQ=255)
                                                          aGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGAt                           <  1:2451042/55‑1 (MQ=255)
                                                            gagaGCGTAAAATGAAAATTCGTTCACAAGTCGGCAt                                           <  1:3417621/37‑1 (MQ=255)
                                                            gagaGCGTAAAATGAAAATTCGTTCACAAGTCGGCAt                                           >  1:2335326/1‑37 (MQ=255)
                                                            gagaGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATaa                         >  1:56236/1‑55 (MQ=255)
                                                            gagaGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa                    >  1:1536728/1‑60 (MQ=255)
                                                            gagaGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa                    >  1:2403855/1‑60 (MQ=255)
                                                             agagCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCcaca        <  1:2470985/71‑1 (MQ=255)
                                                              gagCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGa                            >  1:53865/1‑50 (MQ=255)
                                                                 cGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTg    <  1:456784/71‑1 (MQ=255)
                                                                    tAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGaa                                  >  1:1742712/1‑39 (MQ=255)
                                                                    |                                                                      
GTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTT  >  minE/835212‑835349

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: