Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,262 C→T 12.6% D1240D (GAC→GAT)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,2620CT12.6% 229.6 / 19.5 87D1240D (GAC→GAT) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (36/40);  new base T (5/6);  total (41/46)
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.04e-01

CCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTG  >  minE/835196‑835329
                                                                  |                                                                    
ccGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAg                                                                   <  1:3194256/70‑1 (MQ=255)
 cgcgATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGt                                                                 <  1:3802636/71‑1 (MQ=255)
    gatgatTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACa                                                              <  1:652433/71‑1 (MQ=255)
     aTGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAg                                                             <  1:853548/71‑1 (MQ=255)
       gAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAg                                                             >  1:959944/1‑69 (MQ=255)
        aGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGga                                                           <  1:1839570/70‑1 (MQ=255)
        aGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGga                                                           <  1:3819905/70‑1 (MQ=255)
           ttgttgTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGgaga                                                         >  1:780615/1‑69 (MQ=255)
           ttgttgTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGgaga                                                         >  1:2702740/1‑69 (MQ=255)
            tgttgtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCg                                                      >  1:541158/1‑71 (MQ=255)
             gttgtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCg                                                      >  1:3855325/1‑70 (MQ=255)
               tgtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAg                                                             >  1:3705822/1‑61 (MQ=255)
               tgtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑Ta                                                   >  1:3829155/1‑70 (MQ=255)
                gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAg                                                             >  1:1157715/1‑60 (MQ=255)
                gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑Taa                                                  >  1:3759474/1‑70 (MQ=255)
                gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑Taa                                                  >  1:328648/1‑70 (MQ=255)
                gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑Taa                                                  <  1:2938315/70‑1 (MQ=255)
                 tAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCa                                                                    <  1:1170022/52‑1 (MQ=255)
                 tAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCa                                                                    <  1:3943261/52‑1 (MQ=255)
                 tAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑Taaa                                                 <  1:3301694/70‑1 (MQ=255)
                  aGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑Taaaa                                                <  1:2017284/70‑1 (MQ=255)
                  aGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAAt                                               >  1:3537989/1‑71 (MQ=255)
                    tGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTAc                                                               <  1:504547/54‑1 (MQ=255)
                    tGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTAc                                                               <  1:3371468/54‑1 (MQ=255)
                     gCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGa                                             >  1:2445047/1‑70 (MQ=255)
                      cgcaaGATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAG‑GCGAAAAAATGaaa                                           <  1:1657235/68‑1 (MQ=255)
                      cgcaaGATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAG‑GCGAAAAAATGaaa                                           <  1:2209642/68‑1 (MQ=255)
                       gcaaGATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGga                                                           >  1:765425/3‑55 (MQ=255)
                        tAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGa                                             <  1:3772670/67‑1 (MQ=255)
                        tAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAAt                                         <  1:751456/71‑1 (MQ=255)
                        tAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGAG‑TAAAATGAAAAt                                         <  1:3862461/71‑1 (MQ=255)
                        caaGATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGgagg                                                         >  1:431940/2‑55 (MQ=255)
                        caaGATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAG‑GCGAAAAAATGAAAAt                                         >  1:2576504/2‑71 (MQ=255)
                        caaGATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAG‑GCGAAAAAATGAAAAt                                         >  1:3320715/2‑71 (MQ=255)
                          aGATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAG‑GCGAAAAAATGAAAAtcc                                       >  1:2944020/1‑69 (MQ=255)
                          aGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGga                                                           >  1:2597219/1‑52 (MQ=255)
                          aGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTc                                       >  1:3828491/1‑71 (MQ=255)
                          aGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTc                                       >  1:614527/1‑71 (MQ=255)
                          aGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTc                                       <  1:1001081/71‑1 (MQ=255)
                           gATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAG‑GCGAAAAAATGAAAATCCg                                      <  1:1923858/71‑1 (MQ=255)
                           gATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAG‑GCGAAAAAATGAAAATCCg                                      <  1:148799/71‑1 (MQ=255)
                               aagaaCGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAG‑GCGAAAAAATGAAAATCCGTTca                                  <  1:1018540/71‑1 (MQ=255)
                                agaaCGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAG‑GCGAAAAAATGAAAATCCGTTcac                                 <  1:1445858/71‑1 (MQ=255)
                                 gaaCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTcac                                 <  1:1311226/70‑1 (MQ=255)
                                   aCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTc                                       >  1:2275607/1‑62 (MQ=255)
                                     aTTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGaaaa                                          <  1:3940656/57‑1 (MQ=255)
                                        gACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCgg                          <  1:3132414/70‑1 (MQ=255)
                                        gACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCgg                          <  1:1714093/70‑1 (MQ=255)
                                          cTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCAt                       >  1:1855724/1‑71 (MQ=255)
                                          cTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCAt                       >  1:1350295/1‑71 (MQ=255)
                                           tGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTc                                       >  1:1875490/1‑54 (MQ=255)
                                           tGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGtt                                    >  1:2995808/1‑57 (MQ=255)
                                           tGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATg                      >  1:595237/1‑71 (MQ=255)
                                             gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATgg                     <  1:2762370/70‑1 (MQ=255)
                                             gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATgg                     <  1:1763620/70‑1 (MQ=255)
                                             gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATgg                     <  1:3760649/70‑1 (MQ=255)
                                             gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATg                      <  1:1628791/69‑1 (MQ=255)
                                             gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGt                    >  1:2219106/1‑71 (MQ=255)
                                             gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGt                    <  1:3780801/71‑1 (MQ=255)
                                             gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGt                    >  1:2223794/1‑71 (MQ=255)
                                              ttAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTg                   >  1:3238161/1‑71 (MQ=255)
                                               tAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTcaca                                >  1:2932218/1‑57 (MQ=255)
                                               tAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTcaca                                >  1:1305391/1‑57 (MQ=255)
                                                aGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGaaaa                                          <  1:484874/46‑1 (MQ=255)
                                                aGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCt                 >  1:265481/1‑71 (MQ=255)
                                                 gATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCa                        <  1:1642333/63‑1 (MQ=255)
                                                 gATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCt                 <  1:3464655/70‑1 (MQ=255)
                                                     gCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCgg                          <  1:1154401/57‑1 (MQ=255)
                                                     gCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCa                        <  1:1689228/59‑1 (MQ=255)
                                                     gCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAt             <  1:3189656/70‑1 (MQ=255)
                                                      cGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtct           >  1:3818181/1‑71 (MQ=255)
                                                      cGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtct           <  1:1833539/71‑1 (MQ=255)
                                                       gAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGa                                             <  1:123019/36‑1 (MQ=255)
                                                       gAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCa                        >  1:713048/1‑57 (MQ=255)
                                                        aaGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtctc          >  1:1295489/1‑70 (MQ=255)
                                                         aGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtct           >  1:1676601/1‑68 (MQ=255)
                                                         aGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGa        >  1:632673/1‑71 (MQ=255)
                                                         aGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGa        <  1:665169/71‑1 (MQ=255)
                                                          ggCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGTCATGGTGCTGAAtc            <  1:1098645/66‑1 (MQ=255)
                                                           gCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGt                                     <  1:2951322/40‑1 (MQ=255)
                                                            cAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGAt       >  1:2548940/1‑69 (MQ=255)
                                                            cAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATaa     <  1:2397116/71‑1 (MQ=255)
                                                             aaTGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGtt                                    <  1:1026738/39‑1 (MQ=255)
                                                             aaTGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTcaca                                >  1:3176505/1‑43 (MQ=255)
                                                              aTGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATgg                     <  1:260240/53‑1 (MQ=255)
                                                               tGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTcac                                 >  1:1793006/1‑40 (MQ=255)
                                                               tGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTg  >  1:2816562/1‑71 (MQ=255)
                                                                  |                                                                    
CCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTG  >  minE/835196‑835329

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: ‑