Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,291 T→C 9.2% I3I (ATT→ATC 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,2910TC 9.2% 304.3 / 18.9 109I3I (ATT→ATCnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base T (42/57);  new base C (8/2);  total (50/59)
Fisher's exact test for biased strand distribution p-value = 4.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 110 total aligned reads displayed.

AGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTAC  >  minE/835222‑835356
                                                                      |                                                                 
aGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTc                                                                  >  1:1499077/1‑71 (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:1757746/70‑1 (MQ=255)
     aagaaCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTc                                                              <  1:4019999/70‑1 (MQ=255)
     aagaaCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTc                                                              <  1:2883549/70‑1 (MQ=255)
     aagaaCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTc                                                              <  1:2921813/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)
       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:1300697/66‑1 (MQ=255)
           aTTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTcac                                                            <  1:426998/66‑1 (MQ=255)
            ttGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTc                                                       <  1:2556883/70‑1 (MQ=255)
             tGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCgg                                                     >  1:1462467/1‑71 (MQ=255)
              gACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCgg                                                     <  1:1879892/70‑1 (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:1561820/1‑71 (MQ=255)
                          ggCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAt                                        >  1:1741163/1‑71 (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:3187346/71‑1 (MQ=255)
                                ggCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGAt                                  <  1:538969/71‑1 (MQ=255)
                                 gCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTcac                                                            <  1:2928468/44‑1 (MQ=255)
                                 gCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATg                                                 <  1:315717/55‑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:3291216/1‑71 (MQ=255)
                                       aCCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa                           >  1:3325566/1‑71 (MQ=255)
                                       aCCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa                           <  1:271334/71‑1 (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)
                                         cAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCCCAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATc                         <  1:412281/71‑1 (MQ=255)
                                         cAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCg                                    <  1:1992091/60‑1 (MQ=255)
                                         cAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATc                         <  1:3986898/71‑1 (MQ=255)
                                         cAGGTACAGGAGGCGAAAAAATGAAAATCCGTTcaca                                                           >  1:3893664/1‑37 (MQ=9)
                                          aGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATgg                                                <  1:4096275/47‑1 (MQ=255)
                                          aGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCg                        <  1:3659142/71‑1 (MQ=255)
                                          aGGTACAGGAGGCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATc                         <  1:3340659/70‑1 (MQ=255)
                                           ggTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCg                        >  1:3315964/1‑70 (MQ=255)
                                            gTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCAt                                                  <  1:926956/43‑1 (MQ=255)
                                            gTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCAt                                                  <  1:1538958/43‑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)
                                              aCAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCAt                          >  1:3697308/1‑65 (MQ=255)
                                              aCAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTg                    >  1:650800/1‑71 (MQ=255)
                                              aCAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTg                    >  1:1720348/1‑71 (MQ=255)
                                              aCAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCCGc                                                    >  1:358216/1‑39 (MQ=255)
                                               cAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGaa                                         >  1:3334287/1‑49 (MQ=255)
                                               cAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAt                                        >  1:3197154/1‑50 (MQ=255)
                                               cAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGc                   <  1:1878940/71‑1 (MQ=255)
                                               cAGGAGAGCGTAAAATGAAAATTCGATCACACGTCGTCATGGTGCTGAAtc                                       >  1:2101013/1‑51 (MQ=255)
                                                aGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtctc                                     <  1:2356822/52‑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:742718/1‑50 (MQ=255)
                                                    gagCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGa                                   >  1:53865/1‑50 (MQ=255)
                                                       cGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTg           <  1:456784/71‑1 (MQ=255)
                                                          tAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGaa                                         >  1:1742712/1‑39 (MQ=255)
                                                          tAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGtt         >  1:1315777/1‑71 (MQ=255)
                                                          aaaaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGtt         >  1:1565500/2‑71 (MQ=255)
                                                            aaaTGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCa       <  1:3013980/71‑1 (MQ=255)
                                                            aaaTGAAAATTCGTTCACAAATCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc        <  1:2036912/70‑1 (MQ=255)
                                                             aaTGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCTGCTGCCACACCTGTTCa       <  1:1130553/70‑1 (MQ=255)
                                                               tGAAAATCCGTACACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGt     >  1:760565/1‑70 (MQ=255)
                                                                 aaaaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGTCTGCCACACCTGTTCAGTTAc  >  1:1714717/1‑71 (MQ=255)
                                                                 aaaaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTAc  >  1:525186/1‑71 (MQ=255)
                                                                    aTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATa                                 <  1:2149539/37‑1 (MQ=255)
                                                                     ttCGTTCACAAGTCGGCATGGTGCTGAATCTCGAt                                  <  1:4516/35‑1 (MQ=255)
                                                                      |                                                                 
AGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTAC  >  minE/835222‑835356

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: