Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 835,891 G→A 13.7% 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
*Exported835,8910GA13.7% 238.8 / 22.9 95K203K (AAG→AAA) narHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base G (27/55);  new base A (4/9);  total (31/64)
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

TGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGG  >  Exported/835821‑835959
                                                                       |                                                                    
tGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACATGCCCGAGCGGTGCGATTTACAAg                                                                      >  1:4295982/1‑71 (MQ=255)
  ccGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCg                                                                    >  1:4131440/1‑71 (MQ=255)
  ccGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGT‑GCGACCTGCCCGAGCGGTGCGATTTACAAGCGt                                                                   >  1:4347390/1‑71 (MQ=255)
      gcCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTga                                                                 <  1:1424948/70‑1 (MQ=255)
      gcCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCgaa                                                                <  1:4008683/71‑1 (MQ=255)
         tgtgCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCgaag                                                               <  1:5114094/69‑1 (MQ=255)
          gtgCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaa                                                                >  1:3075944/1‑67 (MQ=255)
           tgCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                                             <  1:5674081/69‑1 (MQ=255)
           tgCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCgaagaa                                                             <  1:1690358/69‑1 (MQ=255)
           tgCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCgaagaa                                                             <  1:5420271/69‑1 (MQ=255)
             cGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCg                                                                    <  1:4920128/60‑1 (MQ=255)
             cGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTg                                                                  >  1:1109417/1‑62 (MQ=255)
              gAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga                                                           >  1:5406577/1‑68 (MQ=255)
               aaCACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATgg                                                        >  1:2365420/1‑70 (MQ=255)
               aaCACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCGAagaagacg                                                         <  1:809713/69‑3 (MQ=255)
                   cTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                                             <  1:4544790/61‑1 (MQ=255)
                   cTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCg                                                   <  1:3433784/71‑1 (MQ=255)
                        tGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATTGCATCGTCCt                                               >  1:854048/1‑70 (MQ=255)
                        tGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATc                                                    >  1:1163641/1‑65 (MQ=255)
                        tGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                                               >  1:5067144/1‑70 (MQ=255)
                          aa‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaag                                                               >  1:3137770/1‑52 (MQ=255)
                          aa‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGGGAAGAAGATGGCATCGTCCTGAt                                            <  1:4859008/71‑1 (MQ=255)
                             cccGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCg                                          <  1:4446315/71‑1 (MQ=255)
                              ccGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                                             <  1:998695/51‑1 (MQ=255)
                               cGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGaagaag                                                            <  1:3814880/51‑1 (MQ=255)
                               cGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAc                                        <  1:2858466/71‑1 (MQ=255)
                                ggCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGaagaag                                                            <  1:2387615/50‑1 (MQ=255)
                                ggCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAcc                                       <  1:4683953/71‑1 (MQ=255)
                                  cATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGaagaag                                                            <  1:4593058/48‑1 (MQ=255)
                                    tgtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                                             <  1:2313650/45‑1 (MQ=255)
                                    tgcgtgGCGACCTGCCCAAGCGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGa                                   >  1:2592635/4‑71 (MQ=255)
                                     gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                                             <  1:1287185/44‑1 (MQ=255)
                                     gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                                             <  1:1390108/44‑1 (MQ=255)
                                     gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                                             <  1:5629327/44‑1 (MQ=255)
                                     gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                                             <  1:543383/44‑1 (MQ=255)
                                     gtgtgGCGACCTGCCCGAGCGGGGCGATTTACAAGCGTgaagaa                                                             <  1:1740031/44‑1 (MQ=255)
                                     gcgtgGCGACCTGCCCAAGCGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGa                                   >  1:4385442/3‑70 (MQ=255)
                                      tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt                                  >  1:456002/1‑70 (MQ=255)
                                      tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt                                  >  1:4396349/1‑70 (MQ=255)
                                      tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt                                  >  1:1775170/1‑70 (MQ=255)
                                      tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt                                  >  1:534635/1‑70 (MQ=255)
                                      tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt                                  >  1:2016248/1‑70 (MQ=255)
                                          gCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa                               >  1:5052181/1‑69 (MQ=255)
                                               cTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTcc                                                >  1:2232410/1‑47 (MQ=255)
                                               cTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtg                        <  1:4317549/71‑1 (MQ=255)
                                                 gCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTg                                              >  1:3817871/1‑47 (MQ=255)
                                                  cccGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtggc                      <  1:512664/70‑1 (MQ=255)
                                                  cccAAGCGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGtggc                      <  1:2476477/70‑1 (MQ=255)
                                                    cGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtggctgg                   >  1:5396429/1‑71 (MQ=255)
                                                      aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                 <  1:4305586/71‑1 (MQ=255)
                                                        cGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATg                             <  1:4212846/57‑1 (MQ=255)
                                                        cGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGtt               <  1:3102966/71‑2 (MQ=255)
                                                         ggTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAt              <  1:2869767/71‑1 (MQ=255)
                                                          gTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAcc                                       <  1:4924439/45‑1 (MQ=255)
                                                          gTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAcc                                       <  1:3894254/45‑1 (MQ=255)
                                                          gTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAcc                                       <  1:3645624/45‑1 (MQ=255)
                                                           tGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATa                                 <  1:4264622/50‑1 (MQ=255)
                                                           tGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtggc                      <  1:2644161/61‑1 (MQ=255)
                                                           tGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTa               <  1:1157081/68‑1 (MQ=255)
                                                           tGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtgt            <  1:3130205/71‑1 (MQ=255)
                                                           tGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATAGACCAGGATa                                 <  1:4733414/50‑1 (MQ=255)
                                                            gCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtgtg           >  1:3749414/1‑71 (MQ=255)
                                                            gCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtgtg           >  1:3141955/1‑71 (MQ=255)
                                                             cGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtgtg           >  1:5115664/1‑70 (MQ=255)
                                                               aTTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCAt        <  1:3329174/71‑1 (MQ=255)
                                                               aTTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCAt        <  1:2886675/71‑1 (MQ=255)
                                                               aTTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCAt        <  1:4883487/71‑1 (MQ=255)
                                                                tttACAAGCGTGAAGAAGATGGCATCGTCCTGATc                                           <  1:936548/35‑1 (MQ=255)
                                                                tttACAAGCGTGAAGAAGATGGCATCGTCCTGATc                                           <  1:821844/35‑1 (MQ=255)
                                                                tttACAAGCGTGAAGAAGATGGCATCGTCCTGATc                                           <  1:1611923/35‑1 (MQ=255)
                                                                tttACAAGCGTGAAGAAGATGGCATCGTCCTGATc                                           <  1:2013520/35‑1 (MQ=255)
                                                                tttACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATa                                 <  1:3165576/45‑1 (MQ=255)
                                                                tttACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATa                                 <  1:262897/45‑1 (MQ=255)
                                                                tttACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATc       <  1:870498/71‑1 (MQ=255)
                                                                tttACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATc       <  1:2778101/71‑1 (MQ=255)
                                                                tttACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATc       <  1:564971/71‑1 (MQ=255)
                                                                tttACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATc       <  1:371163/71‑1 (MQ=255)
                                                                tctACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGtt               <  1:4137480/61‑2 (MQ=255)
                                                                 ctACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCAt        >  1:4756963/2‑69 (MQ=255)
                                                                  tACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATc       <  1:4344729/69‑1 (MQ=255)
                                                                  tACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATc       <  1:2259651/69‑1 (MQ=255)
                                                                  tACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATc       <  1:4979942/69‑1 (MQ=255)
                                                                  tACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATc       <  1:2745206/69‑1 (MQ=255)
                                                                  tACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCa      >  1:4708042/1‑70 (MQ=255)
                                                                  tACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCa      <  1:1514189/70‑1 (MQ=255)
                                                                  tACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCAc     >  1:5361167/1‑71 (MQ=255)
                                                                  tACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCAc     >  1:5483122/1‑71 (MQ=255)
                                                                  tACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCAtaag     >  1:3331763/1‑68 (MQ=255)
                                                                    cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg   <  1:2055133/71‑1 (MQ=255)
                                                                    cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg   <  1:1800628/71‑1 (MQ=255)
                                                                    cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg   <  1:1970781/71‑1 (MQ=255)
                                                                    cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg   <  1:5134126/71‑1 (MQ=255)
                                                                    cAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCAtaag     <  1:1404041/69‑4 (MQ=255)
                                                                     aaGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTgg  <  1:1942584/71‑1 (MQ=255)
                                                                      aGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtg                        >  1:1531544/1‑48 (MQ=255)
                                                                      aGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtg                        >  1:2214513/1‑48 (MQ=255)
                                                                       gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaa                                >  1:755488/1‑39 (MQ=255)
                                                                       gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaa                                >  1:903774/1‑39 (MQ=255)
                                                                       gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtg             >  1:128815/1‑58 (MQ=255)
                                                                       |                                                                    
TGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGG  >  Exported/835821‑835959

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