Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,891 G→A 7.6% K203K (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
*minE835,8910GA 7.6% 292.7 / 10.1 92K203K (AAG→AAAnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base G (41/44);  new base A (3/4);  total (44/48)
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

TGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGG  >  minE/835821‑835959
                                                                       |                                                                    
tGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAg                                                                      >  1:1829469/1‑71 (MQ=255)
tGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAg                                                                      >  1:4045314/1‑71 (MQ=255)
    gcgcCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACgcg                                                                  <  1:2385307/71‑3 (MQ=255)
    gcgcCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACgcg                                                                  <  1:373081/71‑3 (MQ=255)
     cgcCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTg                                                                  >  1:1196786/1‑70 (MQ=255)
       ccTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaag                                                               <  1:2951842/71‑1 (MQ=255)
         tgtgCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaag                                                               >  1:1926426/1‑69 (MQ=255)
         tgtgCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaag                                                               >  1:3680443/1‑69 (MQ=255)
           tgCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                                             <  1:669598/69‑1 (MQ=255)
           tgCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                                             <  1:2590296/69‑1 (MQ=255)
           tgCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGaagaag                                                            <  1:2782525/70‑1 (MQ=255)
             cGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATg                                                         >  1:1010523/1‑71 (MQ=255)
             cGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATg                                                         >  1:3143601/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:4021566/1‑71 (MQ=255)
                      ccTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTcc                                                >  1:911336/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:3533152/1‑51 (MQ=255)
                                ggCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga                                                           >  1:1090911/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:1035948/58‑1 (MQ=255)
                                     gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                                               <  1:3808084/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:1439225/48‑1 (MQ=255)
                                                      aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAc                                        <  1:2819306/48‑1 (MQ=255)
                                                      aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa                               >  1:2295813/1‑57 (MQ=255)
                                                      aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa                               >  1:697619/1‑57 (MQ=255)
                                                      aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtggctg                    >  1:727784/1‑68 (MQ=255)
                                                      aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                 >  1:1429452/1‑71 (MQ=255)
                                                      aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                 <  1:3693529/71‑1 (MQ=255)
                                                      aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                 >  1:1832239/1‑71 (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:3754100/1‑37 (MQ=255)
                                                             cGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGAt                                            >  1:578630/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:1445361/71‑1 (MQ=255)
                                                                    cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg   <  1:2122629/71‑1 (MQ=255)
                                                                    cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg   <  1:2085022/71‑1 (MQ=255)
                                                                    cAAGCGCGAGGAGGACGGCATCGTCCTGATCGACCAGGACAAGTGCCGCGGCTGGCGGATGTGCa         <  1:1434330/65‑1 (MQ=255)
                                                                     aaGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATa                                 <  1:1203018/40‑1 (MQ=255)
                                                                     aaGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTgg  <  1:770933/71‑1 (MQ=255)
                                                                     aaGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTgg  <  1:3258100/71‑1 (MQ=255)
                                                                       gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa                               >  1:3903465/1‑40 (MQ=255)
                                                                       gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa                               >  1:3891022/1‑40 (MQ=255)
                                                                       gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                 <  1:3126887/54‑1 (MQ=255)
                                                                       |                                                                    
TGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGG  >  minE/835821‑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: