Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 836,481 A→G 16.1% K400R (AAA→AGA) ‡ 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,4810AG16.1% 140.3 / 19.9 62K400R (AAA→AGA) ‡narHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base A (27/25);  new base G (5/5);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

CCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGG  >  minE/836414‑836550
                                                                   |                                                                     
ccAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCCGTGATACCAAAc                                                                     <  1:2405206/70‑1 (MQ=255)
 cAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCg                                                                   <  1:426982/71‑1 (MQ=255)
  aGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCgg                                                                  <  1:2976921/71‑1 (MQ=255)
      gTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACt                                                              <  1:3147642/71‑1 (MQ=255)
      atCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTa                                                                >  1:2848576/2‑69 (MQ=255)
       tCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACt                                                              <  1:1755072/70‑1 (MQ=255)
        cGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTgc                                                            >  1:842499/1‑71 (MQ=255)
         gAAAGCCTGCGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTgc                                                            <  1:2414884/70‑1 (MQ=255)
            aGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTgcgcg                                                         >  1:2001257/1‑70 (MQ=255)
              ccTGCGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTgcgggcgc                                                      <  1:711000/71‑3 (MQ=255)
                 gCGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTg                                                    >  1:68702/1‑70 (MQ=255)
                 gCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCgg                                                                  >  1:2983601/1‑56 (MQ=255)
                   gTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTgcggg                                                         >  1:23711/1‑61 (MQ=255)
                     aTTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACg                                               >  1:721657/1‑71 (MQ=255)
                      ttCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGAAACg                                               <  1:3651336/70‑1 (MQ=255)
                      ttCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGAAACg                                               <  1:2141125/70‑1 (MQ=255)
                      ttCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACg                                               <  1:2835664/70‑1 (MQ=255)
                       tCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTa                                             >  1:518274/1‑71 (MQ=255)
                        ccGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTAt                                            >  1:3076874/1‑71 (MQ=255)
                          ggTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGc                                          >  1:123649/1‑71 (MQ=255)
                           gTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCt                                         <  1:2599427/71‑1 (MQ=255)
                             aCAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTgg                                       >  1:3192528/1‑71 (MQ=255)
                             aCAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTgg                                       >  1:721740/1‑71 (MQ=255)
                             aCAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTgg                                       >  1:3529363/1‑71 (MQ=255)
                              cAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGc                                      >  1:1492931/1‑71 (MQ=255)
                              cAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGc                                      <  1:1249856/71‑1 (MQ=255)
                               aGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCg                                     <  1:1119820/71‑1 (MQ=255)
                                 tATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCa                                                       <  1:1286605/51‑1 (MQ=255)
                                 tATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGAAACGGATGATGGCGAt                                   >  1:140326/1‑71 (MQ=255)
                                 tATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGAAACGGATGATGGCGAt                                   >  1:263630/1‑71 (MQ=255)
                                  aTCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGt                                                                 <  1:1456546/40‑1 (MQ=255)
                                   tCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGc                                          <  1:190613/62‑1 (MQ=255)
                                     tGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCa                                                       >  1:535422/1‑47 (MQ=255)
                                              tgctgACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACaaa                      >  1:3777084/1‑71 (MQ=255)
                                                ctgACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACg                    >  1:1953207/1‑71 (MQ=255)
                                                 tgACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAAc                     >  1:1255969/1‑69 (MQ=255)
                                                  gACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACg                                               <  1:3094015/42‑1 (MQ=255)
                                                   aCCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTg                  <  1:862594/70‑1 (MQ=255)
                                                   aCCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGc                 <  1:351318/71‑1 (MQ=255)
                                                    ccgccgGTGATACCAAACCGGTACTGCGCGCACTGAAACGTa                                             >  1:2808009/1‑42 (MQ=255)
                                                     cgccgGTGATACCAAACCGGTACTGCGCGCACTGAAAc                                                >  1:1782418/1‑38 (MQ=255)
                                                      gccgGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGa              <  1:3328494/71‑1 (MQ=255)
                                                       ccgGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGa              >  1:1654898/1‑70 (MQ=255)
                                                         ggTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAAc           >  1:3136597/1‑71 (MQ=255)
                                                         ggTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAAc           >  1:2637478/1‑71 (MQ=255)
                                                             aTACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGtt       <  1:3222740/71‑1 (MQ=255)
                                                             aTACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGtt       <  1:2858762/71‑1 (MQ=255)
                                                             aTACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGtt       <  1:663707/71‑1 (MQ=255)
                                                             aTACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGtt       <  1:2773572/71‑1 (MQ=255)
                                                             aTACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGtt       <  1:3892482/71‑1 (MQ=255)
                                                             aTACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGtt       <  1:2069211/71‑1 (MQ=255)
                                                              tACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGt                   <  1:2990046/58‑1 (MQ=255)
                                                               aCCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGAt                                   >  1:675893/1‑41 (MQ=255)
                                                               aCCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGAt                                   >  1:1998264/1‑41 (MQ=255)
                                                               aCCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCg         <  1:2118596/67‑1 (MQ=255)
                                                                ccAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACg                    >  1:2734011/1‑55 (MQ=255)
                                                                ccAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACg                    >  1:2390871/1‑55 (MQ=255)
                                                                ccAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGa     <  1:3570913/70‑1 (MQ=255)
                                                                ccAAACCGGTACTGCGCGCACTGAAACCTATGCTGGc                                      <  1:1992818/37‑1 (MQ=255)
                                                                 cAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACg   >  1:1338379/1‑71 (MQ=255)
                                                                 cAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACg   >  1:1592283/1‑71 (MQ=255)
                                                                 cAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACg   >  1:2032387/1‑71 (MQ=255)
                                                                  aaaCCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACa                        <  1:3699867/49‑1 (MQ=255)
                                                                  aaaCCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACgg  >  1:1705896/1‑71 (MQ=255)
                                                                   aaCCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGt                               >  1:617399/1‑41 (MQ=255)
                                                                   aaCCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACgg  >  1:3178845/1‑70 (MQ=255)
                                                                   aaCCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACgg  >  1:837261/1‑70 (MQ=255)
                                                                   |                                                                     
CCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGG  >  minE/836414‑836550

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: