Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,278:1 +A 17.1% coding (3736/3744 nt) narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,2781.A17.1% 244.9 / 64.7 82coding (3736/3744 nt)narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base . (31/37);  new base A (5/9);  total (36/46)
Fisher's exact test for biased strand distribution p-value = 5.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTG  >  minE/835208‑835347
                                                                       |                                                                      
tgttgtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCg                                                                         >  1:541158/1‑71 (MQ=255)
 gttgtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCg                                                                         >  1:3855325/1‑70 (MQ=255)
   tgtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTa                                                                     >  1:3829155/1‑70 (MQ=255)
    gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTaa                                                                    <  1:2938315/70‑1 (MQ=255)
    gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTaa                                                                    >  1:328648/1‑70 (MQ=255)
    gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTaa                                                                    >  1:3759474/1‑70 (MQ=255)
     tAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTaaa                                                                   <  1:3301694/70‑1 (MQ=255)
      aGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTaaaa                                                                  <  1:2017284/70‑1 (MQ=255)
      aGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAAt                                                                 >  1:3537989/1‑71 (MQ=255)
         gCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGa                                                               >  1:2445047/1‑70 (MQ=255)
          cgcaaGATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAGGCGAAAAAATGaaa                                                             <  1:2209642/68‑1 (MQ=255)
          cgcaaGATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAGGCGAAAAAATGaaa                                                             <  1:1657235/68‑1 (MQ=255)
            tAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGa                                                               <  1:3772670/67‑1 (MQ=255)
            tAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAAt                                                           <  1:751456/71‑1 (MQ=255)
            tAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGAGTAAAATGAAAAt                                                           <  1:3862461/71‑1 (MQ=255)
            caaGATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAGGCGAAAAAATGAAAAt                                                           >  1:3320715/2‑71 (MQ=255)
            caaGATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAGGCGAAAAAATGAAAAt                                                           >  1:2576504/2‑71 (MQ=255)
              aGATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAGGCGAAAAAATGAAAAtcc                                                         >  1:2944020/1‑69 (MQ=255)
              aGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTc                                                         >  1:3828491/1‑71 (MQ=255)
              aGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTc                                                         >  1:614527/1‑71 (MQ=255)
              aGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTc                                                         <  1:1001081/71‑1 (MQ=255)
               gATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAGGCGAAAAAATGAAAATCCg                                                        <  1:1923858/71‑1 (MQ=255)
               gATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAGGCGAAAAAATGAAAATCCg                                                        <  1:148799/71‑1 (MQ=255)
                   aagaaCGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAGGCGAAAAAATGAAAATCCGTTca                                                    <  1:1018540/71‑1 (MQ=255)
                    agaaCGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAGGCGAAAAAATGAAAATCCGTTcac                                                   <  1:1445858/71‑1 (MQ=255)
                     gaaCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTcac                                                   <  1:1311226/70‑1 (MQ=255)
                       aCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTc                                                         >  1:2275607/1‑62 (MQ=255)
                         aTTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGaaaa                                                            <  1:3940656/57‑1 (MQ=255)
                            gACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCgg                                            <  1:3132414/70‑1 (MQ=255)
                            gACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCgg                                            <  1:1714093/70‑1 (MQ=255)
                              cTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCAt                                         >  1:1855724/1‑71 (MQ=255)
                              cTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCAt                                         >  1:1350295/1‑71 (MQ=255)
                               tGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTc                                                         >  1:1875490/1‑54 (MQ=255)
                               tGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGtt                                                      >  1:2995808/1‑57 (MQ=255)
                               tGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATg                                        >  1:595237/1‑71 (MQ=255)
                                 gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATgg                                       <  1:2762370/70‑1 (MQ=255)
                                 gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATgg                                       <  1:1763620/70‑1 (MQ=255)
                                 gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATgg                                       <  1:3760649/70‑1 (MQ=255)
                                 gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATg                                        <  1:1628791/69‑1 (MQ=255)
                                 gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGt                                      >  1:2219106/1‑71 (MQ=255)
                                 gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGt                                      <  1:3780801/71‑1 (MQ=255)
                                 gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGt                                      >  1:2223794/1‑71 (MQ=255)
                                  ttAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTg                                     >  1:3238161/1‑71 (MQ=255)
                                   tAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTcaca                                                  >  1:2932218/1‑57 (MQ=255)
                                   tAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTcaca                                                  >  1:1305391/1‑57 (MQ=255)
                                    aGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGaaaa                                                            <  1:484874/46‑1 (MQ=255)
                                    aGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCt                                   >  1:265481/1‑71 (MQ=255)
                                     gATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCa                                          <  1:1642333/63‑1 (MQ=255)
                                     gATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCt                                   <  1:3464655/70‑1 (MQ=255)
                                         gCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCgg                                            <  1:1154401/57‑1 (MQ=255)
                                         gCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCa                                          <  1:1689228/59‑1 (MQ=255)
                                         gCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAt                               <  1:3189656/70‑1 (MQ=255)
                                          cGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtct                             >  1:3818181/1‑71 (MQ=255)
                                          cGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtct                             <  1:1833539/71‑1 (MQ=255)
                                           gAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGa                                                               <  1:123019/36‑1 (MQ=255)
                                           gAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCa                                          >  1:713048/1‑57 (MQ=255)
                                            aaGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtctc                            >  1:1295489/1‑70 (MQ=255)
                                             aGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtct                             >  1:1676601/1‑68 (MQ=255)
                                             aGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGa                          <  1:665169/71‑1 (MQ=255)
                                             aGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGa                          >  1:632673/1‑71 (MQ=255)
                                              ggCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGTCATGGTGCTGAAtc                              <  1:1098645/66‑1 (MQ=255)
                                               gCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGt                                                       <  1:2951322/40‑1 (MQ=255)
                                                cAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGAt                         >  1:2548940/1‑69 (MQ=255)
                                                cAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATaa                       <  1:2397116/71‑1 (MQ=255)
                                                 aaTGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGtt                                                      <  1:1026738/39‑1 (MQ=255)
                                                 aaTGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTcaca                                                  >  1:3176505/1‑43 (MQ=255)
                                                  aTGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATgg                                       <  1:260240/53‑1 (MQ=255)
                                                   tGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTcac                                                   >  1:1793006/1‑40 (MQ=255)
                                                   tGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTg                    >  1:2816562/1‑71 (MQ=255)
                                                      tcAGGTACAGGAGGCGAAAAAATGAAAATCCGTTCACaa                                                 <  1:1640459/38‑1 (MQ=9)
                                                      tcAGGTACAGGAGGCGAAAAAATGAAAATCCGTTCACaa                                                 <  1:2193019/38‑1 (MQ=9)
                                                      tcAGGTACAGGAGGCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTAt                 <  1:1038941/70‑1 (MQ=255)
                                                       cAGGTACAGGAGGCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTAt                 >  1:322274/1‑70 (MQ=255)
                                                        aGGTACAGGAGAGTGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCAt                 <  1:3222746/69‑1 (MQ=255)
                                                        aGGTACAGGAGGCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCg               >  1:580056/1‑71 (MQ=255)
                                                         ggTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCgg              <  1:2976711/71‑1 (MQ=255)
                                                          gTACAGGAGGCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGc             >  1:3140096/1‑71 (MQ=255)
                                                            aCAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa                  <  1:3797985/64‑1 (MQ=255)
                                                            aCAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa                  <  1:1988319/64‑1 (MQ=255)
                                                             cAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTg           <  1:2270425/70‑1 (MQ=255)
                                                                gagaGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGcc         >  1:3111585/1‑69 (MQ=255)
                                                                 agagCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGc                   <  1:549269/58‑1 (MQ=255)
                                                                  gagCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATaa                       >  1:681826/1‑53 (MQ=255)
                                                                     cGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTg  <  1:1840511/71‑1 (MQ=255)
                                                                       |                                                                      
TGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTG  >  minE/835208‑835347

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