Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,244 A→G 14.1% L1234L (TTA→TTG) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,2440AG14.1% 208.8 / 23.4 85L1234L (TTA→TTG) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (37/36);  new base G (6/6);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCT  >  minE/835175‑835314
                                                                     |                                                                       
ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTa                                                                         <  1:2382585/70‑1 (MQ=255)
ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTa                                                                         <  1:1907475/70‑1 (MQ=255)
ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTa                                                                         <  1:3220800/70‑1 (MQ=255)
ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTa                                                                         <  1:3456589/70‑1 (MQ=255)
  acggCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGTTAACTGGCTGGAt                                                                      >  1:2152391/3‑71 (MQ=255)
   tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGAt                                                                      <  1:1095397/70‑1 (MQ=255)
        ccGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAGGAACATTGACTGGTTAGATGGCGaa                                                                <  1:487614/71‑1 (MQ=255)
        ccGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGaa                                                                <  1:725893/71‑1 (MQ=255)
           tAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGc                                                             <  1:464788/71‑1 (MQ=255)
               ttCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATg                                                         <  1:1214518/71‑1 (MQ=255)
                tCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATg                                                                     >  1:1144809/1‑58 (MQ=255)
                  tAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAAt                                                          >  1:82640/1‑67 (MQ=255)
                   aaCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAAt                                                          <  1:2957897/66‑1 (MQ=255)
                    aCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATg                                                                     >  1:2595832/1‑54 (MQ=255)
                    aCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATg                                                                     <  1:2070912/54‑1 (MQ=255)
                    aCCGAGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAAt                                                          >  1:2624426/1‑65 (MQ=255)
                     ccGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAAt                                                          >  1:1615579/1‑64 (MQ=255)
                     ccGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAg                                                    <  1:3194256/70‑1 (MQ=255)
                      cgcgATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGt                                                  <  1:3802636/71‑1 (MQ=255)
                         gatgatTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACa                                               <  1:652433/71‑1 (MQ=255)
                         gATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGa                                                                 <  1:3274004/53‑1 (MQ=255)
                         gATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGa                                                                 <  1:370518/53‑1 (MQ=255)
                         gATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGa                                                                 <  1:3853295/53‑1 (MQ=255)
                          aTGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAg                                              <  1:853548/71‑1 (MQ=255)
                            gAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAg                                              >  1:959944/1‑69 (MQ=255)
                             aGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCaa                                                           >  1:3146941/1‑55 (MQ=255)
                             aGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGga                                            <  1:3819905/70‑1 (MQ=255)
                             aGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGga                                            <  1:1839570/70‑1 (MQ=255)
                                ttgttgTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGgaga                                          >  1:780615/1‑69 (MQ=255)
                                ttgttgTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGgaga                                          >  1:2702740/1‑69 (MQ=255)
                                 tgttgtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGaa                                                                >  1:3896716/1‑46 (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)
                                     gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGa                                                                       >  1:294759/1‑35 (MQ=255)
                                     gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGAt                                                                      >  1:79511/1‑36 (MQ=255)
                                     gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAg                                              >  1:1157715/1‑60 (MQ=255)
                                     gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑Taa                                   >  1:328648/1‑70 (MQ=255)
                                     gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑Taa                                   >  1:3759474/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:2209642/68‑1 (MQ=255)
                                           cgcaaGATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAG‑GCGAAAAAATGaaa                            <  1:1657235/68‑1 (MQ=255)
                                           cGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAg                                                               >  1:2462320/1‑37 (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)
                                               aGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATg                                                         >  1:138884/1‑39 (MQ=255)
                                               aGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGga                                            >  1:2597219/1‑52 (MQ=255)
                                               aGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTc                        >  1:614527/1‑71 (MQ=255)
                                               aGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTc                        <  1:1001081/71‑1 (MQ=255)
                                               aGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTc                        >  1:3828491/1‑71 (MQ=255)
                                                gATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAG‑GCGAAAAAATGAAAATCCg                       <  1:148799/71‑1 (MQ=255)
                                                gATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAG‑GCGAAAAAATGAAAATCCg                       <  1:1923858/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:1350295/1‑71 (MQ=255)
                                                               cTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCAt        >  1:1855724/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:3760649/70‑1 (MQ=255)
                                                                  gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATgg      <  1:2762370/70‑1 (MQ=255)
                                                                  gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATgg      <  1:1763620/70‑1 (MQ=255)
                                                                  gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATg       <  1:1628791/69‑1 (MQ=255)
                                                                  gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGt     >  1:2223794/1‑71 (MQ=255)
                                                                  gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGt     >  1:2219106/1‑71 (MQ=255)
                                                                  gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGt     <  1:3780801/71‑1 (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)
                                                                     |                                                                       
CTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCT  >  minE/835175‑835314

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