Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,315 G→T 10.9% L11L (CTG→CTT 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,3150GT10.9% 295.9 / 12.2 91L11L (CTG→CTTnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base G (35/46);  new base T (5/5);  total (40/51)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAG  >  minE/835245‑835385
                                                                       |                                                                      
gATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTg                                                                        <  1:1145921/71‑1 (MQ=255)
   ggCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAt                                                                     >  1:1561820/1‑71 (MQ=255)
   ggCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAt                                                                     >  1:1741163/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:3187346/71‑1 (MQ=255)
         ggCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGAt                                                               <  1:538969/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:3325566/1‑71 (MQ=255)
                aCCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa                                                        >  1:3291216/1‑71 (MQ=255)
                aCCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa                                                        <  1:271334/71‑1 (MQ=255)
                aCCAGGTAAAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa                                                        >  1:3377333/1‑71 (MQ=255)
                 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)
                   aGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCg                                                     <  1:3659142/71‑1 (MQ=255)
                   aGGTACAGGAGGCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATc                                                      <  1:3340659/70‑1 (MQ=255)
                    ggTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCg                                                     >  1:3315964/1‑70 (MQ=255)
                     gTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTg                                                                        <  1:1424238/50‑1 (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)
                        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)
                           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)
                                       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:1209399/71‑1 (MQ=255)
                                                 gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa                        <  1:584646/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:1258608/55‑1 (MQ=255)
                                                    cacaAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                     <  1:2091827/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:3414327/71‑1 (MQ=255)
                                                         gTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTg                <  1:142995/71‑1 (MQ=255)
                                                         gTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTg                >  1:2017075/1‑71 (MQ=255)
                                                          tCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                     <  1:3542954/49‑1 (MQ=255)
                                                          tCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                     <  1:555095/49‑1 (MQ=255)
                                                          tCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                     <  1:2673847/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:2638249/37‑1 (MQ=255)
                                                                      tGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                     <  1:883458/37‑1 (MQ=255)
                                                                      tGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCa                                    <  1:2111351/38‑1 (MQ=255)
                                                                      tGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCa                                    <  1:3554844/38‑1 (MQ=255)
                                                                       gAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAg  >  1:2275248/1‑71 (MQ=255)
                                                                       |                                                                      
GATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAG  >  minE/835245‑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: