Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 836,005 T→A 7.3% S241S (TCT→TCA narH → nitrate reductase 1, beta (Fe‑S) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE836,0050TA 7.3% 193.1 / 6.8 68S241S (TCT→TCAnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base T (31/32);  new base A (2/3);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

TGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGT  >  minE/835936‑836073
                                                                     |                                                                    
tggctggcGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTg                                                                     <  1:2854133/71‑1 (MQ=255)
  gctggcGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTgag                                                                   <  1:1387852/71‑1 (MQ=255)
         tATGTGCATCACTGTATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCa                                                            <  1:1657722/71‑1 (MQ=255)
         tATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCa                                                            <  1:1182694/71‑1 (MQ=255)
          aTGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCa                                                            >  1:1844103/1‑70 (MQ=255)
           tgtgCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCAtc                                                          >  1:1746518/1‑71 (MQ=255)
             tgCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAtctt                                                        <  1:1243368/71‑1 (MQ=255)
               cATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAtcttct                                                      <  1:1351853/71‑1 (MQ=255)
                 tCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGc                                                    >  1:1081727/1‑71 (MQ=255)
                  cACTGGATGCCCGAACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAt                                                           >  1:3821169/1‑63 (MQ=255)
                   aCTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCt                                                   <  1:760497/70‑1 (MQ=255)
                   aCTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTa                                                  >  1:1308476/1‑71 (MQ=255)
                   aCTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTa                                                  >  1:1611484/1‑71 (MQ=255)
                       gATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATcc                                               <  1:3891101/70‑1 (MQ=255)
                       gATGCCAGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa                                                                 >  1:540667/1‑52 (MQ=255)
                         tGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCcgcg                                            <  1:2440273/71‑1 (MQ=255)
                             cGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTgaga                                                                  <  1:1620997/45‑1 (MQ=255)
                              gTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTgag                                                                   <  1:1735364/43‑1 (MQ=255)
                               tACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTg                                                                     <  1:3202829/40‑1 (MQ=255)
                               tACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGa                                      <  1:975383/71‑1 (MQ=255)
                                aCAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTg                                                                     <  1:2922962/39‑1 (MQ=255)
                                aCAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa                                                                 <  1:1913945/43‑1 (MQ=255)
                                aCAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa                                                                 <  1:3135431/43‑1 (MQ=255)
                                 cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTgaga                                                                  >  1:992809/1‑41 (MQ=255)
                                 cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTga                                                                    <  1:3459567/39‑1 (MQ=255)
                                 cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTga                                                                    >  1:3792626/1‑39 (MQ=255)
                                 cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTga                                                                    >  1:1926320/1‑39 (MQ=255)
                                 cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa                                                                 >  1:913501/1‑42 (MQ=255)
                                 cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa                                                                 >  1:4024130/1‑42 (MQ=255)
                                 cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa                                                                 <  1:2943451/42‑1 (MQ=255)
                                 cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa                                                                 >  1:2902138/1‑42 (MQ=255)
                                 cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAg                                    >  1:206487/1‑71 (MQ=255)
                                 cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAg                                    >  1:3513274/1‑71 (MQ=255)
                                   aaaaaaTCTACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGagtc                                   <  1:3543412/70‑4 (MQ=255)
                                   aaaaaaTCTACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGagtc                                   <  1:1046301/70‑4 (MQ=255)
                                    aaaaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAtct                                                         <  1:1759857/47‑1 (MQ=255)
                                    aaaaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAt                                                           >  1:1905121/1‑45 (MQ=255)
                                    aaaaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAt                                                           >  1:591480/1‑45 (MQ=255)
                                     aaaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCggg                                <  1:2924591/71‑1 (MQ=255)
                                       aaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa                                                                 >  1:3983243/1‑36 (MQ=255)
                                        aTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa                                                                 >  1:2057001/1‑35 (MQ=255)
                                        aTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGaa                                     <  1:1196987/63‑1 (MQ=255)
                                        aTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCa                             <  1:1158317/71‑1 (MQ=255)
                                        aTCTACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGa                                      >  1:3249947/1‑62 (MQ=255)
                                        aTCTACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCa                             <  1:3802249/71‑1 (MQ=255)
                                           tACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTg                                       >  1:3051319/1‑58 (MQ=255)
                                             cTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAg                            <  1:3433341/67‑1 (MQ=255)
                                             cTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAg                            <  1:1827491/67‑1 (MQ=255)
                                             cTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGa                        >  1:3220856/1‑71 (MQ=255)
                                                cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAtctt                                                        >  1:3751946/1‑36 (MQ=255)
                                                cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCg                                  >  1:2901892/1‑58 (MQ=255)
                                                cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg                     >  1:1894805/1‑71 (MQ=255)
                                                 aaCTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAg                                    >  1:922271/1‑55 (MQ=255)
                                                 aaCTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGAcc                      <  1:3493792/69‑1 (MQ=255)
                                                 aaCTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg                     <  1:1790614/70‑1 (MQ=255)
                                                  aCTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCgtg                   <  1:3979604/71‑1 (MQ=255)
                                                   cTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTa                                                  <  1:2609985/39‑1 (MQ=255)
                                                   cTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCgt                    >  1:1707050/1‑69 (MQ=255)
                                                   cTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCgt                    >  1:3961544/1‑69 (MQ=255)
                                                     ggAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCcg                                              >  1:1260889/1‑41 (MQ=255)
                                                      gAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGc                <  1:2191634/70‑1 (MQ=255)
                                                       aaGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCt               >  1:3657457/1‑70 (MQ=255)
                                                        agagCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg                     <  1:1445157/63‑1 (MQ=255)
                                                         gagCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCt               >  1:754678/1‑68 (MQ=255)
                                                             ggTAAGGCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAAc        <  1:3216790/71‑1 (MQ=255)
                                                                aaGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCg                                  >  1:3771073/1‑42 (MQ=255)
                                                                aaGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCtg     >  1:1992359/1‑71 (MQ=255)
                                                                   tCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTgtgt  <  1:4030856/71‑1 (MQ=255)
                                                                     |                                                                    
TGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGT  >  minE/835936‑836073

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: