Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 836,011 G→A 9.4% K243K (AAG→AAA 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,0110GA 9.4% 160.1 / 6.1 53K243K (AAG→AAAnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base G (22/26);  new base A (2/3);  total (24/29)
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.84e-01

TATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCG  >  minE/835945‑836079
                                                                  |                                                                    
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)
                tGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCcgcg                                                  <  1:2440273/71‑1 (MQ=255)
                      tACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGa                                            <  1:975383/71‑1 (MQ=255)
                        cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAg                                          >  1:3513274/1‑71 (MQ=255)
                        cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAg                                          >  1:206487/1‑71 (MQ=255)
                          aaaaaaTCTACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGagtc                                         <  1:1046301/70‑4 (MQ=255)
                          aaaaaaTCTACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGagtc                                         <  1:3543412/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)
                               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:1827491/67‑1 (MQ=255)
                                    cTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAg                                  <  1:3433341/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:3961544/1‑69 (MQ=255)
                                          cTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCgt                          >  1:1707050/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)
                                                              agaAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTgtcggt    <  1:2091282/71‑1 (MQ=255)
                                                              agaAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTgtcgg     <  1:588547/70‑1 (MQ=255)
                                                                aaGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTgtcggtcg  <  1:2220681/71‑1 (MQ=255)
                                                                aaGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTgtcggtcg  <  1:3940566/71‑1 (MQ=255)
                                                                  gTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGaa                <  1:609223/55‑1 (MQ=255)
                                                                  |                                                                    
TATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCG  >  minE/835945‑836079

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: