Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,318 T→C 9.1% N12N (AAT→AAC 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,3180TC 9.1% 211.8 / 11.5 77N12N (AAT→AACnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base T (29/41);  new base C (2/5);  total (31/46)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

GGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAG  >  minE/835248‑835385
                                                                       |                                                                   
ggCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAt                                                                     >  1:1741163/1‑71 (MQ=255)
ggCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAt                                                                     >  1:1561820/1‑71 (MQ=255)
 gCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtc                                                                    <  1:2442591/71‑1 (MQ=255)
  cGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtct                                                                   <  1:1533646/71‑1 (MQ=255)
    aaGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtc                                                                    <  1:2869373/68‑1 (MQ=255)
      ggCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGAt                                                               <  1:538969/71‑1 (MQ=255)
      ggCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGAt                                                               <  1:3187346/71‑1 (MQ=255)
         aaTGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAt                                                                     <  1:3007371/62‑1 (MQ=255)
          aTGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGt                                                           <  1:3008300/71‑1 (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)
               cAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCCCAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATc                                                      <  1:412281/71‑1 (MQ=255)
               cAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCg                                                                 <  1:1992091/60‑1 (MQ=255)
               cAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATc                                                      <  1:3986898/71‑1 (MQ=255)
                aGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCg                                                     <  1:3659142/71‑1 (MQ=255)
                aGGTACAGGAGGCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATc                                                      <  1:3340659/70‑1 (MQ=255)
                 ggTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCg                                                     >  1:3315964/1‑70 (MQ=255)
                  gTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTg                                                          >  1:1295622/1‑64 (MQ=255)
                   tACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtct                                                                   >  1:2213868/1‑54 (MQ=255)
                    aCAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCAt                                                       >  1:3697308/1‑65 (MQ=255)
                    aCAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTg                                                 >  1:1720348/1‑71 (MQ=255)
                    aCAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTg                                                 >  1:650800/1‑71 (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)
                        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)
                          gagCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGa                                                                >  1:742718/1‑50 (MQ=255)
                             cGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTg                                        <  1:456784/71‑1 (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)
                                    aTGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCAt                                                       >  1:1942934/1‑50 (MQ=255)
                                     tGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGtt                                 <  1:2563524/71‑1 (MQ=255)
                                     tGAAAATCCGTACACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGt                                  >  1:760565/1‑70 (MQ=255)
                                       aaaaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGTCTGCCACACCTGTTCAGTTAc                               >  1:1714717/1‑71 (MQ=255)
                                       aaaaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTAc                               >  1:525186/1‑71 (MQ=255)
                                          atccGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTc                                     <  1:1220197/59‑1 (MQ=255)
                                          aTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATa                                                              <  1:2149539/37‑1 (MQ=255)
                                           ttCGTTCACAAGTCGGCATGGTGCTGAATCTCGAt                                                               <  1:4516/35‑1 (MQ=255)
                                              gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa                        <  1:3949171/71‑1 (MQ=255)
                                              gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa                        <  1:2609513/71‑1 (MQ=255)
                                              gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa                        <  1:584646/71‑1 (MQ=255)
                                              gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa                        <  1:1209399/71‑1 (MQ=255)
                                               ttCACAAGTCGGCATGGTGCTGAATCTCGATAAGt                                                           <  1:3827467/35‑1 (MQ=255)
                                               ttCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa                        <  1:2759346/70‑1 (MQ=255)
                                                tCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCAt                                                       >  1:4119775/1‑38 (MQ=255)
                                                tCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGt                                  >  1:959615/1‑59 (MQ=255)
                                                 cacaAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                     <  1:2091827/55‑1 (MQ=255)
                                                 cacaAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                     <  1:1258608/55‑1 (MQ=255)
                                                 cacaAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAAc                     >  1:218643/1‑71 (MQ=255)
                                                  acaAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAAc                     <  1:1398341/70‑1 (MQ=255)
                                                  acaAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAAc                     <  1:1864425/70‑1 (MQ=255)
                                                   caAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGt                   <  1:2737055/71‑1 (MQ=255)
                                                    aaGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCt                             >  1:3289879/1‑60 (MQ=255)
                                                      gTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTg                >  1:2017075/1‑71 (MQ=255)
                                                      gTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTg                <  1:3414327/71‑1 (MQ=255)
                                                      gTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTg                <  1:142995/71‑1 (MQ=255)
                                                       tCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                     <  1:2673847/49‑1 (MQ=255)
                                                       tCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                     <  1:3542954/49‑1 (MQ=255)
                                                       tCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                     <  1:555095/49‑1 (MQ=255)
                                                       tCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTgg               <  1:1611548/71‑1 (MQ=255)
                                                        cGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTgg               >  1:435241/1‑70 (MQ=255)
                                                         ggCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGAc             <  1:3251052/71‑1 (MQ=255)
                                                             tGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                     <  1:1698773/43‑1 (MQ=255)
                                                               gTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGcc        <  1:377194/70‑1 (MQ=255)
                                                               gTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCg       >  1:829858/1‑71 (MQ=255)
                                                                tGCTGAATCTCGATAAGTGCATCGGCTGCCAAACCTg                                        >  1:620591/1‑37 (MQ=255)
                                                                 gCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTaaa                        <  1:2487693/52‑1 (MQ=255)
                                                                  cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg     <  1:2688763/70‑1 (MQ=255)
                                                                  cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg     >  1:696405/1‑70 (MQ=255)
                                                                  cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg     >  1:701179/1‑70 (MQ=255)
                                                                   tGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                     <  1:883458/37‑1 (MQ=255)
                                                                   tGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                     <  1:2638249/37‑1 (MQ=255)
                                                                   tGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCa                                    <  1:3554844/38‑1 (MQ=255)
                                                                   tGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCa                                    <  1:2111351/38‑1 (MQ=255)
                                                                    gAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAg  >  1:2275248/1‑71 (MQ=255)
                                                                       tctcGATAAGTGCATCGGCTGCCACACCTGTTCAGTTAc                               <  1:1352567/39‑1 (MQ=255)
                                                                       |                                                                   
GGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAG  >  minE/835248‑835385

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: