Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,894 T→C 6.3% R204R (CGT→CGC narH → nitrate reductase 1, beta (Fe‑S) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,8940TC 6.3% 272.8 / 7.5 95R204R (CGT→CGCnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base T (43/46);  new base C (3/3);  total (46/49)
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

CGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGG  >  minE/835826‑835959
                                                                    |                                                                 
cgcCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTg                                                                  >  1:1196786/1‑70 (MQ=255)
  ccTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaag                                                               <  1:2951842/71‑1 (MQ=255)
    tgtgCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaag                                                               >  1:3680443/1‑69 (MQ=255)
    tgtgCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaag                                                               >  1:1926426/1‑69 (MQ=255)
      tgCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                                             <  1:2590296/69‑1 (MQ=255)
      tgCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                                             <  1:669598/69‑1 (MQ=255)
      tgCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGaagaag                                                            <  1:2782525/70‑1 (MQ=255)
        cGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATg                                                         >  1:3143601/1‑71 (MQ=255)
        cGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATg                                                         >  1:1010523/1‑71 (MQ=255)
          aaCACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATgg                                                        >  1:1651483/1‑70 (MQ=255)
          aaCACTGCCTCAATCCCAGATGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCGAAGAAGACgg                                                        >  1:2350301/1‑70 (MQ=255)
           acacTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                                             >  1:2775066/1‑64 (MQ=255)
           acacTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaa                                                                >  1:3369353/1‑61 (MQ=255)
           acacTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCa                                                      <  1:71421/71‑1 (MQ=255)
             acTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaag                                                               <  1:3717996/60‑1 (MQ=255)
               tGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                                             <  1:294023/60‑1 (MQ=255)
               tGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGt                                                  >  1:3699780/1‑71 (MQ=255)
               tGCCTGAACCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                                             <  1:39978/59‑1 (MQ=255)
                 ccTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTcc                                                >  1:911336/1‑71 (MQ=255)
                 ccTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTcc                                                >  1:4021566/1‑71 (MQ=255)
                  cTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                                               <  1:1308139/71‑1 (MQ=255)
                  cTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                                               <  1:2695292/71‑1 (MQ=255)
                     aaCCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                                               >  1:757825/1‑68 (MQ=255)
                      aCCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga                                                           >  1:1386404/1‑55 (MQ=255)
                      aCCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATc                                           >  1:3624603/1‑71 (MQ=255)
                      aCCCGGCATGTGTGGCGACATGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga                                                           >  1:3363334/1‑55 (MQ=255)
                       cccGGCATGTGTGGCGACCTGCGCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATc                                           >  1:4018316/1‑70 (MQ=255)
                       cccGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATc                                           >  1:3070801/1‑70 (MQ=255)
                        ccGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                                             <  1:3014250/51‑1 (MQ=255)
                         cGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAc                                        <  1:4018509/71‑1 (MQ=255)
                          ggCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga                                                           >  1:1090911/1‑51 (MQ=255)
                          ggCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga                                                           >  1:3533152/1‑51 (MQ=255)
                              tgtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                                             <  1:3252888/45‑1 (MQ=255)
                              tgtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                                             <  1:1350903/45‑1 (MQ=255)
                              tgtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGc                                                       <  1:1662162/51‑1 (MQ=255)
                              tgtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTg                                              <  1:2435676/60‑1 (MQ=255)
                              tgtgtgGCGACATGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGa                                   <  1:1152876/71‑1 (MQ=255)
                               gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                                               <  1:35998/58‑1 (MQ=255)
                               gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                                               <  1:3808084/58‑1 (MQ=255)
                               gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                                               <  1:1035948/58‑1 (MQ=255)
                                tgtgGCGACCTGCCCGAGCGGTTCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATa                                 <  1:3843982/71‑1 (MQ=255)
                                tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                                             >  1:1343375/1‑43 (MQ=255)
                                tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt                                  <  1:3039922/70‑1 (MQ=255)
                                cgtgGCGACCTGCCCAAGCGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGa                                   <  1:201995/68‑1 (MQ=255)
                                   ggCGCCCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTcgt                                               >  1:2612787/1‑52 (MQ=255)
                                    gCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAAt                              <  1:1605920/70‑1 (MQ=255)
                                      gACCTGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga                                                           >  1:2066572/1‑39 (MQ=255)
                                      gACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCg                                          <  1:2757191/56‑1 (MQ=255)
                                       aCCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                                             <  1:2988911/36‑1 (MQ=255)
                                        ccTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGa                                         >  1:1070637/1‑55 (MQ=255)
                                        ccTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAc                                        >  1:2224999/1‑56 (MQ=255)
                                         cTGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga                                                           >  1:110676/1‑36 (MQ=255)
                                         cTGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga                                                           >  1:2541091/1‑36 (MQ=255)
                                         cTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGa                                         >  1:1070854/1‑54 (MQ=255)
                                         cTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtg                        >  1:2985414/1‑71 (MQ=255)
                                             ccGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtggctg                    <  1:1017280/71‑1 (MQ=255)
                                              cGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtggctgg                   >  1:853566/1‑71 (MQ=255)
                                               gAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTg                                              <  1:325348/43‑1 (MQ=255)
                                                aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAc                                        <  1:2819306/48‑1 (MQ=255)
                                                aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAc                                        <  1:1439225/48‑1 (MQ=255)
                                                aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa                               >  1:697619/1‑57 (MQ=255)
                                                aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa                               >  1:2295813/1‑57 (MQ=255)
                                                aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtggctg                    >  1:727784/1‑68 (MQ=255)
                                                aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                 >  1:1429452/1‑71 (MQ=255)
                                                aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                 >  1:1832239/1‑71 (MQ=255)
                                                aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                 <  1:3693529/71‑1 (MQ=255)
                                                  cGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTa               <  1:1031246/71‑1 (MQ=255)
                                                  cGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGtt               <  1:3105730/71‑2 (MQ=255)
                                                     tGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATa                                 <  1:2986992/50‑1 (MQ=255)
                                                     tGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                 <  1:3485686/66‑1 (MQ=255)
                                                      gCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtgtg           >  1:3744958/1‑71 (MQ=255)
                                                       cGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGAt                                            >  1:578630/1‑37 (MQ=255)
                                                       cGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGAt                                            >  1:3754100/1‑37 (MQ=255)
                                                       cGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtgtg           <  1:3394863/70‑1 (MQ=255)
                                                          tttACAAGCGTGAAGAAGATGGCATCGTCCTGATGGACCAGGATAAATGCCGt                         <  1:4032851/53‑1 (MQ=255)
                                                          tttACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtg                        <  1:2983292/54‑1 (MQ=255)
                                                          tttACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATc       <  1:1977388/71‑1 (MQ=255)
                                                            tACAAGCGTGAAGAAGATGGCATCGTCCTGATCGa                                         >  1:4123587/1‑35 (MQ=255)
                                                            tACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCAc     >  1:1857657/1‑71 (MQ=255)
                                                            tACAAGAGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCa      >  1:539242/1‑70 (MQ=255)
                                                            tACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGtggc                      >  1:2352836/1‑54 (MQ=255)
                                                            tACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGtggc                      >  1:1667466/1‑54 (MQ=255)
                                                             aCAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACt    >  1:2751332/1‑71 (MQ=255)
                                                              cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg   <  1:2122629/71‑1 (MQ=255)
                                                              cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg   <  1:2085022/71‑1 (MQ=255)
                                                              cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg   <  1:1445361/71‑1 (MQ=255)
                                                              cAAGCGCGAGGAGGACGGCATCGTCCTGATCGACCAGGACAAGTGCCGCGGCTGGCGGATGTGCa         <  1:1434330/65‑1 (MQ=255)
                                                               aaGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATa                                 <  1:1203018/40‑1 (MQ=255)
                                                               aaGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTgg  <  1:3258100/71‑1 (MQ=255)
                                                               aaGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTgg  <  1:770933/71‑1 (MQ=255)
                                                                 gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa                               >  1:3891022/1‑40 (MQ=255)
                                                                 gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa                               >  1:3903465/1‑40 (MQ=255)
                                                                 gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                 <  1:3126887/54‑1 (MQ=255)
                                                                  cGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtggc                      <  1:1464355/48‑1 (MQ=255)
                                                                  cGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                 <  1:3506649/53‑1 (MQ=255)
                                                                    tGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGt                <  1:917889/52‑1 (MQ=255)
                                                                    |                                                                 
CGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGG  >  minE/835826‑835959

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: